Minimizers of the dynamical Boulatov model

We study the Euler–Lagrange equation of the dynamical Boulatov model which is a simplicial model for 3d Euclidean quantum gravity augmented by a Laplace–Beltrami operator. We provide all its solutions on the space of left and right invariant functions that render the interaction of the model an equilateral tetrahedron. Surprisingly, for a non-linear equation of motion, the solution space forms a vector space. This space distinguishes three classes of solutions: saddle points, global and local minima of the action. Our analysis shows that there exists one parameter region of coupling constants for which the action admits degenerate global minima.


Introduction
In three dimensions, general relativity can be formulated as a B F-theory [1]. Its functional integral quantization discretized over simplicial complexes leads to the Ponzano-Regge model [2,3], which can be regarded as a quantum gravity model of discrete geometry. A cornerstone of this approach then is to recover continuum geometries with all desired requirements and properties of a three-dimensional spacetime. Such a description, however, as well as a mechanism, which could successfully sort it, remains an open problem in background independent approaches to quantum gravity.
The Boulatov model of group field theory (GFT) [4,5] provides one way to address this issue for Euclidean signature. The model is formally defined by the generating functional, a e-mail: jobengeloun@lipn.univ-paris13.fr b e-mail: kegeles@aei.mpg.de c e-mail: andreas.pithis@kcl.ac.uk where S(ϕ) denotes the Boulatov action [6]. The striking fact about this generating functional is that its Feynman graphs correspond to simplicial complexes and its Feynman amplitudes coincide with Ponzano-Regge spin foam amplitudes [2,3]. One concludes that a perturbative expansion of Eq.
(1) provides a discrete model of quantum gravity and that a description of continuum geometries will require a non-perturbative understanding of Eq. (1). The construction of a full non-perturbative quantum field theory is rarely possible, but often it is already enough to construct a perturbation theory around a non-perturbative vacuum [7]. Moreover, if quantum fluctuations are not too strong, a non-perturbative vacuum can be reasonably well approximated by the minimum of the classical action S, called the minimizer. In that case, the mean-field approximation around the minimizer prompts an effective field theory that will capture the non-perturbative regime of the model. For that reason, a study of minimizers of the Boulatov action is an important step towards a better understanding of continuous quantum geometries.
Despite their importance, however, the extrema of the Boulatov action are poorly understood in the literature. This is mostly because the Euler-Lagrange equations of the Boulatov action are non-linear differential equations that also involve integrals. Such equations are called integrodifferential equations; generally, they are notoriously difficult to solve. In the Boulatov model, these integro-differential equations can be formulated as integral equations with an integral kernel given by the Wigner 6J -symbol. A solution of the extremal equations then requires full control of the zeros of the 6J -symbol, which remains an open problem despite decades of research [8][9][10][11]. This makes the complete analysis of the problem out of reach.
In addition to this, there seems to be no consensus on the signs of the coupling constants in GFT models. For instance, the convention used in renormalization analyses [12] is opposite to the one used in the context of the GFT condensate cosmology investigations [13][14][15][16][17][18][19][20][21]. Despite this ambiguity in the sign convention both analyses rely on the existence of global or at least local minimizers and for that reason require a good understanding of the extrema in GFT.
In this work, we address the minimizers of the Boulatov action augmented by a Laplace-Beltrami operator, hereafter called dynamical Boulatov action [22]. To make the problem tractable, we look for minimizers in the space of left and right invariant fields corresponding to equilateral triangles. Section 2 gives the definition of the model and the space of functions considered in this article. On this space the Euler-Lagrange equations of the action become solvable, allowing us to provide a full characterization of solutions in Sect. 3.1. We then identify the parameter regimes in which the action admits minima and characterize the minimizers in Sect. 3.2. Our main result regarding the extrema is presented in Theorem 1 and the subsequent discussion. The characterization of minimizers is provided in Theorem 2. Implications of our results on the quantum theory are discussed in Sect. 4. A closing appendix gathers useful identities and the proofs of some statements in the text.

