How robust are particle physics predictions in asymptotic safety?

The framework of trans-Planckian asymptotic safety has been shown to generate phenomenological predictions in the Standard Model and in some of its simple new physics extensions. A heuristic approach is often adopted, which bypasses the functional renormalization group by relying on a parametric description of quantum gravity with universal coefficients that are eventually obtained from low-energy observations. Within this approach a few simplifying approximations are typically introduced, including the computation of matter renormalization group equations at 1 loop, an arbitrary definition of the position of the Planck scale at 1019GeV,\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$10^{19}\,\textrm{GeV},$$\end{document} and an instantaneous decoupling of gravitational interactions below the Planck scale. In this work we systematically investigate, both analytically and numerically, the impact of dropping each of those approximations on the predictions for certain particle physics scenarios. In particular we study two extensions of the Standard Model, the gauged B-L\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B-L$$\end{document} model and the leptoquark S3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$S_3$$\end{document} model, for which we determine a set of irrelevant gauge and Yukawa couplings. In each model, we present numerical and analytical estimates of the uncertainties associated with the predictions from asymptotic safety.

The framework of trans-Planckian asymptotic safety has been shown to generate phenomenological predictions in the Standard Model and in some of its simple new physics extensions.A heuristic approach is often adopted, which bypasses the functional renormalization group by relying on a parametric description of quantum gravity with universal coefficients that are eventually obtained from lowenergy observations.Within this approach a few simplifying approximations are typically introduced, including the computation of matter renormalization group equations at 1 loop, an arbitrary definition of the position of the Planck scale at 10 19 GeV, and an instantaneous decoupling of gravitational interactions below the Planck scale.In this work we systematically investigate, both analytically and numerically, the impact of dropping each of those approximations on the predictions for certain particle physics scenarios.In particular we study two extensions of the Standard Model, the gauged B − L model and the leptoquark S 3 model, for which we determine a set of irrelevant gauge and Yukawa couplings.In each model, we present numerical and analytical estimates of the uncertainties associated with the predictions from asymptotic safety.

Introduction
Asymptotic safety (AS) is the property of a quantum field theory to develop ultraviolet (UV) fixed points of the renormalization group (RG) flow of the action [1].Following the development of functional renormalization group (FRG) techniques three decades ago [2,3], it was shown in numerous papers that AS may arise quite naturally in quantum gravity and provide the key ingredient for the non-perturbative renormalizability of the theory.Fixed points for the rescaled Newton coupling and the cosmological constant were identified initially in the Einstein-Hilbert truncation of the effective action [4][5][6], based on two operators, and later confirmed in the presence of gravitational operators of increasing mass dimension [7][8][9][10][11][12][13][14][15], and of matter-field operators [16][17][18].
From the point of view of particle physics in four space-time dimensions, a particularly exciting possibility is that not only the gravitational action but the full system of gravity and matter may feature UV fixed points in the energy regime where gravitational interactions become strong [19][20][21][22][23][24][25][26][27][28][29][30].A trans-Planckian fixed point may provide in that case specific boundary conditions for some of the a priori free couplings of the matter Lagrangian, as long as they correspond to irrelevant directions in theory space.In this context, indications of AS emerging in U(1) gauge theories have led to discovering a gravity-driven solution to the triviality problem of the Standard Model (SM) hypercharge coupling [31][32][33].Early achievements of embedding realistic matter systems in AS include the ballpark prediction for the value of the Higgs mass (more precisely, of the quartic coupling of the Higgs potential) obtained a few years ahead of its discovery [34] (see also Refs.[35][36][37]), and the retroactive "postdiction" of the top-mass value [38].
In order to properly complete a matter system with trans-Planckian AS, one should consistently calculate gravitational corrections to the matter beta function using the formalism of the FRG.The ensuing high-scale boundary conditions for the matter couplings, obtained this way from first principles, should be then compared with observations (cross sections, decay rates, etc.) after running the couplings to the low scale.It has been long known, however, that calculations based on the FRG can be subject to large theoretical uncertainties, stemming from a variety of sources -from the choice of truncation in the gravity action [7,12,[39][40][41], to cutoff-scheme dependence [6,42], to the backreaction of matter, which may introduce a dependence of the gravitational fixed point on the specifics of the matter sector [43][44][45].Various results can differ by up to several times [45] so that, often, a precise quantitative determination of the high-scale value of the matter couplings is not possible.Moreover, it has been recently pointed out [46,47] that some conceptual questions regarding the appropriateness of the FRG itself as the right tool for investigating/enforcing AS in quantum gravity remain open.
On the other hand, because of the universal nature of trans-Planckian interactions, and thanks to the "rigidity" of the dynamics in the vicinity of an interactive fixed point -by virtue of which all quantities can be predicted except for a handful of relevant parameters that will have to be determined experimentally -a first-principle calculation of the gravitational contribution to the matter couplings is not necessarily needed to prove the consistency of certain low-energy predictions with quantum gravity.Often one is content with establishing a heuristic framework in which the trans-Planckian interactions are parameterized by coefficients that are eventually obtained from low-energy observations [48][49][50][51][52][53][54][55][56][57][58].
In such an effective trans-Planckian embedding, one generally introduces parametric corrections to the renormalization group equations (RGEs) of the renormalizable matter couplings, which take the form ) where t = ln µ (renormalization scale), g i and y j (with i, j = 1, 2, . . . ) indicate, respectively, the set of gauge and Yukawa couplings of the theory, and β matter i,j are the beta functions of the matter theory, which can be evaluated at 1 loop in dimensional regularization (DREG). 1 The two "gravitational" coefficients f g and f y multiply linearly all matter couplings of the same kind.They are thus universal, in the sense that gravity is not expected to be affected by the internal degrees of freedom of the matter system.They appear only in the regime where the gravitational action develops an interactive fixed point, at µ > M Pl , and serve the purpose of inducing trans-Planckian zeros on the matter beta functions.If some of the emerging fixedpoint coupling values correspond to irrelevant directions of the RG flow, one can estimate f g and f y by requiring that the irrelevant fixed points should be connected, along a unique RG trajectory, to quantities measured in experiments at the low scale.
Such heuristic embedding of a gauge-Yukawa system in trans-Planckian AS has been used in the SM to attempt a prediction of the top/bottom mass ratio [48], the Cabibbo-Kobayashi-Maskawa [51], and Pontecorvo-Maki-Nakagawa-Sakata [54] matrix elements.In the context of neutrino mass generation, this effective approach was employed to theorize the interesting possibility that the tiny Yukawa couplings of Dirac-type right-handed neutrinos may find their origin in the dynamics of the trans-Planckian flow [54,58].New physics (NP) predictions were extracted for neutrino masses [49], leptoquarks [52], vector-like fermions [53], and γ/Z ′ kinetic mixing [55].An expression similar to Eqs. ( 1) and ( 2) potentially applies to the quartic couplings of the scalar potential as well, to which beta functions one could in principle assign a universal correction factor f λ .Following this, predictions were made for the relic abundance of dark matter [50,60], baryon number [56,57], as well as axion models [61]. 2hile in the absence of a fully developed theory of quantum gravity the heuristic approach described above has proven to be extremely fruitful for phenomenological studies in particle physics, it is also important to be aware that it is based on several simplifying approximations: • The DREG matter beta functions are typically computed at 1-loop • The Planck scale is set arbitrarily at M Pl = 10 19 GeV • The scale dependence of f g and f y , which should parameterize the cross-over from the interactive to non-interactive regime of quantum gravity, is neglected.f g and f y are treated as constants above the Planck scale and are set to zero below.In other words, gravity decouples instantaneously at M Pl = 10 19 GeV.
The question then naturally arises as to how robust the predictions derived in this way can be considered, and to what extent dropping any of the approximations listed above may affect a potential observational strategy to test these predictions at the low scale.
In this study we attempt to address the issue in a systematic way.We analyze the effects of discarding one by one the approximations of the minimal parametric setup.Specifically we consider • the inclusion of higher-order corrections in the matter sector • changing the position of the Planck scale by a few orders of magnitude • the non-trivial functional dependence of the running gravitational couplings, f g,y (t), resulting in the non-instantaneous decoupling of the trans-Planckian UV completion.
We limit ourselves to the study of gauge-Yukawa systems of the type (1) and ( 2).As we have mentioned above, the dimensionless parameters of the scalar potential lie on a slightly different footing.Besides the issue of whether the gravitational correction is multiplicative or not, it has also been shown that in the SM [36] and some models of NP [55] it is difficult to obtain the precise value of the Higgs mass if the flow originates from a UV irrelevant fixed point, so that the predictivity of AS for the scalar potential is often in question and a model-dependent issue.On the other hand, from the point of view of the RG flow, the parameters of the scalar potential result somewhat "decoupled" from the gauge-Yukawa system, as they can only affect Eq. ( 1) from the third-loop level up, and Eq. ( 2) from the second loop.We shall quantify, when necessary, the size of the corrections that unknown parameters of the scalar potential induce on the gauge-Yukawa system via higher-order contributions.
Our main finding is that the predictions in the gauge sector are extremely robust.Each and every one of the effects listed above induce an uncertainty that does not exceed the one percent level.The situation is slightly more intricate in the Yukawa sector, where the uncertainty strongly depends on the size of the predicted coupling itself.We find that, while Yukawas as large as the SM hypercharge coupling are subject to uncertainties not exceeding the five percent level, if they turn out to be much smaller -as a result of fine tuning -their uncertainties increase significantly.
The paper is organized as follows.In Sec. 2 we introduce some general notions about the heuristic approach to AS, which we then use to derive benchmark predictions for two NP models.In Sec. 3 we investigate the impact of abandoning one by one the simplifying approximations listed above and we derive analytical and numerical estimates of the resulting uncertainties.We present our conclusions in Sec. 4. We give the explicit form of the 1-and 2-loop RGEs used in this work in Appendix A.

