Scheme dependence of asymptotically free solutions

Recent studies have provided evidence for the existence of new asymptotically free trajectories in non-Abelian particle models without asymptotic symmetry in the high-energy limit. We extend these results to a general SU(NL)×SU(Nc)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm {SU}(N_{\mathrm {L}})\times \mathrm {SU}(N_{\mathrm {c}})$$\end{document} Higgs-Yukawa model that includes the non-Abelian sector of the standard model, finding further confirmation for such scenarios for a wide class of regularizations that account for threshold behavior persisting to highest energies. We construct these asymptotically free trajectories within conventional MS¯\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\overline{\text {MS}}$$\end{document} schemes and systematic weak coupling expansions. The existence of these solutions is argued to be a scheme-independent phenomenon, as demonstrated for mass-dependent schemes based on general momentum-space infrared regularizations. A change of scheme induces a map of the theory’s coupling space onto itself, which in the present case also translates into a reparametrization of the space of asymptotically free solutions.


Introduction
Universality in physics characterizes the fact that long-range effective properties of a system can be largely independent of its microscopic details. In particle physics, where microscopic details, say at the Planck scale, are neither known nor currently experimentally accessible, universality is often quantified in terms of observables which should be independent of the choice of the regularization and the renormalization scheme.
On a technical level, universality can also be visible in properties of renormalization group (RG) functions such as β functions specifying the behavior of couplings under a a e-mail: holger.gies@uni-jena.de b e-mail: rene.sondenheimer@uni-graz.at c e-mail: alessandro.ugolotti@uni-jena.de d e-mail: luca.zambelli@uni-jena.de change of scale. A standard textbook result is the scheme independence of the perturbative one-loop β function coefficient; in a mass-independent scheme, also the two-loop coefficient is universal. These results form the basis of classifying theories according to their weak-coupling behavior with a prominent example being asymptotic freedom (AF) towards high energies for the case of a negative one-loop coefficient [1][2][3][4][5][6][7][8][9][10][11][12][13].
Though being universal, the one-loop β function coefficent does not necessarily provide a reasonable measure for the physical scale dependence of couplings. A simple example is the running of the QED fine-structure constant at, say, nano-electron-Volt scales: here the standard one-loop coefficient still assumes its standard value, whereas the coupling (as, for instance, measured by Thomson scattering) does not run at all, because the electron fluctuations decouple below the electron mass threshold.
The reason for this apparent mismatch is that standard β function definitions make implicit use of the deep Euclidean region (DER), where all physical mass scales or external momenta are assumed to be small with respect to the loop momenta of the fluctuations. By contrast, definitions of RG functions that take mass or momentum thresholds explicitly into account lead to β functions that describe the decoupling adequately. A famous example is given by RG functions defined by the Callan-Symanzik equation [14][15][16].
The price to be paid for including physical threshold phenomena in an RG description is that the corresponding β functions become scheme dependent even at one-loop order. This is natural, as this dependence parametrizes the details of the physical decoupling of massive modes; of course, such a scheme dependence cancels in physical observables such as cross sections.
While threshold phenomena in RG functions are wellknown and controlled by standard procedures [17][18][19][20][21][22], their potential role towards higher energies has been studied very little. Here, the analysis in the DER seems only natural, as highest momentum fluctuations are assumed to always exceed any mass scale. In the case of mass generation through spontaneous symmetry breaking this expectation is summarized as "asymptotic symmetry" [23].
By contrast, new RG trajectories have recently been discovered in non-Abelian gauge theories with various matter content that invalidate the assumption of asymptotic symmetry [24]. Most importantly, these trajectories give rise to new routes to AF and thus ultraviolet (UV) complete scenarios in non-Abelian Higgs models [24,25] as well as gauged Yukawa models [26], with large classes of models remaining to be explored and used for model building. In fact, AF theories still enjoy an unabated interest for the construction of UV complete models in particle physics [27][28][29][30][31][32][33][34][35][36].
Since the occurrence of symmetry-breaking-induced thresholds on all scales is an essential ingredient for the corresponding RG flows, the standard reasoning used for the DER and implying one-loop universality is no longer applicable. This raises naturally the question of scheme dependence: is the existence of these new AF UV completions an universal statement? Can it be verified in a scheme-independent fashion? Answering these questions is a goal of the present work.
For this, we first generalize previous studies to a Yukawa model with an SU(N L ) × SU(N c ) gauge symmetry, covering the non-Abelian part of the Standard Model (SM); also, the previously considered Z 2 -Yukawa-QCD and non-Abelian Higgs models represent limiting cases. In order to make contact with the most widely used MS scheme of standard perturbation theory, we elucidate the construction of AF trajectories on the basis of one-loop β functions obtained from dimensional regularization. The new UV-complete trajectories become visible from these RG functions upon inclusion of a running expectation value and higher dimensional operators, as is familiar from an effective-field theory (EFT) approach.
A functional approach for the full Higgs potential can also be set up within the MS scheme. We present several approaches to analyze the resulting β functional also including its global stability features towards the UV limit. While the MS scheme -though widely used -is a rather particular projection scheme, a more comprehensive analysis can be performed on the basis of general mass-dependent schemes with momentum-space regularization, as featured, e.g., by the functional RG (FRG). Here, we provide further evidence for the existence of these AF trajectories for all admissible regulator functions.
In agreement with earlier findings [24][25][26], the new RG trajectories occur as quasi-fixed points (QFPs) of the β functions in the matter sector. These QFPs are driven by the AF gauge couplings to the non-interacting Gaußian fixed point (FP) towards higher energies. The presence of an AF gauge sector -potentially also beyond the DER -hence forms a crucial ingredient in our construction. The important point, however, is that this feature of AF can fully extend to further sectors of the system which may not seem to be AF in the conventional perturbative analysis restricted to the DER. For future work, an analysis going beyond the DER may also be worthwhile for asymptotically safe particle-physics scenarios [37][38][39][40][41] which have recently attracted substantial attention for concrete model building [42][43][44][45][46][47][48][49][50][51][52].
The paper is organized as follows: In Sect. 2, we first introduce the class of models featuring a local SU(N L ) × SU(N c ) gauge symmetry, highlight several relevant limiting cases, and review the standard perturbative analysis also establishing our notation. In order to transcend the limitations of the DER, Sect. 3 presents a perturbative weak-coupling analysis within an EFT approach, allowing for an inclusion of higher-dimensional operators. Here, the analysis primarily relies on the most widely used MS scheme based on dimensional regularization, elucidating how the new AF trajectories become visible in the most conventional scheme using standard methods. Section 4 is devoted to a functional analysis of the flow of the effective potential, still using the MS scheme. We construct various functional approximations to the QFP trajectories; this includes also a controlled weakcoupling expansion, illustrating how the perturbative EFT emerge in the functional picture. In order to discuss general classes of RG schemes, we set up the FRG equations for the models in Sect. 5, employing a derivative expansion of the effective action and also accounting for threshold effects in the gauge sectors. The scheme-independent existence of the new AF trajectories is then demonstrated in Sect. 6 using a weak-coupling analysis of the FRG equations. As expected, a change of regularization scheme induces a map of the coupling space onto itself, thereby rearranging the space of initial conditions used for specifying AF trajectories. As a new ingredient, this theory space also includes rescaling parameters that distinguish between different AF trajectories.

Asymptotic freedom within perturbative renormalizability
Ultimately aiming at the SM, we base our concrete studies in this work on a toy model which comprises both non-Abelian sectors, the SU(N L ) gauge group as part of the electroweak interaction coupled to scalars and fermions, described in [53], and the strong-interaction-type SU(N c ) gauge group coupled only to fermions as studied in [26]. Whenever feasible, we specialize to the SM matter content including its flavor and generation substructure. We implicitly assume the existence of further sectors -to be ignored for the purpose of this work -that cancel perturbative gauge anomalies or a global Wit-ten anomaly possibly occurring for certain N L and fermion content. More explicitly, let us consider a complex scalar for arbitrary gauge functions α i (x), α sI (x). Here T I are the generators of the su(N c ) Lie algebra with I = {1, . . . , N 2 c − 1} andḡ s is the strong gauge coupling.
The essential part of the classical action that we address in four-dimensional Euclidean spacetime reads where the scalar field amplitudeρ is the SU(N L ) invariant ρ = φ †a φ a . The indices a, b, c, · · · and A, B, C, · · · starting at the beginning of the alphabet are associated to the fundamental representations of SU(N L ) and SU(N c ), respectively, a = {1, . . . , N L } and A = {1, . . . , N c }. We explicitly introduce only one Yukawa coupling. For the limiting case of the SM, this Yukawa coupling will play the role of the top Yukawa coupling which is quantitatively the most relevant Yukawa coupling for the running of the Higgs potential. For N L = 2 it is also straightforward to introduce a bottom-like Yukawa coupling for the second component via the charge conjugated Higgs field. We suppresse possible lepton terms as well as generation indices which are implicitly understood and will be included in the counting of degrees of freedom whenever relevant. The right-handed quarks are coupled to the gluons G μ I through the covariant derivative The covariant derivative acting on the left-handed components involves also the gauge boson vector fields W μ i , The complex scalar is coupled only to the W bosons, The classical parameter space of this model is spanned by five bare couplings: the weak gauge couplingḡ, the strong gauge couplingḡ s , the Yukawa coupling for the top quarkh, the scalar mass parameterm and the scalar quartic couplinḡ λ. While the mass parameter is power-counting relevant, all other couplings are marginal.
In the remainder of this section, we review the standard perturbative analysis for the above model at one loop and only for perturbatively renormalizable couplings in the DER. In the latter approximation, we set any propagator masses to zero since they are supposed to be negligible with respect to the RG scale in the UV limit. In particular we do not consider any contributions coming from the scalar mass parameterm as well as from a nontrivial vacuum expectation value in the case where the scalar potential is in the spontanously symmetry-broken (SSB) regime. Moreover, we focus on the UV behavior of this toy model and look for totally AF trajectories. In order to address this point, we need to study the RG flow equations for the renormalized dimensionless couplings g, g s , h and λ. Their definitions in terms of the bare couplings and wave-function renormalizations are detailed later on in Sect. 5.

Gauge sector
Let us start by analyzing the RG flow equation for the gauge couplings. The RG equation for the gauge coupling of the SU(N L ) group is [3] ∂ t g 2 = η W g 2 , Here, d L γ denotes the dimension of the Clifford-algebra representation (with d L γ = 2 for the left-handed Weyl spinors of the SM). We also introduced N L f as the number of fermionic N L -tuples. In the SM, we have 3 doublets for the leptons and 9 doublets for the quarks, accounting for their SU(3) color, therefore N L f = 12. The number of scalar N L -tuples is counted by N sc , with N sc = 1 for the SM. The RG equation for the strong gauge coupling g s reads where d c γ denotes the dimension of the combined left-and right-handed Clifford algebra, i.e, d c γ = 4, and N c f is the number of quark flavors. For the SM, we have in summary In this case, both one-loop β functions are negative such that g and g s approach the AF Gaußian FP in the UV limit. In the present work, we use these β functions for various specific models differing by their matter and gauge-symmetry content.

Yukawa sector
In the present section, we retain the general N c , N L dependence as well as generic fermionic matter content specified by N L f and N c f , while we set N sc = 1. The standard one-loop RG flow equation for the top-Yukawa coupling h 2 in the DER reads where the anomalous dimensions for the scalar, the left-and right-handed Weyl spinors are