The dynamical Boulatov model
This section reviews the construction of the Boulatov model and sets up our notations. We assume the reader to be familiar with the harmonic analysis on SU(2) (Appendix A however, reports useful notions on this topic).
Let C ∞ (M) be the space of smooth, real-valued functions defined on the compact Lie group M = SU(2) ×3 . The components of elements of M are denoted by a subindex such that Define the space S of right and cyclic invariant functions. That is functions f in C ∞ (M) that satisfy right invariance: for any R ∈ SU(2) and any x ∈ M, The dynamical Boulatov action is a functional S m,λ on S, given by the integral where m 2 and λ are real, possibly negative, coupling constants, dx is the Haar measure on M, is the Laplace-

Fig. 1 Combinatorics of a tetrahedron
Beltrami operator on M with the canonical metric, 1 and the integral kernel Tet is given by This kernel encodes the combinatorics of a tetrahedron ( Fig. 1) and is symmetric under cyclic permutations of its arguments.
To address the variational problem, we topologize S by the family of semi-norms f n . = sup x∈M | n f (x)|, with the neighborhood base given by semi-balls [23], N ,n (0) = { f n < | f ∈ S}, for n ∈ N and > 0.
Leading the analysis further, we will restrict the space S by requiring (1) left invariance: for any L ∈ SU(2), . By the Peter-Weyl theorem, every f ∈ S can then be written as , Cycl denotes cyclic permutations of the set {1,2,3} and χ j i denotes the character of an SU(2) representation of dimension d j = 2 j + 1 for j ∈ N 2 . By [23, theorem 3] the sequence of coefficients ( f J ) J ∈J is a rapidly decreasing sequence of real numbers and the equality is understood such that the right hand side of the equation converges to f (x) in the aforementioned topology.
Furthermore, for such functions we require (2) the equilateral condition: let f ∈ S, then f is an equilateral function Peter-Weyl mode Mode with left invariance Mode with left and right invariance We denote the restriction of S to left invariant equilateral functions by S EL and the space of equilateral triples by J EL = {( j, j, j) | j ∈ N}. Note, that J EL contains only integer multiindices, since for any half-integer j the matrix coefficients vanish: In the following we will sometimes use the notation f ∈ S (EL) and f J with J ∈ J (E L) to signal that the statement holds equally for S and S EL and, correspondingly, with a set of indices belonging to J or to J EL . For clearer notation we also define the square of the triple J as J 2 . = j 1 ( j 1 + 1) + j 2 ( j 2 +1)+ j 3 ( j 3 +1), and its modulus as |J | . = j 1 + j 2 + j 3 .

