Discovering and quantifying nontrivial fixed points in multi-field models.

We use the functional renormalization group and the [Formula: see text]-expansion concertedly to explore multicritical universality classes for coupled [Formula: see text] vector-field models in three Euclidean dimensions. Exploiting the complementary strengths of these two methods we show how to make progress in theories with large numbers of interactions, and a large number of possible symmetry-breaking patterns. For the three- and four-field models we find a new fixed point that arises from the mutual interaction between different field sectors, and we establish the absence of infrared-stable fixed-point solutions for the regime of small [Formula: see text]. Moreover, we explore these systems as toy models for theories that are both asymptotically safe and infrared complete. In particular, we show that these models exhibit complete renormalization group trajectories that begin and end at nontrivial fixed points.

Introduction

The O(N) Wilson–Fisher fixed point appears in a large variety of systems where it controls the universal critical behavior in the infrared (IR) scaling regime [1-5]. Generalizations of this universality class appear in the context of coupled-field models, e.g., for the two-field model [6-11]. Depending on the number of field components and dimension d, one finds that different fixed points (FP) govern the IR behavior of the model. Two of these, the decoupled (DFP) and isotropic fixed point (IFP), can be deduced from the existence of the Wilson–Fisher fixed point. While the DFP is characterized by a complete decoupling of the fields and therefore can effectively be regarded as a model for two independent vector fields, the IFP displays a complete symmetry enhancement to an rotational symmetry. However, the two-field model features another so-called biconical fixed point (BFP), which emerges due to the nontrivial interactions between the two field sectors. The BFP is fully coupled, i.e., mixed interactions are nonvanishing, but it does not show an enhanced symmetry, as the IFP does; see, e.g., Refs. [10-12]. Interestingly, it turns out that the model in dimensions exhibits exactly one IR-stable fixed point (with no more than two relevant directions) for any pair of values and . In this work, we address the question whether generalizations of the two-field model to the case of three and four fields allow for further unprecedented fixed-point solutions that are relevant for the IR scaling behavior of the respective model. Moreover, we look for additional confirmation of our previous study of fields in three dimensions [13], where no stable fixed point (with no more than three relevant directions) was found for small values of – in contrast to the two-field model.1 In our previous study of this system [13], we searched for fixed points using the nonperturbative functional renormalization group (RG) [14-19]. Within this scheme, the -functions are non-polynomial functions of the couplings and it is therefore challenging to make sure that numerical fixed-point searches do indeed uncover all stable fixed points of the system. To address this problem, we match the solutions of the renormalization group -functions derived within the framework of the functional RG to those obtained with the Wilsonian momentum-shell RG by employing an expansion in . The -expansion features -functions that are polynomials of the couplings. A comprehensive study of all fixed points that are continuously connected to the Gaussian FP at , is therefore straightforward. On the other hand, the full nonlinear -functions of the functional RG yield reasonable estimates of the stability of fixed points already at low orders in the approximation. Indeed, comparing fixed points in both RG schemes, we can convincingly identify stable fixed points and determine their stability regions in the space spanned by the values of the field components in arbitrary dimensions.

More recently, interacting fixed points have become an active field of research in four-dimensional models, which are explored in the context of an ultraviolet (UV) completion for gravity [20, 21] as well as QFTs including matter fields [22-24]. In this setting, an interacting fixed point provides a well-defined microscopic starting point from which a fundamental quantum field theory (QFT) valid on all scales, can be defined. Here, we add another example to the collection of toy models for asymptotic safety that apply to QFTs in low dimensions (see, e.g., Refs. [25, 26]). In our example, we focus on the question how both the UV and the IR limit of the RG trajectory are determined by interacting fixed points with different degrees of symmetry.

Effective action functional for multi-field models

Building on previous work [13], we derive the functional RG -functions from the nonperturbative flow equation for the scale-dependent effective action functional in d-dimensional Euclidean space [14] (see, e.g., [15-19] for reviews). Our starting point is an ansatz for to leading order in the derivative expansionwhere the summation runs over n distinct field degrees of freedom (defined in the -dimensional vector space representation of the symmetry group). The renormalization factors and the effective potential U are both assumed to be scale-dependent.2