Z 2 -Yukawa-QCD model
In order to obtain a better understanding of the RG trajectories towards the UV in the three dimensional space (h 2 , g 2 , g 2 s ), let us start with the flow within the (h 2 , g 2 s ) plane. This corresponds to setting g 2 = 0 inside the RG flow equation for the top-Yukawa coupling. Moreover we choose N L = 2 for illustration as in the SM. In this case we expect a similar behavior as for the Z 2 -Yukawa-QCD model analyzed in Ref. [26], for which there is an AF region, bounded by a special AF trajectory along which h 2 is proportional to g 2 s . This behavior can be characterized in terms of a rescaled Yukawa coupling, Its β function is where The β function in Eq. (15) has only one nontrivial zero for g 2 s = 0 which isĥ 2 * = χ 2 s . A partial fixed point for a ratio of AF couplings such as in Eq. (14) has been called quasifixed point (QFP) in Ref. [25]. It is a defining condition for AF scaling solutions and a useful tool to search for such trajectories [3,7,13,26,27]. In Ref. [26] we observed also that AF requires the matter content parameter N c f to stay within a finite window for fixed N c . The upper bound of this window is given by the requirement that η G < 0, while the lower bound can be derived from Eq. (16) by demanding χ 2 s > 0. Thus, we obtain the criterion which is fulfilled by the SM parameters of Eq. (10). In this case, the ratio in Eq. (16) attains the value

Non-Abelian Higgs model
A similar analysis can be performed also by projecting the flow in Eq. (11) onto the (h 2 , g 2 ) plane, corresponding to taking the g 2 s → 0 limit. Setting N c = 3 for illustration, we can search for AF trajectories along which the Yukawa coupling becomes proportional to g 2 . Namely, we are interested in a QFP for the rescaled coupling The corresponding RG flow equation is in this case where χ 2 g reads The constraint on the matter content in order to have AF for both couplings g 2 and h 2 then is However, the lower bound is not fulfilled for the SM, resulting in a negative value for χ 2 g , We can therefore conclude that nontrivial solutions for the QFP equation ∂ th 2 = 0 do not exist in the positive part of the (h 2 , g 2 ) plane, where the only possible solution for g 2 = 0 is the trivial one h 2 = 0.
The two scenarios are illustrated for the SM case in Fig. 1, where the flows projected onto the (h 2 , g 2 s ) (left) and (h 2 , g 2 ) (right) planes are depicted. On the left panel, the special AF trajectory expressed in Eq. (18) in the g 2 → 0 limit is highlighted by a green line. On the right panel, the flow in the g 2 s → 0 limit is shown, where the trivial solution corresponding to the axis h 2 = 0 is highlighted by a purple line. It is clear from the left panel that the trajectoryĥ 2 * = χ 2 s represents an upper bound for TAF. In fact, if at some initializing RG scaleĥ 2 0 > χ 2 s , the Yukawa coupling hits a Landau pole at some finite energy scale towards the UV within this one-loop approximation. On the other hand the Yukawa coupling becomes AF for those initial values which fulfill the constraintĥ 2 0 ≤ χ 2 s . For the same reasons, there are no AF trajectories for the top-Yukawa coupling due to the negative plane. The UV repulsive QFP trajectory highlighted by the green line corresponds to the solution in Eq. (18). Right panel: the corresponding RG flow projected onto the (h 2 , g 2 ) plane. The UV repulsive QFP trajectory highlighted by the purple line corresponds to the trivial solution h 2 = 0. The arrows for both flows are pointing towards the UV value of χ 2 g in the physical quadrant of the plane (h 2 , g 2 ) > 0. The only possible way to have TAF is to set h 2 equal to the trivial null value. The QFP nature of these two special trajectories can be better understood by looking at the RG flow for the rescaled couplingsĥ 2 orh 2 as a function of g 2 s or g 2 , respectively. They correspond to IR attractive trajectories, which govern the low-energy behavior of the model, enhancing its predictive power.

SU(2) L × SU(3) c model
Let us next address the running of the Yukawa coupling in the presence of both gauge couplings. It is possible to analytically integrate the RG β function of h 2 in Eq. (11) together with the β functions for g 2 and g 2 s . As explained in Appendix 1 the matter content parameters N L f and N c f must fulfill the following necessary but not sufficient condition in order to feature total AF, which generalizes the two lower bounds χ 2 s > 0 and χ 2 g > 0 previously obtained for N L = 2 and N c = 3. In the SM case, the inequality (24) is satisfied, since χ 2 = 227/266. If the condition (24) holds, we can identify a critical surface parametrized by a function h 2 = (g 2 s , g 2 ) which represents the upper bound for TAF. In other words, for any initial condition such that h 2 0 ≤ (g 2 s0 , g 2 0 ) the top-Yukawa coupling becomes AF and approaches the Gaußian fixed point in the UV limit. As detailed in Appendix 1, this surface is a UVrepulsive surface along its normal directions, while all the trajectories on the surface itself are in the UV limit attracted towards a special one where the top-Yukawa coupling and the gauge couplings are proportional to each other. In order to find the corresponding equation for this trajectory, let us use the QFP criteria and consider first the ratio of the two gauge couplingŝ The flow equation forĝ 2 is then which has a QFP solution for g 2 = 0 at The solution g 2 =ĝ 2 * g 2 s identifies a plane in the three dimensional space of parameters (g 2 , g 2 s , h 2 ) whose intersection with the critical surface is a trajectory along which the top-Yukawa coupling is proportional to both gauge couplings. Assuming g 2 =ĝ 2 * g 2 s , we perform the same rescaling as in Eq. (14), arriving at the β function for the rescaled top-Yukawa coupling which has a nontrivial QFP solution at Alternatively, Eq. (19) could be used in the same way. Equations (27) and (29) are positive for the SM set of parameters and attain the QFP valueŝ The RG flow in the three dimensional space of couplings is plotted in Fig. 2 exhibiting the critical surface (g 2 s , g 2 ) for the SM case and the RG flow on top of it. Since this surface represents the upper bound for TAF, the directions normal to it are UV repulsive. On the surface itself however, all the trajectories are attracted towards the special one described in Fig. 2 The upper critical surface h 2 = (g 2 s , g 2 ) of total AF for the perturbatively renormalizable model in the DER and for the SM set of parameters summarized in Eq. (10). The special trajectory in Eq. (30) along which h 2 is proportional to g 2 as well as g 2 s is highlighted by a red line. It is a UV attractive (repulsive) trajectory along the directions tangent (orthogonal) to the critical surface. The intersection of with the g 2 = 0 plane is highlighted by a green line with a slope given by Eq. (18); the intersection of with the g 2 s = 0 plane is shown as a purple line satisfying h 2 = g 2 s = 0. The arrows of the RG stream flow on the top of the critical surface are pointing towards the UV Eq. (30) in the UV limit, as is highlighted by a red line in Fig. 2. In the same plot, the two trajectories in the g 2 = 0 (green) and g 2 s = 0 (purple) planes are also highlighted, corresponding to those of Fig. 1.

Scalar sector
Now we investigate the scalar sector and thus include also the running of the quartic scalar coupling λ. Its β function at one loop in the DER for our model defined by Eq. (4) is where η φ is given by Eq. (12). Since we are interested in the special trajectory described by Eqs. (27) and (29), along which the top-Yukawa and the gauge couplings are proportional, we can express h 2 and g 2 as a function of g 2 s . Thus the beta function ∂ t λ turns out to be just a function of λ and g s . Any AF solution must correspond to a particular scaling of the quartic coupling with respect to the gauge coupling. The latter is best revealed by inspecting the flow for the ratiô Here, the positive power P is either fixed by the QFP condition forλ 2 at nonvanishing g 2 s , or remains a free parameter. The β function for this rescaled Higgs coupling then receives an extra contribution coming from the running of g s . Indeed Here, we have introduced the rescaled gluon and scalar anomalous dimensionŝ which assume constant values on the QFP. Close inspection reveals that a nontrivial finite QFP solution forλ 2 in the UV limit requires P = 1/2. By choosing the SM set of parameters, see Eq. (10), we find the two rootŝ The stability properties of these two QFPs can be deduced by plotting the RG flow of λ 2 orλ 2 as a function of g 2 s , as shown in Fig. 3. The positive root (red line) corresponds to a UV-repulsive trajectory. By contrast, the negative root (green line) characterizes a UV-attractive trajectory. For any initial condition withλ 2 <λ + 2 , the rescaled quartic scalar coupling is attracted towards the negative root in the UV, and the perturbative potential appears to become unstable. On the other hand for an initial value bigger thanλ + 2 , the scalar coupling hits a Landau pole at some finite energy scale towards the UV. Thereforeλ 2 =λ + 2 corresponds to the only trajectory along which the theory is UV-complete. The perturbatively renormalizable potential is automatically stable then. This trajectory is IR attractive and the low-energy behavior is governed by the QFP valueλ + 2 which means that the theory exhibits an high degree of predictivity.
Comparing our toy-model flow to that of the SM, current data suggests that the SM flow is governed by its vicinity to the analogue of the critical surface , with the gauge couplings, the top-Yukawa coupling h and the scalar coupling λ all exhibiting a flow to smaller values above the Fermi scale. As the strong coupling g s is larger than the weak coupling g, the gauge sector has not yet reached its QFP (27). Also, the top-Yukawa coupling is below its QFP value (29), and behaves AF, cf. Fig. 1 (left panel). The scalar coupling appears to be near critical [54][55][56][57][58][59], withλ 2 being slightly below (the analogue of)λ + 2 , such that λ appears to approach zero or potentially drop below zero towards higher scales, cf. Fig. 3 (left panel). Of course, the contribution of the hypercharge U(1) group that would dominate the flow far above the Planck scale are ignored in the present discussion.
As a last remark of this section, we observe from the oneloop β functions for the Yukawa coupling and the quartic scalar coupling in the DER, Eqs. (11)- (13) and (31) that it is possible to recover the corresponding β functions for various limiting models. For instance, recovering the non-Abelian Higgs model from the general case is straightforwardly possible in the DER, by setting h 2 → 0 and g 2 s → 0. Naively, the flow equation ∂ t h 2 reduces to the one for the Z 2 -Yukawa-QCD case [26] by taking the limits g 2 → 0 and N L → 2. Whereas the flow ∂ t λ would reduce to the Z 2 -Yukawa-QCD model in the limits g 2 → 0 and N L → 1/2. This seeming contradiction can be resolved by taking the unitary-gauge limit before approaching the DER. In this way the Goldstone modes decouple from the theory and do not propagate because of their infinite mass. Thus, the Z 2 -Yukawa-QCD model simply corresponds to the limit g 2 → 0 of the present model in the unitary gauge. More details on this reduction are given in Sect. 5.

Effective field theory approach to the scalar potential in MS
Let us generalize the previously outlined construction to the inclusion of perturbatively nonrenormalizable interactions which will result in the description of new AF models. In adding higher-dimensional operators to the scalar potential in Eq. (4), we follow the EFT paradigm and start with the most widely used MS scheme for concreteness. Here, we concentrate on momentum-independent scalar self-interactions which form the effective potential. As detailed in the next sections, the consistency of these solutions requires an infinite number of higher-dimensional operators. The class of point-like scalar self-interactions is such an infinite set that becomes manageable by functional methods, as discussed in the following. The goal of the present section is to explain how to reveal these solutions and to properly account for some of their properties in a parameterization where first only a finite number of couplings with higher dimension is included. These steps then generalize to the inclusion of all interactions up to some given dimensionality in the effective Lagrangian. Still, the crucial ingredient in the construction is a treatment of the β functions of these operators that slightly differs from the standard EFT one: the scale dependence of one coupling or Wilson coefficient in the EFT expansion has to be treated as free. The subsequent sections then demonstrate that this additional freedom has to be present in any rigorous definition of the RG flow of the model, because of the infinite dimensionality of the theory space. It plays the role of a boundary condition in a functional representation of the quantum dynamics.
In the DER, where all mass parameters are neglected, it is a well known fact that higher-dimensional scalar self-interactions do not influence the running of the lower dimensional ones. This is because divergences giving rise to powers of the renormalization scale are replaced by corresponding powers of the masses in the MS scheme. However, the DER does not exhaust all possible asymptotic behaviors of a quantum field theory, as we show in the following.
We begin with a systematic polynomial expansion of the scalar effective potential. For convenience, we now switch to dimensionless renormalized quantities which are obtained by rescaling the dimensionful ones with suitable powers of the RG scale and wave function renormalizations. Precise definitions will be given below in Sect. 5. Let us call the dimensionless effective potential u. We expand the potential about the scale-dependent minimum κ, which is vanishing in the symmetric (SYM) regime, and positive in the spontaneoussymmetry-broken (SSB) regime. For instance, in the latter case a polynomial approximation of the effective potential reads where ρ is the dimensionless renormalized analog of the squared scalar field amplitudeρ = φ † φ. Generically, we expect all couplings to be generated by fluctuations. Truncating the sum at some finite N p corresponds to a polynomial approximation of the potential. For reasons of clarity, we first study the simpler limiting models, the Z 2 -Yukawa-QCD and the non-Abelian Higgs model, separately. Both models represent well-defined limiting cases of the general model with SU(N L ) × SU(N c ) gauge symmetry. In either case, we choose the remaining matter content as in the SM, cf. Eq. (10), for illustration, and perform the analysis in the massless MS scheme. Most of our results will be generalized to the full SU(N L ) × SU(N c ) model and to more general RG schemes in Sect. 6. However, already in Sect. 4 we unveil novel AF solutions for the general SU(N L ) × SU(N c ) model in the MS scheme.

