Discovering and quantifying nontrivial fixed points in multi-field models

We use the functional renormalization group and the ϵ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\epsilon $$\end{document}-expansion concertedly to explore multicritical universality classes for coupled ⨁iO(Ni)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\bigoplus _i O(N_i)$$\end{document} 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 Ni\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$N_i$$\end{document}. 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][2][3][4][5]. Generalizations of this universality class appear in the context of coupledfield models, e.g., for the O(N 1 ) ⊕ O(N 2 ) two-field model [6][7][8][9][10][11]. Depending on the number of field components N i 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 a e-mail: a.eichhorn@imperial.ac.uk b e-mail: mesterh@itp.unibe.ch c e-mail: scherer@thphys.uni-heidelberg.de fields and therefore can effectively be regarded as a model for two independent vector fields, the IFP displays a complete symmetry enhancement to an O(N 1 + N 2 ) 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][11][12]. Interestingly, it turns out that the O(N 1 ) ⊕ O(N 2 ) model in d = 3 dimensions exhibits exactly one IR-stable fixed point (with no more than two relevant directions) for any pair of values N 1 and N 2 . 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 n = 3 fields in three dimensions [13], where no stable fixed point (with no more than three relevant directions) was found for small values of N i -in contrast to the twofield model. 1 In our previous study of this system [13], we searched for fixed points using the nonperturbative functional renormalization group (RG) [14][15][16][17][18][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 = 4 − d. 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 d = 4, 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 N i 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][23][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 = [{φ i }] in d-dimensional Euclidean space [14] (see, e.g., [15][16][17][18][19] for reviews). Our starting point is an ansatz for to leading order in the derivative expansion where the summation runs over n distinct field degrees of freedom φ i (defined in the N i -dimensional vector space representation of the O(N i ) symmetry group). The renormalization factors Z i 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 ρ i = φ 2 i /2. Furthermore, we choose to expand the effective potential to some finite order M ≥ 2 around a possibly nonvanishing scale- 2 Here, we include terms to O(∂ 2 ) and neglect a possible field dependence of the scale-dependent renormalization factors Z i (that are evaluated at the minimum of the effective potential U ). Terms of the type ∼ ∂φ 2 i 2 , which in principle contribute at the same order in the derivative expansion, are not taken into account.
where λ m 1 ··· m n are the scale-dependent couplings. The above expansion (2) effectively introduces a large number of couplings (at each value of the RG scale parameter, 0 ≤ k < ), but only a few of them appear in the bare action S = (k = ) 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 d = 3 dimensions with three different field degrees of freedom, φ 1 , φ 2 , and φ 3 , with N 1 , N 2 , and N 3 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.
We solve for the zeros of the β-functions, Eqs. (5)-(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 β M ≡ β m 1 ··· m n and λ M ≡ λ m 1 ··· m n . We refer to a FP as IR-stable, if all eigenvalues (11) are negative. Masslike 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.
Within the -expansion, we confirm our previous finding that the three-field models in d = 3 dimensions exhibit regions in the space of field components N i where no IRstable FP exists; cf. Fig. 1. Specifically, this implies that particular three-field models with a given set of (N 1 , N 2 , N 3 ) 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]).

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

Asymmetrically coupled FP
Our main result is the discovery of a new FP, which is completely coupled, i.e., λ 101 = 0, λ 110 = 0, and λ 011 = 0, 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 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 N 2 and N 3 (white circle) 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 λ m 1 m 2 m 3 and the critical exponents θ 4 , θ 5 , . . . at selected points in the Table 1. Mass parameters do not appear in the β-functions to the given order of the -expansion. Therefore, the three relevant scaling exponents θ 1 , θ 2 , and θ 3 are not provided in the following.
We follow the fully coupled asymmetric FP along the N ≡ N 2 = N 3 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 (N = 10, N 1 = 1). That is, at the collision point the exponent that decides about the stability properties of the scaling solution, θ 4 , 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 N = 8, N 1 = 1, together with another fully coupled FP which is always unstable -both FPs become inaccessible for small values of N i . From these results one might conclude that the fully coupled asymmetric FP will not be of any significance experimentally: The one-loopexpansion seems to suggest that there is a threshold value N i 5, for all i = 1, 2, 3, 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 N i .
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 O(N 1 ) ⊕ O(N 2 ) models. In the case of two coupled order parameters, we may associate exactly one IR-stable FP to each pair of values (N 1 , N 2 ). Such a scaling solution can then be continued to all values of N i , 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 N 1 and N 2 via additional fine tuning of the parameters. This is in sharp contrast to the three-and four-field models (see Sect. 4).  . 2 Varying the number of field components N ≡ N 2 = N 3 , while keeping N 1 = 1 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 θ 4 , which determines the stability of the considered scaling solution, we conclude that the ACFP is IR stable until N = 10 where it collides with the DFP (green upright triangles). From there the DFP takes over stability