Definition 1
A local minimizer of the action S m,λ on S EL is a field ϕ ∈ S EL , that for some n ∈ N and > 0 satisfies for any φ ∈ N ,n (ϕ) ∩ S EL . If condition (5) is satisfied on the whole space S EL we call the minimizer global.
In the following we will characterize all minimizers of the action S m,λ on S EL for the four different parameter regions For each of the parameter regions, we will characterize all extrema of the action S m,λ on S EL and identity, which of the extrema are minimizers.
We now briefly motivate the restrictions made in our analysis and point out the geometrical considerations behind the use of the space S EL .
The space S By the Peter-Weyl theorem we can decompose any smooth field f on M in modes such that 3 , are the Wigner-matrix coefficients for the product representation of M. To gain intuition on the construction, we depict the Peter-Weyl coefficients by stranded lines, emanating from a single point (Fig. 2a). Then the right invariance of f ensures a closure of the dual edges to form a triangle (Fig. 2b). Hence, the right invariance is necessary to give a geometric interpretation to the fields and it is thus crucial for the connection between the Boulatov group field theory and the Ponzano-Regge spin-foam model [1,2,6]. In addition, the invariance under cyclic relabeling of the field arguments ensures that the ordering of the field arguments has no physical meaning. 2 The space S EL To enforce rotational symmetry of the triangles, one requires left invariance of the field [15,16], such that for any h ∈ SU(2) the field f satisfies f (hx 1 (Fig. 2d). Applications of GFT to quantum cosmology demonstrate that this symmetry is needed to identify the domain space of the fields with the superspace of homogeneous spatial geometries [15].
The equilateral condition on functions ensures that their modes correspond to equilateral triangles. This condition can relate to isotropy in quantum cosmology studies of GFT, reflecting that we need to set all edges of the triangle to equal length in order to ensure equality in all directions (Fig. 2d).
It is crucial for the recovery of a Friedmann-like dynamics from GFT [16,17,20,21]. 3,4 Besides these arguments, there is also an algebraic reason to consider the restricted space S EL . In GFT the action S m,λ defines statistical weights of a generating functional using a functional integral (1). It has been shown, however, that on S the action S m,λ is not bounded from below, regardless of the parameter region [32,33]. For this reason, the above integral is dominated by those field configurations that make the action S m,λ arbitrarily negative making Eq. (1) ill-defined. As we will show below, this problem gets resolved on S EL , where global minimizers of the action exist (at least for some parameter regions). This allows us to define the generating functional perturbatively, and could lead to a well-defined statistical theory. 5

Extrema and minimizers
Let I ⊂ R denote an interval containing zero; for t ∈ I and ϕ, f ∈ S EL a necessary condition for ϕ to be a local minimizer on S EL is given by for any f ∈ S EL . In the following we will investigate the extremal condition (6) for the model (2). We will then check if some solutions are minimal and thus fulfill (7) and the condition in Definition 1.

Proposition 1 ϕ ∈ S (EL) is an extremum of S if and only if
the Peter-Weyl coefficients of ϕ -denoted ϕ J -satisfy for any J ∈ J (EL) , 3 Notice that for the subsequent analysis of extrema on S EL the cyclicity property mentioned above has no impact and could in principle be lifted from the outset. 4 The restriction to equilateral configurations bears strong resemblance to what is done in the closely related contexts of dynamical triangulations [28,29] and tensor models for quantum gravity [30,31] where the use of standardized building blocks -by universality arguments -is believed not to affect the continuum results. 5 One can bound the Boulatov action by adding a so-called pillow term to the action [32]. We leave the impact of such a modification onto the ensuing analysis to future investigations.
Proof See Appendix B1.
The extremal condition (8) is a non-linear tensor equation with an integral kernel given by the 6J -symbol squared. To this issue adds the fact that the non-trivial zeros of the 6Jsymbol are still under investigation, making (8) inherently difficult to solve in full generality. Some specific solutions for the case without the Laplace-Beltrami operator and λ < 0 have been introduced in Ref. [34], however, a systematic analysis of extrema was not performed therein.
Although the extremal condition (8) is difficult to solve on S, it turns out to be solvable on S EL , because in this case the 6J -symbol significantly simplifies.

Extrema
In the following we will denote the Wigner 6J -symbol for In order to characterize the extrema of the action, we define the space of extremal sequences. Let C = (C J ) J ∈J (EL) denote the sequence of (possibly complex) numbers such that for and for J ∈ J (EL) /J S (EL) Since J 2 > 0, the first case in (10) can happen only when m 2 is negative and for J ∈ J EL , m 2 has to be an even integer. For simplicity, we will exclude this case in the following analysis, because it requires a fine-tuning on the parameter m 2 . It is convenient to define the length of the sequence C such that with the convention sgn(0) = 0.

Definition 2
We define the space of extremal sequences as where the coefficients of each sequence are of the form (9).
This space of course depends on the values of m 2 and λ, since different choices of these parameters may violate the reality condition C J ∈ R. E m,λ fully characterizes the space of extrema of the action as states the following theorem.
is an extremum of the action S m,λ . Moreover, every equilateral extremum of S m,λ is of the above form.
Proof See Appendix B2.
We denote the space of extremal functions byẼ m,λ . It is worth mentioning that, in spite of the non-linearity of the Euler-Lagrange equations, its solutions form a vector space over (Z 3 , +, ·).
Corollary The spaceẼ m,λ is a vector space over the discrete algebraic field (Z 3 , +, ·).
Proof Denote the space of sequences with finitely many nonzero elements over Z 3 by c 00 (Z 3 ). It is a vector space over Z 3 . Consider the map with the convention sgn(0) = 0. I is one-to-one on its image, however, it may not be onto c 00 (Z 3 ) simply because the non-trivial zeros of the 6J -symbol are not fully characterized. Nevertheless, the image of I is algebraically closed and forms a subspace of c 00 (Z 3 ). For any s = (s 0 , s 1 , . . .) ∈ I(Ẽ m,λ ), the inverse mapping is given by Since there are only finitely many non-zero coefficients, s j = 0, the sum trivially converges in S EL . Since I is linear it is an isomorphism betweenẼ m,λ and I(c 00 (Z 3 )).
We define the sum onẼ m,λ by We now discuss the space of extremal sequences according to different parameter regions, whose major difference is captured by the sign of the radicand in (9). We obtain the four cases: (a) m 2 < 0, λ < 0 : the radicand is positive only if which is the case when j satisfies where · denotes the ceiling function. The space of extremal sequences contains infinitely many sequences of the form 0, . . . , 0, C J min , C J min +1 , . . . , where we used the notation J min + n . = ( j min + n, j min + n, j min + n) for n ∈ N, with finitely many non-zero elements C J .
or equivalently for j satisfying, where · denotes the floor function. In this case E m,λ contains finitely many sequences of the form At this point, a few comments are in order: according to the geometrical interpretation in the previous section, each Fourier mode represents a triangle with the edge length j and the area proportional to J 2 . In the parameter regime (d) relation (17) provides an upper bound on the possible j's for the extrema of the action. Hence, in this case, |m 2 | can be interpreted as the bound on the area of the triangles determined by the extremal solutions. This is an interesting geometrical fact that deserves further investigation.
A second remark is that the method of resolution restricting to equilateral configurations used to tackle (8) certainly exports to GFT models on higher dimensional manifolds M = G ×D with G = SU(2), SO(4) and D ∈ N. We expect that a similar result as in (9) will hold if we replace the 6Jsymbol by the appropriate Wigner symbol and replace the square root by the D − 2 root. However, the search of minimizers for these theories as performed in the subsequent analysis might be different.

Minimizers
We now seek the minimizers of the action and show that only two parameter regions admit global minimizers.
First, notice that in the case, m 2 < 0, λ > 0, the value of |m 2 | can determine, whether or not the action S m,λ is bounded from below. To agree with this, assume the first non-trivial zero of the 6J -symbol to be at J 0 ∈ J EL and choose a function Hence, the action can become arbitrarily negative and thus is unbounded from below. On the other hand, for |m 2 | < J 2 0 the action has a global minimum as we will show in the following.
In order to give a general classification of solutions, we need to exclude cases when the 6J -symbol vanishes. A quick numerical analysis shows that for |m 2 | ≤ 10 9 , the space of non-trivial zeros of the 6J -symbol with J 2 ≤ |m 2 | is empty, Therefore, Theorem 2 captures all possible solutions up to this order. In fact, we conjecture that for equilateral configurations, J EL /J S EL = ∅, and our theorem holds for any value of |m 2 |. Theorem 2 Let |m 2 | be such that for j ≤ j max every J ∈ J S EL and such that there is no J ∈ J EL /J S EL such that J 2 − |m 2 | = 0. Then the equilateral extrema of the dynamical Boulatov action are of the following type: Proof of Theorem 2 In the following, let ϕ(x) denote an extremum and let f ∈ S EL be a generic function with the Peter-Weyl decomposition given by f (x) = J ∈J EL f J X J (x). We remind here that a necessary condition for an extremum ϕ(x) to be a minimizer (maximizer, respectively) is given by for any f ∈ S E L . In the Peter-Weyl decomposition the second variation recasts as where ϕ K is the Peter-Weyl coefficient of the extremum ϕ. The above condition is necessary but not sufficient. Nevertheless, it turns out to be useful to exclude some extrema. Case (a) (m 2 ≤ 0, λ ≤ 0): By Theorem 1, extremal solutions contain only finitely many non-zero Fourier coefficients. Therefore it is possible to find J > ∈ J E L such that which violates the maximizer condition. To see that the minimizer condition is also violated, choose f < (x) . = f J < X J < (x) such that J 2 < − |m 2 | < 0. Then the second variation is written as Hence, each extremum in this parameter region violates the minimizer and the maximizer condition and therefore is a saddle point. Case (b) (m 2 ≥ 0, λ ≤ 0): For the non-trivial minimizer the above argument can also be applied in this case. Choosing the functions f > (x) and f < (x) as above we find Hence, non-trivial extrema are saddle points. For the trivial extremum the second variation of S m,λ reads for any f ∈ S EL and the necessary condition is satisfied. Indeed, the trivial extremum is a local minimum. To prove this we first notice that the Peter-Weyl transform is a topological isomorphism from S EL to the space of rapidly decreasing sequences S(N) with topology given by the family of semi-norms [23, theorem 4], The action evaluated at f becomes Since the Wigner-6J -symbol is upper-bounded by 1, we can estimate Since Peter-Weyl transform is a topological isomorphism, we get for any f ∈ S EL with f 0 ≤ 4!m 2 |λ| , an estimate on the Fourier coefficients Inserting this bound in (27) we have for any f ∈ S EL Hence, S m,λ (0) = 0 ≤ S m,λ ( f ), ∀ f ∈ S EL . We obtain a global minimizer, since the minimal condition is satisfied on the whole S EL .
Case (d) (m 2 < 0, λ > 0): For any f ∈ S EL the action evaluated at f gives Splitting Hence, verifying the minimizer condition, it is enough to show that The space of functions of the form f − is finite-dimensional and we can use the usual minimization procedure for functions. More specifically, let s J : R → R be a function such that The action S m,λ ( f − ) is smallest when each s J is minimal on R for each J ≤ J max . Taking the first and second derivative of s J we see that the minimum is achieved by the coefficients C J from (9). Hence, an extremum given by an extremal sequence of maximal length is a global minimizer on the whole S EL . If ϕ is given by an extremal sequence C of length (C) < j max , then there exists a X J 0 with J 0 ≤ J max and ϕ J 0 = 0. For δ ∈ R define the function g(x) = ϕ(x) + δ · X J 0 (x). Inserting g into the action we get If δ 2 is in the range 0 < δ < 2C J 0 the square bracket is negative and it follows S m,λ (g) ≤ S m,λ (ϕ). Moreover, for any > 0 and δ < J 2n 0 we have since the characters are bounded by one, |X J 0 (x)| ≤ 1. Hence, g ∈ N ,n (ϕ). For any > 0 choosing δ < min J 2n 0 , C J 0 we get S m,λ ( f ) < S m,λ (ϕ). This shows that we can find a function g in any neighborhood of ϕ that decreases the value of the action, and hence, ϕ is not a minimizer.

Concluding remarks
We investigated the minimizers of the dynamical Boulatov action in four different parameter regions of the coupling constants. Our analysis is restricted to the space of smooth, equilateral, left and right invariant functions, also invariant under cyclic permutations of its variables, S EL . This restriction ensures that the action is bounded from below for some parameter regions. Moreover, it is motivated by quantum cosmology studies on GFT. It appears that the very same restrictions allow us to solve the Euler-Lagrange equations for the dynamical Boulatov action and lead to a complete characterization of minimizers on the restricted space. Our result characterizes the space of solutions by extremal sequences of finite length and shows that it forms a vector space over Z 3 , which is surprising for the set of solutions to a nonlinear integro-differential equation. Furthermore, in the most interesting parameter region (d), the non-vanishing Fourier modes of extremal solutions are bounded by the coupling constant m 2 , which suggests a connection between m 2 and the area of the triangle of the largest Peter-Weyl mode of the GFT field.
Our analysis shows that the region (a) does not have any minimizers on S EL , which makes this parameter region perhaps the least suitable for the definition of the statistical measure in (1). For the parameter regions (b) and (c) there is a single (local respectively global) minimizer given by the trivial extremum, ϕ = 0. Finally, in the region (d) the action has 2 j max degenerate global minimizers, where j max is a function of the coupling constant m 2 . The rich structure of global minima makes this region most interesting for further investigations, especially for the statistical theory.
On the space of equilateral functions only two possible parameter regions (c) and (d) allow for the presence of global minimizers and hence could induce a meaningful definition of a non-perturbative statistical measure.
Case (c) admits a single global minimizer ϕ = 0. Perturbation theory around this minimizer defines the perturbation theory in the coupling constant λ and is used in the GFT literature to draw a connection to spin-foam models. Hence, our analysis would suggest that this regime is suitable for such a relation.
Case (d), on the other hand, may suggest more structure for the quantum theory: a degenerate global minimum could lead to instantons or symmetry breaking in the corresponding statistical field theory in the following sense: Instantons The full non-perturbative formulation of a model is given by the minimizer of its quantum effective action. The latter is commonly assumed to be convex [35] and therefore admits a single, unique minimizer. Hence, the difference between the minimizers of the classical and the quantum effective action becomes apparent, especially in the case when the classical action admits degenerate minimiz-ers. In this case, a perturbative description around any of the minimizers of the classical action does not capture the non-perturbative effects of the theory. In quantum field theory, these non-perturbative effects can be understood as "tunneling" between the perturbative vacua, where the instanton action describes the tunneling probability. Thus, the degenerate structure of global minimizers in our case, suggests the necessity of instantons in the statistical formulation of GFT at least for the parameter region (d) (for a similar result see Ref. [19]).
Symmetry breaking This mechanism happens when the classical action admits degenerate global minimizersrelated by a symmetry of the classical action -but the tunneling probability between them vanishes. As we already mentioned, the tunneling probability is described by the instanton action, which in ordinary field theory is often proportional to the volume of the base manifold. On a manifold with a finite volume, the tunneling probability is therefore finite. This often pertains to the statement that spontaneous symmetry breaking cannot occur in quantum field theories in a box. This realization, however, contains further assumptions that are satisfied in ordinary field theories but do not hold for GFT. It has been recently shown that even on the compact base manifold, M = SU (2) d the tunneling between different perturbative minima can vanish [36], leading to a similar phenomenon of symmetry breaking. In order to talk about symmetry breaking, we need to identify the symmetry, which in our case, is given by a flip of the sign of at least one of the modes in the Peter-Weyl decomposition of the minimizer (this can be modeled as a Z 2 -symmetry). Since the action is of even power in the fields, such a flip will not affect the value of the action and will correspond to a discrete symmetry. For this reason, it is possible that the global minimizers of the action provoke the breaking of sign-flip symmetry. This needs to be investigated more rigorously in future work.
For ordinary local quantum field theories, a symmetry breaking mechanism can sometimes be related to a phase transition and the formation of a condensate. In particular, this could be the signal of a Bose-Einstein condensation just as expected for quantum cosmology studies in GFT. A closer look at the solutions found for sector (d) shows that these might bear intriguing perspectives. Indeed, the 'particle' number, used in the condensate cosmology context, is computable in terms of the L 2 -norm of the minimizer. In the present situation, that very number proves to be bounded by the parameter m 2 : . For |m 2 | 1 we can approximate j max further as j max ≤ 2|m 2 | and obtain a simpler bound on the L 2 -norm of the minimizers The coupling constant m 2 (or |m 4 |/λ 1) could be large but that itself is not enough to ensure N = ∞. Nevertheless, starting from our solutions, a divergent parameter m 2 is a necessary condition for the divergent L 2 -norm. A large particle number would be desirable for the condensate cosmology approach because such configurations could then be interpreted as to define non-trivial homogeneous and isotropic background geometries in 3d with Euclidean signature. This point deserves further investigations.
We should mention here that our analysis does not capture minimizers with a divergent L 2 -norm (dealing only with integrable functions), and some modifications will be in order to consider these cases. One necessary modification would be to relax the smoothness condition of the minimizers and use the space of tempered distributions instead. This could be particularly interesting for GFT models without the Laplace-Beltrami operator, which correspond to a topological BFtheory. Due to the distributional nature of minimizers their L 2 -norm will sometimes diverge making them potentially interesting for quantum cosmological studies [19,36] and spin-foam models [34]. The solutions to these GFT models must be addressed differently but certainly deserve further attention.
There are several models using tensor fields (with interesting properties such as perturbative renormalizability) which do not impose strong symmetry conditions on the fields. These models' interactions could also be radically different from that of Boulatov [37]. Their corresponding Euler-Lagrange equation (without 6J -symbols) still involves a nonlinear tensor like equation, and it remains a difficult task to solve them. In this case, an approach to circumvent the nonlinearity and to obtain solution fields which are more general than equilateral configurations is to consider symmetric tensor fields and to decompose the field into its traceless part and the rest, namely vector-like components [38]. Such a decomposition could help to solve the extremal conditions on S which might find applications in GFT studies of inhomogeneous and anisotropic quantum cosmologies.
On the other hand, the existence of global minima on S EL suggests that we can define a self-consistent statistical theory using only this space. This theory could potentially be welldefined due to the bound of the action on S EL and may have implications for quantum cosmological studies of GFT.
We briefly recall the most important properties of the Peter-Weyl transform and Wigner matrices, needed for the harmonic analysis on SU(2). Let C ∞ (SU (2)) be the space of smooth functions f on SU(2) which is equipped with the topology given by semi-norms with the Laplace-Beltrami operator and n ∈ N.
If we call the space of rapidly decreasing sequences S(N), then (A3) defines a family of semi-norms on S(N) and in the corresponding topology it becomes a Fréchet space. Then the Peter-Weyl transform F : C ∞ (SU(2)) → S(N) is a topological isomorphism between the space of smooth functions and the space of rapidly decreasing sequences [23]. In our work, we deal with functions on three copies of SU (2). For this reason, we introduce M = SU(2) ×3 as a Lie For the "if" direction we observe the following: since the set {X J } J ∈J (EL) is dense in S (EL) , for any f ∈ S (EL) there exists a family of real numbers { f J } J ∈J (EL) such that the sequence of functions given for all N ∈ N as converges to f . Then c = sup x∈M sup N ∈N | f N (x)|, exists and dominates each f N such that, | f N | ≤ c. Moreover, c, seen as a constant function on M, is integrable since M is compact.
For any f ∈ S (EL) the extremal condition for the action S m,λ reads as Using the Peter-Weyl decomposition for f , we can interchange the limit and the integral by the dominant convergence theorem (using the bound c) and obtain for any f ∈ S (E L) , from which the statement follows.

Corollary ϕ ∈ S (EL) is an extremum of S if and only if the
Peter-Weyl coefficients of ϕ -denoted by ϕ J -satisfy for any J ∈ J (EL) , Proof From Lemma 1 the extremal condition is given by the variation in the basis direction X J for any J ∈ J (EL) . Inserting the Peter-Weyl decomposition of ϕ in the action S m,λ (ϕ), interchanging the limit with the integral by the dominant convergence theorem and using the relation in (A10) we obtain the desired statement.

B.2. Proof of Theorem 1
Theorem For any C ∈ E m,λ the field ϕ ∈ S EL is an extremum of the action S m,λ . Moreover, every equilateral extremum of S m,λ in S EL is of the above form.
Proof To show that ϕ solves the extremal condition we need to show, by Proposition 1, that each C J satisfies (8), which follows by direct calculation. Conversely, every equilateral function can be written as with (A J ) J ∈J EL being a rapidly decreasing sequence [23]. Using Proposition 1, we find that the extremal solutions have coefficients A J which satisfy or A J ∈ R for J ∈ J EL /J S EL with J 2 + m 2 = 0. If A J is not trivial we can estimate its growth using the asymptotic behavior of 6J -symbols [40] as However, for (A J ) J ∈J EL to be a rapidly decreasing sequence, the coefficients have to satisfy for any n ∈ N, lim j→∞ | j| n |A J | → 0. (B9)