Z 2 -Yukawa-QCD model
For the Z 2 -Yukawa-QCD model, the flow equations for the nontrivial minimum κ and the quartic scalar coupling λ 2 , obtained by dimensional regularization in the MS scheme, are The β function for the minimum κ involves only the couplings λ 2 and λ 3 , while the β function of a general selfinteraction λ j depends on all λ n up to n = j + 2. These β functions follow straightforwardly from the functional MS flow of the effective potential discussed in the next section.
Compared to the standard one-loop flow in the DER, which is contained in Eq. (39) in the limit κ → 0, it appears that nonvanishing values of κλ 3 or κ 2 λ 4 can considerably influence the flow of the quartic coupling. In fact, this has implications for the construction of AF trajectories. As in Sect. 2.2, we look for AF scaling solutions by means of a QFP condition forλ 2 , as defined in Eq. (32). Beyond the restriction to perturbatively renormalizable couplings, and in the parameterization of Eq. (37), similar conditions can be imposed on suitably defined rescaled couplingŝ with P 2 = 2P, cf. Eq. (32). Also the coupling κ may scale asymptotically as a definite power of g 2 s , where the real power Q is a priori arbitrary. Recursive solutions to the QFP condition can be constructed by keeping one coupling of the scalar potential as a free parameter. See Sect. V of Ref. [26] for a general description of this recursive problem. Various ways to search for scaling solutions and for performing the recursive procedure are possible. In practice, we find it useful, to express allλ n as a function of g s andκ, cf. [25]. For definiteness, we concentrate in this work on solutions exhibiting the property thatλ 2 = 0 at the QFP (though this might be a scheme-dependent statement). We now illustrate this process by considering N p = 2. At this order, we set λ 4 = 0, such that the beta functions for the two ratios of Eqs. (32) and (41) become functions ofκ,λ 2 , λ 3 , h 2 and g 2 s . The dependence on the Yukawa coupling can be eliminated by considering the special trajectory along which h 2 = χ 2 s g 2 s where the QFP value for χ 2 s is given by Eq. (16). For the SM parameters for the remaining matter content, the finite ratio χ 2 s takes the value as in Eq. (18). Thus the RG flow equations forκ andλ 2 within the MS renormalization scheme read The two terms proportional to the rescaled powers P and Q are the contributions coming from the running of the strong gauge coupling g 2 s . Its β function in the MS scheme equals the flow equation within the DER, namely Eq. (9). Then we look for QFPs with nonnegativeλ 2 andκ in the g 2 s → 0 limit, while leaving λ 3 as a free parameter. It is straightforward to rediscover the Cheng-Eichten-Li (CEL) solution, for P = 1/2,κ = 0, and λ 3 = 0. In this case, we find two QFP solutions [4,26] The resulting flow structure is similar to Eq. (36), and the same conclusions as outlined below Eq. (36) apply. Let us generalize our discussion by considering the case P > 1/2 and focus on Eq. (43). For small g 2 s , we can neglect the terms proportional to g 4P s and g 2 s in comparison with g 4(1−P) s , and we retain also λ 3 . Solving the QFP equation ∂ tλ2 = 0 for λ 3 , in the g 2 s → 0 limit, we obtain Inserting Eq. (45) into Eq. (42) and again keeping only the leading g 2 s dependence gives the QFP and a freeκ. Thus, there is a two-parameter family of AF solutions labeled byκ > 0 and Q > −2. It is identified by Eq. (46) and For completeness, let us discuss the case P < 1/2, still in the N p = 2 truncation. By the same process we first solve the equation ∂ tλ2 = 0 forλ 3 , and investigate which g 2 s scaling of this coupling might produce QFPs forλ 2 . Then we input such scaling with an arbitrary coefficientλ 3 and search for QFPs forκ andλ 2 where such couplings are finite and nonnegative. It turns out that, for any P < 1/2, there is no acceptable solution.
After having worked out the problem at order N p = 2, one might increase N p and check the stability of the known solution upon inclusion of more couplings. However, this returns the QFPs described by Eqs. (46) and (47), whereκ is free and anyλ n is a function of it, which separately arises as solutions of the equation ∂ tλn−1 = 0. For instance, at N p = 3, one finds again the same solution as before, that is Eq. (47), complemented bŷ The simplest way to address the result of this recursive problem is by considering all theλ n 's at once. As we will show in the following Sect. 4, this can be done by a functional approach where the full scalar potential u(ρ) is accounted for.

Non-Abelian Higgs model
Let us apply the same strategy as before to the non-Abelian Higgs model, setting g s → 0 and h → 0; for simplicity, we work with N L = 2. We expand the dimensionless potential around a nontrivial minimum κ as in Eq. (37). Let us choose the polynomial expansion parameter N p = 2 and retain the nonperturbatively renormalizable coupling λ 3 as a free parameter. This leads to the RG flow equation for κ and λ 2 obtained by dimensional regularization in the MS scheme and for the SU(2) L gauge group, Since we are interested in looking for AF trajectories, we rescale the couplings similar to those in Eqs. (32), (40), and (41), where the strong gauge coupling g s is replaced by the gauge coupling g, The corresponding RG flow equations forκ andλ 2 read Apart from the present use of the MS scheme, these equations generalize the ones discussed in [25] by an independent Q rescaling of the minimum κ. The two terms proportional to the rescaled powers P and Q are the contributions coming from the running of the weak gauge coupling g 2 . Its β function in the MS scheme equals the flow equation within the DER, namely Eq. (8). Since we want to construct QFP solutions whereκ approaches a finite value in the UV limit, it is possible to see from Eq. (52) that the following three values for the rescaled powers are allowed: P = 1, or P 3 = Q or In the first case where P = 1, only the first two terms in Eq. (52) contribute to the QFP equation ∂ tκ = 0 at leading order in g 2 , providing a constant solution forλ 2 . Substituting P = 1 in Eq. (53), the value of the rescaled power P 3 is fixed by the relation P 3 = Q + 2 in order to have a finite g 2 → 0 limit for the β function ∂ tλ2 . To summarize this first possible solution, we havȇ whereλ 3 and Q ≥ −2 remain two free parameters. However, this solution has to be rejected since it is not compatible with our assumption to expand the potential at its minimum, i.e., λ 2 > 0, in order to interpret the coefficients as couplings and mass parameters during the flow towards the UV. Analogous considerations can be performed also for the second possibility where P 3 = Q. In this case, the QFP solution is which admits a suitable solution with a positive value forλ 2 .
In addition, the presence of a nontrivial minimum requires thatλ 3 < 0. For completeness we stress that the third possibility with P 3 = Q +2−2P does not lead to any real solution since the QFP equation ∂ tλ2 = 0 admits only complex roots at leading order in g 2 . The construction generalizes to higher orders in the polynomial expansion N p . For instance, the P = 1 solution, c.f. Eq. (54), survives and we have, for example, for N p = 4 This still represents a two-parameter family of solutions with couplingsλ n>2 having alternating signs forκ > 0. The solution in Eq. (55) acquires a different QFP value forλ 2 , as its β function receives leading-order contributions both fromλ 3 as well as fromλ 4 . For example for N p = 4, the solution reads whileλ 3 and P 3 are still given by Eq. (55). This solution has again alternating signs for the higher order couplings if the potential is in the SSB regime. These findings motivate a full functional analysis beyond the polynomial expansion of the potential.

Full effective potential in MS
The existence of nonpolynomial structures in the functional RG flow of the scalar potential can already be anticipated from classic results of one-loop computations with fielddependent thresholds [60,61]. In this section, we stick to evaluating the loop integrals in dimensional regularization. For examples of this procedure, see [62][63][64]. According to the MS prescription, the β function equals the residue of the (d − 4) −1 poles of these integrals, which can be singled out by taking RG time derivatives followed by the d → 4 limit. Recent applications of these flow equations have shown several advantages of dealing with functional perturbative beta functions, see [65,66].
Let us begin with the general SU(N L ) × SU(N c ) model with SM matter content, i.e., N L = 2 and N c = 3. The functional flow equation for the dimensionless scalar potential at one loop in the MS scheme is where we have used the Landau gauge, and ω H and ω θ are the bosonic thresholds associated to the radial Higgs fluctuation and the three Goldstone fluctuations, The arguments associated to the gauge boson and fermionic threshold contributions are defined as The scalar anomalous dimension η φ is given by Eq. (12).
In the limiting case of the Z 2 -Yukawa-QCD model, the scalar field φ is real, therefore only the physical Higgs excitation contributes to the scalar threshold function. Also the degrees of freedom associated to the weak gauge bosons do not occur. Then, the RG flow equation for u(ρ) in this model reads On the other hand, the non-Abelian Higgs model is recovered simply by ignoring the quantum effects arising from the fermions. The β function for the dimensionless potential then reads We already know from the previous sections that AF trajectories in the theory space can be detected by simply looking for QFPs of the flow for rescaled couplings. To implement this condition in a functional approach we define a rescaled field variable x and its potential f (x) as The field amplitude ρ is then multiplied by an appropriate power of an AF gauge coupling, which is either the weak gauge coupling for the non-Abelian Higgs model or the strong gauge coupling in the Z 2 -Yukawa-QCD model. In the general model, we also use the strong gauge coupling for the rescaling. Denoting the nontrivial minimum by x 0 , we have The arbitrary rescaling power P has to be chosen as the P value corresponding to the scaling of the quartic scalar coupling such that ξ 2 =λ 2 , because we specifically look for QFPs whereλ 2 approaches a finite value in the UV limit. Notice that the relation between ξ n andλ n (and between x 0 andκ) at finite value of g 2 is a simple rescaling, but in the {g 2 , g 2 s } → 0 limit these couplings might attain different fixed-point values. Thus, the rescaling of Eq. (63) is expected to be useful as long as the quartic scalar coupling is the leading term in the approach of the scalar potential to flatness. According to the rescaling in Eq. (63), the functional RG flow equation for f (x) is thus depending on whether we use g s or g to rescale the field amplitude ρ. The anomalous dimensions η G/W , as well as η φ and ∂ t h 2 , in the MS renormalization scheme are the same as in the DER, cf. Eqs. (8), (9), (11) and (12), due to the vanishing of power-like divergent diagrams.