Generalized BFP and regions without IR stable FP
In general, the space of renormalized couplings λ m 1 ... m n 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 N i , 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 N i 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, N i ) 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 N 1 = 1, N 2 = N 3 = 11, 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 N i 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.

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 O(N 1 ) ⊕ O(N 2 ) ⊕ O(N 3 )-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)- (10) to the case where one additional field degree of freedom with O(N 4 ) 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 N i , 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.
The new FP collides with the DFP at N i = 10 and becomes unstable. Moreover, our results indicate that no FP is stable, e.g., at the point N 1 = 1, N 2 = N 3 = N 4 = 8. Together with the results in Table 2, this suggests that a structurally similar picture to the three-field case carries over to the fourfield case: The IFP will be stable for very small values of the N i , before it is destabilized. Keeping N 1 = 1 fixed and increasing N 2 = N 3 = N 4 , we pass through a regime without a stable FP, i.e., with nonuniversal behavior only. At N i = 9, 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 N i . Based on our findings in the three-field model, we expect that our O( ) estimates for the N i , 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. Table 2 Stable, fully coupled FP in the four-field case. There are four relevant directions; the corresponding exponents θ 1 , …, θ 4 are not provided

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 N i , 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., η i ≡ −k∂ ln Z i /∂k = 0 (see Eq. (1) for the definition of the parameters and couplings in the given model). The LPAtype truncation of the functional RG provides a quantitatively more precise estimate than the one obtained from the oneloop -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., η i = 0 (which we will refer to as LPA 4 + η in the following), we confirm the qualitative behavior of theexpansion: Fixing N 1 = 1 and increasing N ≡ N 2 = N 3 , the IFP becomes unstable around N 1.85. 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 N 2.8. Finally, at N 3 the asymmetrically coupled FP exchanges its stability with the DFP, which remains IR stable for all N 3. 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 We find the same picture as in the -expansion: The ACFP is IR stable in a small region 2.8 N 3 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 ( = 1): 8 ≤ N ≤ 10 (cf. Fig. 2) 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 4 + η 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.
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 d = 4, 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 N i while minimizing the computational effort.

Conclusions
With this study we identify a new fully coupled FP in the d = 3 dimensional three-and four-field models. While we find that this FP is indeed IR stable for some values of N i , it does not lie at real fixed-point values for the couplings at other values of the N i . This is in stark contrast to the O(N ) theory or the class of O(N 1 ) ⊕ O(N 2 ) models [6][7][8][9][10][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][7][8][9][10][11] different FPs exchange stability only through a collision of two fixed points at real values of N i , 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 d = 4 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][36][37] to determine the exponent θ 4 from critical exponents of the O(N ) Wilson-Fisher FP [1][2][3][4][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 2 < d < 4 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 N i 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 N i , 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 d = 4 dimensions. We observe that the asymmetrically coupled FP can be found in the -expansion, i.e., it emerges from the Gaussian FP at d < 4. Thus we conclude that no nontrivial FP exists for these models in d = 4, 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 (n ≥ 2) 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 d > 4. Following the methods discussed in [40][41][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.     6 We show the IFP (blue inverted triangles), the DIFP (purple filled circles), and the coupled isotropic BIFP (orange multiple sign) as a function of N 2 = N 3 for fixed N 1 = 1. The BIFP collides with the IFP at N 2 = N 3 = 1.5 and renders the IFP unstable without becoming stable itself. At larger values of N 2 = N 3 it collides with the DIFP. The coupled BIFP is shown only between the collision points to illustrate its role as a mediator between the IFP and DIFP their critical exponents and selected fixed-point coordinates in Fig. 5. The collision that destabilizes the IFP leaves no stable FP. The coupled IFP then moves on to collide with the DIFP which features a similar symmetry enhancement to an O(N 1 + N 2 ) ⊕ O(N 3 ) symmetry, cf. Fig. 6. The existence of new, symmetry-enhanced FPs is thus responsible for the early destabilization of the IFP (Fig. 6).