In the following, we exploit the symmetry of the model and write U in terms of the field invariants . Furthermore, we choose to expand the effective potential to some finite order around a possibly nonvanishing scale-dependent minimum , where and . We havewhere are the scale-dependent couplings. The above expansion (2) effectively introduces a large number of couplings (at each value of the RG scale parameter, ), but only a few of them appear in the bare action defined at the scale .

The -functions to one-loop order in the -expansion may also be obtained from the nonperturbative RG flow equation. This is achieved by employing an expansion around the upper critical dimension and restricting the functional space to those operators that appear in the bare action.3 In the case of two fields these -functions agree with those given in Ref. [11], as expected by one-loop universality. In the following, we consider a model in dimensions with three different field degrees of freedom, , and , with , , and field components, respectively, and compare our results [13] explicitly with the Wilsonian momentum-shell RG to one-loop order in the -expansion. As outlined in Sect. 1, our main goal is to complement the functional RG with the -expansion to identify and characterize all possible multicritical scaling solutions relevant in the IR scaling regime. A similar strategy was also chosen in Ref. [27], where the Polchinski version of the nonperturbative RG [28] was contrasted to the -expansion to investigate multicritical points for a scalar theory with a single order parameter.

RG fixed points to in the three-field model

Employing the Wilsonian momentum-shell RG to one-loop order in the -expansion, we obtain the following set of -functions:4

Let us briefly highlight the difference of the above -functions to those derived within the functional RG approach [13]: In the latter case higher-order interactions (generated by the RG flow toward the IR) are explicitly taken into account. Therefore, the -functions for the quartic couplings receive contributions from higher-order couplings (their scale-dependence being characterized by additional -functions that are determined within this approach). Moreover, the functional RG represents a massive renormalization scheme, and accordingly, mass parameters explicitly enter all -functions, resulting in their non-polynomial structure. Details can be found in Refs. [12, 13].

We solve for the zeros of the -functions, Eqs. (5 6 7 8 9)–(10) to obtain the FPs of the RG flow. Their stability is captured by the (critical) scaling spectrum, defined in terms of the eigenvalues of the stability matrix at the FP:where and . We refer to a FP as IR-stable, if all eigenvalues (11) are negative. Mass-like perturbations in the bare action are always relevant. It is the parameters corresponding to these mass-like operators that need to be tuned in order to reach the IR scaling solution (assuming that the microscopic parameters of the model are in the domain of attraction of that particular FP).5 If no stable FP is found, additional fine tuning might be necessary to observe a continuous phase transition with universal scaling exponents. However, usually this is experimentally unfeasible, and we therefore conjecture that the corresponding systems will not feature multicritical behavior. In such a case the RG flow trajectories diverge and we expect that if a phase transition is observed it will be of first order. Nevertheless, the divergence of RG trajectories might manifest itself only deep in the IR. E.g., the RG trajectory of a theory for which the couplings at the UV scale are chosen to be close to a particular symmetry-enhanced subspace (where another fixed point with four relevant directions exists) might display a weak scale-dependence, a very slow walking, over a wide range of scales. For all practical purposes, such a scenario is hard to distinguish from conformal scaling behavior.

Within the one-loop -expansion our study uncovers the following expected scaling solutions:These FPs were previously identified in the framework of the functional RG in Ref. [13] where their scaling and stability properties were discussed in detail.

Isotropic fixed point (IFP), featuring a symmetry enhancement to an symmetry and consequently has coordinates ;

Decoupled fixed point (DFP) is characterized by a complete decoupling of all sectors, i.e., ;

Decoupled isotropic fixed points (DIFP) are characterized by a partial decoupling of the sectors as well as partial symmetry enhancement. One representative in this class is given by , while , and features an enhanced symmetry;