Reference model predictions at 1 loop
In order to quantify the theoretical uncertainties originating from the approximations listed in Sec. 1, one needs a set of "standard candle" predictions to observe how they are modified when those approximations are dropped.In this section we introduce two popular renormalizable NP scenarios as reference models and we derive, in the approximations of Sec. 1, the AS-driven predictions for the couplings of their Lagrangian.In Sec. 3, we will drop the approximations one by one and quantify the relative change in the predicted values.

General notions
We consider typical extensions of the SM characterized by a set of new abelian gauge couplings and new Yukawa couplings.We compute at 1 loop the beta functions of the full, extended gauge-Yukawa system, and we correct the RGEs at µ ≥ M Pl linearly in the couplings, like in Eqs. ( 1) and (2).
Besides the fundamental assumptions of universality and linearity in the matter couplings, the gravity parameters f g and f y must satisfy some loose requirements to be in agreement with the current knowledge of explicit FRG calculations.In particular, it has long been known that f g is likely not negative, irrespective of the chosen RG scheme [28].More specifically, imposing f g > 0 will enforce asymptotic freedom in the non-abelian gauge sector of the SM.Conversely, no essential constraints apply to the leading-order gravitational term f y .The gravitational contribution to the Yukawa coupling was investigated in a set of simplified models [16,24,25,29], but no general results and definite conclusions regarding the size and sign of f y are available.
A fixed point of the system is defined by any set {g * i , y * j }, generically denoted with an asterisk, for which the beta functions develop a zero: The RGEs of couplings {α k } ≡ {g i , y j } are subsequently linearized around the fixed point to derive the stability matrix M ij , which is defined as Eigenvalues of the stability matrix define the opposite of critical exponents θ i , which characterize the power-law evolution of the couplings in the vicinity of the fixed point.If θ i is positive, the corresponding eigendirection is dubbed as relevant and UV-attractive.All RG trajectories along this direction will asymptotically reach the fixed point and, as a consequence, a deviation of a relevant coupling from the fixed point introduces a free parameter in the theory (this freedom can be used to adjust the coupling at the high scale so that it matches an eventual measurement at the low scale).If θ i is negative, the corresponding eigendirection is dubbed as irrelevant and UV-repulsive.As we have stated already, there exists in this case only one trajectory that the coupling's flow can follow in its run to the low scale, thus potentially providing a clear prediction for its value at the experimentally accessible scale.Finally, θ i = 0 corresponds to a marginal eigendirection.The RG flow along this direction is logarithmically slow and one ought to go beyond the linear approximation provided by the stability matrix to decide whether a fixed point is UV-attractive or UV-repulsive.While the approximations listed in Sec. 1 can be abandoned to evaluate the precision of the obtained predictions, the whole heuristic approach relies on a few assumptions which cannot be removed without compromising critically the predictive power of the fixed point analysis: 1.The trans-Planckian UV completion (whether it be gravity of else) should be responsible for the disappearance of the Landau pole in the hypercharge gauge coupling of the SM and for the appearance of an irrelevant fixed point in its place 2. The trans-Planckian UV completion should be responsible for the appearance of an irrelevant fixed point in at least one of the SM Yukawa couplings 3. The RGEs of the NP model completed in the trans-Planckian UV should not be altered by large couplings below the Planck scale.Additional light or heavy states can actually exist, but they ought to be characterized by feeble interactions with the visible sectorsee, e.g., Appendix D in Ref. [55].
Point 1 and 3 allow one to extract the value of f g uniquely, by connecting the hypercharge gauge coupling flow from the Planck scale down to the electroweak symmetry breaking (EWSB) scale.
Point 1, together with the fact that f g > 0, implies that the non-abelian gauge couplings of the SM, g 3 and g 2 , remain relevant and asymptotically free.Point 2 implies that one can connect an irrelevant SM direction from the Planck to the EWSB scale and extract uniquely the trans-Planckian value of f y .All the NP Yukawa couplings for which we seek a prediction should correspond to irrelevant directions, while the SM Yukawa couplings that are not employed for the determination of f y , as well as other possible couplings (point 3), will vanish in the deep UV and correspond to relevant directions in the coupling space.We emphasize that the predictions from the trans-Planckian fixed point analysis are extracted for specific NP models characterized by specific symmetries and particle content, whose RGEs should not be altered significantly below the Planck scale (point 3).As the obtained predictions are typically rather precise, this assumption, of quantum gravity being the only allowed strongly interacting UV completion to the chosen model, can be experimentally falsified.Any other strongly (or weakly) interactive UV completion defined below the Planck scale would have to be treated as a different model altogether, giving rise to different predictions and experimental signatures.
We introduce in the rest of this section two models: the B − L model, for which we seek to predict the dark abelian gauge coupling, the kinetic mixing, and the NP Yukawa couplings associated with the right-handed neutrino sector; and a leptoquark (LQ) model, for which we seek to predict the NP Yukawa interaction of the color-charged leptoquark with a SM quark and a lepton.For each of these models we derive the RGEs using PyR@TE 3 [62,63]. 3he explicit forms of the beta functions for the models considered in this study are given in Appendix A.