φ 4 -dominance approximation
In order to get closer to a full functional description, we first use a simple approximation of the β function ∂ t f (x) by asserting that the scalar fluctuations are dominated by the marginal quartic coupling in the UV limit. More precisely, we assume that the scalar potential appearing in the threshold functions takes the form u(ρ) = λ 2 ρ 2 /2. Nevertheless, we still retain the full u(ρ) dependence in the scaling term and on the left-hand side of ∂ t u(ρ) as an unknown arbitrary function of ρ. This assumption leads to the following approximation for the radial Higgs excitation and Goldstone fluctuations: We start with the general model, specifically considering the trajectories described by Eq. (30) along which the top-Yukawa and weak gauge coupling become proportional to g s in the UV limit. This yields a β function for the rescaled scalar potential f (x) that depends only on the AF strong gauge coupling g s , where the scaling dimension d x of the rescaled field includes also a contribution from the running of the strong gauge coupling, in fact where η G is given by Eq. (9). The QFP solutions for the ratioŝ g 2 * andĥ 2 * are given by Eq. (30). The QFP equation, which is obtained by the requirement that the left-hand side of Eq. (68) is vanishing, is solved by where C f is a free integration constant, parameterizing the general solution for the associated homogeneous equation. Setting C f = 0 and requiring the consistency condition f (0) = ξ 2 singles out the same solution with P = 1/2 and ξ 2 =λ + 2 as it was found in Sect. 2, cf. Eq. (36). For any nonvanishing C f the QFP potential behaves as a nonrational power of x at the origin; therefore, its second order derivative at x = 0 is singular for η x > 0. If the system is in the SYM regime, the anomalous dimension for the rescaled field is indeed positive in the Z 2 -Yukawa-QCD model for all values of P. Hence, the singularity would affect large classes of correlation functions expanded about the symmetric ground state, such that we consider such solutions as unphysical. On the other hand, in the general model and the non-Abelian Higgs model, η x can be negative for small enough values of P, because of the negative gauge-loop contribution entering in η φ .
The problematic singular behavior at the origin might be avoided in all models if there is at least one nontrivial minimum for f (x), in the spirit of the Coleman-Weinberg mechanism [60]. In fact, the system of two equations that arises by setting n = 2 in Eq. (65) can be solved for C f and ξ 2 as functions of x 0 . The additional requirement that ξ 2 is finite and positive in the g 2 s → 0 limit can be fulfilled only when P = 1. The expressions for C f and ξ 2 at leading order in g 2 If x 0 attains a finite value in the g 2 s → 0 limit, this corresponds to a potential that has a finite minimum as well as finite derivatives at this minimum, which are given by We can thus construct a family of solutions parametrized by the nontrivial minimum x 0 with the desired property that the rescaled quartic coupling at x 0 is finite in the UV limit. This is in fact a two-parameter family of solutions, as Eq. (72) is compatible with an arbitrary asymptotic scale dependence of x 0 of the form The appearance of the additional parameter Q occurs as in the EFT analysis of the previous section. More details are provided in the following for the specific case of the Z 2 -Yukawa-QCD model. Evidence for the global stability of the scalar potential f (x) can be obtained by studying the asymptotic behavior for large amplitudes x. In fact, the flow equation allows to study two different asymptotic limits, both corresponding to large amplitudes and small gauge coupling, but differing by the product g 2P s x being either small or large. The former asymptotic region is addressed by taking first the g 2 s → 0 limit and then the x → ∞ limit, where we find the following asymptotic behavior The latter asymptotic regime is obtained by the opposite order, yielding In both regimes, we find a stable potential, providing evidence for global stability.

Z 2 -Yukawa-QCD model
Within the φ 4 -dominance approximation, we can address the limiting case of the Z 2 -Yukawa-QCD model by substituting the expressions in Eq. (67) into the RG flow equation (61) for the scalar potential where the QFP solution for the rescaled top-Yukawa coupling assumes the value as in Eq. (18). The QFP equation where C f is again a free integration constant, parameterizing the general solution for the associated homogeneous equation. Setting C f = 0 and requiring the consistency condition f (0) = ξ 2 singles out the CEL solution with P = 1/2 and ξ 2 =λ ± 2 as for Eq. (44). As discussed in the general model, the potential has a log-type singularity in the second derivative at the origin for any C f = 0, as η x is always positive in this model. This problem can be avoided if f (x) admits a nontrivial minimum x 0 . Keeping x 0 as a parameter, it is possible to solve for C f and ξ 2 , which are to leading order in g 2 s . As for the general model, the rescaled quartic coupling ξ 2 can be finite only for P = 1. Moreover, the higher-order couplings at the nontrivial minimum are We can thus construct a one-parameter family of solutions which enjoy all the desired properties usually expected for a QFP potential. Their singular behavior at vanishing field values makes them invisible in an expansion for small field amplitudes. To recover the two-parameter family of solutions observed in Sect. 3.1, it is sufficient to notice that Eq. (79) still holds if x 0 scales as in Eq. (73). Inserting the latter scaling into Eq. (79), we would find precisely the results shown in Eqs. (47) and (48), as well as the predictions for all higherorder couplingŝ which can be verified within the EFT approach. Furthermore, the large-field behavior for any g 2 s > 0 is whereas in the other asymptotic regime where the limit g 2 s → 0 is taken before considering the x → ∞ limit, the large-field behavior reads In both cases, the potential appears stable.

Non-Abelian Higgs model
The limiting case of the non-Abelian Higgs model can be recovered from Eq. (68) by the substitutions The RG flow equation for the rescaled potential f (x) then becomes The quantum dimension d x includes a contribution from the anomalous dimension of the gauge vector fields, since the scalar amplitude ρ is rescaled with the weak gauge coupling g 2 : The QFP solution of Eq. (84) is that features a log-type singularity in the second derivative at the origin f (0), as long as the integration constant C f is different from zero and η x positive.
In contrast to the general SU(2) L × SU(3) c model or the Z 2 -Yukawa-QCD model, there is no real solution compatible with the consistency condition f (0) = ξ 2 for C f = 0. This reflects the conventional conclusion of triviality as seemingly evidenced by Landau-pole singularities in perturbation theory. A different situation occurs for C f = 0. In this case indeed, the second derivative at the nontrivial minimum is finite only for P = 1 and takes the value to leading order in g 2 . As this is negative, it contradicts one of our selection criteria. We can moreover find a recursive formula for all the higher-order couplings which is This is in agreement with the solution found within the EFT approximation in the MS scheme. In fact, if we express the latter equation in terms of the finite rescaled couplingsλ n andκ as in Eq. (51), we would find which coincide with Eqs. (54) and (56).
Working out the behavior of the potential in the two asymptotic regions, we find in the intermediate asymptotic regime, taking first the g 2 → 0 and then the x → ∞ limit while in the opposite order yields the large-field asymptotics Both asymptotic regions reveal that the potential is not stable. We conclude this subsection on the non-Abelian Higgs model by comparing the present φ 4 -dominance approximation with the EFT analysis: Within both approximations we discovered the P = 1 solution, which has to be rejected as it involves a negative quartic coupling which violates our assumptions. Only the EFT analysis can reveal the acceptable solution associated to P = 1/2, c.f. Eq. (57). In the latter case, we have indeed observed that the contributions from the higher-order couplingsλ 3 andλ 4 were crucial for finding the QFP value forλ 2 . We expect that neglecting the presence of these interaction terms within the loops results in the impossibility to reveal these additional AF solutions in the present context. This illustrates the limitations of the non-systematic but nevertheless useful φ 4 -dominance approximation.

Weak coupling expansion
In order to abandon the assumption about the dominance of the local four-point interaction, we now perform a parametrically controlled functional weak-coupling analysis. For this, we neglect the subleading corrections in powers of the gauge couplings to the UV-asymptotic behavior of f (x), and expand its β function for weak coupling. We first write the full flow equations of f (x) at one loop in the MS renormalization scheme for all the models under investigation. In the general case, the RG flow equation for the rescaled scalar potential reads, cf. Eq. (58), where the arguments of the threshold functions can be obtained from Eqs. (59) and (60) after having rescaled the field amplitude ρ according to Eq. (63). We thus have where as usual we consider the special trajectories along which h 2 and g 2 become proportional to g 2 s . These are characterized by the QFP values in Eq. (30). The full quantum dimension of x is given by Eq. (69).
In the Z 2 -Yukawa-QCD model, the degrees of freedom associated to the gauge bosons and the Goldstone mode are not present thus the beta function of f (x) becomes Also the QFP forĥ 2 * has to be changed and takes the value as in Eq. (18).
The non-Abelian Higgs model can be recovered from the general case by simply performing the substitutions in Eq. (83), thus ∂ t f reduces to where d x is given by Eq. (85).
Since the scalar, fermion and gauge boson loops appear with different powers of the gauge couplings, we distinguish several cases corresponding to the classification of leadingorder terms in the UV limit where {g 2 s , g 2 } → 0 but x, f (x), and f (x) stay finite.
If P > 1/2, the bosonic contributions arising from the radial and/or Goldstone fluctuations are negligible with respect to the fermionic and/or gauge boson fluctuations. Also the {g 2 s , g 2 } dependence of d x is negligible which can nevertheless be easily accounted for. Under these approximations, the QFP equations ∂ t f = 0 for the different models reduce to the flow equations obtained within the φ 4 -dominance approximation if we set ξ 2 = 0. As a consequence, the weak-coupling expansions of Eqs. (92), (94) and (95) for P > 1/2 agree with the approximation made in Sect. 4.1 as far as the UV limits {g 2 s , g 2 } → 0 are concerned.
For P < 1/2, the effects from the fermion and/or gauge boson loops are negligible and the flow of the scalar potential is the same as the flow of a purely scalar quantum field theory, where the dimension of the field is externally driven towards the canonical one as the classical sources given by g 2 s or g 2 vanish. In this case, as well as in the P = 1/2 case, the QFP equation remains nonlinear and of second order, and it does not offer straightforward analytical solutions. As such, the present approximation is not helpful, and does not offer better perspectives with respect to the results of the EFT-like analysis in Sect. 4.
This concludes our first analysis of a possible existence of further AF trajectories within the well-known perturbative MS scheme. To summarize for example the Z 2 -Yukawa-QCD model: in addition to the well-known CEL solution, we have found evidence for a family of further AF trajectories. Let us substantiate these findings by a more comprehensive analysis also addressing the question of scheme dependence further.