Decoupled biconical fixed points (DBFP) are partially decoupled but feature no symmetry enhancement, e.g., we might have and .

Within the -expansion, we confirm our previous finding that the three-field models in dimensions exhibit regions in the space of field components where no IR-stable FP exists; cf. Fig. 1. Specifically, this implies that particular three-field models with a given set of do not feature multicritical behavior without additional fine tuning. A similar absence of IR-stable multicritical FPs was observed in Ref. [29] where the effect of competing order was investigated on fermionic quantum criticality (see also Ref. [30]).

Fig. 1

IR-stable and marginally stable FPs for the three-field model in the one-loop -expansion () in the -plane, where (left) and (right). We identify the IFP (blue inverted triangle), the DFP (green upright triangles), the DIFP (purple filled circles), and the DBFP (blue filled squares). In addition to these scaling solutions, we find another IR-stable fully coupled FP, the ACFP (pink filled diamonds), while no stable FP is found for small values of and (white circle)

Fully coupled FPs

To uncover additional IR-stable FPs, we inspect the scaling solutions as a function of the parameters ; see Fig. 1.

Asymmetrically coupled FP

Our main result is the discovery of a new FP, which is completely coupled, i.e., , , and , but which does not feature any symmetry enhancement. In the following, we will refer to this scaling solution as the asymmetrically coupled fixed point (ACFP). It defines a genuine new universality class that cannot be obtained as a generalization of the Wilson–Fisher FP, and occurs for the first time in the three-field model. This new universality class relies crucially on the presence of three competing orders, and cannot occur in systems with a smaller number of order parameters.

To illustrate its properties, we give the corresponding values of the dimensionless, renormalized couplings and the critical exponents at selected points in the -plane, for ; see Table 1. Mass parameters do not appear in the -functions to the given order of the -expansion. Therefore, the three relevant scaling exponents , , and are not provided in the following.

Table 1

As we decrease the number of field components and pass through the point , (), the fully coupled FP ceases to be IR stable – the scaling exponent becomes non-negative as the value of is lowered

N_1	N_2	N_3	λ_200	λ_020	λ_002	λ_110	λ_101	λ_011
1	8	8	0.105	0.062	0.062	0.020	0.020	-0.002
1	8	9	0.109	0.062	0.059	0.013	0.008	-4×10^-4
1	8	10	0.110	0.063	0.056	0.012	0.002	-7×10^-5
N_1	N_2	N_3	θ_4	θ_5	θ_6	θ_7	θ_8	θ_9
1	8	8	0	-0.096	-0.255	-0.996	-0.996	-1.000
1	8	9	-0.019	-0.046	-0.274	-0.977	-0.999	-1.000
1	8	10	-0.007	-0.037	-0.294	-0.985	-1.000	-1.000

We follow the fully coupled asymmetric FP along the direction, where we expect that it should collide with the DFP at some critical value of N; cf. Fig. 2. In fact, we find that the two FPs exchange their stability properties at (, ). That is, at the collision point the exponent that decides about the stability properties of the scaling solution, , changes its sign for each of the two solutions. If we attempt to continue the asymmetric FP to smaller values of N, we observe that it disappears into the complex plane at , , together with another fully coupled FP which is always unstable – both FPs become inaccessible for small values of . From these results one might conclude that the fully coupled asymmetric FP will not be of any significance experimentally: The one-loop -expansion seems to suggest that there is a threshold value , for all , below which the ACFP disappears completely (cf. Fig. 1). We show in Sect. 5 that the functional RG provides a quantitatively more reliable estimate for the critical values of .

Fig. 2

Varying the number of field components , while keeping fixed, we observe that the asymmetric fully coupled FP (pink filled diamonds) appears together with a second fully coupled (unstable) FP (rose diamond). By inspecting the sign of the exponent , which determines the stability of the considered scaling solution, we conclude that the ACFP is IR stable until where it collides with the DFP (green upright triangles). From there the DFP takes over stability