Gauged B − L
We derive predictions for the NP sector of the gauged B−L model [64,65].The SM symmetry is extended by an abelian gauge group U(1) B−L , with gauge coupling g B−L .The particle content of the SM is extended by one or more right-handed neutrinos and a complex scalar field S, whose vacuum expectation value (vev) spontaneously breaks U(1) B−L .The abelian charges of the SM and NP fields can be found, e.g., in Refs.[64,65].
The Yukawa part of the SM Lagrangian is extended by terms where we have used two component spinor notation and ε = iσ 2 .Spinor and SU (2) indices are understood to be contracted trivially, following matrix multiplication rules.H and L are the SM Higgs and lepton doublets, respectively, and N is the SM-singlet right-handed neutrino.Y ν and Y N are typically 3-by-3 matrices in flavor space.Here, for simplicity, we will only focus on a single entry of each matrix: y ν ≡ (Y ν ) 33 and y N ≡ (Y N ) 33 .When the scalar S acquires a vev, the right-handed neutrino develops a Majorana mass term.
The abelian gauge part of the Lagrangian takes the form where we indicate with Bµ and Xµ the gauge bosons of U(1) Y and U(1) B−L , respectively, and B µν and X µν are the corresponding field strength tensors.SM fermions f transform under both symmetry factors with charges Q Y , Q B−L so that kinetic mixing ϵ is generated between the two abelian groups.
It is convenient to work in a basis in which the gauge fields are canonically normalized.This can be achieved by a rotation [66,67] which parameterizes the gauge interaction vertices of Lagrangian (5) in terms of a "visible" gauge boson V µ and a "dark" gauge boson D µ : The elements g Y , g d , and g ϵ are related to the original couplings as The value of the gravitational parameter f g is determined by equating the 1-loop RG flow of the hypercharge gauge coupling onto the low-scale MS value in the SM, in agreement with assumption 1 in Sec.2.1.We take g SM,MS Y (M t = 173.1 GeV) = 0.36 [68].We extract f y by equating the 1-loop flow of the top Yukawa coupling to its MS value, y SM,MS t (M t ) = 0.95 [69], in agreement with assumption 2 in Sec.2.1.We obtain f g (1 loop) = 0.0097 , f y (1 loop) = 0.0020 , (9) which lead to the following predictions at the irrelevant fixed point: The subsequent low-scale predictions for the dark sector coupling and kinetic mixing, at the scale M t , are For practicality reasons, we choose to present the low-scale predictions at the EWSB scale as reference.In a realistic neutrino-mass model one should decouple the NP at the scale of the vev of the scalar field S. Note, however, that the dimensionful parameters of the Lagrangian are canonically relevant and as such they do not emerge as predictions of the fixed-point analysis.
In Fig. 1(a) we show for illustration the RG flow of the three gauge couplings from their UV fixed point, cf.Eq. ( 10), down to the EWSB scale, cf.Eq. (12).Gravity decouples at the scale M Pl = 10 19 GeV, marked as a vertical gray line.The RG flow of the three Yukawa couplings is presented in Fig. 1(b).Note the deviation of y t from its UV fixed-point value, due to the presence of relevant directions associated with the gauge couplings g 3 and g 2 .

Leptoquark S 3
The SM is extended in this case by a complex scalar that carries both lepton and baryon number [70,71].Its SU(3) c ×SU(2) L ×U(1) Y quantum numbers are Following the notation of Ref. [52], we introduce the Yukawa-type interaction of S 3 with SM fermions as follows where, again, spinor and SU(2) indices are contracted trivially, Q and L are the SM quark and lepton doublets, respectively, and Y LQ is a 3-by-3 matrix in flavor space.For simplicity we shall focus only on the 3rd generation, denoting y LQ ≡ (Y LQ ) 33 .
As before, the gravitational parameters f g and f y are determined by equating the 1-loop RG flow of the hypercharge gauge coupling and the top Yukawa coupling onto their low-scale MS values.The resulting gravitational parameters read leading to the following fixed-point values for the irrelevant couplings: Finally, the low-scale prediction for the LQ coupling y LQ reads

Estimation of the uncertainties
The predictions from AS for a specific NP model find perhaps their greatest usefulness when deriving phenomenological bounds on the relevant parameters of the model in question (e.g. its masses) [50,54,60] or, alternatively, when confronting the model with an experimental anomaly [52,53,55].In particle physics phenomenology it is natural to interpret a new experimental bound (or, in alternative, a deviation from the SM emerging in the precision measurement of a certain process) as an estimate of the NP contribution to the Wilson coefficient of a higher-dimensional operator in the effective field theory (EFT).While the size of a Wilson coefficient generally points to a UV scale of interest, this corresponds to the actual mass of the particles of the UV completion only if the latter generate the Wilson coefficient at the tree level with couplings of order one.In practice the picture is more obscure, due to our complete ignorance of the actual type and strength of UV interactions.More specifically, the generic EFT matching of a Wilson coefficient C NP to a UV Lagrangian with couplings c i of unknown strength depends possibly on several factors: where Λ sets the scale of the UV physics that is integrated out and the last factor, containing the Higgs vev v H , is meant to indicate effects that depend explicitly on the breaking of electroweak symmetry, e.g., the generation of a neutrino mass in see-saw models, or the presence of a chiral enhancement in dipole operators.As was discussed in Sec. 2, the assumption of AS can lead to predictions for the dimensionless parameters of the Lagrangian that are not constrained otherwise.Thus, once the l.h.s. of Eq. ( 19) is determined following a measurement at the low scale, m NP remains the only unknown variable on the r.h.s., which can be extracted straightforwardly.Following this procedure, under the three assumptions of Sec.2.1, in Ref. [52] a fairly precise prediction for the mass of the S 3 LQ was obtained from a global fit of the Wilson coefficients of the Weak EFT; in Ref. [53] predictions for the mass of two vector-like fermions were obtained from the measurement of leptonic dipole operators; and in Ref. [54] a prediction for the Majorana mass of sterile neutrinos was obtained from the see-saw mechanism.
In this context, pinning down the theoretical uncertainties that mar the predicted couplings c i , which are then fed into Eq.( 19) to obtain m NP , becomes of crucial importance.Ideally, the error should remain significantly below the experimental uncertainty on the determination of the Wilson coefficient C NP .To get a quantitative sense of the precision we expect from the asymptotically safe prediction, we point out that observing, for example, an hypothetical deviation from the SM at the p σ level in some experiment, would induce a relative uncertainty on a related Wilson coefficient C NP around its central value, at the level of Thus, in order for this approach to make sense, we expect the theoretical uncertainty on the coupling prediction to remain below the few percent level.