Renormalization in mass-dependent IR schemes
The RG equations for action functionals as obtained from a masslike scale-dependent deformation of the Gaußian part of the action have been known for a long time [67,68]. For the purpose of extending the analysis of the scheme dependence of AF solutions, we use the exact RG flow equation for the one-particle irreducible effective average action k [69] at an RG scale k, given by the Wetterich equation [69][70][71][72] where t = log k is the RG time. Here, k [ ] denotes the second functional derivative with respect to the collective field variable . The function R k encodes a general IR regularization of momentum modes near the scale k. The derivative ∂ t R k in the numerator provides for a UV regularization. The detailed form of R k therefore defines a regularization scheme within the FRG approach. Results that hold for any physically admissible regulator therefore provide evidence for scheme independence. A solution k to the Eq. (96) interpolates between the initial condition at some UV scale , k= = S cl in the form of a classical action, and the effective action k=0 = generating the 1PI correlation functions of the full quantum theory, see [73][74][75][76][77][78] for reviews.
In the following, we focus on the beta functional for a general scalar potential. For this, we solve Eq. (96) on a projected theory space spanned by the truncated action: All couplings, wave function renormalizations Z , and the effective potential U are k dependent. This truncated theory space can be viewed as a leading-order derivative expansion of the action in terms of local operators which has been proven useful, e.g., in the analysis of the RG flow of the Higgs potential [24,25,53,[79][80][81][82][83][84][85][86][87][88][89][90][91].
For simplicity, we refer to the N L = 2 case for the remainder of this section, but we will consider N c and the spacetime dimension d as arbitrary parameters. We use a gauge-fixing Lagrangian L gf of the general form where F i and F sI are the gauge-fixing conditions for the weak and the strong gauge group, respectively. The corresponding gauge-fixing parameters are ζ and ζ s ; below, we mostly quote results obtained in the Landau gauge, ζ, ζ s → 0. Moreover, the ghost Lagrangian L gh where c i ,c i , b I , andb I are the ghost fields, encodes the determinants of the Faddeev-Popov operators where α j and α sJ are the local parameters for the finite gauge transformations in Eqs. (1) and (3). In order to take into account also the threshold effects coming from the SSB regime, we decompose the scalar field into the bare vevv and the fluctuations around it. Without loss of generality we choose the radial mode in the first real component such that [25,53]: where the radial fluctuation H corresponds to the Higgs excitation and the Goldstones form a triplet. We choose the gaugefixing functional for the SU(2) L gauge group such that no mixing terms between the Goldstone modes and the gauge bosons appear in the propagators, As usual, the Higgs excitation H is not included, thus the gauge-fixing condition involves only the Goldstone bosons and not the radial mode. From the definition in Eq. (100), we identify the Faddeev-Popov operator In the gluon sector, we use standard Lorenz gauge, such that the gauge-fixing functional for the SU(N c ) gauge group and the corresponding Faddeev-Popov operator read Let us introduce also the mass parameters for the elementary fields of the Lagrangian in the SSB regime. The unrenormalized mass matrix for the gauge boson fields is where {·, ·} denotes the anticommutator. The generators for the SU(2) L gauge group are t i = σ i /2, therefore all gauge bosons acquire the same mass, Introducing the decomposition as in Eq. (101), we recover the following formulas for the mass of the scalar fluctuations, In the SYM regime, where the minimum of the potential is v = 0, all the scalar fluctuations acquire the same mass. By contrast, in the SSB regime wherev = 0 and by definition U (v 2 /2) = 0, only the radial fluctuation H becomes massive which corresponds precisely to the Higgs excitation. The angular fluctuations are instead massless in the Landau gauge and correspond to the Goldstone modes, Furthermore the unrenormalized mass for the top quark is given bȳ At this point we would like to emphasize that these parameters of the elementary fields do not necessarily have to coincide with observables of the theory in the IR. Also the term spontaneous symmetry breaking is misleading, although often used in this context, as a local gauge symmetry cannot be spontaneously broken [92]. Moreover, the vev of the Higgs field is not a reliable order parameter [93,94] as it depends on the gauge choice even if the potential has a Mexican hat-type form [95]. To formulate the spectrum of a theory with a Brout-Englert-Higgs (BEH) effect in a gaugeinvariant manner is cumbersome on a nonperturbative level due to the Gribov-Singer ambiguity [96][97][98][99][100][101]. This ambiguity states that commonly used gauge-fixing conditions like (102) and (104) are insufficient to fully fix the gauge. However, how this problem affects a BEH theory is still under investigation [102][103][104][105][106].
That the elementary fields are not observable quantities is a consequence of the Gribov-Singer problem. Nevertheless, gauge-invariant approaches have been developed to formulate the spectrum appropriately [107][108][109]. Describing the observables in terms of gauge-invariant bound states, proposed by Fröhlich, Morchio, and Strocchi, is a useful procedure [110,111]. First, it can easily be generalized to other gauge groups [112][113][114][115][116]. Second, it explains why the perturbative description of the spectrum of the weak sector of the standard model is so successful by using a one-to-one mapping of the symmetry structures among the bound states to the weak symmetry group. Thus, we will stick with the standard nomenclature (SYM, SSB, mass, ...) throughout this paper as we will concentrate on N L = 2, keeping in mind that actually not the vev breaks the gauge symmetry but the gauge fixing term and that the gauge-variant objects can be used to describe gauge-invariant observables with high precision for the weak sector of the standard model.
Since we are interested in FPs where the model asymptotically features a self-similar behavior, we study the RG flow for the renormalized dimensionless quantities. Let us introduce therefore the dimensionless renormalized U(N L )invariant scalar field amplitude and the dimensionless renormalized couplings Inserting our truncation of the effective average action Eq. (97) into the Wetterich equation (96) and projecting onto the scalar sector allows to extract the RG flow equation for the dimensionless potential In a similar manner, the β function for the dimensionless renormalized top-Yukawa coupling can be extracted. Similarly, we can obtain the anomalous dimensions for the fields which are defined as encoding the running of the scale-dependent wave function renormalizations. The functional flow equation for the dimensionless renormalized potential in the Landau gauge is given by [26,53] where v −1 d = 2 d+1 π d/2 (d/2) and the arguments of the threshold functions ω H , ω θ , ω W as well as ω F have already been defined in Eqs. (59) and (60). Let us remark here that additional contributions coming from the ghost loop, the gluon loop, and the bottom-quark loop contribute only to the running of the ρ-independent vacuum energy and thus can be ignored for our present purpose. From Eq. (115), we can extract also the flow equation for the nontrivial minimum κ in the SSB regime, The threshold functions l ( ) 0 (ω), with ∈ {H, θ, F, W}, carry the dependence on the momentum-space regularization of loop integrals specified by the form of R k . Physically, they quantify how massive modes decouple from the flow, once the RG scale crosses the mass threshold. For their general definitions see the discussion in Appendix 1. From the β functional (115) for the scalar potential, the RG flow for the scalar self-couplings can straightforwardly be derived to any order by polynomial expansion. More generally, Eq. (115) encodes the flow of the global properties of the Higgs potential to be studied below.
Similar FRG flow equations for the Yukawa coupling and for the gauge coupling, as well as FRG expressions for the anomalous dimensions of the fields, are presented in Appendix 1. In addition, we address also the two limiting cases of the Z 2 -Yukawa-QCD and the non-Abelian Higgs models.
In the SU(2) L ×SU(3) c model, we have decided to rescale the field amplitude ρ with the strong gauge coupling, i.e., x = g 2P s ρ. Therefore the functional flow equation for the rescaled potential is where the anomalous scaling dimension for the rescaled field d x is given by Eq. (69) and the arguments of the threshold functions are given in Eq. (93). The QFP solutionsĝ 2 * andĥ 2 * take the same values as in Eq. (30). Let us remind the reader here that these values have been calculated by assuming that any mass contributions induced by a nontrivial minimum in the scalar potential are negligible in the flow equations for the gauge couplings and the top-Yukawa coupling. In other words we have considered the latter beta functions in the DER, where the arguments z θ , z H , z W and z F are assumed to go to zero in the UV limit, as in Ref. [26]. The consistency of this assumption has to be tested once a scaling solution for the Higgs potential is found.
In the Z 2 -Yukawa-QCD model, the flow equation for the rescaled scalar potential takes the form Within this case, the QFP solutionĥ 2 * takes the value given in Eq. (18).
In the non-Abelian Higgs model with SU(2) L gauge group the β function for f (x) reads where the quantum dimension d x is given by Eq. (85). Since only the weak gauge coupling is involved in this model, the scalar field is rescaled via an appropriate power of g 2 , namely x = g 2P ρ.
Although in Eqs. (117)- (119) we have used the same notation for the threshold functions as in the FRG case of Sect. 5, in this section we generalize their scope and we interpret them as threshold functions in a generic scheme. In other words, the threshold functions in Eqs. (117)-(119) represent the loop contributions of the scalars, fermions, and gauge bosons in any arbitrary regularization and renormalization scheme. To each of the loop-momentum integrals we can associate generic regularization schemes which can even be different for each field. As an example, the MS scheme discussed in Sects. 3 and 4, without RG improvement, i.e., suppressing the anomalous dimensions in the threshold functions, would correspond to in d = 4. By Taylor expanding for small gauge couplings, these loop integrals to first order take the form where the first term is the β function in the UV limit where {g 2 s , g 2 } → 0 and the second term is the leading {g 2 s , g 2 } contribution upon expanding the loops and the anomalous dimension η x . This last point requires an important comment: since η x depends on the properties of the Higgs potential at the nontrivial minimum κ which could be a general function of the gauge couplings, a self-consistency check of the Taylor expansion has to be performed, once the analytic QFP solution for f (x) is computed.
According to the rescaling in Eq. (63), quantum fluctuations can contribute to the zeroth-order term [β f ] 0 only for P = 1. For example, in the two limiting models we have for both models, and for P = 1, (123) in the Z 2 -Yukawa-QCD model, and in the non-Abelian Higgs model. Thus, for P < 1 the zeroth order in the gauge couplings is trivial since no quantum fluctuations are retained. On the other hand for P = 1, we need a more detailed specification for the regulator in order to address explicit properties of the QFP solutions, as discussed below.
Aiming at the leading {g 2 s , g 2 } corrections, the values of P = 1 are simpler to address in a generic scheme, since the vertices of the theory, by assumption, scale like positive powers of g 2 s or g 2 , depending on the model under consideration. The leading contribution to δβ f is produced by Taylor expanding the threshold functions to first order in the gauge couplings. This gives rise to several coefficients which account for all the scheme dependence of the QFPs, namely with ∈ {H, θ, F, W} labeling the fluctuation modes.
Within an FRG scheme these coefficients can also be written as Here, the operator∂ t denotes differentiation with respect to t = log k acting only on the regulators, and P is the inverse regularized propagator of each field. In the FRG formalism P depends on the regularization kernel R k in the Wetterich equation, cf. Eq. (96). For more details and explicit expression see Appendix 1. The framework of Eq. (96) has been derived with k as a 1PI effective action in the presence of an IR regularization at the scale k, rather then a UV one [71]. As such, the requirement by which the shape functions provide a physical coarse-graining is that they should diverge for k → +∞ for a given p 2 , and vanish for k → 0. For monotonic shape functions, the RG time derivative is always positive, therefore A > 0. For example the piecewise linear regulator [117,118], discussed in Appendix 1, leads to A = 1/(32π 2 ).
Let us now address the case with P ≤ 1, in the general SU(2) L × SU(3) c model and also in the two limits of the Z 2 -Yukawa-QCD and non-Abelian Higgs models. 6.1 P ∈ (0, 1/2) In this window of P values and for all models under consideration, only the scalar loops contribute to the first correction in the β function for f (x). In the general model, the leading correction scales as g 2P s , where we have distinguished the contributions coming from the Goldstone or radial modes with the labels θ and H, respectively. For the Z 2 -Yukawa-QCD model, the leading correction δβ f can be recovered by simply setting A θ = 0. On the other hand, for the non-Abelian Higgs model δβ f is the same as in Eq. (127) with the substitution g 2 s ↔ g 2 . By including the leading-order correction in g 2 s , the QFP equation for f (x) becomes a second order ordinary differential equation (ODE) which can be analytically solved and leads to two different solutions. The first one is given by a special case of the Kummer function which reduces to a quadratic polynomial, The second one grows exponentially for large field amplitudes. However, we are only interested in solutions that obey power-like scaling for x → ∞, since a scalar product can then be defined on the space of eigenperturbations of these solutions [65,119,120]. Thus, we set the second integration constant to zero. By imposing the defining properties for the nontrivial minimum and the rescaled quartic scalar coupling, namely f (x 0 ) = 0 and f (x 0 ) = ξ 2 respectively, we find Therefore we can infer that, in the general model as well as in the non-Abelian Higgs model, the condition for having a nontrivial positive minimum is whereas in the Z 2 -Yukawa-QCD model the latter expression becomes simply As these conditions are satisfied for all admissible FRG regularization schemes, within the latter framework the existence of these solutions is a scheme-independent result. A particular limiting case is the one of MS, where these solutions are not present as A H = A θ = 0. We provide an interpretation of this fact at the end of this section. Let us perform the consistency check mentioned above, testing if the gauge and Yukawa coupling QFP values remain unaffected by threshold effects. The solution in Eq. (128) has been found by expanding the beta function in Eq. (117) for small g 2 s while keeping x, f (x), and its derivatives finite. By inserting the QFP solution for f (x) into the expression for the anomalous dimension η x , we find that there is no contribution to δβ f coming from η x for any P < 1/2. In fact, the leading terms in η x scale as either g 2 s or g 8P s , and thus are negligible with respect to the scalar contributions in Eq. (127) scaling as g 2P s . As a matter of fact, it is also true that contributions proportional to g 8P s in the scalar anomalous dimension modify also the asymptotic UV behavior of the top-Yukawa coupling for P ≤ 1/4, eventually leading to a different QFP value for h 2 . A more detailed explanation of this fact can be found in Appendix 1. Thus, the consistency test is passed only by the 1/4 < P < 1/2 scaling solutions. For P = 1/4 this does not require a change of the QFP solution for f (x) we have just discussed, as only the QFP value ofĥ 2 * changes. For P < 1/4 instead the current solutions for f (x) are no longer valid for the Z 2 -Yukawa-QCD model and for the general model, and they survive only in the non-Abelian Higgs model, where the Yukawa coupling is not present. A different QFP for P < 1/4 might still be possible in Yukawa models, if κ and h 2 exhibit asymptotic scaling powers different from the ones discussed in this section. This behavior might require strong threshold phenomena and decoupling of some degrees of freedom. We leave this analysis for future investigations.

P = 1/2
For this value of P, the contributions from the scalar loops mix with the contribution from the fermionic loop and/or the gauge boson loop. Let us consider first the general SU(2) L × SU(3) c model. In this case the first leading correction to β f is where the QFP valuesĝ 2 * andĥ 2 * are given in Eq. (30). The QFP equation [β f ] 0 + δβ f = 0 is a second-order inhomogeneous linear ODE. Its solution consists of the general solution of the homogeneous part, c.f. Eq. (128), plus a particular solution of the nonhomogeneous equation, which is a quadratic polynomial in x. We find therefore the QFP solution where the rescaled quartic coupling ξ 2 remains a free parameter. The rescaled potential has a nontrivial minimum at whose positivity requires Let us now investigate the two limiting cases. In the non-Abelian Higgs model, we have implying the QFP potential which has a nontrivial minimum at The condition for having a positive x 0 implies Also the results for the Z 2 -Yukawa-QCD model are easily recovered from the general model by setting where the QFP value ofĥ 2 * is exactly the one in Eq. (18). For example, the expression for the nontrivial minimum and the condition for its positivity become We conclude that, in all models under consideration and for P = 1/2, the QFP equation β f = 0 admits scaling solutions which are in the SSB regime. For any FRG scheme, the scheme-dependent coefficients A assume positive values and thus define a range of values for the rescaled quartic coupling ξ 2 that parametrizes the new AF solutions. The changes of this range for different FRG schemes, corresponds to the expected mapping of the coupling space onto itself induced by a scheme change. Again MS appears special. The vanishing of the A's signals the fact that the leading order contribution in a weak coupling expansion is in fact quadratic in f , such that at the present order only the canonical scaling term contributes. The latter is expected to be purely classical and the corresponding QFP solutions unphysical. Hence, also for P = 1/2 and within the leading weak-coupling approximation we find no new AF scaling solution in MS. This can be seen, as already stated in Sect. 4.2, as a deficiency of the present approximation strategy in the MS case, since at leading order it suppresses all interaction effects and at next-to-leading order it returns the full nonlinear secondorder ODE for the QFP potential. Let us recall however, that analyzing the latter equation with an EFT-like approximation in Sect. 3, we did observe novel AF trajectories for P = 1/2 in the non-Abelian Higgs model, within the MS scheme.
Also for P = 1/2, the contribution due to the anomalous dimension η x is subleading with respect to the correction in Eq. (133). In addition, these solutions pass the consistency check that the β function for the top-Yukawa coupling still has the same QFP solution as given in Eq.
where the QFP valuesĝ 2 * andĥ 2 * are given in Eq. (30). The FP equation [β f ] 0 + δβ f = 0 remains a first-order ODE whose analytical solution is where the rescaled quartic scalar coupling ξ 2 remains an arbitrary integration constant. The QFP potential has a nontrivial positive minimum at If one adopts the same scheme for both loops, such that A F = A W , the latter condition is not satisfied in the SM case, where Let us then consider the two limiting models separately. For the Z 2 -Yukawa-QCD model obtained from the substitutions in Eq. (141), the expression for the nontrivial minimum becomes which would be negative for any positive value of ξ 2 and general FRG scheme coefficient A F > 0. The opposite situation occurs in the non-Abelian Higgs model. From the substitution in Eq. (137), the nontrivial minimum now reduces to For any positive value of A W and ξ 2 the scalar potential for the non-Abelian Higgs model is in the broken regime and features a nontrivial minimum. Also for these values of P, we need to perform the consistency check that the QFP of the top-Yukawa coupling remains unaffected: we observe that the contribution to δβ f from the anomalous dimension η x is subleading with respect to the fermionic and gauge boson contributions and that the top-Yukawa coupling scales as g 2 s in the UV limit. It thus has the same QFP solutions as in Eq. (30).
We emphasize that for all values of P < 1 the QFP solutions we have obtained are analytic in x in the present weak-coupling approximation. In the previous Sect. 3, this was implemented by construction, since we have projected the functional flow equation onto a polynomial ansatz. In the present analysis, this happens because the contributions to β f producing non-analyticities are accompanied by subleading powers of g 2 s or g 2 for P < 1. Indeed, both the appearance of the anomalous dimension η x in the scaling term (2 + η x )x f (x), and the contributions from the threshold functions proportional to x 2 would produce a singularity of f (0) for any {g 2 s , g 2 } = 0, as discussed before in Sect. 4, see also below. Knowing about the presence of this singularity for any P at {g 2 s , g 2 } = 0, we can accept the previous solutions only if x 0 > 0. This appears to be possible in all models under investigation for P ≤ 1/2. In addition, it is possible in the non-Abelian Higgs model, and in the general SU(2) L × SU(3) c model for P ∈ (1/2, 1), in the family of FRG schemes where an IR regularization is provided.
Moreover, we want to stress that all the solutions obtained for P < 1 are consistent with the assumptions made at the beginning: the arguments of the threshold functions z θ , z H , z F and z W go to zero in the UV limit, and the flow equations for the gauge couplings and the top-Yukawa coupling can be treated as in the DER.
A general observation that could be raised against the validity of the present approximation, in the FRG framework, is that neglecting the nonlinear contributions to the flow equation of f (x) might miss crucial terms and produce spurious or unphysical QFP solutions. In particular, part of the universal one-loop contribution, the one which arises from the Taylor expansion of the threshold functions to second order in z θ , z H , z F and z W , is not accounted for in the present discussion of the P < 1 scaling solutions. However, the inclusion of the latter contributions as well as of further nonlinearities, up to the full complexity of the FRG flow equations of Sect. 5, has been performed, with specific regulator choices, for the non-Abelian Higgs model [24,25] as well as for the Z 2 -Yukawa-QCD model [26] in fact confirming the results of this leading-order weak coupling expansion.
As for the fate of these solutions in MS, a similar conclusion as for the P ≤ 1/2 cases holds. The vanishing of the quadratically divergent coefficients A results in the vanishing of the interaction effects retained at the leading order of the present approximation scheme. Thus, no new AF trajectory is visible in the scheme for P < 1. However, this does not exclude the existence of new AF scaling solutions in MS altogether. In fact, as discussed in Sects. 3 and 4 as well as in the following, a two-parameter family of these solutions occur in the general SU(2) L × SU(3) c model and in the Z 2 -Yukawa-QCD model, at P = 1.

P = 1
For P ≥ 1 the β functional to zeroth order in the gauge couplings [β f ] 0 accounts for the full nonlinearity of the gauge and the fermion loops. In fact, let us consider first the general SU(2) L × SU(3) c model. The leading RG flow equation for the rescaled scalar potential reads For P = 1, the QFP solutionsĝ 2 * andĥ 2 * will not have the same values as in Eq. (30), since it is no longer true that the RG flow equations for the gauge couplings and the top-Yukawa coupling can be treated as being in the DER. Indeed, the arguments z W and z F are finite and do not approach zero in the UV limit. As a consequence, we expect that g 2 and h 2 are still proportional to g 2 s in the UV limit, but with QFP solutions which depend nontrivially on x 0 . We can nevertheless considerĝ 2 * andĥ 2 * as finite ratios. The corresponding QFP equation [β f ] 0 = 0 can be solved analytically and leads to the integral solution where c is an arbitrary integration constant. For instance, within the FRG framework the piecewise linear regulator of Appendix 1 would give a Coleman-Weinberg-like potential, which has a log-type singularity at the origin in the second derivative, given by the term ∼ x 2 log x. This remains true in any scheme. In fact by Taylor expanding the threshold functions l (F) 0 (y) and l (W) 0 (y) in Eq. (150) around y = 0, a logarithmic divergence of the integral arises precisely from the quadratic term of this expansion. In other words, the appearance of this singular behavior is as universal as the one-loop beta function of the marginal couplings. Thus we expect that this feature survives also in the full g 2 s -dependent solution. The freedom of choosing the parameter c allows us to construct QFP solutions, which circumvent the problem of nonanalytic structures at the minimum of the potential by developing a nontrivial minimum away from the origin. The defining equation for the minimum f (x 0 ) = 0, involves an integral of two arbitrary threshold functions and might be hard to solve analytically for x 0 . Still, it can straightforwardly be used to express c as a function of x 0 . From the point of view where the latter is the free parameter labeling the QFP solutions, the natural question then is as to whether it can be chosen such that f (x 0 ) = ξ 2 is positive and finite in the g 2 s → 0 limit. As for P < 1 the answer to this question involves some scheme dependence, which is encoded in the coefficients where ∈ {F, W}. For all schemes which can be embedded into the FRG, the coefficients A F,W (x 0 ) of Eq. (152) are similar to the ones defined in Eq. (125). The evaluation of these loop integrals at nonvanishing values of z F,W simply results in adding k 2 z F,W to P F,W in the denominator of the integrand function, whose sign is positive for any IR regularization scheme, since the shape function is monotonically increasing.
For any threshold functions, the rescaled quartic scalar coupling simplifies to and the condition that ξ 2 be positive requires For small values of x 0 this condition reduces to the one in Eq. (145). In the FRG framework, we already observed in the previous section that the latter is not fulfilled in the SM case, if the same regulator is chosen for both fields. We expect this conclusion holds for any value of x 0 , and for generic FRG schemes. In fact, in the opposite limiting case of large x 0 , A (x 0 ) reduces to an x 0 -independent number, and Eq. (155) again is equivalent to Eq. (145). For the special case of the piecewise linear regulator of Appendix 1 the coefficients of Eq. (153) read such that Eq. (155) is not fulfilled for all x 0 ≥ 0. Let us now consider the two limiting cases. For the Z 2 -Yukawa-QCD model where the gauge-boson loop is absent, the inequality in Eq. (155) becomes which cannot be fulfilled for any admissible IR regularization scheme. Conversely, in the non-Abelian Higgs model the fermion loop is absent, resulting in the condition which is satisfied by any admissible regulator in the FRG framework.
In MS these conclusions get twisted. In fact, Eq. (120) results in which does indeed hold for the SM, as was already observed in Sect. 4.1, e.g., c.f. Eq. (71).

Summary
Let us summarize the results obtained so far. In most cases, a simple leading-order weak-coupling expansion of the flow equations suffices to unveil a two-parameters family of novel AF trajectories. In a general FRG scheme based on an IR regulator, these solutions are effectively labeled by the position of the nontrivial minimum x 0 of the rescaled scalar potential f , or equivalently by its finite asymptotic ratioκ, and by the power P defining the approach of the quartic coupling λ 2 to the Gaussian FP.
In the general SU(2) L × SU(3) c model, for the specific case of SM matter content and under the assumption that the same kind of regulator is chosen for all fields, we have successfully constructed solutions with P ∈ [1/4, 1/2]. For P < 1/4 we do not exclude the possible existence of solutions for which h 2 exhibits an asymptotic scaling different from the one supported within the DER. The same conclusions can be drawn for the Z 2 -Yukawa-QCD model. In the non-Abelian Higgs model, we have recovered the P ∈ [0, 1] solutions already described in [24,25]. For completeness, the P > 1 case is discussed in Appendix 1, where we recover the known solutions for the non-Abelian Higgs model and we find no AF trajectories for the Yukawa models.
Different IR regulators result in a change of range of the attainable values for the finite ratioλ 2 = ξ 2 . Special regulator choices which assign different regularizations to different degrees of freedom can also change the range of attainable values of the asymptotic power P. For instance, they can make the P ∈ (1/2, 1] solutions available also in the SU(2) L × SU(3) c model, c.f. Eqs. (145) and (155).
The MS scheme appears to correspond to a peculiar limit of this remapping of allowed parameter ranges. In fact, the two-parameter family of new AF solutions in this case is labeled by the finite ratioκ and by the power Q, which describes the asymptotic scaling of the running dimensionless vev κ. The power P instead is fixed to P = 1 in the general SU(2) L × SU(3) c model as well as in the Z 2 -Yukawa-QCD model, while the non-Abelian Higgs model exhibits these solutions for P = 1/2.
These results are schematically summarized in Table 1. Here, for each model under consideration, we list the attainable values of the asymptotic power P related to the scaling of the quartic coupling λ 2 . Furthermore, we recall the two remaining continuous parameters that play the role of coordinates on the space of AF scaling solutions in the different RG schemes: Q andκ in the MS scheme, P and ξ 2 in the FRG schemes. Table 1 Summary of the family of new AF solutions constructed in this work, for the models described in Sect. 2, and in the RG schemes we have analyzed. For each solution, we provide the value(s) of P for which they occur, and the (remaining) variables that parametrize the space of solutions. In the FRG schemes, we refer to the setup where the same regulator is used for all fields, as detailed in Sect. 6.1. For the most general model we specify the SM matter content in the flow of the gauge and Yukawa couplings SU(2) L × SU(3) c Z 2 -Yukawa-QCD Non-Abelian Higgs

Conclusions
The recently discovered class of new asymptotically free RG trajectories in various non-Abelian particle models [24][25][26] calls for a critical assessment of their scheme independence.
In contrast to standard perturbative investigations, these new solutions become visible beyond the deep Euclidean region, because threshold effects can play an important role on all scales. As one-loop universality is no longer guaranteed beyond the deep Euclidean region, we have investigated in this work whether the existence of these UV complete trajectories is universal.
We have specifically studied an SU(N L ) × SU(N c ) gauge theory coupled to fermions which are either charged under both gauge groups or only with respect to the SU(N c ) group. Moreover, we consider a scalar field being in the fundamental representation of SU(N L ) but a singlet of SU(N c ). Thereby, we cover the non-Abelian subsector of the standard model, and generalize previous investigations of Z 2 -Yukawa-QCD and non-Abelian Higgs models. The latter models represent particular limits of our general SU(N L ) × SU(N c ) model.
In order to make contact with the most widely used RG scheme, we have started with the one-loop β functions in the MS scheme. The new solutions become visible by inclusion of higher-dimensional operators as well as a nonvanishing expectation value for the scalar field inducing mass thresholds for various fields. This becomes even more transparent from the flow equation for the entire scalar potential. The β functional for the potential u(ρ) or its rescaled form f (x) allows to include arbitrary generalizations of the scalar selfcoupling as well as to address the global stability properties of the fixed-point potential.
The functional formalism is well suited for the investigation of the present scenario of UV completion, as it allows to monitor the effect of higher-dimensional operators on the perturbatively renormalizable couplings which is crucial to trigger asymptotic freedom. From a naive effective-fieldtheory perspective, higher-dimensional operators would be suspected to have stronger effects in the UV, thus suggesting a loss of control on the high-energy asymptotics. By contrast, the inclusion of infinitely many such operators in the functional renormalization group equations provides access to the novel asymptotically free solutions. Here, higher-order operators also feature asymptotic scaling and are controlled by progressively higher powers of the gauge coupling, thus resulting in a stronger suppression at high energies. In this way, higher-order couplings evade the naive power-counting argument that associates such operators to nonrenormalizable theories in perturbation theory. Also, our solutions do not exhibit a high-energy enhancement suggested by powercounting, but remain controlled by the asymptotically free gauge coupling. For the same reason, we do not observe violations of perturbative unitarity. The failure of perturbativeunitarity-violation arguments is also well known for theories which are non-perturbatively renormalizable, i.e. possessing a fixed point with a finite number of relevant directions [121]. Indeed functional descriptions involving higherdimensional operators have been provided for several unitary models, e.g. minimal conformal [122] or superconformal models [123] in two dimensions, and their Wilson-Fisher relatives in three Euclidean dimensions with [124][125][126] and without supersymmetry [127][128][129][130], by means of the methods we have discussed in Sects. 5 and 6.
In the case of the Higgs-Yukawa models addressed above, the functional analysis reveals that the family of new asymptotically free trajectories can be classified by rescaling parameters (such as P) and the position of the rescaled potential minimum (such asκ or x 0 ). These parameters can in turn be related to the boundary conditions specified for the quasi fixed-point potential, and thus the boundary conditions for correlation functions [25]. For instance, the standard perturbative solution of the Z 2 -Yukawa-QCD model, i.e., the Cheng-Eichten-Li solution [4], can be recovered as a particular boundary condition, which is implicitly assumed within perturbative computations.
As a probe of scheme independence, we have investigated the flow equations also using the functional-renormalizationgroup framework, that facilitates to perform the momentum space regularization with a general class of regulator functions. We demonstrated that the new solutions exist for any physically admissible functional-renormalization-group regulator which encodes the details of the decoupling of massive modes. On a line of constant physics, a change of the regularization scheme induces a map of the coupling space of initial conditions, e.g., given in terms of bare couplings at an initial UV scale , onto itself. Our results for the new trajectories show that such a mapping in theory space also includes the boundary conditions, i.e., parameters such as the rescaling power P, which classify the new asymptotically free trajectories.
The present work exemplifies common features of those models studied so far which allow for the construction of new asymptotically free trajectories: (i) a model needs an asymptotically free sector with a coupling that can serve as a rescaling parameter. While the rescaling is merely a technical step, it helps introducing the quasi-fixed-point concept that can be extended to full functions such as quasi-fixed-point scalar potentials. (ii) the (asymptotic) boundary conditions for such full quasi-fixed-point functions need to be specified explicitly or implicitly, e.g., through a suitable choice of rescaling parameters or the scaling of higher-dimensional operators. (iii) threshold effects need to be accounted for as they can invalidate the conventional naive analysis in the deep Euclidean region. Whereas we have demonstrated that the new asymptotically free trajectories are clearly visible in a weak-coupling analysis, the conventional perturbative treatment does not fully account for (ii) and (iii) but is confined to implicit trivial choices. We emphasize that ingredient (i) does not mean that the full perturbatively renormalizable model needs to be asymptotically free; in fact, the non-Abelian Higgs model is a counterexample, as generic perturbatively renormalizable models are not completely asymptotically free due to scalar triviality. Our construction demonstrates that the asymptotic freedom of the gauge sector can be sufficient to seed asymptotic freedom in the complete model using ingredients (ii) and (iii).
This last point inspires to use the new trajectories for more realistic model building and, eventually, phenomenology. Asymptotic freedom can render a particle physics model high-energy complete. Our construction thus has the potential to solve the consistency problems of the standard model, namely the triviality problem being manifest in the Higgs and the U(1) sector of the standard model at high energies. Our results for the non-Abelian models with semi-simple gauge groups may already be useful for the construction of asymptotically free grand unified models. If the construction can also be extended to a U(1) hypercharge sector, a more direct completion of the standard model could exist as a fundamental particle physics theory, offering new routes to model building. This is also attractive in the light of our observation that quasi-fixed-point trajectories can exist with a higher degree of predictivity than their conventional perturbatively renormalizable counterparts. The first requirement is necessary in order to provide for the existence of AF trajectories in the space of parameters (h 2 , g 2 , g 2 s ). Incidentally, this criterion is fulfilled for the SM set of parameters. In fact, if this condition is not satisfied the integral I (∞) diverges; consequently, there will be a finite RG time t critic at which the right-hand side of Eq. (A4) is zero and the top-Yukawa coupling hits a Landau pole, i.e., I (t critic ) = (h 2 0 a h ) −1 . On the other hand, whenever the integral I (∞) converges, AF is guaranteed as long as the righthand side of Eq. (A4) remains positive. This is precisely the second inequality in Eq. (A6).
By denoting (g 2 s0 , g 2 0 ) = (a h I (∞)) −1 , the equation h 2 0 = (g 2 s0 , g 2 0 ) represents an upper critical surface for TAF in the (h 2 , g 2 , g 2 s ) space. In fact, if h 2 0 > (g 2 s0 , g 2 0 ), the negative one-loop contributions to ∂ t h 2 involving the gauge boson fluctuations are suppressed in the UV limit, and the dominant positive scalar fluctuations drive the top-Yukawa coupling towards a Landau pole. By contrast, all couplings approach the Gaußian FP towards the UV for all initial conditions such that h 2 0 ≤ (g 2 s0 , g 2 0 ). Let us concentrate in more detail on the RG flow on the surface : By virtue of its critical nature, represents a UV repulsive surface along its orthogonal directions. Only for those initial conditions such that h 2 0 = (g 2 s0 , g 2 0 ), the integrated RG trajectories will remain on the critical surface itself for all t > 0 and will approach the Gaußian FP in the UV limit t → ∞. Next, we observe that the one-loop β functions for the gauge couplings are independent of the Yukawa-coupling in the DER. As a consequence, the gauge-coupling flow on each slice of constant h 2 looks the same as within the (g 2 s , g 2 ) plane, see Fig. 4 (left panel). Therefore, there will be a special AF trajectory also on the surface along which the two gauge couplings are proportional to each other, representing a UV attractive trajectory. This special RG trajectory can be characterized more explicitly by introducing "rotated" gauge couplings G 2 s = cos θ g 2 s + sin θ g 2 , G 2 = − sin θ g 2 s + cos θ g 2 . (A7) Choosing tan θ =ĝ 2 * , i.e., corresponding to the proportionality factor of the gauge couplings, it suffices to study the flow of h 2 for G 2 = 0. Here, the β functions for h 2 and G 2 exhibiting the phase diagram and RG flow as plotted in Fig. 4 (right panel). Here the highlighted red line represents the QFP trajectory along which the top-Yukawa coupling is proportional to G 2 s , and the corresponding QFP value satisfies We conclude that the special trajectory in Fig. 4 along which all three perturbatively renormalizable couplings are proportional to each other corresponds exactly to the QFP solution characterized in Eq. (30).

Appendix B: RG flow equations in mass-dependent IR schemes
This appendix complements Sect. 5 and presents the FRG equations for the Yukawa coupling, the fields anomalous dimensions, and the beta functions of the gauge couplings.

RG flow equations for the matter content
Within the truncation of Eq. (97), different choices are possible for projecting on the scalar anomalous dimension, the anomalous dimension for the left-handed fermions, as well as for the Yukawa coupling [53]. For the present purpose, we compute ∂ t h 2 and η φ by projecting the flow onto the radial scalar operator H in the SSB regime which represents the physical Higgs excitation. Moreover, we concentrate on the fermionic wave function renormalizations associated with the massive top quark. Accordingly, the flow equation for the Yukawa coupling in the Landau gauge reads [26] coupling at one loop is then determined by the wave function renormalization. For the SU(2) L gauge group and for one complex scalar as well as N L f left-handed fermion doublets, η W reads The threshold functions L W,F,φ guarantee the decoupling across mass thresholds in the SSB regime. They are normalized such that L W,F,φ (0) = 1, implying that we obtain the standard results in the DER. In the SSB regime, the renormalized and dimensionless mass parameters are proportional to the vev κ of the scalar potential, Here, we have also allowed for a bottom-type Yukawa coupling h 2 b j in addition to the top-type Yukawas h 2 t j associated to the j-th generation in order to model the decoupling of all quark mass thresholds. Ignoring higher-loop resummations, the threshold functions read (B10) For the SU(3) c gauge group, the wave function renormalization η G for the gluon field is where in this case the fermionic threshold function L F takes the form where μ 2 Q j = h 2 Q j κ is the mass for the j-th quark where j has to be understood as a multiindex, labeling the position within the left-handed doublet as well as possible generation copies.
Let us briefly comment on the case of a rescaling power P > 1 where further simplifications arise: as discussed in Appendix 1, the gauge boson and fermion fluctuations decouple from the dynamics for P > 1, since their masses diverge in the UV limit. In this regime, all loop contributions from the massive modes drop out of the gauge coupling flows. Depending on the nature of the neutrinos, they either decouple as well if they are Dirac neutrinos with a mass term generated via a Yukawa coupling to the vacuum expectation value. Or as Majorana neutrinos, they could essentially behave as nearly massless particles in the DER and thus would not decouple from η W . Counting the essentially massless neutrinos by n ν , we obtain In this case, the ratio of the two gauge couplings, defined in Eq. (25), takes the QFP valuê On the other hand, if we treat the neutrinos as Dirac particles, their contribution decouples from η W and the latter QFP value changes intô

Appendix C: Threshold functions
For the application of the functional-RG equation (96), we choose a regulator R k which is diagonal in field space. We keep the freedom to have different regulators, specified by corresponding sub-or superscripts, for the Higgs scalar (H) and for the Goldstone bosons (θ ), as this is possible in the SSB regime. Notice however that we do not distinguish between different runnings of the wave function renormalization for the radial excitation Z H and the Goldstone modes Z θ in the SSB regime at the present level of our truncation. Thus we use Z φ as a collective wave function renormalization. Its flow is stored in the anomalous dimension η φ which coincides with the projection rule for the radial excitation. Similarly, due to the choice of covariant gauges, we can have different regulators for the transverse gluons (GT) or W bosons (WT) and for the longitudinal gluons (GL) or W bosons (WL). In the Landau gauge used in this work, only transverse gauge bosons propagate, hence we can avoid the further specifications L and T, and simply write G or W. Finally, we account for independent regularizations of the left-handed (L) and right-handed (R) spinors. One of the lefthanded Weyl spinors together with its right-handed partner becomes massive in the SSB regime. The contributions of the corresponding Dirac field are denoted with F. The regularized kinetic (or squared kinetic) terms are The definition of analogous terms for θ , G, and W are identical to the one for H. Their RG time derivative is defined through the operator Recall however that we approximate The loop momentum integrals appearing on the r.h.s. of Eq. (96) are classified by defining the corresponding threshold functions. Most of the threshold functions used in this work can be found in Appendix A of Ref. [53]. However, the present abbreviations differ from the ones adopted there. In the latter reference, any scalar contribution carries the label B, and the letter G was used for the gauge bosons, corresponding to the W bosons in our work. This applies to the different versions of l d 0 appearing in Eq. (115); to those of l d . Two kinds of threshold functions require a more detailed discussion. One threshold function is called a d 3 in Ref. [53] and m d 1,2 in the present work, as well as in Ref. [88]. Also, the function called a d 1 in Ref. [53] is correspondingly renamed m d 2 in this work. For clarity, we provide here the explicit definitions of these two kinds of threshold functions Here, the operator∂ t denotes differentiation with respect to t = log k acting only on the regulators.
As an example, the piecewise linear regulator [117,118], where x = q 2 /k 2 , yields the following results for these threshold functions For any regulator, we note that ω 1 = 0 renders m ( becomes independent of η L ). Correspondingly, we can drop the associated η argument, as has been used in Eqs. (B4) and (B5).

Appendix D: Weak gauge-coupling expansion for P < 1/4
At the end of Sect. 6.1 we have found indications that the UV behavior of the top-Yukawa coupling might change for P ≤ 1/4 due to leading terms proportional to g 8P s in the scalar anomalous dimension. Even though we demonstrated that this does not modify the QFP solution for the rescaled potential f (x), the persistence of the QFP value of the top-Yukawa coupling is an important consistency check of our construction.
We start with QFP solution for f (x), as given in Eq. (128) for P < 1/2, which has a nontrivial minimum x 0 and substitute it into the RG flow equation for h 2 , c.f. Eq. (B1). Within a weak-coupling expansion, we recover the same expression as for the DER, c.f. Eq. (11), plus an extra term proportional to g 8P s , The extra term arises from the scalar loop threshold function m For P = 1/4, the extra term is of the same order as the DER limit for ∂ t h 2 . Therefore the QFP value for the ratio among the top-Yukawa coupling and the strong-gauge coupling becomes x 0 -dependent and readŝ whereĝ 2 * is given by Eq. (30). Moreover, the QFP value stays positive as long as The situation is different for P < 1/4 where the last term in Eq. (D1) becomes leading. In order to capture the (in)existence of possible different scaling solutions for the top-Yukawa coupling with respect to the strong-gauge coupling, we look for QFP solutions for the ratiô with E > 0. With such a rescaling, the only possible QFP value iŝ In view of the condition in Eq. (131), this solution is, however, negative and thus unphysical. In other words, the presence of a nontrivial minimum for the scalar potential prevents the existence of scaling solutions for the top-Yukawa coupling for all P < 1/4. Nevertheless scaling solutions for f (x) do exist also for all P < 1/4 and do not depend on the asymptotic behavior of the top-Yukawa coupling.

Appendix E: Weak gauge-coupling expansion for P > 1
It is worthwhile to study the possibility of new AF trajectories for the case P > 1 in the general model. Previous studies found that they exist in the non-Abelian Higgs model [24,25], whereas no valid solutions have been found in the Z 2 -Yukawa-QCD model [26]. As the general model interpolates between the two limiting cases, a search for their scheme-independent (in)existence is particularly instructive. Our result confirms their existence in the non-Abelian Higgs model as a special limiting case, whereas the general model does not feature the same mechanism. For P > 1, the arguments in the fermionic loop z F and in the gauge boson loop z W defined in Eq. (93) diverge in the {g 2 s , g 2 } → 0 limit. Therefore, in order to capture the correct UV behavior, we are led to Taylor expand the threshold functions l with ∈ {F, W}, such that B > 0 for general RG schemes providing an IR regularization. For instance, the piecewise linear regulator yields the positive value B = 1/(32π 2 ). On the other hand, the bosonic thresholds associated to the radial Higgs fluctuation and the three Goldstone fluctuations are always subleading in the UV for P > 1. Moreover, the anomalous dimension η x can contribute to leading order in the β function for f (x) for these values of P. For example, it has been observed in the non-Abelian Higgs model that η W becomes leading for P > 2 since it is proportional to g 2 [25]. The same conclusion holds also for the Z 2 -Yukawa-QCD model: the anomalous dimension η G , being proportional to g 2 s , contributes to leading order for P > 2. In order to discuss also the possibility to have an anomalous dimension for the rescaled field, we therefore solve a QFP differential equation where η x is retained as a parameter which has to go to zero in the UV limit.
Let us start by investigating the general SU(2) L × SU(3) c model. By keeping the terms linear in z −1 F,W , the β function for the rescaled potential becomes x . (E2) The presence of a singular term at the origin induces a corresponding pole in the QFP solution which is obtained from integrating the QFP condition β f = 0 (at fixed g s ), where c is the integration constant of the first-order ODE. Additionally, there is also a log-type singularity in the second derivative at the origin. In fact by Taylor expanding the scaling term for small η x , we get a contribution proportional to x 2 log x. This singularity can be avoided if the potential admits a nontrivial minimum x 0 such that f (x 0 ) = 0. The latter condition can be solved for c, yielding a function c(g 2 s , x 0 ), and by substituting it into the definition of the rescaled quartic scalar coupling ξ 2 = f (x 0 ). This provides the following expression We observe that the QFP potential has a nontrivial minimum for positive ξ 2 only if The anomalous dimension η x depends nontrivially on x 0 , but the consistency criterion holds that η x → 0 in the UV limit. From this property we can infer that for any finite value of ξ 2 the behavior of the nontrivial minimum as a function of the strong gauge coupling is x 0 ∼ g 2(P−1)/3 s in the g 2 s → 0 limit. By substituting this scaling inside the definitions for z W,F we observe that z W,F → ∞ in the UV limit, as stated above.
For the non-Abelian Higgs model, there are no fermion fluctuations, so the right-hand side of Eq. (E5) vanishes and the criterion is satisfied in any scheme. The evidence found in [24,25] for the existence of new AF trajectories is thus confirmed in a scheme-independent manner. By contrast, there are no weak gauge contributions in the Z 2 -Yukawa-QCD model, implying that the left-hand side of Eq. (E5) is zero in this limiting case. Hence, the criterion cannot be satisfied.
For the SU(2) L × SU(3) c model, a diagonal choice of regulators B W = B F inside Eq. (E4) would result in ξ 2 < 0 for SM matter content. It seems that the Eq. (E5) still leaves room for legitimate models in the general case. However, this is not the case as detailed in the following.
In writing Eq. (E2), we have also assumed that h 2 and g 2 scale with g 2 s . This is true in the DER but has to be verified outside this regime. For P > 1, the arguments z F and z W diverge in the UV limit, corresponding to a divergence of the gauge-boson and the top-quark thresholds. Physically, this means that they decouple from the theory and do not propagate. In Sect. 1, we have seen that the anomalous dimensions η G and η W within the decoupled regime are still proportional to g 2 s and g 2 , respectively. Therefore, the scaling g 2 ∼ g 2 s is still valid also outside the DER but the constant of proportionality depends on the number of decoupled degrees of freedom. For the SM case, the QFP value forĝ 2 is given by Eq. (B14). By contrast, the β function for the top-Yukawa coupling changes drastically beyond the DER. Thus the scaling h 2 ∼ g 2 s might no longer be valid outside the DER. This leaves a loophole in the argument presented in Ref. [26], where we had assumed that the approximation of the β function ∂ t h 2 in the DER holds also for P > 1. This loophole will be closed by the following analysis: as an example, let us assume that {z W , z F } → ∞ and z H → 0. Since we are looking for solutions with a nontrivial minimum, we can set z θ = 0. By expanding Eq. (B1) and retaining only the leading terms in g 2 s , the β function for the rescaled top-Yukawa coupling, defined in Eq. (D4), reduces to Let us assume the nontrivial minimum to scale with a power (P − Q) of the strong gauge coupling, x 0 = g 2(P−Q) sκ , with P ≤ Q or P > Q.
A careful analysis among all the possible combinations between the scaling powers P > 1, E > 0, and Q, leads to the conclusion that one or more of the assumptions above are violated for any combination which allows to have a QFP solution forĥ 2 .
As an example, let us consider the case where P = E and P = Q. The rescaled top-Yukawa coupling has a QFP solution depending on the nontrivial minimum which readŝ However this solution is not compatible with our assumptions, since z F would stay finite and does not diverge in the UV limit. All other cases can be analyzed analogously.
We conclude that the general model does not feature new AF trajectories for P > 1 similarly to the Z 2 -Yukawa-QCD model, whereas they do exist in the non-Abelian Higgs model in any scheme covered by our analysis.