It is interesting to note that the spontaneous creation/annihilation of two FPs, at least one of which is IR stable and therefore might be relevant for the multicritical scaling behavior in the IR, does not appear in the models. In the case of two coupled order parameters, we may associate exactly one IR-stable FP to each pair of values . Such a scaling solution can then be continued to all values of , but in that process it might lose its stability to another FP. That is, in principle each of the possible multicritical universality classes is accessible for all values of field components and via additional fine tuning of the parameters. This is in sharp contrast to the three- and four-field models (see Sect. 4).

Generalized BFP and regions without IR stable FP

In general, the space of renormalized couplings features closed subspaces that are characterized by enhanced symmetries and the decoupling phenomenon: Whenever one of the sectors decouples, and the couplings between sectors vanish, fluctuations cannot regenerate the mixed couplings, and therefore the RG flow stays within that space, making it an RG-invariant subspace. With the discovery of the new asymmetrically coupled FP we may complete this picture in the following way: We may state that each of these subspaces (excluding its symmetry-enhanced or decoupled subspaces) contains at least one FP. In fact, from our analysis we find that almost all of these subspaces will feature an IR-stable FP for a particular set of values , with one notable exception, the BIFP, cf. Appendix 1. This scaling solution is associated to the partially symmetry-enhanced subspace and is nowhere stable.

While in principle such a universality class exists, it would require a higher degree of fine tuning to reach it. Thus, the associated pattern of symmetry-breaking is not expected to be relevant experimentally. The BIFP would be a natural candidate to take over stability from the IFP as soon as it becomes unstable, just as the BFP takes over stability from the IFP in the two-field model. The additional relevant directions of the BIFP prevent this scenario from being realized, and imply that the three-field case features a region in the space of the that is devoid of stable FPs.

Multi-field theories as toy models for asymptotic safety and IR-completeness

In order for a QFT to provide a viable description of a set of degrees of freedom and their interactions on all scales, i.e., in order for the theory to be fundamental, it must reach a renormalization group FP in the UV and IR, respectively.6 Here, we provide a set of models that feature a large number of complete trajectories that run into nontrivial FPs both in the UV and IR.

In this context, it is important to realize that in principle a given FP can be reached asymptotically in either one of the two limits, if it features at least one critical exponent that differs in sign from the others. If a FP should be reached in the UV, all irrelevant couplings need to be tuned in such a way that the RG trajectory lies within the UV-critical hypersurface of the FP. In the context of high-energy physics, this implies that the values of all irrelevant couplings correspond to predictions of the model, i.e., for the model to be asymptotically safe, there is exactly one possible value for each irrelevant coupling. On the other hand, if the FP is reached in the IR, the renormalization group flow is automatically drawn toward it along the irrelevant directions, and it is the relevant directions that require tuning.

The large number of interacting FPs in our model provides a variety of different complete trajectories, connecting pairs of nontrivial FPs, subject to global properties of the flow. Furthermore, due to the possibility of symmetry enhancement at FPs, multi-field models are also of interest from the point of view of fundamental physics. For instance, it has been conjectured that quantum gravity should exhibit a violation of Lorentz symmetry, connected to anisotropic scaling in the ultraviolet [31, 32]. As violations of Lorentz symmetry are strongly constrained in the IR, such a setting requires a rather precise restoration of Lorentz symmetry at small scales. Here, we identify a set of models (defined by the number of field components in the different field sectors, ) where a symmetry enhancement – in our case an enhanced rotational symmetry in field space – requires additional tuning. That is, the symmetry enhancement scenario is thus considered “unnatural”. Interestingly, we also observe a number of realizations of these theories, in which an enhancement of symmetry is the most natural IR-endpoint of a trajectory, as all other existing FPs require a higher degree of tuning in order to reach them. Note that some degree of tuning is always required, as no fixed point comes with only IR-attractive directions. Demanding that the trajectory ends in an IR fixed point thus requires tuning at least three parameters. To avoid symmetry enhancement in these cases requires additional tuning.