Impact of higher-order corrections to the beta functions
In this subsection, we describe how the predictions we obtained in Secs.2.2 and 2.3 are modified when we drop the first approximation listed in Sec. 1.In other words, we consider the DREG beta functions in the matter sector at perturbation orders higher than 1.Note, however, that in order to ensure maximal predictivity assumptions 1, 2, and 3 in Sec.2.1 must hold throughout this analysis.4 Comment on threshold corrections Since we impose initial conditions on the solutions of 2-loop RGEs for the SM couplings, consistency requires that we include threshold corrections at the 1-loop order.Such corrections would also reduce the uncertainty associated with the unphysical matching scale.
In the case of an unbroken gauge symmetry, a generic expression for threshold correction exists.For example, in the S 3 LQ model the matching condition for g 3 yields (see for example [72]) Assuming, as we do in this work, that the matching scale is set to µ = M t , this correction is small even for leptoquarks with masses of a few TeV.If, for example, m TeV, the relative shift due to Eq. ( 21) is smaller than 1%.In the case of a broken gauge symmetry the expression of threshold corrections is more involved, but the final result is similar to Eq. ( 21) in its analytical form, being a sum of terms involving logs of the NP particle masses.Therefore, as long as the NP scale is not much different from the matching scale, threshold corrections remain small.Matching conditions for the SM Yukawa couplings are obtained by matching the SM fermion masses, which are measured directly.For example, for the top Yukawa coupling one obtains [72] where M W is the measured mass of the W boson and δ denotes the finite part of the full self energy.To quantify the impact of threshold corrections on y t we have created an S 3 LQ model in FlexibleSUSY [73][74][75][76][77][78][79].For a parameter point with m 2 S 3 = 1 TeV, λ S 3 = λ HS 3 = 0.1 (see Eq. ( 64) in Appendix A.2 for the definition of those parameters) and y LQ = 1 we obtain y NP t = 0.92, as opposed to 0.95 which is the value we use in this work.The correct matching of the NP theory to the SM is required to make precise statements about the values of its parameters.This however depends on the mass spectrum of the model and possibly on mixing matrices and cannot be done without discussing also the dimensionful parameters of the matter theory, which are canonically relevant at the fixed point and therefore cannot be predicted.Since none of the statements made in this work depend on the treatment of threshold corrections, we do not lose generality by neglecting them.In the numerical analysis we will therefore equate the NP equivalents of the SM parameters to their SM values (as was done in Sec. 2) also at the 2-loop level.

Gauge couplings
Given assumptions 1 and 3 in Sec.2.1, one can estimate the numerical value of the gravitational correction f g at the fixed point.Once this is extracted, the property of universality in gravitational interactions can be invoked to find the irrelevant fixed points of other gauge couplings.As was discussed above, the non-abelian gauge couplings remain relevant, so that we can limit our analysis to the generic system of abelian couplings.In agreement with the B − L example, let us consider U(1) Y ×U(1) B−L .It is straightforward to extend our conclusions to any pair of abelian symmetries.As can be seen in Appendix A, the pertinent RGEs take the following parametric form, common to all models with abelian mixing: Equations ( 23)-( 25 expressed in terms of loop coefficients α The c l parameters indicate the gauge and Yukawa couplings of the theory (at perturbation order n = 2), or the gauge, Yukawa, and (quartic) 1/2 couplings (at order n > 2).We require that Eqs. ( 23)-( 25) develop a zero above the Planck scale.
Assuming for the moment that the fixed point for the gauge couplings is developed sharply at a M Pl = 10 19 GeV, one can express the value of f g in terms of the (known) U(1) Y trans-Planckian fixed point g * Y , with an accuracy that increases at each successive order, where the asterisk refers to all couplings being set at their UV fixed-point value.The predicted ratios of the gauge couplings do not depend explicitly on the value of f g .One can define where we have adopted a simplified notation, bi In order to obtain some quantitative estimates, let us retain for simplicity only the 2-loop corrections, Π (i) 2 , and quantify the uncertainty on the ratios r * g,i(=d,ϵ) by calculating  26) are by far those corresponding to the product of two gauge couplings.Retaining, thus, only the abelian gauge-coupling contributions in the 2-loop RGEs of Appendix A, and assuming in the first approximation that all gauge couplings acquire the same fixed-point value, one can factor out the sums of coefficients α (i) lk and compute them separately keeping track of the relative coupling signs.The percent uncertainty is presented in Fig. 2. It is very similar in the two cases i = d, ϵ.It is shown as a function of the generic product of abelian gauge couplings |g * l g * k | at the fixed point.We note that the error remains at the percent level, far below the uncertainty typically associated with Wilson coefficients of higher-dimensional operators of the EFT, see Eq. (20).Similar levels of uncertainty apply to the UV predictions of gauge couplings in other models with abelian gauge mixing, e.g., those analyzed in Ref. [55].

Yukawa couplings
The discussion can be repeated with small modifications for the Yukawa couplings.Let us parameterize the full set of Yukawa-coupling RGEs of a generic SM+NP theory in the following way: where y r(j)=1,2,... label the set of Yukawa couplings in the Lagrangian and g l(k) label the gauge couplings of the theory.In Eq. (32) we have defined the generic multiplicative n-loop contri-bution as and the generic additive n-loop piece as in terms of n-loop coefficients α , and of all gauge, Yukawa, and (quartic) 1/2 couplings c l i entering the Yukawa RGEs at order n ≥ 2 .
Neglecting for the moment the additive part of the beta functions, let us assume we solve the system (32) when all the beta functions are equal to zero.Repeating the steps that led to the 1-loop results in Sec. 2, we express the unknown gravitational contribution f y in terms of one known "reference" fixed-point Yukawa coupling y * 1 , which is in general the top quark's.5Let us retain, for simplicity, only the 2-loop correction.One gets where F 0 , G 0 , and H 0 are rational functions of the a (r) j coefficients, with the latter two being linear in lk a 2 , respectively.Their analytic form is model-dependent and is not needed for the discussion.Note that the reference fixed-point value y * 1 (2 loops) is obtained by equating the 2-loop RGE flow to the MS values of the couplings at the EWSB scale.As such, it can be significantly shifted with respect to y * 1 (1 loop), mostly due to the contributions of relevant parameters like g 3 and g 2 to the 2-loops RGEs.
Let us thus define the shift at the fixed point due to relevant parameters in the flow, and derive the parametric expression for the fixed point of a second Yukawa coupling of the irrelevant type.This is our prediction, which we indicate as y * 2 .Its value must be a function of f y .By using Eq. ( 35) we trade f y for y * 1 and plug into the fixed-point solution for y * 2 : where the F 1 , G 1 , are some other rational functions of coefficients of order one, and the H 1 function carries the 2-loop terms.Using Eq. ( 33), one can see that the higher-loop correction H 1 only depends, up to order-one coefficients, on sums of fixed-point products 1 and the G 1 function.As a consequence, the uncertainty due to δy * 2 1 and H 1 on the final result would grow.Conversely, obtaining y * 2 ≈ y * 1 , g * k at 1 loop would imply that such prediction can be trusted to a very good approximation.
For a quantitative example, let us choose a fixed point with one irrelevant gauge coupling (l = k = 1) and two Yukawa couplings (r, j = 1, 2).Equation (37) becomes (38) This is the case, for example, of the S 3 LQ.With respect to the notation of Sec.2.3, here g 1 ≡ g Y , y 1 ≡ y t , and y 2 ≡ y LQ .All other couplings are chosen relevant and zero at the fixed point.The 1-loop coefficients in this model can be found in Appendix A and read a We note that in this example the uncertainty is not negligible.This was expected, since at  1.The discussion can be extended, finally, to cases in which an additive contribution to the Yukawa coupling RGEs -the fourth and fifth addends in Eq. ( 32) -is present.It is not difficult to realize that in those cases the Yukawa system only admits a set of real interactive solutions y * r ̸ = 0 of the same order of magnitude as the irrelevant gauge-coupling fixed point.Since there is no fine tuning involved, such solutions are not very sensitive to the n-loop corrections.

Numerical examples
We provide quantitative estimates of the higher-loop effects discussed in this section in the two NP models introduced in Sec. 2.
In Table 1 we present the 2-loop determinations of the gravity parameters f g and f y , the fixed-point values of the model couplings, as well as their percent uncertainties at the fixed point and at the low-energy scale, µ = M t .The latter are derived under the assumption that the low-scale value of the reference SM coupling (g  1 confirm our discussion.The impact of higher-order corrections on the fixed-point value of a coupling is almost negligible for the gauge couplings, while in the Yukawa sector it depends on the relative size of the reference and the predicted couplings.In the B − L model, the fixed-point values of all three couplings are comparable in size.Hence, the impact of 2-loop corrections is small and the uncertainty on the NP couplings is of the order of a few percent.Conversely, in the S 3 LQ model, y LQ is smaller by a factor of 2 than y t .The impact of the 2-loop corrections is then more pronounced, and the resulting fixed-point uncertainty increases. Closer inspection of the last column of Table 1 reveals that the low-scale uncertainties of the NP Yukawa couplings are smaller than those at the fixed-point.This is not a coincidence and it is intrinsically related to modifications of the RG flow due to higher-order corrections.To get a better understanding of the underlying mechanism, we show in Fig. 3(a) the 1-loop (dashed) and 2-loop (solid) running of y t (red), y ν (blue) and y N (green) in the B − L model.The analogous flow of the top (red) and LQ (blue) Yukawa couplings of the S 3 model is presented in Fig. 3(b).Due to the negative contribution of the gauge coupling g 3 to the beta functions of the Yukawa couplings of particles carrying the color charge, the running value of y t at 2 loops is generically smaller than at 1 loop.As was discussed below Eq. ( 37), it follows that the fixed-point values of the NP Yukawa couplings also tend to be shifted downwards.The effect is stronger for the couplings that depend on y t already at 1 loop (y ν and y LQ ).
On the other hand, since the same reference low-scale value y t (M t ) is used in the fixedpoint analysis at any loop order, the RG trajectories of the NP Yukawa couplings have the tendency to focus towards their 1-loop value in the infrared, which results in a reduction of the uncertainty with respect to the prediction at the fixed point.A generic rule of thumb thus applies: the uncertainty in the determination of the fixed-point value of a NP Yukawa coupling provides an upper bound on the uncertainty of the same prediction at the low scale.

Dependence on the position of the Planck scale
In this subsection, we drop the second of the approximations listed in Sec. 1, i.e., we consider what happens if gravity decouples from the matter RGEs sharply at a scale that differs from 10 19 GeV by a few orders of magnitude.

Gauge couplings
When assessing the impact of the Planck-scale position on the predictions for gauge couplings one should note that an uncertainty on the position of the Planck scale is effectively equivalent to an uncertainty on the fixed-point value of the hypercharge gauge coupling g * Y , hence on f g .On the other hand, since the f g dependence cancels out from Eqs. ( 28) and ( 29), moving the Planck scale back and forth does not affect the predicted ratios r * g,i at the 1-loop order.This feature is not preserved at higher orders in the perturbative expansion.However, the impact of the Planck scale position remains negligibly small, at about 0.01%.

Yukawa couplings
Neglecting for this discussion the 2-loop contribution, the percent uncertainty on the predicted Yukawa-coupling ratios propagates in Eq. ( 37) as It is straightforward to re-frame the two cases described in Sec.3.1 in light of Eq. ( 41).40).
If the predicted Yukawa fixed point is of comparable size to the irrelevant gauge couplings, r * y,2 ≈ r * g,l(k) , the uncertainty is dominated by the shifting of the gauge coupling fixed point: Conversely, if the predicted coupling is small, r * y,2 ≪ r * g,l(k) , the relative uncertainty can grow substantially and we lose control of the prediction, As a quantitative example, let us consider again the simple case of two Yukawa couplings, y * 1 , y * 2 , and one gauge coupling g * 1 (hypercharge), all irrelevant at the fixed point.Equation ( 41) becomes in this case δr * y,2 This is the case, e.g., of the S 3 LQ model, whose 1-loop coefficients were given below Eq. (38).One can show numerically (1-loop RGEs will suffice) that an eventual shift in the value of the Planck scale will induce a variation in the fixed point values g * 1 , y * 1 , and, consequently, r * g,1 ≡ g * 1 /y * 1 .For example, let us move the Planck scale to M Pl = 10 16 GeV in one case, and to M Pl = 10 20 GeV in another.This results in the changing of g * 1 = g * Y from the value in Eq. ( 17) to 0.45 and 0.49, respectively.One finds δr * g,1 /r * g,1 (10 16 ) ≈ −10% and δr * g,1 /r * g,1 (10 20 ) ≈ 4% .We can thus use Eq.(44)

Numerical examples
In Table 2 we summarize the findings of this subsection for our benchmark models.For two different values of the Planck scale, M Pl = 10 20 GeV and M Pl = 10 16 GeV, we show determinations of the gravity parameters f g and f y , the fixed-point values of the model couplings, their percent uncertainties at the fixed point and the uncertainties at the scale µ = M t .One can immediately see that all the fixed-point values of the gauge couplings are rescaled by the same amount, so the low-scale predictions are not altered.
The RG-invariance of the gauge-coupling ratios r * g,i is made transparent in Fig. 4(a), where we show the difference in the flow of the gauge couplings of the B − L model for different choices of the Planck scale.The shift in the flow of the Yukawa sector of the B − L models is presented in Fig. 4(b), whereas the difference in the flow of the S 3 LQ model is shown in Fig. 4(c).
As was discussed in Sec.3.1, we observe focusing of the RG trajectories of the Yukawa couplings below the Planck scale.This indicates that the uncertainty calculated at the fixed point according to Eq. ( 41) is the maximal uncertainty that can be observed in a given NP model.