In particular, we will focus on two examples: The first involving the ACFP as a UV fixed point, thus defining a toy model for an asymptotically safe model. Here we pick , where the ACFP has one IR-relevant direction; this triggers a flow to the DFP in the IR. As a second example, we consider a region of where, e.g., the IFP is stable (with three relevant directions) it is a natural candidate FP for RG trajectories in the IR; see Fig. 3. Thus, a model that has been rendered asymptotically safe, e.g., by defining it at the DFP, can only be infrared complete when at least three directions are tuned and it is the symmetry-enhanced IFP which provides the lowest number of relevant directions. In this case, IR FPs with a lower degree of symmetry will typically require a higher degree of fine tuning.

Fig. 3

Upper panel We plot RG flow trajectories connecting the DFP in the UV (large positive s) with the symmetry-enhanced IFP in the IR (large negative s). We show the couplings (solid red line) and (dotted blue line) for as a function of the RG scale . In this case the IFP is the stable FP and the enhanced symmetry emerges naturally as a result of the flow toward the IR. Lower panel We show a trajectory that connects the ACFP in the UV to the DFP in the IR for the case , for which the DFP is IR stable, and the ACFP has one IR-relevant direction, which triggers the flow to the IR. We show the couplings (solid red line), (dotted dark red line), (dashed blue line) and (dot-dashed dark blue line)

Four-field model to one-loop order in the -expansion

We proceed in an analogous manner for the four-field model. Our main goal here is to confirm that the two novel features of the class of -field models – the possible existence of theories without an IR-stable FP and the existence of a new fully coupled FP – carry over to the case of larger numbers of fields. The -functions are given by the obvious generalization of Eqs. (5 6 7 8 9)–(10) to the case where one additional field degree of freedom with symmetry is added. Determining the zeros of the beta functions, we find that a new FP which is fully coupled and does not feature any symmetry enhancement, exists and is stable at selected points in the space of the , i.e., a FP that appears for the first time in the four-field case similar to the role of the ACFP in the three-field case, cf. Table 2. In this context, stability is of course defined as the existence of no more than four relevant directions.

Table 2

Stable, fully coupled FP in the four-field case. There are four relevant directions; the corresponding exponents , ..., are not provided

N_1	N_2	N_3	N_4	λ_2000	λ_0200	λ_0020	λ_0002	λ_1100	λ_1010	λ_0110	λ_1001	λ_0101	λ_0011
1	9	9	9	0.110	0.059	0.059	0.059	0.007	0.007	-1×10^-4	0.007	-1×10^-4	-1×10^-4
1	10	9	9	0.110	0.056	0.059	0.059	0.001	0.006	-2×10^-5	0.006	-2×10^-5	-1×10^-4
1	10	9	8	0.109	0.055	0.059	0.062	0.003	0.008	-8×10^-5	0.013	-1×10^-5	-3×10^-4
N_1	N_2	N_3	N_4	θ_5	θ_6	θ_7	θ_8	θ_9	θ_10	θ_11	θ_12	θ_13	θ_14
1	9	9	9	-0.012	-0.028	-0.028	-0.295	-0.295	-0.295	-0.986	-1.000	-1.000	-1.000
1	10	9	9	-0.004	-0.016	-0.024	-0.295	-0.314	-0.314	-0.993	-1.000	-1.000	-1.000
1	10	9	8	-0.009	-0.022	-0.046	-0.294	-0.306	-0.314	-0.975	-0.999	-1.000	-1.000

The new FP collides with the DFP at and becomes unstable. Moreover, our results indicate that no FP is stable, e.g., at the point , . Together with the results in Table 2, this suggests that a structurally similar picture to the three-field case carries over to the four-field case: The IFP will be stable for very small values of the , before it is destabilized. Keeping fixed and increasing , we pass through a regime without a stable FP, i.e., with nonuniversal behavior only. At , the new FP then appears from the complex plane, and is stable until it collides with the DFP, that takes over stability for all larger values of the . Based on our findings in the three-field model, we expect that our estimates for the , at which FPs are stable, are considerably larger than the correct values. As can be tested within, e.g., the LPA 4, the functional RG is more reliable when it comes to quantitative estimates.