Scale-dependence of the gravitational corrections
Finally, in this subsection we are going to drop the third of the approximations listed in Sec. 1, i.e., we include in the analysis the scale dependence of the gravitational parameters.The actual functional form of f g and f y is subject to technical assumptions within the FRG framework.For example, in the Einstein-Hilbert truncation and in the Landau-gauge limit the relevant formulae can be found in Ref. [38].In the following analysis we will avoid sticking to any particular form of f g (t) and f y (t) to be able to draw conclusions that remain as generic as possible.
Gauge couplings Consider Eq. ( 1), applied e.g. to the three (irrelevant) gauge couplings of the B − L model, i = Y, d, ϵ.Due to its universality, the explicit gravitational contribution to the matter beta function f g entirely factors out of the running ratios: where the matter beta functions β matter i are given in parametric form in Eqs. ( 23)-( 25) and we have formally defined slope functions F d and F ϵ which do not depend on f g .We are going to show that the 1-loop ratios g d /g Y and g ϵ /g Y are exact invariants of the RG flow, whereas, at order n ≥ 2, RG invariance is respected up to a very good approximation.
To understand the underlying mechanism, let us apply the definition of total derivative to Eq. ( 45) and focus on a sequence of infinitesimal scale intervals, ... t 2 < t 1 < t 0 , with the system lying at the fixed point at t 0 .Moving backwards in t, one gets where r * g,d was defined in Eq. ( 28) for REGs at arbitrary loop order.The same steps can be repeated for Eq.(46).One can check, by directly imposing boundary condition (28) into Eq.( 45) at t 0 , that F d vanishes at the fixed point: F d (g * Y , ...) = 0. Equivalently, one can show that F ϵ (g * Y , ...) = 0. Thus, g d (t 1 )/g Y (t 1 ) = r * g,d and g ϵ (t 1 )/g Y (t 1 ) = r * g,ϵ .For the next time interval one can plug Eqs. ( 28), (29) into Eqs.( 23)-( 25) and express the slope functions in terms of fixed-point ratios: where, in agreement with Eq. ( 30), we have defined bi (t) away from the fixed point.At order n = 1, bi (t) = b i .Boundary conditions ( 28) and ( 29) ensure that F (n=1) d (t 1 ) = 0 and F (n=1) ϵ (t 1 ) = 0, independently of the actual values of the running gauge couplings.As an immediate consequence, g d (t)/g Y (t) = r * g,d and g ϵ (t)/g Y (t) = r * g,ϵ along the entire RG flow.Given that the ratios r * g,d and r * g,ϵ are, at 1 loop, uniquely determined by the gauge quantum numbers, we can conclude that the predictions of AS for the gauge couplings do not depend on the particular functional form of the gravitational contribution f g (t).
On the other hand, the scale-invariance of the gauge coupling ratios is washed out by higherorder loop contributions, which introduce a t-dependent deviation from the 1-loop prediction due to appearance in Eq. ( 51) of couplings (relevant gauge couplings, Yukawa and quartic couplings) that lie outside of the system ( 23)- (25).While these deviations from scale-invariance could potentially build up over scales spanning many orders of magnitude, they remain in practice quite limited.
For illustration, we show in Fig. 5 the RG flow of the three gauge couplings of the B − L model: g Y (red), g d (blue), and −g ϵ (green).Dotted lines correspond to the benchmark scenario with f g and f y constant above the Planck scale (M Pl = 10 19 GeV).Darker solid lines indicate two random parametrizations of the functional dependence f g (t), where the gravity parameter is allow to vary by a factor 10 in the range between 10 16 GeV and 10 20 GeV.One can see that, in spite of different fixed-point values in each case, the RG invariance of the coupling ratios leads to unchanged low-scale predictions.Finally, the dashed lines correspond to the f g (t) parametrization based on the FRG results of Ref. [38].Note that in that framework it is reasonable to expect only very moderate changes in the gravity parameters.

Yukawa couplings
In the case of the Yukawa coupling RGEs given in Eq. ( 2) one can consider the ratio of two irrelevant Yukawa couplings y 1 and y 2 : where the parametric form of the matter beta functions is given in Eq. ( 32) and the slope function G does not depend explicitly on f y .Unlike the gauge sector, which features scaleinvariant ratios at 1 loop, for the Yukawa couplings the t-dependence of f g and f y can impact the running of the y 2 /y 1 ratio even at 1 loop.  .Dashed lines correspond to the FRG results of Ref. [38].
In order to see this, let us neglect for simplicity the additive contribution to the Yukawa beta functions in Eq. (32).In analogy to the gauge coupling case, one can then express where r * y,2 was defined before Eq.(40).Note that the non-abelian gauge couplings and some of the Yukawa couplings contributing to the sums within square brackets in Eq. ( 55) correspond to relevant directions of the parameter space.Thus, even at 1 loop, if any of the couplings deviates from its fixed point (due to the change of the gravity parameters or the presence of relevant directions in the coupling space) the contribution from the first line of Eq. ( 55) does not vanish exactly at t 1 and the ratio y 2 (t)/y 1 (t) starts to flow.
This effect is illustrated in Fig. 6(a), where we show the actual size of the function |G(t)| for different parametrizations of f g (t) and f y (t) in the S 3 LQ model.Gray dashed line indicates our benchmark scenario with both gravitational parameters constant above the Planck scale (M Pl = 10 19 GeV).Red solid line corresponds to a situation in which f g is kept constant, but f y is allowed to vary by a factor 10 in the range between 10 16 GeV and 10 20 GeV.Once can see that there is no difference with the previous case, confirming the fact that the f y dependence factors out from the RG running of the Yukawa-coupling ratio.Finally, in green and blue we show two arbitrary parametrizations of f g (t), f y (t).Despite broad fluctuations, |G(t)| remains small in size.We can conclude that the flow of the ratio y 2 (t)/y 1 (t) remains fairly stable throughout.
As a matter of fact, the dominant source of uncertainty for the prediction of a Yukawa coupling is not given by the changing of y 2 (t)/y 1 (t) along the RG flow, but rather by the fact that the fixed-point ratio r * y,2 itself becomes unknown once we factor in the t-dependence of f g and f y .As was shown in several papers [51][52][53], at 1 loop in the matter RGEs the fixed-point ratio of two Yukawa couplings takes the parametric form Numerical coefficients A, A ′ , B, B ′ are combinations of the 1-loop coefficients of the matter beta functions.They can be finite or zero.Referring, again, to the simple case of two irrelevant Yukawa couplings, y 1 , y 2 , and one gauge coupling, g 1 , one can write where a (r) lk are coefficients of the Yukawa-coupling beta functions, cf.Eq. (32), and b 1 is the 1-loop coefficient of the g 1 beta function.
In the limiting case of one of the two gravitational parameters being strongly dominant with respect to the other, the ratio becomes known, albeit possibly extreme (even zero or infinity).For example, in the S 3 LQ model one gets which are not far from the prediction obtained from Eq. ( 17), r * y,2 ≈ 0.48.In general, however, cancellations in the numerator/denominator of Eq. ( 56) might take place, so that there is no real handle on the prediction.This is illustrated in Fig. 6(b), where we show the RG flow of the top Yukawa (red), LQ Yukawa (blue) and hypercharge (green) gauge coupling of the S 3 model with different functional forms of f g (t) and f y (t).The color code is the same as in Fig. 5(a).While it was shown above that even extreme fluctuations do not modify the low-scale predictions for the gauge couplings, we can see in Fig. 6(b) that the trans-Planckian behavior of the running y LQ (blue) and y t (red) couplings is drastically different in the cases in solid vs. those in dotted/dashed.
An assessment from first principles of the uncertainties on the action was attempted in Ref. [45].Several choices of the truncation, and different matter content, were shown to induce significant shifts in the fixed-point values of the cosmological and Newton constant.On the other hand, when those shifts are plugged into typical FRG computations of f g and f y -e.g., those in Ref. [38] -the actual percent uncertainty on f g (t) and f y (t) rarely exceeds the O(1) level along the entire trans-Planckian flow.In Fig. 6(b), the parametrization based on the explicit FRG results of Ref. [38] is depicted in dashed lines.It introduces a small uncertainty on the prediction of y LQ , of around 5%.
Incidentally, in the sub-Planckian regime of Fig. 6(b) we observe focusing due to the presence of a large coupling g 3 , in agreement with the discussion at the very end of Sec.3.1.The resulting low-scale uncertainty on y LQ reduces to 4% for the FRG-inspired functional dependence of the gravitational parameters.For the case in solid, which corresponds to f g (t) ≪ |f y (t)|, the low-scale uncertainty on y LQ reads ∼ 19%.