Based on our findings in the three- and four-field case, we therefore conjecture that models with larger numbers of competing orders will not feature multicritical behavior without additional fine tuning.

Results from the functional RG

Details on the derivation of the functional RG equations for the model can be found in Ref. [13]. As an important difference to the results of the -expansion we note that the critical values of , at which the DFP becomes stable, are found to lie at much lower values even in the simplest possible truncation of the functional RG equations, the local potential approximation (LPA). It is defined by an expansion of the scale-dependent effective potential (2) in terms of point-like interactions, without taking into account the scale-dependence of the renormalization factors, i.e., (see Eq. (1) for the definition of the parameters and couplings in the given model). The LPA-type truncation of the functional RG provides a quantitatively more precise estimate than the one obtained from the one-loop -expansion, which may be confirmed independently by employing nonperturbative scaling relations (see Ref. [33]).

Employing an LPA truncation to fourth order in the fields and including the scale-dependence of the renormalization factors, i.e., (which we will refer to as LPA in the following), we confirm the qualitative behavior of the -expansion: Fixing and increasing , the IFP becomes unstable around . For larger values of N there is no IR-stable FP, until the new ACFP appears and becomes stable; cf. Fig. 4. In contrast to the -expansion, this already happens at . Finally, at the asymmetrically coupled FP exchanges its stability with the DFP, which remains IR stable for all . The mechanism by which the new ACFP appears is completely analogous to the situation observed in the one-loop -expansion: It appears from the complex plane together with another FP and immediately takes over stability. Note, however, that the region of values N where the ACFP is stable is shifted to significantly smaller values of N bringing it into the reach of physically interesting models. Thus, the asymmetrically coupled FP might actually be of interest for efforts to establish the phase diagram of strongly correlated many-body systems either experimentally or via lattice Monte Carlo techniques, for an overview, see, e.g., Ref. [34]. Using scaling relations to estimate the stability regime for the DFP, we find that the LPA slightly overestimates the width of the region where the ACFP is stable. We expect that the region where the ACFP is stable becomes even smaller at higher orders of the LPA. In fact, this might account for the fact that it was not discovered in our previous analysis [13] based on a LPA to eighth order in the fields. We generically expect that a truncation of eighth will be sufficient to provide quantitatively reasonable estimates for the critical exponents.

Fig. 4

We show the values for the critical exponent obtained from the functional RG (LPA 4 + ) as a function of (, ), for the IFP (blue inverted triangles), the DFP (green upright triangles), and the new coupled asymmetric FP (pink filled diamonds). We find the same picture as in the -expansion: The ACFP is IR stable in a small region before the DFP takes over stability. The upperand lower boundary values to this region are significantly smaller than those obtained from the one-loop -expansion (): (cf. Fig. 2)

Our present results clearly highlight the strength of a combination of the functional RG with the -expansion: The latter allows us to compile a complete list of all FPs that can be continuously connected to the Gaussian FP at , whereas the former provides us with a quantitatively more reliable estimate of the stability regions of the different FPs. Combined, these methods allow us to arrive at a complete picture of stable FPs in the space of the while minimizing the computational effort.

Conclusions

With this study we identify a new fully coupled FP in the dimensional three- and four-field models. While we find that this FP is indeed IR stable for some values of , it does not lie at real fixed-point values for the couplings at other values of the . This is in stark contrast to the O(N) theory or the class of models [6-11] where the relevant scaling solution(s) are either IR stable or can be reached via additional fine tuning. In addition to identifying a new FP, we confirm our previous finding [13] that for certain multi-field models there is no IR-stable multicritical scaling solution. This behavior is directly tied to the properties of the new fully coupled fixed points, the ACFP and the BIFP. Thus, this study has further clarified the reason for the absence of multicritical scaling solutions: While in the two-field model [6-11] different FPs exchange stability only through a collision of two fixed points at real values of , the three- and higher-field models feature the additional possibility that FPs emerge from the complex plane.

Our study plays out the strengths of two methods: The -expansion allows for a straightforward identification of all FPs that can be continuously connected to the Gaussian FP in dimensions, as the -functions are polynomial in the couplings. In contrast, the fixed-point search is more involved with the functional RG due to the non-polynomial nature of the -functions. However, the functional RG provides better quantitative results already at low orders of the LPA. This can be seen clearly for the example of the DFP. To estimate its stability regime, we may apply an exact scaling relation [8, 9, 35–37] to determine the exponent from critical exponents of the O(N) Wilson–Fisher FP [1-5]. By doing so, we find that the result from the LPA at fourth order in the fields provides a quantitatively more reliable estimate for the scaling exponent than the -expansion at one-loop order. Taken together, the two methods thus allow for an efficient identification of all existing FPs, using the -expansion at low orders, followed by a leading order determination of the stability regimes and critical exponents with the functional RG.

The identification of distinct interacting FPs in three- and four-field models also allows us to explore RG trajectories that define both UV- and IR-complete QFTs in dimensions. In general, multi-field theories provide a large number of such trajectories and typically feature two distinct regimes when it comes to the question of symmetry enhancement in the IR: For values of where a symmetry-enhanced FP is IR stable, all other FPs require a higher degree of fine tuning to reach them in the IR. Thus, IR symmetry enhancement appears as a “natural” possibility that requires the least amount of fine tuning. In contrast, for other values of , the same symmetry-enhanced FP will feature additional relevant directions, giving rise to the familiar notion that an enhancement of symmetry typically requires additional fine tuning.

Our findings might have implications for possible UV completions of coupled scalar models in dimensions. We observe that the asymmetrically coupled FP can be found in the -expansion, i.e., it emerges from the Gaussian FP at . Thus we conclude that no nontrivial FP exists for these models in , unless it lies within a strongly nonperturbative regime at very large values of the couplings. This implies that, e.g., inflationary models with several scalar fields are not UV complete, but instead they most probably feature Landau poles at finite scales. Interestingly, a coupling to gravity could facilitate a UV completion in the context of asymptotically safe models. Studies suggest that a gravitational FP persists when the effects of several minimally coupled scalars are taken into account [38]. It is of course interesting to understand whether a similar statement applies to interacting matter models. In particular, the new universality classes that we discuss in this paper and which are inherent to n field models () could potentially survive an extension to 4 dimensions, when gravitational effects are added, as these generically seem to shift Gaussian FPs to interacting FPs [39]. Thus gravity might extend the upper critical dimension for this interacting FP to . Following the methods discussed in [40-42], an assessment of this scenario could be possible. In the context of scalar dark-matter models, where the coupling to other matter fields is less relevant, the existence of such scalar-gravity FPs could provide a predictive UV completion.

Table 3

Fixed-point values and critical exponents for the symmetry-enhanced BIFP. A positive exponent indicates that this FP is not IR stable

N_1	N_2	N_3	λ_200	λ_020	λ_002	λ_110	λ_101	λ_011
1	1	1	0.084	0.084	0.068	0.084	0.115	0.115
1	2	2	0.085	0.085	0.091	0.085	0.054	0.054
N_1	N_2	N_3	θ_4	θ_5	θ_6	θ_7	θ_8	θ_9
1	1	1	0.130	-0.011	-0.101	-0.491	-0.708	-1
1	2	2	0.097	-0.017	-0.053	-0.667	-0.836	-1

Table 4

For comparison, we show fixed-point values and critical exponents for the biconical FP in the two-field case. There are only two relevant directions in this model and it is therefore the exponent , which determines the IR stability of the scaling solution

N_1	N_2	λ_20	λ_02	λ_11
1	2	0.084	0.067	0.115
3	2	0.085	0.091	0.054
N_1	N_2	θ_3	θ_4	θ_5
1	2	0.130	-0.491	-1
3	2	-0.053	-0.836	-1