Relevant abelian gauge couplings
The systematic analysis featured in the three subsections presented above relies on the implicit choice of a maximally predictive fixed point for the particle physics plus gravity system.This was guaranteed by enforcing assumptions 1, 2, and 3 in Sec.2.1.
As was mentioned at the beginning of Sec.3.1, however, often one may be interested in obtaining predictions for irrelevant Yukawa couplings in a scenario where all the gauge couplings, including the abelian ones, correspond to relevant directions around a Gaussian fixed point -g * Y = 0, etc.All gauge couplings become, in other words, asymptotically free like the nonabelian ones.Such situation emerges quite naturally if f g is large enough (larger, for example, than the values given in Eq. ( 9) and Eq. ( 16) for the B − L and S 3 LQ models, respectively) and it corresponds, effectively, to relaxing assumption 1 in Sec.2.1.
We discuss in this subsection what uncertainties we expect in this case for the predicted Yukawa couplings.When g * Y = 0, the irrelevant fixed points of the Yukawa-coupling system do not depend on the value of f g .Consequently, the dependence on f y factors out of the ratio r * y,2 , which then becomes fully determined by the 1-loop coefficients of the beta functions, see Eq. (56).That does not mean, however, that predictivity in the Yukawa sector is fully restored.In fact, the presence of relevant parameters typically affects the RGEs of different Yukawa couplings unequally.
Let us consider two arbitrary values f g , f ′ g > f AS g , where f AS g would be the value making the fixed point in g * Y irrelevant -for example, the value given in Eq. ( 9) or the one in Eq. ( 16).Unlike f AS g , f g and f ′ g cannot be predicted from low-energy measurements and are thus unknown.Given the beta-function system of Eq. ( 32) (at 1-loop order) and following the steps that led to Eq. ( 37), we can derive the shift on the Yukawa-coupling fixed-point prediction computed with f ′ g with respect to f g : where Note that, unlike in Eq. ( 36), the fixed-point values cannot be here controlled without explicit knowledge of the quantum gravity parameters.
To get a quantitative sense of the shift observed in predictions, let us focus on the S 3 LQ model and consider two arbitrary choices leading to relevant fixed points in the gauge-coupling system: f g = 0.0107 and f ′ g = 0.0212.Recalling that y * 2 ≡ y * LQ , the relative uncertainty reads approximately of the size of the relative uncertainty on the parameter f g itself.We observe similar numerical values for the predictions of the B − L model.On the other hand, one should keep in mind that, notwithstanding what the actual numerical value of f g > f AS g may be, it should always lead the running gauge couplings g 3 , g 2 , and g Y to match their low-scale determination (which is independent of f g ).As such, we expect the predictions in the Yukawa sector to be subject to significant focusing in their flow from the fixed point down to the EWSB scale.In fact, in the case of the S 3 LQ we observe a significant reduction of the uncertainty: As was the case in Sec.3.1, we can state a general rule of thumb, that the uncertainty at the fixed point is always larger that the actual one at the low-energy scale.

Conclusions
In this paper, we have investigated in detail the issue of the theoretical uncertainties associated with predictions for irrelevant Lagrangian couplings emerging from trans-Planckian AS.
We have chosen a complementary approach with respect to Ref. [45]: while that article evaluated the uncertainties of the quantum-gravity sector in the context of calculations within the FRG, in this work we have decided to bypass strictly gravitational aspects following a popular phenomenological approach which relies on a parametric description of quantum gravity with universal coefficients that are eventually obtained from low-energy observations.Our focus lies squarely on the uncertainties pertaining to the matter sector, which takes the form of a system of gauge-and Yukawa-coupling RGEs.
We have quantified the impact of relaxing several approximations commonly used in the literature to extract phenomenological predictions.In particular we have considered: the effect of including higher-order corrections in the matter RGEs of the gauge-Yukawa system; the impact of selecting a value different from 10 19 GeV for the (somewhat arbitrary) position of the Planck scale; and the effects of using scale-dependent parametrizations of the gravitational UV completion in the matter RGEs, resulting in a non-instantaneous decoupling of gravity from matter around the Planck scale.
The main findings of our study are the following.First, in the gauge sector the uncertainty induced by relaxing any of the simplifying assumptions never exceeds the 1% level.We can conclude that fixed-point predictions for the irrelevant gauge couplings of the SM and/or NP models are extremely robust, even when they are obtained in a heuristic, simplified approach to AS that is based on some approximations.
A similar conclusion can be drawn for the Yukawa sector of the NP theory, if the predicted Yukawa couplings are of comparable size to the irrelevant gauge couplings.The uncertainties remain at bay, not exceeding ∼ 10% at the fixed point if higher-order corrections are included, or if the Planck scale is moved to, e.g., 10 16 GeV.The situation is additionally helped by focusing of the RG trajectories in the sub-Planckian regime.
Potentially more dangerous uncertainties could stem from considering the non-trivial scale dependence of the gravitational contributions to the matter beta functions, parameterized by functions f g (t) and f y (t), as in this case we lose the ability of determining the actual value of the Yukawa couplings at the fixed point.However, we have argued, based on both an analytical and numerical discussion, that in the range of variability of the gravitation parameters that can be realistically expected in the framework of the FRG, the resulting uncertainty is moderate.
Finally, we have identified one situation in which the AS-based predictions cannot be trusted as any modification of the setup would lead to a drastic change on the predicted value of a NP coupling.This happens if the predicted Yukawa coupling is much smaller than the irrelevant gauge couplings, since in this case it is a result of an accidental, precise cancellation of two quantities of comparable size.This is a manifestation of fine tuning, exactly analogous to some of the cases described in, e.g., Ref. [54], in which the gravitational parameter f y happened to lie unnaturally close to a critical value f crit Z,XY which depends exclusively on the gauge quantum numbers of the matter sector.Since the prediction then finds its origin on a fortuitous cancellation of unrelated quantities it is very subject to the approximations employed for its derivation.

A Renormalization group equations
In this appendix we presents the DREG matter RGEs used in this work.They are listed order-by-order according to the following notation: We do not write explicitly the contributions from the down-type quark and charged-lepton Yukawa matrices which do not affect our numerical analysis.
β (1)  β (1)  We have checked that these couplings do not influence predictions for the gauge and Yukawa sector in a significant manner.As such, they are set to 0 in our numerical analysis.
Note that there was an error in the normalization of Y LQ in the RGEs given in Ref. [52].Below we present the correct expressions.

Figure 1 :
Figure 1: RG flow of the (a) gauge and (b) Yukawa couplings in the gauged B − L model from trans-Planckian energies down to the EWSB scale, for a scenario characterized by the UV fixed points of Eq. (10) and Eq.(11).The gray vertical line indicates the fixed position of the Planck scale, M Pl = 10 19 GeV.
) are given in terms of the one-loop coefficients b Y , b d , and b ϵ , and generic n-loop contributions,

Figure 2 :
Figure 2: Estimated percent uncertainty on the predicted fixed-point ratios of abelian gauge couplings |δr * g,i /r * g,i | in the B − L model due to missing higher-loop contributions.The uncertainty is given as a function of the product of abelian gauge couplings at the fixed point.

2 ≈ − 3 ×
The shift in the top Yukawa fixed point can be computed numerically: δy * 2 1 ≈ −9 × 10 −3 .This should be summed to the explicit 2-loop contribution to the beta functions, Π can be obtained from Appendix A. The latter turns out to be approximately one order of magnitude smaller than δy * 2 1 when the quartic couplings of the scalar potential are neglected.If the quartic couplings are assumed to be of order one, we get instead Π 10 −2 .The first of these determinations leads to the largest uncertainty on the prediction.By plugging the 1-loop coefficients into Eq.(38), together with the 1-loop fixed-point values given in Sec.2.3, one can easily compute the percent deviation,

Figure 3 :
Figure 3: (a) RG flow of the top (red), Dirac neutrino (blue) and sterile neutrino (green) Yukawa couplings of the gauged B − L model at the 1-loop (dashed) and 2-loop (solid) orders.The gray vertical line indicates the fixed position of the Planck scale, M Pl = 10 19 GeV.(b) Same as in (a), but for the top (red) and LQ (blue) Yukawa couplings of model S 3 .

Figure 4 :
Figure 4: (a) 1-loop RG flow of the hypercharge (red), dark gauge (blue), and kinetic mixing (green) couplings of the B − L model with the Planck scale set at 10 19 GeV (dashed) and 10 16 GeV (solid).The gray vertical line indicates the fixed positions of the Planck scale.(b) Same as in (a) for the top Yukawa (red), the Dirac neutrino Yukawa (blue) and sterile neutrino Yukawa coupling (green) of the B − L model.(c) Same as in (a) for the top Yukawa (red) and LQ Yukawa coupling (blue) of the S 3 model.

Figure 5 :
Figure 5: RG flow of the hypercharge (red), dark gauge (blue), and kinetic mixing (green) couplings of the B − L model at 1 loop.Dotted lines correspond to the benchmark scenario with f g constant above the Planck scale.Darker solid lines indicate two arbitrary parametrizations of f g (t) .Dashed lines correspond to the FRG results of Ref. [38].

Figure 6 :
Figure 6: (a) Function |G(t)| in the S 3 LQ model for different parametrizations of f g (t) and f y (t).Gray dashed line indicates the benchmark scenario with both gravitational parameters constant above the Planck scale (M Pl = 10 19 GeV).Red solid line corresponds to f g constant and f y allowed to vary by a factor 10 in the range between 10 16 GeV and 10 20 GeV.Blue and green lines show two arbitrary parametrizations of f g (t) and f y (t).(b) 1-loop RG flow of the top Yukawa (red), LQ Yukawa (blue) and hypercharge gauge (green) couplings of the S 3 model for different parametrizations of f g (t) and f y (t).Dotted lines correspond to the benchmark 1-loop scenario with f g constant above the Planck scale.Dashed lines indicate a parametrization based on the FRG results of Ref. [38].Solid lines show a parametrization resulting in f g (t) ≪ |f y (t)| in the trans-Planckian regime.
It is thus very similar in size for all the predicted Yukawa couplingsy * 2 , y * 3 , etc.The shift uncertainty δy * 2 1 also enters in the same way in all predicted Yukawa couplings.However, the relative uncertainties affecting the predictionsδy * 2 /y * 2 , δy * 3 /y * 3 , etc. -differ for the different Yukawa couplings because they are very sensitive to the actual size of the coupling itself.If one, for example, had obtained at 1 loop y * 2 ≪ y * 1 , g * k , that would have implied a precise cancellation (fine tuning) between F 1 (a

Table 1 :
(39)op determination of the gravity parameters f g and f y , fixed-point values of the reference and the to-be-predicted couplings, percent uncertainty at the fixed point, and percent uncertainty at the low scale for the models introduced in Sec. 2. The uncertainties are defined w.r.t. the 1-loop results of Sec. 2, cf.Eq.(39).Yukawa sector of the B − L model and in the S 3 LQ model) is not varied when performing the fixed point analysis at 2 loops.The results shown in Table (37)he other hand, since Eq.(37)depends explicitly on the fixed-point values of the couplings, changing the position of the Planck scale will alter the prediction for the Yukawa couplings, as it changes the fixed-point values of the abelian gauge couplings, which we indicate collectively with g * k=1,2,... , and of the reference Yukawa coupling y *

Table 2 :
Impact of the position of the Planck scale on the determination of the gravity parameters f g and f y , fixed-point values of the reference and the to-be-predicted couplings, percent uncertainty at the fixed point, and percent uncertainty at the low scale for the models introduced in Sec. 2. The uncertainties are defined w.r.t. the 1-loop results of Sec. 2, cf.Eq. ( In this model the 2-loop Yukawa RGEs depend also on the dimensionless BSM Higgs potential parameters λ S 3 and λ HS 3 , which we define as 2 S 3 Y LQ + λ 2 HS 3 Y LQ +