Luca Lionni
email: luca.lionni@th.u-psud.fr

Johannes Thürigen
email: johannes.thuerigen@physik.hu-berlin.de

Multi-critical behaviour of 4-dimensional tensor models up to order 6

Tensor models generalize the matrix-model approach to 2-dimensional quantum gravity to higher dimensions. Some models allowing a 1/N expansion have been explored, most of them generating branched-polymer geometries. Recently, enhancements yielding an additional 2d pure gravity (planar) phase and an intermediate regime of proliferating baby universes have been found. It remains an open issue to find models escaping these lower-dimensionality universality classes.

Here we analyze the dominant regime and critical behaviour of a range of new models which are candidates for such effective geometries, in particular interactions based on the utility graph K3,3. We find that, upon proper enhancement, the two-phase structure of a branched-polymer and a 2d gravity regime is the common case in U (N )-invariant rank D = 4 tensor models of small orders. Not only the well known so-called necklace interactions but also K3,3-type interactions turn out as the source for the planar regime. We give a systematic account of the enhancement scaling, the counting of leading-order diagrams and the multi-critical behaviour of a wide range of interactions, in particular for all order-6 interactions of rank 3 and 4. These findings support the claim of universality of such mixtures of branched-polymer and planar diagrams at criticality. In particular, this hints at the necessity to consider new ingredients, or interactions of higher order and rank, in order to obtain higher dimensional continuum geometry from tensor models.

INTRODUCTION

Tensor models [START_REF] Gurau | Random Tensors[END_REF][START_REF] Gurau | Invitation to Random Tensors[END_REF][START_REF] Gurau | Colored Tensor Models -a Review[END_REF], which generalize matrix models [START_REF] Di Francesco | 2D gravity and random matrices[END_REF] to higher dimensions, were introduced as a non-perturbative approach to quantum gravity, and an analytical tool to explore random geometries in dimension higher than two. These geometries are generated as Feynman diagrams in the perturbative expansion of a tensor path integral, and the amplitudes associated to discrete geometries are discretizations of the quantum path integral of general relativity. Restricting the tensor action to be U (N ) D -invariant such that interactions can be represented as D-coloured graphs [START_REF] Bonzom | Random tensor models in the large N limit: Uncoloring the colored tensor models[END_REF] allows for a 1/N expansion [START_REF] Gurau | The 1/N expansion of colored tensor models in arbitrary dimension[END_REF][START_REF] Gurau | The Complete 1/N Expansion of Colored Tensor Models in Arbitrary Dimension[END_REF]. A continuum limit can be obtained by tuning correlation functions to criticality in the coupling parameters [START_REF] Bonzom | Random tensor models in the large N limit: Uncoloring the colored tensor models[END_REF][START_REF] Bonzom | Critical behavior of colored tensor models in the large N limit[END_REF][START_REF] Bonzom | Multi-critical tensor models and hard dimers on spherical random lattices[END_REF].

So far, only continuous geometries which come as limits of one or two dimensional structures have been identified. Tensors of rank D = 1, i.e. vector models, generate plane trees in the large-N limit, which converge towards the continuous random tree [START_REF] Aldous | The Continuum Random Tree. I[END_REF][START_REF] Aldous | The continuum random tree. II. An overview[END_REF][START_REF] Aldous | The Continuum Random Tree III[END_REF] in the continuum limit, also called branched polymers, and are characterized by the entropy critical exponent γ = 1/2, Hausdorff dimension d h = 2 and spectral dimension d S = 4/3. In the D = 2 case of matrix models, amplitudes scale as N 2-2g leading to planar geometries (genus g = 0) in the large-N limit, which converge at criticality towards the Brownian sphere [START_REF] Marckert | Limit of normalized quadrangulations: The Brownian map[END_REF][START_REF] Gall | Scaling limits of bipartite planar maps are homeomorphic to the 2-sphere[END_REF][START_REF] Gall | The topological structure of scaling limits of large planar maps[END_REF], with γ = -1/2, d H = 4, and d S = 2, which has recently been proven to be equivalent to 2d Liouville quantum gravity [START_REF] Polyakov | Quantum geometry of bosonic strings[END_REF][START_REF] Miller | Liouville quantum gravity and the Brownian map I[END_REF][START_REF] Miller | Liouville quantum gravity and the Brownian map II: geodesics and continuity of the embedding[END_REF][START_REF] Miller | Liouville quantum gravity and the Brownian map III: the conformal structure is determined[END_REF][START_REF] David | Liouville Quantum Gravity on the Riemann Sphere[END_REF]. Tensor models of rank D ≥ 3 are generically dominated in the large-N limit by so-called melonic diagrams [START_REF] Bonzom | Critical behavior of colored tensor models in the large N limit[END_REF] encoding discrete D-dimensional spherical manifolds. Surprisingly, in the continuum limit, they converge towards branched polymers [START_REF] Gurau | Melons are branched polymers[END_REF]. Nevertheless, it has recently been shown that non-melonic interactions coined "necklaces" can be enhanced to obtain continuum 2d geometry in this limit [START_REF] Bonzom | Enhancing non-melonic triangulations: A tensor model mixing melonic and planar maps[END_REF] using an intermediate-field (IF) representation [START_REF] Gurau | The 1/ N Expansion of Tensor Models Beyond Perturbation Theory[END_REF][START_REF] Delepouve | Borel summability and the non perturbative 1/N expansion of arbitrary quartic tensor models[END_REF].

The question is thus whether it is possible to find an enhancement of certain interactions such that the large-N continuum limit has properties suited for a quantum space-time of dimension higher than two. Such new continuum geometries should have spectral dimension d S > 2 and are expected to be characterized also by new critical exponents γ. While there are no guiding lines from mathematics where no such geometries are known so far, tensor models generating D-dimensional discrete topologies could, in principle, contain such new universality classes in the continuum limit.

The motivation of the current work is that there are classes of tensor models not investigated before which, on first sight, have the potential to describe such new universality classes of higher-dimensional continuum geometries. The naive idea is that such geometries might arise when there are sufficiently large non-planarities, e.g. crossing edges on the diagrams such that the diagrams cannot be embedded on the 2-sphere. Intuitively, a space obtained as a limit of discrete random spaces with diverging number of holes competing with the diverging volume should have interesting fractal properties. Natural candidates are non-local interactions with a combinatorial structure based on the utility graph K 3,3 . According to Kuratowski's theorem [START_REF] Kuratowski | Sur le problème des courbes gauches en Topologie[END_REF], finite graphs are not planar if and only if they contain (a subdivision of) either the utility graph K 3,3 or the complete graph K 5 . As bipartite graphs, U (N ) D -invariant tensor interactions cannot contain K 5 . It is thus natural to expect new quantum geometries to result from sufficiently dominant interactions which are extensions of K 3,3 .

To test this hypothesis we have explored the full space of sixth-order tensor-invariant interactions at rank D = 3 and 4. There are preliminary results that the K 3,3 -interaction for D = 3 leads only to small crossings [START_REF] Bonzom | Colored triangulations of arbitrary dimensions are stuffed Walsh maps[END_REF], which should not be dominant in the continuum limit, as we indeed prove in this paper. At D = 4, it is possible to have K 3,3 -type interactions which lead, on first sight, to large crossings. Nevertheless, studying the dominant diagrams and the critical behaviour in detail we find that these models are still in the universality class of branched-polymer and 2d-gravity geometries. We find an explanation for this in terms of a bijection between the Feynman diagrams and embedded diagrams, that is, graphs drawn on surfaces.

Our systematic study of non-melonic tensor models covers a number of interesting aspects. First of all, it is the first study of interactions of order higher than 4, apart from a specific rank-3 model of order 8 leading to branched polymers in the continuum at large N [START_REF] Bonzom | Counting gluings of octahedra[END_REF]. So far, the focus in the literature has been on quartic interactions [START_REF] Delepouve | Borel summability and the non perturbative 1/N expansion of arbitrary quartic tensor models[END_REF][START_REF] Dartois | Double scaling in tensor models with a quartic interaction[END_REF][START_REF] Delepouve | Phase transition in tensor models[END_REF][START_REF] Nguyen | An analysis of the intermediate field theory of t4 tensor model[END_REF][START_REF] Delepouve | Universality and Borel Summability of Arbitrary Quartic Tensor Models[END_REF][START_REF] Delepouve | Quartic Tensor Models[END_REF] and straightforward extensions. We are able to extend these studies to higher order using the recent generalization of intermediate-field representation to arbitrary order [START_REF] Bonzom | Colored triangulations of arbitrary dimensions are stuffed Walsh maps[END_REF]. This method relies on a bijection between stranded Feynman diagrams and embedded diagrams. It is a powerful tool for finding enhancements of less dominant interactions and for identifying the resulting diagrams at leading order in the 1/N expansion. We provide a systematic analysis of the enhancement scalings of all sextic interactions.

The bijection to simpler embedded diagrams provides means to determine amplitudes from counting diagrams. Tensor models, just like matrix models, are purely combinatorial in that amplitudes in the perturbative expansion are given, up to factors in couplings and N , simply by the number of diagrams contributing at that order. Accordingly, leading-order correlation functions are generating functions in the combinatorics sense. In this way we obtain explicit Dyson-Schwinger-type equations for two-point functions and find solutions solving them. On these grounds we are able to determine the critical behaviour in many cases by performing explicitly the Newton-Puiseux expansion in fractional powers at critical loci. The critical behaviour of correlations, characterized by the entropy exponent γ, provides a first hint of what the continuum limit of the model is.

As a consequence, we are for the first time able to analyse models generating topologically non-trivial discrete geometries. Indeed, the tensor interactions are dual to D-dimensional building blocks, and the Feynman diagrams are dual to discrete spaces obtained by gluing such building blocks together. So far, only models involving interactions with boundaries of spherical topology, thus generating discrete geometries obtained from gluing balls, have been investigated. The interactions we are considering in this paper are dual to building blocks whose boundaries are spherical, toroidal, or may even be singular (pseudo)-manifolds. Still, we find that they all belong to the universality classes with continuum limits of branched-polymer or planar geometries.

That is, all models considered here belong to the universality classes of branched polymers (γ = 1/2), of 2d pure gravity (planar diagrams, with γ = -1/2), or to an intermediate transitional regime of "proliferating baby universes", with critical exponent γ = 1/3, through various combinations of enhanced interactions. This last regime was discovered in matrix models by considering modified Einstein-Hilbert actions with higher curvature terms, resulting into multi-trace interactions [START_REF] Das | New Critical Behavior in d = 0 Large-N Matrix Models[END_REF][START_REF] Korchemsky | Matrix model perturbed by higher order curvature terms[END_REF][START_REF] Alvarez-Gaumé | A Proposal for strings at D > 1[END_REF][START_REF] Ambjorn | Baby universes in 2-d quantum gravity[END_REF][START_REF] Gubser | Scaling functions for baby universes in two-dimensional quantum gravity[END_REF]. On the contrary, we recover these regimes as different phases of tensor models with connected interactions, corresponding to discrete theories of gravity without higher-order curvature terms. In combinatorics, these three regimes appear also in the context of (non-embedded) planar graphs with given 2-or 3-connected components [START_REF] Giménez | Graph classes with given 3-connected components: Asymptotic enumeration and random graphs[END_REF].

A peculiar feature of tensor models is that the branched-polymer-planar phase diagram can even be obtained from a single type of interaction. Indeed we show here that it occurs already for two quartic necklace interactions of different colouring, extending and completing previous work [START_REF] Bonzom | Enhancing non-melonic triangulations: A tensor model mixing melonic and planar maps[END_REF][START_REF] Bonzom | Large N limits in tensor models: Towards more universality classes of colored triangulations in dimension d ≥ 2[END_REF]. In the present model, the branched-polymer phase is not obtained from a planar regime by coupling it to a polymerization term, as done in [START_REF] Bonzom | Enhancing non-melonic triangulations: A tensor model mixing melonic and planar maps[END_REF][START_REF] Das | New Critical Behavior in d = 0 Large-N Matrix Models[END_REF][START_REF] Korchemsky | Matrix model perturbed by higher order curvature terms[END_REF][START_REF] Alvarez-Gaumé | A Proposal for strings at D > 1[END_REF][START_REF] Bonzom | Large N limits in tensor models: Towards more universality classes of colored triangulations in dimension d ≥ 2[END_REF], but rather as an interference between planar regimes, which instead of building up interesting behaviour, reduces to a branched polymer regime for a small region of the phase space. As expected, at the boundary of this critical region we observe the intermediate regime of proliferation of baby universes.

The multi-critical structure of the models turns out to be more involved than anticipated in the literature [START_REF] Bonzom | Enhancing non-melonic triangulations: A tensor model mixing melonic and planar maps[END_REF][START_REF] Das | New Critical Behavior in d = 0 Large-N Matrix Models[END_REF][START_REF] Alvarez-Gaumé | A Proposal for strings at D > 1[END_REF]. It is well known that with several couplings one can have sequences of multi-critical subspaces with critical exponents γ = 1/2, 2/3, 3/4, ... in the branched-polymer regime [START_REF] Ambjorn | Summing over all genera for d ¿ 1: a toy model[END_REF], and the Kazakov exponents γ = -1/2, -1/3, -1/4, ... in the planar regime [START_REF] Kazakov | The Appearance of Matter Fields from Quantum Fluctuations of 2D-GRAVITY[END_REF] which in the continuum comes to coupling non-unitary conformal matter to 2d gravity. If such multi-critical subspaces are incident with the baby-universe transition regime, special multi-critical values could be found. In the case of matrix models it is stated that such cases occur [START_REF] Korchemsky | Matrix model perturbed by higher order curvature terms[END_REF][START_REF] Alvarez-Gaumé | A Proposal for strings at D > 1[END_REF]. However, our explicit calculations of critical loci in the tensor models considered here indicate that branched-polymer-regime multi-critical subspaces are not incident with the transition regime. Instead, we find that the presence of couplings relevant for the planar regime can enrich also the multi-critical structure of the branched-polymer regime.

Apart from these new details in the multi-critical structure, our main result is the confirmation that order-6 tensor models of rank 3 and 4 belong still to the same universality classes as order-4 models. Initial hope to find new universality classes for K 3,3 -type interactions could not be fulfilled. Indeed, the bijection with embedded diagrams shows clearly that the influence of non-planarities in the interactions on the structure of leading order diagrams is not strong enough by itself to change their large scale properties, and therefore the corresponding continuum limit. This could be taken as an argument that U (N ) D -invariant tensor models do not cover new regimes.

Conversely, the absence of new regimes in the models considered might suggest that the structure of the embedded interaction, after the intermediate-field bijection, should be a more pertinent criterion than the non-planarity of its stranded version. Indeed, all order 4 and 6 tensor model diagrams of rank 3 and 4 allow for some bijection to twodimensional combinatorial structures. The non-planarity becomes a condition on colourings of edges of ribbon graphs which at best has no effect with respect to their usual planar large-N behaviour, or tends to prevent large loops from appearing instead of favouring them, as was initially hoped for. However, this argument is not applicable in general for higher-order or higher-rank models since the bijection maps interactions to "stuffed" embedded vertices with a non-trivial internal structure in general.

The structure of the paper is the following: In the first section, we introduce rank-D enhanced random tensor models, describe their relevance as discrete formulations of gravity in dimension D, and comment on known critical behaviour and continuum limits. We introduce some of the models we solve in the following. The exhaustive list of order-3 and order-4 U (N ) D -invariant interactions can be found in Tables I and II, as well as the scalings and critical exponents computed in this paper. Section II is dedicated to the characterization of leading order Feynman diagrams throughout the bijection introduced in [START_REF] Bonzom | Colored triangulations of arbitrary dimensions are stuffed Walsh maps[END_REF] with embedded diagrams. Finally, in Section III, we analyse the critical and multi-critical behaviour of the models.

I. TENSOR MODELS WITH BRANCHED-POLYMER-PLANAR-TRANSITION CRITICAL BEHAVIOUR

A. Tensor models

We consider tensor models as defined by a partition function with a finite number of U (N ) ⊗D -invariant interactions indexed by B in the set of interactions B defining the model,

Z B {t B }, N = C D exp N D-1 B∈B t B N s B Tr B (T, T ) dµ 0 (T, T ), (1)
where s B ∈ N are the scalings, and dµ 0 (T, T ) is the Gaussian measure

dµ 0 (T, T ) = 1 Z 0 e -N D-1 T • T N a1,...,a D =1 dT a1•••a D d Ta1•••a D , (2)
in which we denoted

T • T = a1,•••a D T a1•••a D Ta1•••a D . The partition function Z 0 is normalized such that Z B (0, N ) = 1.
More specifically, Tr B (T, T ) are polynomials of T and T invariant under unitary transformations in

U (N ) ⊗D = U (N ) ⊗ • • • ⊗ U (N )
. Therefore they must have the same number of copies of T and T , and indices may only be contracted between a T and a T . An example in D = 4 is the following, (3)

Tr B 12
An interaction B, usually called a bubble, is represented graphically by drawing a line labelled with the index c between a T and a T if the index c is summed between the two corresponding tensors, as pictured below on the left for the previous example. (4)

For a given choice of set B of interactions, the perturbative expansion of the free energy

F B {t B }, N = ln Z B ({t B }, N ) (5)
is labelled by connected Feynman diagrams which are generalized ribbon graphs, also called stranded graphs. The propagators contain D strands, one for each index of T , as shown above on the right of (4), and the interactions B ∈ B have a non-trivial internal structure. 6) is an example of a contribution to the 2-point function of a quartic rank-4 model. The corresponding diagrams Γ contain a certain number F c (Γ) of closed curves for each colour c, which we name faces. We denote

F (Γ) = D c=1 F c (Γ) (7)
the total number of faces of a diagram Γ. In the following, n B is half the order of the interaction B. Using the fact that there are 2n B half-propagators per interaction B, each contributing with a factor N -(D-1) , the amplitude associated to a Feynman diagram Γ is

A(Γ) = 1 c(Γ) B∈B t B B (Γ) B N B∈B B B (Γ) s B -(n B -1)(D-1) +F (Γ) , (8)
where B B (Γ) is the number of interaction bubbles of type B in Γ, and c(Γ) is a symmetry factor. The 1/N expansion of the free-energy of such a random tensor model can be understood as a discrete version of the formal Einstein-Hilbert partition function of D-dimensional pure Euclidean gravity [START_REF] Gurau | Invitation to Random Tensors[END_REF],

Z EH = topologies dg e -S EH , S EH = 1 16πG n d D x |g|(2Λ -R), ( 9
)
where Λ is the cosmological constant, G n is Newton's constant, and R is the Ricci scalar of the metric g. Indeed, every discrete manifold can be equipped with a geometry by taking all the edges to have a fixed length. Therefore, the integration over geometries can be replaced with a sum over discretized manifolds of a discrete version of the Einstein-Hilbert action. Stranded tensor diagrams Γ are dual to pseudo-manifolds such that for some set B of interactions

Z B EH = F B {t B }, N = Γ e B∈B B B (Γ) κ D-2 (Gn)Fint(B)-κ D (Gn,Λ)2n B +κ D-2 (Gn)F (Γ) , (10)
where κ D (G, Λ) and κ D-2 (G) do not depend on B and are obtained from Regge's action [START_REF] Regge | General relativity without coordinates[END_REF] for fixed lengths, and the number of internal faces F int (B) depends only on the stranded structure of the interaction B. We recover the amplitudes 1 Eq. ( 8) by setting N = e κ D-2 (Gn) , and t B = e κ D-2 (Gn) Fint(B)-s B +(n B -1)(D-1) -2κ D (Gn,Λ)n B . [START_REF] Aldous | The continuum random tree. II. An overview[END_REF] This explains the physical meaning of parameters N and t B in tensor-model quantum gravity, and justifies that enhanced tensor models can be interpreted as discrete Einstein-Hilbert theories of gravity without higher order curvature terms, something which has not been stressed before. In particular, the large-N limit corresponds to a weak-coupling limit in quantum gravity as κ D-2 (G n ) is inversely proportional to Newton's constant G n [START_REF] Gurau | Invitation to Random Tensors[END_REF].

Continuum limit and critical behaviour

In the large N limit, only contributions which maximize the exponent of N in Eq. ( 8) survive. This corresponds to the physical limit of small Newton's constant. The large-N free energy,

F LO {t B } = lim N →∞ 1 N D F {t B }, N , (12)
is a formal power series in the coupling constants. There are regions of the phase space for which F LO becomes singular (the inverse functions develop extrema). When reaching a singularity, the correlation functions do not depend anymore on the details of the Feynman diagrams but only on the way diagrams with large number of interactions scale. This corresponds to a continuum limit, for which the tensor model could, in some relevant cases, describe a non-perturbative theory of random fractal continuum spaces, a short-scale version of the Einstein-Hilbert action for pure gravity.

The relevance of such a continuum limit is partially encoded in the critical behaviour of the correlation functions when reaching singularities, that is, by the corresponding entropy critical exponent γ, also known in matrix models as the string susceptibility [START_REF] Di Francesco | 2D gravity and random matrices[END_REF]. It is obtained from the smallest non-integer exponent of the power expansion of the correlation functions around a singular point. For instance, in the simpler case of a single coupling t with dominant singularity t cr , the critical exponent γ is given by

F LO (t) = k-1 i=0 α i (t cr -t) ai + α k (t cr -t) 2-γ + o(t cr -t) 2-γ , (13)
where the α i may be null, all a i are integers and 2 -γ is not. In the case of a 2n-point function the first non-integer exponent is 2 -n -γ. In practice, it is easier to compute the critical exponent from the 2-point function since then the symmetry factors c(Γ) of the amplitudes in Eq. ( 8) are trivial. We interpret the critical exponent γ as follows.

• If γ = 1/2, the continuum limit is very likely to describe branched polymers [START_REF] Gurau | Melons are branched polymers[END_REF], also known as the continuous random tree [START_REF] Aldous | The Continuum Random Tree. I[END_REF][START_REF] Aldous | The continuum random tree. II. An overview[END_REF][START_REF] Aldous | The Continuum Random Tree III[END_REF]. This is the limit reached in random vector models or coloured random tensor models [START_REF] Bonzom | Critical behavior of colored tensor models in the large N limit[END_REF], and it is not a relevant limit for quantum gravity in our context 2 . We expect branched polymers in the continuum when the Feynman diagrams are in one-to-one correspondence with families of embedded trees which do not carry too much additional information3 .

• If the critical exponent is γ = -1/2, we are in the universality class of pure 2d gravity, and we expect the continuum limit to be the Liouville theory of 2d quantum gravity [START_REF] Polyakov | Quantum geometry of bosonic strings[END_REF][START_REF] Miller | Liouville quantum gravity and the Brownian map I[END_REF][START_REF] Miller | Liouville quantum gravity and the Brownian map II: geodesics and continuity of the embedding[END_REF][START_REF] Miller | Liouville quantum gravity and the Brownian map III: the conformal structure is determined[END_REF][START_REF] David | Liouville Quantum Gravity on the Riemann Sphere[END_REF]. Intuitively, this behaviour is expected when there is a one-to-one correspondence between the Feynman diagrams and families of planar diagrams which do not carry too much information, and which are not restricted by criteria so strong they change the universality class, e.g. by restricting the number of possible loops (embedded trees are planar).

• If the critical exponent is γ = 1/3, and if the characterization of contributing diagrams suggests so, we expect baby universes to proliferate [START_REF] Das | New Critical Behavior in d = 0 Large-N Matrix Models[END_REF][START_REF] Korchemsky | Matrix model perturbed by higher order curvature terms[END_REF][START_REF] Alvarez-Gaumé | A Proposal for strings at D > 1[END_REF]. This balanced regime interpolates between branched polymers and pure 2d gravity. It is obtained in the context of matrix models by adding a non-local 4 Tr(M 2 ) 2 polymerization perturbation to the potential. Such a multi-trace term encodes higher order curvature perturbations to the Einstein-Hilbert action and generates nodal points at which the discrete surfaces generated by one-matrix models "touch" each other without interacting. At leading order, the incidence relations between these surfaces are tree-like. If the coupling of the non-local term is too weak, the perturbation has no effect and the universality class is that of pure 2d gravity. If the polymerization is too strong, the nodal points dominate, and branched polymers are recovered in the continuum.

There is an intermediate value of the coupling for which the proliferation of nodal points competes with the diverging areas of the surfaces, leading to an intermediate regime for which baby universes proliferate, and the entropy exponent γ = 1/3 is positive. We stress however that, for matrix models, this regime is reached at the cost of introducing non-local interactions corresponding to higher order curvature terms. In this paper, we study various 4-dimensional models with local interactions 5 . Furthermore, they correspond to D-dimensional Einstein-Hilbert Euclidean partition functions without higher order curvature terms throughout relations Eq. ( 11). Still, they entail the same three regimes.

Multi-critical behaviour

For some regions of the phase space, the correlation functions exhibit more singular behaviour. For instance, this occurs when higher derivatives of the inverse functions vanish. To reach such regions, one typically has to fine-tune the couplings. Higher order critical exponents appear, which depend on how singular the correlation functions get, but also on the critical exponents when the couplings are not fine-tuned. For instance, when the generic critical behaviour is that of branched polymers γ = 1/2, the multi-critical exponent of order m is [START_REF] Ambjorn | Summing over all genera for d ¿ 1: a toy model[END_REF]

γ = m m + 1 . ( 14
)
When the generic critical behaviour is that of 2d pure gravity, the order-n multi-critical exponent is [START_REF] Kazakov | The Appearance of Matter Fields from Quantum Fluctuations of 2D-GRAVITY[END_REF]

γ = - 1 n + 1 . (15)
To obtain a multi-critical exponent of order q = m (branched polymer) or q = n (planar), one typically needs to consider at least q interactions of different (polynomial) order, and tune the couplings to reach the right subspace of codimension q in the coupling phase space. It is possible to interpret the interactions of higher order as certain harddimer configurations on the random lattices obtained from only the lowest-order interactions. In the case of matrix models and 2d gravity, coupling random planar diagrams to hard-dimers leads to the Yang-Lee edge singularity in the continuum limit, a non-unitary conformal matter coupled to gravity [START_REF] Staudacher | The Yang-Lee edge singularity on a dynamical planar random surface[END_REF]. This motivates the study of multi-critical behaviour in the context of enhanced random tensor models. However, the entropy exponent γ is in general just a first indicator of what the continuum theory can be. For instance, an exponent γ = -1/3 does not necessarily mean that the continuum theory is the Yang-Lee singularity, as the same exponent is recovered for the Ising model on random lattices [START_REF] Kazakov | The Appearance of Matter Fields from Quantum Fluctuations of 2D-GRAVITY[END_REF][START_REF] Kazakov | Ising model on a dynamical planar random lattice: Exact solution[END_REF], leading to unitary conformal matter coupled to 2d gravity in the continuum.

In the transition regime of proliferating baby-universes multi-criticality is a particular subtle issue. In [START_REF] Korchemsky | Matrix model perturbed by higher order curvature terms[END_REF] it is argued that, tuning a matrix model to order-n multi-criticality and adding a strong-enough polymerization (nonlocal) term Tr(M 2 ) 2 , the 1/3 exponent corresponding to proliferating baby-universes could be changed to a

γ = 1 n + 2 (16)
multi-critical exponent. It is further argued in [START_REF] Alvarez-Gaumé | A Proposal for strings at D > 1[END_REF] that the exponent

γ = m n + m + 1 ( 17
)
could be reached using non-local (multi-trace) interactions of higher order. We will show that the latter is not realized in the tensor models considered here.

The scalings s B

In order for the 1/N expansion of 2n-point functions to be defined, the scalings s B in Eq. ( 1) have to be chosen such that the dependence of the amplitude of Feynman diagrams in N is bounded by N f (n) , where f (n) only depends on n. For the free-energy, we expect f (n) = D, D being the rank of the tensor T , as N D is the number of degrees of freedom. However, any choice leading to bounded amplitudes is not considered a good choice: if there are finitely many contributions to the leading order, the large N correlation functions do not exhibit any singularity, and the theory is trivial in the continuum. Gurau proved [START_REF] Gurau | Universality for random tensors[END_REF] that the choice s B = -δ(B) ≤ 0, δ(B) being the Gurau degree of the bubble (see [START_REF] Gurau | Universality for random tensors[END_REF]), leads to a defined 1/N expansion. The interaction bubbles with vanishing Gurau degree have an almost vector-like structure, and are called melonic (e.g. interactions labelled 2 M , 3 M , 3 M in Table I, II). We do not describe them precisely, but for our purposes, the important properties we need are:

-the only scaling leading to a defined and non-trivial 1/N expansion for a melonic interaction is s melonic = 0, -the leading order of melonic random tensor models is tree-like, the 2-point function has a square-root singularity, so that the critical exponent is γ = 1/2, and they lead to branched polymers in the continuum [START_REF] Gurau | Melons are branched polymers[END_REF].

A non-trivial consequence of a result from [START_REF] Gurau | Regular colored graphs of positive degree[END_REF] is that for every non-melonic case, the choice s B = -δ(B) leads to a finite number of Feynman diagrams contributing to each order of the 1/N expansion, and in particular to a trivial leading order. For known cases of bubbles B there is only one particular integer value of s B satisfying the two conditions above, and one has to look at each model individually to determine it. This value is positive for non-melonic interactions. Models with a choice of scaling leading to a defined and non-trivial 1/N expansion are known as enhanced tensor models in the literature [START_REF] Bonzom | Enhancing non-melonic triangulations: A tensor model mixing melonic and planar maps[END_REF][START_REF] Bonzom | Large N limits in tensor models: Towards more universality classes of colored triangulations in dimension d ≥ 2[END_REF].

The aim of this paper is to compute the appropriate scaling s B and the critical exponent γ for every type of quartic and sextic interaction B in dimension 3 and 4. The results we obtain are combined in Tables I and II. As will be clear in the following section, some interactions only differ by exchanging colours. In the tables, we only represent one interaction of each type, but we provide in the second column the number of interactions which only differ by an exchange of colours. In this way, the number of sextic interactions of rank 3 and 4 sums up to 7 and 41 respectively, in agreement with [START_REF] Ben Geloun | Counting Tensor Model Observables and Branched Covers of the 2-Sphere[END_REF].

We recall that the amplitude of a vacuum diagram Γ is given by ( 8), up to a symmetry factor. The symmetry factor is one for the two-point function. The accurate scaling s B leading to a defined and non-trivial 1/N expansion is deduced from the number of faces [START_REF] Gurau | The Complete 1/N Expansion of Colored Tensor Models in Arbitrary Dimension[END_REF] of leading-order diagrams, and the critical exponent γ is given by the critical behaviour of the resulting leading-order two-point function, expected to be an indicator for the continuum limit as explained.

In this paper, we have chosen to use the stranded representation of interactions leading to diagrams which generalize ribbon graphs. Another graphical representation is commonly used in the tensor models literature, in which tensors T are represented as white vertices, tensors T are represented as black vertices, and a line of colour c is drawn between a black vertex and a white vertex if a strand of colour c goes between the corresponding tensors T and T in the stranded interaction. We represent these graphs as well in the third column of the tables. Finally, we represent the corresponding intermediate-field vertex which will be explained in detail in Section II.

B. Phase transition from a single type of interaction: quartic necklaces

The leading order of matrix models with single-trace interactions of orders {k 1 , k 2 , • • • } involves planar ribbon graphs with vertices that have {k 1 , k 2 , • • • } incident propagators. The critical behaviour of the correlation functions is given by the exponent γ = -1/2, or, if the couplings are fine-tuned, by the multi-critical exponents -1/m. It is still possible to obtain other critical regimes, such as that of proliferating universes, characterized by the exponent γ = 1/3, but one needs to consider multi-trace interactions [START_REF] Das | New Critical Behavior in d = 0 Large-N Matrix Models[END_REF][START_REF] Alvarez-Gaumé | A Proposal for strings at D > 1[END_REF], graphically represented by non-connected stranded interactions. Here, we introduce a one-tensor model with connected stranded interactions, which presents a branchedpolymer phase, a 2d planar phase, and a proliferation of baby universes at the phase transition between them. It happens to be equivalent to a multi-matrix N 2 × N 2 hermitian model (with usual trace). This tensor model has three interactions, that previously introduced in (4) and two obtained by exchanging colours.

B 12 2 N
These interactions have been called (quartic) necklaces in [START_REF] Bonzom | Enhancing non-melonic triangulations: A tensor model mixing melonic and planar maps[END_REF] because of their representation in graphs (see the second and fifth rows of Table II). We conserve this denomination here. The corresponding partition function is dµ 0 (T, T ), [START_REF] Miller | Liouville quantum gravity and the Brownian map III: the conformal structure is determined[END_REF] with coupling g i = t 1i 2 N corresponding to the interactions B 1i 2 N . We know from [START_REF] Bonzom | Enhancing non-melonic triangulations: A tensor model mixing melonic and planar maps[END_REF] that the appropriate scaling is s B = 1. The amplitude (8) rewrites as

Z 4 (g 2 , g 3 , g 4 , N ) = C D
A(Γ) ∼ g B2(Γ) 2 g B3(Γ) 3 g B4(Γ) 4 N F (Γ)-2B(Γ) , (20)
where B i denotes the number of interaction bubbles of type B 1i 2 N , and B = B 2 + B 3 + B 4 . A simple argument shows that the random tensor model for only a single non-vanishing coupling, e.g. g 3 = 0, g 4 = 0, is described by a one-matrix model: by grouping together the colours 1,2, and 3,4 in (4), one obtains ribbon graphs for a matrix M (12)(34) = T 1234 in which strands carry two colours instead of one.

M M M M M M (21)
Such ribbon graphs label the perturbative expansion of a quartic N 2 × N 2 matrix model, with partition function

Z 4 (g 2 , N ) = dM d M π N 2 e -TrM M † -g2N -2 TrM M † M M † . (22)
In the case where not all three couplings vanish, the tensor model can also be rewritten as a matrix model. Indeed, it was shown in [START_REF] Bonzom | Enhancing non-melonic triangulations: A tensor model mixing melonic and planar maps[END_REF] that the Feynman diagrams of this theory are one-to-one with embedded diagrams, or ribbon graphs, with vertices of any valency (the valency is the number of incident propagators), and three kinds of propagators, those carrying colours 1 and 2 and coming with a factor g 2 , those carrying colours 1 and 3 and coming with a factor g 3 , and those carrying colours 1 and 4 and comming with a factor g 4 . Embedded diagrams label the perturbative expansion of matrix models, and the model turns out to be described by the following three-matrix model

Z 4 (g 2 , g 3 , g 4 , N ) = e -Tr ln 1 ⊗4 -i √ 2g2M12⊗134-i √ 2g3M13⊗124-i √ 2g4M14⊗123 dν(M 12 , M 13 , M 14 ), (23)
where M 12 , M 13 , and M 14 are N 2 ×N 2 hermitian matrices that propagate their subscript colours, 1 ij is the identity in colours i and j, and the measure is dν(M ) , Z 0 being chosen so that Z 4 (0, 0, 0, N ) = 1. This is the intermediate-field representation of the tensor model. Denoting E 1i the number of propagators of type 1i, the amplitude of a ribbon graph M is

A(M ) ∼ g E12(M ) 2 g E13(M ) 3 g E14(M ) 4 N F (M )-2E(M ) , (24)
where now

F (M ) = 4 i=1 F i (M )
is the sum of the faces of colours 1, 2, 3 and 4 of the embedded diagram, and E is the total number of propagators. Leading-order diagrams are planar, but also satisfy an additional constraint on the possible loops, which depends on the colours of the propagators.

The usual three regimes also appear in this case, as we show in Subsection III C: for two vanishing couplings, this reduces to a usual matrix model, in the universality class of pure 2d gravity (γ = -1/2); for g 2 = g 3 = g 4 , there is a destructive interference between the three kind of matrices which reduces the critical behaviour to a square root singularity, with exponent γ = 1/2. Some intermediary values of the coupling constant reproduce the regime of proliferating baby universes, with exponent γ = 1/3.

The only other type of quartic interactions in dimension 4 is melonic (see Table II). The quartic melonic model is the most studied uncoloured tensor model. As all melonic models, it has a tree-like leading order characterized by the critical exponent γ = 1/2 and converging towards branched polymers in the continuum limit. Coupling a quartic melonic interaction and an interaction of the kind in Eq. ( 18), the same three regimes appear, as detailed in [START_REF] Bonzom | Enhancing non-melonic triangulations: A tensor model mixing melonic and planar maps[END_REF].

C. The enhanced K3,3 tensor model Before considering all sextic interactions for rank D = 4 tensors, it is useful to understand the sextic interactions at rank D = 3. There are 3 different kind of sextic interactions in dimension 3. Two of them are melonic (type 3 M and 3 M interactions in Table I), and therefore are not to be enhanced (s = 0). Their critical behaviour is known to be that of branched polymers. We therefore focus on the 3 dimensional theory with the following non-melonic interaction Tr K3,3 (T, T ) = N a,b,c,d,e,f,g,h,i=1 T abc Tade T f dg Thbg T hie Tfic .

and corresponding partition function

Z K3,3 (t 3k , N ) = C D e t3 k N 3 Tr K 3,3 (T, T ) dµ 0 (T, T ). ( 26
)
It is represented below, as well as the corresponding propagator.

It is referred to as the K 3,3 interaction in the random tensor literature, as one obtains the complete bipartite K 3,3 graph by replacing each tensor by a vertex. It is a natural model to consider as it is symmetric upon exchange of colours and no non-planar interactions have been considered before (you can not draw it without a crossing). It has been shown in [START_REF] Bonzom | Colored triangulations of arbitrary dimensions are stuffed Walsh maps[END_REF] that the number of faces Eq. ( 7) of leading-order diagrams is

F LO = 3(1 + B), ( 28
)
where B is the number of interaction bubbles. From Eq. ( 8), this implies that the strongest scaling one can choose is s = 1. Knowing this, the amplitude of Feynman diagrams (8) rewrites as

A(Γ) ∼ t B(Γ) 3k N F (Γ)-3B(Γ) . (29)
In the large N limit, only diagrams which maximize F (Γ) -3V (Γ) survive. They have been identified in [START_REF] Bonzom | Colored triangulations of arbitrary dimensions are stuffed Walsh maps[END_REF], however, although a γ = 1/2 exponent was expected, it had not been computed before. They scale in O(N 3 ) and we show in Section II B that the leading-order two-point function satisfies the following self-consistency equation,

G 2 (t 3k ) = 1 + 3t 3k G 2 t 3k ) 3 + 3t 2 3k G 2 (t 3k ) 6 . ( 30
)
Using standard results in singularity analysis of tree generating functions [56, ch. VII.4] we find that the solution to this equation has a square-root expansion

G 2 (t 3k ) = G 2 (t cr ) -c 1 -t 3k /t cr + o 1 -t 3k /t cr , (31)
indeed leading to a γ = 1/2 exponent, and where t cr > 0, G 2 (t cr ) > 0, and c > 0 can be computed numerically. To summarize, we present all interactions of order 6 for rank D = 3 in table I.

D. Order six theory in four dimensions

In dimension 4, besides melonic interactions 3 M and 3 M (see Table II), there are necklace interactions 3 N , and the interactions of type 3 MN , which all generate leading orders similar to that of quartic necklaces, there is the interactions 3 N , which, as we will explain in Subsection II C, can be understood as an enrichment of the quartic necklaces, and there are two kinds of interactions, 3 K and 3 K , which are obtained from the K 3,3 interaction by adding an additional colour. As detailed in Subsection II C, it is very easy to show that interaction of the type 3 K have a tree-like leading order with exponent γ = 1/2. We therefore describe the 3 K tensor model here. The interaction is obtained from the K 3,3 interaction by adding strands of colour 4 parallel to the strands of colour i, where i = 1, 2, or 3 (and by possibly exchanging colours). An interaction of this type is b,c,d,e,f,g,h,i,j,k,l=1 T abcd Taefg T heij Tkbij T klf g Thlcd .

Tr B3 K (T, T ) = N a,
(

) 32
We will study the critical regimes for an order six theory that also includes the following type 3 M melonic interaction, T abcd Taefg T hef g Thijk T lijk Tlbcd .

The corresponding partition function is

Z 6 (λ, N ) = C D e t3 k N 3+s 3 K Tr3 K (T, T )+t3 M N 3 Tr B 3 M (T, T ) dµ 0 (T, T ). ( 34
)
Because the interaction in Eq. ( 33) is melonic, we should have a square-root critical behaviour in t 3 M , when t 3k = 0. We will show in Subsection II C that the right scaling is s 3K = 2. We will also show that the three critical regimes with exponents 1/2, -1/2, and 1/3 are obtained for this model.

II. LEADING-ORDER DIAGRAMS

In this section we explain the derivation of the Dyson-Schwinger equations for the tensor models under consideration by counting their Feynman diagrams. The evaluation of the perturbative sum of the 2-point function can be performed most efficiently transforming the tensor-model state sum to an intermediate-field representation involving embedded diagrams [START_REF] Gurau | The 1/ N Expansion of Tensor Models Beyond Perturbation Theory[END_REF][START_REF] Delepouve | Borel summability and the non perturbative 1/N expansion of arbitrary quartic tensor models[END_REF][START_REF] Bonzom | Colored triangulations of arbitrary dimensions are stuffed Walsh maps[END_REF]. This was illustrated in section I B for the quartic tensor model with a single interaction of type B 1i 2 N , for which the Feynman diagrams are in one-to-one correspondence with four-valent ribbon graphs, and in which closed strands carry a weight N 2 . For other models, there generally is no one-to-one map as simple but one can develop more subtle transformations.

There is a one-to-one correspondence between Feynman diagrams and embedded diagrams equivalent to ribbon graphs, introduced in [START_REF] Bonzom | Colored triangulations of arbitrary dimensions are stuffed Walsh maps[END_REF], in which interactions are replaced by white vertices of a specific kind. White vertices can only be linked to black vertices, and conversely. Propagators, which we will call edges in this section, carry various colour-sets I ⊂ {1, • • • , D} that depend on the interactions, and black vertices have no restriction on possible incident edges. Importantly, embedded diagrams, also sometimes called fatgraphs, are graphs for which the cyclic ordering of edges around vertices matters: two diagrams which only differ by the ordering of edges around vertices are different diagrams. Consequently, two half-edges incident to the same vertex and following each other are separated by a well-defined corner. A face in an embedded diagrams is a closed walk: a cyclic succession of corners that are on the same side of the same edge (they correspond to the closed strands of the equivalent ribbon graph). The propagators in the stranded representation are mapped to corners in the new intermediate-field representation, and crucially, the faces of colour i of the stranded diagrams (Eq. ( 7)) are mapped to faces around the embedded diagrams obtained by keeping all the black vertices, and only the edges whose colour-set contains colour i, as can be seen on the example of [START_REF] Ambjorn | Baby universes in 2-d quantum gravity[END_REF]. We can therefore apply ribbon graphs tools to solve tensor models. We refer the reader to the Appendix for a more detailed explaination of this correspondence.

To obtain the white vertex corresponding to a given interaction, one has to choose a way to pair every tensor T with a T . The intermediate-field vertices corresponding to each interaction are given in the fifth column of Tables I andII. The choices of pairings on order-6 interactions made to obtain these embedded vertices are given by the following rules,

Given an embedded diagram M, we denote L(M) its number of independent loops, given in terms of its number of edges E, vertices V , and connected components K by

L(M) = E(M) -V (M) + K(M). (37)
A tree is a diagram with no loops. The genus g of an embedded diagram is given in terms of its number of faces F by

2K(M) -2g(M) = V (M) -E(M) + F (M). ( 38
)
It is the minimal number of holes of a surfaces on which the diagram can be drawn without crossings. If M is an intermediate-field diagram, we may consider the colour-i sub-diagram obtained by keeping all the black vertices and all the edges whose colour-set contains colour i. We respectively denote L i (M), F i (M), and g i (M), the number of loops, faces, and the genus of the corresponding colour-i sub-diagram. If M is an intermediate-field diagram and T is an intermediate-field tree with the same number B B (M) of interaction bubbles of each type, then the following identity holds, which was shown in [START_REF] Bonzom | Colored triangulations of arbitrary dimensions are stuffed Walsh maps[END_REF],

F (T ) -F (M) = DL(M) + 2 D i=1 g i (M) -L i (M) . (39)
We will mainly use this formula to determine leading-order diagrams of the tensor models. The number of faces of trees is always linear in the number of interaction bubbles of each kind,

F (T ) = D + B∈B a B B B (T ), (40)
where a B only depends on the interaction bubble and the way the tensors T and T have been paired (it is the number of faces when putting the propagators between the pairs (T, T ), minus the rank of the tensor) [START_REF] Bonzom | Colored triangulations of arbitrary dimensions are stuffed Walsh maps[END_REF]. If trees belong to the leading order, then a B gives the strongest scaling s B one can choose to have a defined 1/N expansion

s B = (n B -1)(D -1) -a B . (41)
All interactions of order 4 and 6 in dimensions 3 and 4 are such that, for the choices of intermediate-field vertices presented in Tables I and II, trees belong to the leading order, and we can therefore deduce the appropriate scalings from the number of faces of trees. With this choice of scaling, the amplitudes of intermediate-field diagrams scale as

A(M) ∼ N D 1-L(M) +2 D i=1 Li(M)-gi(M) = O(N D ). (42)
A. Intermediate-field representation of rank D = 4 quartic tensor models

For quartic theories with interactions of the type 2 N (18) or of the type 2 M (first row of Table II), the transformation can be done directly on the correlation functions using the Hubbard-Stratonovich transformation [START_REF] Stratonovich | On a Method of Calculating Quantum Distribution Functions[END_REF][START_REF] Hubbard | Calculation of Partition Functions[END_REF]. This can be seen e.g. in [START_REF] Delepouve | Quartic Tensor Models[END_REF]. However we here focus on describing this correspondence at the diagram level. We choose to pair a tensor T and a T , and we see the quartic interactions as new kind of edges between these two pairs6 , as illustrated in [START_REF] Kitaev | A simple model of quantum holography[END_REF]. The fact that interactions become edges is why this new representation is called the intermediate-field representation.

1j 1 ↔ 1j 1 (43)
Of course there are several choices of pairings. If the interaction is of type B i 2 M , we always choose to pair the tensors that are linked by 3 strands, and if the interaction is of type B 1i

2 N , we always choose to pair tensors that are not linked by a strand of colour 1. The colours which propagate between the chosen pairs now label the corresponding edges in the new representation. The vertices in the new representation correspond, in the stranded representation, to cycles alternating pairs of tensors not linked by colour 1 and propagators. In the quartic case, one has the following theorem, first proven in [START_REF] Bonzom | Enhancing non-melonic triangulations: A tensor model mixing melonic and planar maps[END_REF]: there is a one-to-one correspondence between Feynman diagrams of enhanced four dimensional quartic tensor models and embedded diagrams with edges carrying either one colour {i}, where i = 1, 2, 3, 4, or a set of two colours {1, 2}, {1, 3} or {1, 4}. Crucially, the faces of colour i of the original Feynman diagrams are mapped to the faces of the embedded colour-i sub-diagram. The total number of faces of an embedded diagram M is

F (M) = 4 i=1 F i (M). ( 44
)
Note that if there are no edge with label-set containing colour i, then F i (M) = V (M). Similarly, if colour 1 appears on every edge, as it is the case if all interactions are of the B 1i 2 N type, then F 1 (M) is just the total number of faces of the diagram M. The dependence in N of the amplitude of a diagram is

A(M) ∼ N F (M)-3 4 i=1 Ei(M)+2 4 i=2 E1i(M) , (45)
where E i are edges carrying colour-i only, and E 1i are edges labelled with colours 1 and i. Consequently, the embedded diagrams contributing to the leading order of the tensor model are those which maximize F at fixed number of edges of each kind. Leading-order contributions are identified by first looking at the behaviour of trees, which are connected embedded diagrams M such that L(M) = 0, and then studying the way F changes when deleting edges. One can easily prove recursively on the number of edges [START_REF] Bonzom | Colored triangulations of arbitrary dimensions are stuffed Walsh maps[END_REF] that a tree T is such that

F (T ) = 4 + 3 4 i=1 E i (T ) + 2 4 i=1 E 1i (T ), (46)
corresponding to Eq. ( 40), and resulting in amplitudes scaling as

A(T ) ∼ N 4 . ( 47
)
Furthermore, one can easily show that if a diagram M contains an edge labelled with a single colour and which is not a cut-edge (a cut-edge is such that the number of connected components increases by 1 when being deleted), then there exist another diagram M such that

A(M) ∼ A(M )N -α , (48)
with α > 0. Thus, M cannot contribute to the leading order. To identify leading-order diagrams we now use Eq. ( 39).

Here, L 1 = L and g 1 = g, so that if M and a T have the same number of edges of each kind and if colour-i edges in M are all cut-edges, the scaling in N of the amplitude is

4 -F (T ) -F (M) = 4 -2L m (M) -2 g(M) + 4 i=2 g i (M) , (49)
where all quantities have been defined before, except

L m (M) = L(M) - 4 i=2 L i (M) ≥ 0, ( 50
)
the number of independent loops in the graph that contain edges of at least two sets of colours 1i and 1j with i = j. Consequently, a diagram contributing to the leading order is planar (g(M) = g i (M) = 0), and only has loops of edges that carry the same two colours 1 and i (L m = 0) for any i. A schematic example is shown in [START_REF] Staudacher | The Yang-Lee edge singularity on a dynamical planar random surface[END_REF] (each blob corresponds to a planar component with only edges carrying the same colour-set {1i}). distinguished edge, so that all symmetries break down. Distinguishing an edge goes back to choosing a particular pair (T, T ). On each side, given a T , there are three different ways of choosing a T . Consequently, the self-consistency equation of the 2-point function

G 2 is G 2 (t 3k ) = 1 + 3t 3k G 2 t 3k ) 3 + 3t 2 3k G 2 (t 3k ) 6 . ( 58
)
We stress that the self-consistency equation for the 2-point function of all sextic melonic interactions would be the same, but without the order 6 term. The effect of the non-planarity is therefore to generate loops at leading order, but which are too small to change the large-scale structure of leading-order diagrams.

For this model, we introduce a new correspondence between Feynman stranded diagrams and embedded diagrams, subtler than that described at the beginning of Sec. II and in the Appendix, which renders the tree-structure of leading-order diagrams explicit at the cost of introducing extra colours on the (white) intermediate-field vertices corresponding to the interaction bubbles. This additional colouring takes into account the way the tensors T and T have been paired to obtain the patterns of (57). 7 Leading-order stranded diagrams are mapped to embedded diagrams with coloured vertices and black vertices such that edges -which carry colour sets {12}, {13}, {23}, or {123} -can only link coloured vertices to black vertices. The black vertices are of any valency and are similar to the vertices introduced in Section II A. The other vertices are of valency 3 (colours a, b and c), or 6 (colours d, e, and f ), and are rigid in the sense that they have a fixed ordering of edges around them. The cyclic counterclockwise ordering of edges around valency-3 vertices a, b and c, is given by their colour-sets according to {12}, {13}, {23}, as shown in [START_REF] Goulden | Combinatorial Enumeration[END_REF]. As before, the faces of colour i of the stranded diagram in the original picture are now mapped to the faces of the corresponding colour-i sub-diagram. Leading-order diagrams are mapped to trees in the new representation, and other diagrams are mapped to embedded diagrams as just described, which have at least one loop and have a few forbidden small sub-diagrams. 8 The number of faces of a tree in the new representation is recursively proven to be [START_REF] Flajolet | Analytic Combinatorics[END_REF], which corresponds to equation Eq. ( 40) and provides a 3K = 3, and therefore s 3K = 1 from Eq. [START_REF] Kazakov | The Appearance of Matter Fields from Quantum Fluctuations of 2D-GRAVITY[END_REF]. As in the previous subsection, we define L 12 ≥ 0, the number of independent loops of M which are not counted in L 1 and L 2 , as

L 12 = L -L 1 -L 2 . (61)
In particular, we will use the two following properties:

L ≥ L i + L j (62)
for any i = j, and

L ij (M) > 0 ⇒
There exist some edge in M which is a cut-edge in Mi and Mj but is not a cut-edge in M .

Using Relation [START_REF] Bonzom | Large N limits in tensor models: Towards more universality classes of colored triangulations in dimension d ≥ 2[END_REF] which compares the number of faces of intermediate-field diagrams and trees, we see that a diagram M satisfies

F (M) = 3(1 + B(M )) -∆ 3 (M) = 0, (64)
where ∆ 3 (M) is the quantity

∆ 3 (M) = L 12 + L 13 + L 23 + 2 3 i=1 g i (M) ≥ 0, (65)
which is non-negative as a sum of positive quantities, and vanishes iff all terms vanish. The amplitude of diagrams in the new representation is

A(M) ∼ t B(M) K N 3-∆3(M) , (66)
where the number of bubbles B(M) = V col (M) is the number of coloured vertices and we can now restate the condition a diagram needs to satisfy in order to contribute at leading order,

M is L.O. ⇔ L 12 = L 13 = L 23 = 0 g 1 = g 2 = g 3 = 0 . ( 67
)
This is the characterization used in [START_REF] Bonzom | Colored triangulations of arbitrary dimensions are stuffed Walsh maps[END_REF][START_REF] Lionni | Colored discrete spaces: higher dimensional combinatorial maps and quantum gravity[END_REF] to prove that the leading order diagrams are as shown in Eq. ( 57). The proof is quite lengthy, so we do not repeat it here, but just sketch it: The leading order embedded diagrams are characterized recursively on the number of white vertices. The leading order diagrams with one and two vertices are identified to be the embedded diagrams corresponding to the diagrams of Eq. ( 57) but without the G 2 blobs. Then it is recursively shown that if ∆ 3 (M) = 0, then M contains a subgraph as in Eq. ( 57), and the blobs still satisfy ∆ 3 (M) = 0, so the recursion hypothesis applies to them. For more details, we refer the reader to [START_REF] Lionni | Colored discrete spaces: higher dimensional combinatorial maps and quantum gravity[END_REF].

C. The sextic case in dimension four

We explain now the leading-order behaviour of all D = 4 sextic interactions. As mentioned in Subsection I D, the interactions of type 3 M and 3 M in Table II are melonic, and therefore their critical behaviour and continuum limit is known. The tensor model built with necklace interactions of type 3 N , as well as melonic interactions 3 M and M is similar to the full quartic rank-4 tensor model in all aspects, and the proof follows the exact same steps. Trees belong to the leading order and behave as Eq. ( 40) with number of faces

F 3N LO = 4 + 4B 3N + 6(B 3M B 3 M ) (68)
such that a 3N = 4, and therefore s 3N = 2 from Eq. ( 41). Denoting P

ns (z) the generating function of planar ribbon graphs with white vertices of valency three and black vertices of any valency, and such that no face visits twice or more the same black vertex, counted according to their number of black vertices9 , the self-consistency equation for the leading-order 2-point function of this model is

G 2 = 1 + 4 i=2 P (3) ns (t 1i 3 N G 3 2 ) -1 + G 3 2 4 i=1 t i 3 M + 4 i<j=1 t ij 3 M (69)
as G 2 is inserted on every corner around a black vertex, and there are as many such corners as edges in the diagram. We expect the same critical regimes as in the quartic case.

Interactions of type 3 MN are obtained from a quartic necklace by "opening" a strand of colour i and inserting two tensors linked by three new strands of colours different from i. These two tensors are paired, and, depending on whether colour i is colour 1 or parallel to colour 1, or not, we have the following intermediate-field vertices (with the pairing rule (35)),

and their equivalents. The edge corresponding to the insertion always carries a single colour. As in the quartic case, such an edge is necessarily a cut-edge in every leading-order diagram. Indeed, unhooking an edge e as shown in (71) (here, we cannot delete an edge because of the rigid structure of white vertices),

e → e (71)
it is easily seen that, if e carries a single colour, the number of faces necessarily increases at least by two. Therefore, if unhooking the edge does not increase the number of connected components, the diagram is suppressed at large N . Leading-order diagrams of a type 3 MN tensor model are obtained from leading-order diagrams from necklace-type models by performing the following operations on the edges 1i L.O.

→ 1i G 2 1 or i 1i or 1i G 2 j = 1, i 1ij , (72)
where G 2 is a 2-point leading-order insertion from the current order-6 model. Consequently, trees belong to the leading order, and it is seen recursively that their number of faces behave as Eq. ( 40) with

F 3MN LO = 4 + 5B, (73)
yielding a 3MN = 5, and therefore s 3MN = 1 from Eq. ( 41). A distinguished edge in a leading-order diagram is always incident to a non-separable planar diagram with all edges of colour 1i, with one full G 2 insertion on every corner, and another G 2 inserted on each edge, with 3 ways of deciding which one is the distinguished edge. For each edge, one has to decide which one of the four colours c of the quartic necklace B 1i 2 N has received an insertion, contributing with a coupling t ci 3MN . The self-consistency equation for the 2-point function of the full type-3 MN tensor model is

G 2 = 1 + 4 i=2 P ns 3G 3 2 4 c=1 t ci 3MN -1 . (74)
In the case where a single interaction B ci 3MN , with c ∈ {1, 2, 3, 4} and i ∈ {2, 3, 4}, is allowed it reduces to the simpler equation

G 2 = P ns 3G 3 2 t ci 3MN . (75)
We now consider an interaction of type 3 K , and show that the leading-order intermediate-field diagrams are precisely trees. Suppose that 4 is the colour which belongs to every double-strands, as in the diagram of Table II. We apply relation Eq. [START_REF] Bonzom | Large N limits in tensor models: Towards more universality classes of colored triangulations in dimension d ≥ 2[END_REF], which compares faces of trees T and diagrams M with the same number of interactions:

F (T ) -F (M) = L(M) -2L 4 (M) + 2g 4 (M) + ∆ 3 (M) (76) ≥ L(M) -2L 4 (M) + 2g 4 (M), ( 77
)
where ∆ 3 has been defined in Eq. (65). As colour 4 does not appear on any edge, L 4 (M) = g 4 (M) = 0 such that

F (T ) -F (M) = L(M). (78)
This vanishes if M is a tree and is positive otherwise. The self-consistency equation is a simple cubic tree-equation, and, as for 3 MN type interactions, a 3 K = 5, and s 3 K = 1. This argument holds for a single type 3 K colouring, however we also expect a tree-behaviour for the full 3 K tensor model.

We consider a single interaction of type 3 K , with the colouring as in Table II. The number of faces of trees is Eq. ( 40)

F (M) = 4(1 + B(M)). (79)
Moreover, diagrams are such that L 3 = L 4 and (61, 62, 63) are true for i = j < 4, so that (39) rewrites as

∆ 4 (M) = F (T ) -F (M) = 2(L 13 (M) + L 23 (M)) + 2 4 i=1 g i (M), (80)
so that trees belong to the leading order, from which we deduce that, as a 3K = 4, s 3K = 2 from Eq. ( 41), and other leading-order contributions are such that L 13 = L 23 = 0, and ∀i, g i = 0. (81)

There are two cases. If in addition L 12 = 0, then we know from the previous section that in our representation, the diagram is a tree. If not, L 12 (M) > 0, and from (63) we know that there exist an edge e which is not a cut-edge, but is so in both M 1 and M 2 . Unhooking e as in (71), we obtain a new diagram M such that,

L = L -1, L 1 = L 1 , and L 2 = L 2 . ( 82
)
Using ( 62) and (81),

0 ≤ L -L 1 -L 3 = L -1 -L 1 -L 3 = L 3 -L 3 -1. ( 83
)
Furthermore, L 3 -L 3 is 0 or 1, but the first value contradicts (83), so that

L 3 = L 3 -1, (84)
and therefore M is obtained from M by hooking an edge of colour i34 to create a loop such that M 3 = M 4 stays planar and such that it is a cut-edge in M i . We still have L ≥ L 1 + L 2 = L 1 + L 2 , so there are two cases. Either L > L 1 + L 2 , in which case we repeat the same steps, or L = L 1 + L 2 , in which case M is a tree. From this tree, M is obtained by hooking a certain number of edges i34 so that every edge of colour 12 stays a cut-edge and such that the overall diagram remains planar.

Leading-order diagrams are quite similar to those obtained from a single type 3 MN interaction, with the difference that here, the tree part has three valency-3 vertices and three valency-6 vertices. We distinguish an edge with colour set {134} or {234}. To recover the leading-order 2-point function, one simply replaces z with 3/2z. There are two cases. In the first case, the distinguished edge is either a cut-edge of colour set {134} or {234}, and it belongs either to one of the three 3-valent vertices or to one of the three 6-valent vertices. The contribution is 6t 3k G 2 (t 3k ) 3 + 6t 2 3k G 2 (t 3k ) 6 . In the other case, the distinguished edge is not a cut-edge, in which case it belongs to a non-separable planar component Eq. [START_REF] Gurau | Universality for random tensors[END_REF] for which all edges contain the colours 34, and we need to insert one 2-point function per black corner of the non-separable component and one 2-point function per edge -corresponding to the colour-{12} cut-edge -that can be added in two possible ways. Indeed, a white vertex is added on every edge, and a colour-{12} bridge is added between that white square and a leading order sub-diagram, and there are two ways of deciding which one of the two colour-{34} segments has colour-set {234}, and which one has colour-set {134}. This second possibility contributes as P ns 2zG 3 2 . It generates a tree part 2zG 3 2 which has to be subtracted from the overall tree part. The 2-point function G 2 therefore satisfies

G 2 = G 2 (t 3k ) = P ns 2t 3k G 3 2 + 4t 3k G 3 2 + 6t 2 3k G 6 2 . ( 85
)
As for interactions of type 3 N , we notice that a propagator can be inserted on every interaction bubble,

transforming every bubble into two quartic necklaces with different colouring and a distinguished propagator between them, without changing the number of faces. A type 3 N diagram belongs to the leading order if the diagram obtained by inserting a propagator for every bubble is a type-2 N leading-order diagram. Because the inserted propagators are distinguished, leading-order diagrams of type 3 N correspond to leading-order diagrams of type 2 N for which there exist a way of choosing propagators such that every bubble is incident to a distinguished propagator, and such that every distinguished propagator is between two quartic necklaces with different colourings.

the additional parameter u (see Eq. ( 53)) and one of the couplings. The advantage is that all critical values including G cr 2 are directly determined.

A second way to find the critical loci is to reduce the parametric equations to a single one and demand its discriminant to be vanishing. For Eq. ( 88), elimination of the extra variable of the parametric definition of P ns Eq. ( 53) in terms of the resultant yields the equation

0 = V 2 (V + 1) -2(8 + 9V )gG k 2 -27(gG k 2 ) 2 . ( 90
)
Up to some monomial factors, its discriminant has the form of a product P bp • (P pl ) 3 of a polynomial P bp determining the surfaces at which G 2 obeys the critical behaviour of branched polymers and P pl yielding planar critical behaviour. The transition line between the two regimes is then determined by the common solution P bp = P pl = 0. A third method is to determine critical loci of G 2 (g, t j ) as extrema of the inverse function t k (G 2 , g, t j ) for some t k ∈ {t j }. While this is straightforward when equations are linear in the couplings, for polynomial equation like in our case a clever rescaling of couplings has to be guessed in order to find an explicit inversion, if possible at all. Indeed, Eq. ( 88) can be appropriately inverted after rescaling t j = λ j t and g = λt. Then V (G 2 , t j ) = t Ṽ (G 2 , λ j ) and the equation takes the form

0 = Ṽ 3 t 2 + (9λG k 2 -Ṽ ) 2 t -108(λG k 2 ) 2 t -16λG k 2 . ( 91
)
such that

t ± = t ± [G 2 , λ, Ṽ (G 2 , λ j )] = 1 2 Ṽ 3 (9λG k 2 -Ṽ ) 2 -108(λG k 2 ) 2 ± (λG k 2 -Ṽ ) 1 2 (9λG k 2 -Ṽ ) 3 2
(92)

where t -is the physical branch with t -(G 2 = 1, λ, Ṽ ) = 0 . The third method is particularly useful to determine multi-critical loci with exponents

δ ≡ 1 -γ = 1 q + 1 → γ = q q + 1 (93) from vanishing derivatives 0 = ∂t ∂ Ṽ = ∂ 2 t ∂ Ṽ 2 = ... = ∂ q t ∂ Ṽ q . ( 94
)
In contrast to the case with only melonic interactions [START_REF] Bonzom | Random tensor models in the large N limit: Uncoloring the colored tensor models[END_REF] where all equations are polynomial, in this case it is not obvious whether solutions to these systems of equations exist up to q being the number of couplings. Nevertheless, in explicit examples we do find multi-critical loci up to q = m + 1 for m ≥ 2 melonic couplings, that is, even including the planar coupling g. This means that planar couplings enrich also the multi-critical structure of the branched-polymer regime which is a new insight to the best of our knowledge.

In the phase space of the models considered here, there are no intersections of multi-critical lines with transition lines. Accordingly, the case of critical exponents

γ ≡ 1 -δ = p m + n + 1 (95)
as argued in [START_REF] Alvarez-Gaumé | A Proposal for strings at D > 1[END_REF] for matrix models and applied to tensor models [START_REF] Bonzom | Enhancing non-melonic triangulations: A tensor model mixing melonic and planar maps[END_REF] does not appear (in this relation, m is the number of fine-tuned branched-polymer-type couplings, n the number of fine-tuned planar couplings, and p ≤ m).

B. Coupling melons and a necklace

In this section we illustrate the phase structure of models with a planar coupling and a number of melonic couplings Eq. (88) for a particular example. At least three couplings are necessary to investigate both the transition from a planar to a branched-polymer regime as well as interference with multi-criticality in the latter one. The simplest non-trivial case is the tensor model with a quartic and sextic melonic and a quartic necklace interaction,

V = -G 2 + t 2 G 2 2 + t 3 G 3 2 . ( 96
)
On the surface

t pl 3 = t pl 3 (g, t 2 ) = 3 4 9g -2 3g(9g + t 2 ) 2 , G pl 2 = G pl 2 (g) = 2 3 √ 3g , (97)
The surprising finding is that in this example the γ = 3/4 point lies close to but not on the phase-transition line Eq. (99). Moreover, we have checked that the transition line and the multi-critical line Eq. ( 104) have no intersection. In particular, there is no critical point on the transition line with exponent γ = 2/4 = 1/2 as would follow from Eq. (95). Therefore, the exponents argued for with matrix-model methods on the grounds of intersections of a multicritical line with the line of phase transition [START_REF] Alvarez-Gaumé | A Proposal for strings at D > 1[END_REF] are not recovered here. For such a result it is necessary to further prove that such intersections actually exist. Our example shows that there are cases where they do not exist.

The properties of the two-point function G 2 detailed in this example are common to solutions of the general type of Eq. (88). We have checked combinations of various melonic interactions and, in particular, also type 3 K interactions instead of quartic necklaces. A peculiar aspect of the type 3 K interaction B 3 K is that it effectively contains already both planar and branched-polymer behaviour since the equations for the two-point function Eq. (85) are of the form

G 2 = P ns gG 3 2 + t 3 G 3 2 + t 2 6 G 6 2 , (105)
only with fixed balance of coupling g = 2t 3k , t 3 = 4t 3k and t 6 = 6t 3k . Nevertheless, the tree-like part dominates in this case and G 2 has critical exponent γ = 1/2 at the critical point. This can be easily freed adding melonic interactions of type 3 M and 6 M such that then t 3 = 4t 3k + t 3 M and t 6 = 6t 3k + t 6 M . Then there is again the full phase-space with γ = 1/3 phase transition between planar and branched-polymer regime. according to Eq. ( 105). The branched-polymer critical curve (magenta) extends here to arbitrary high λ, as does the planar critical curve for arbitrary small λ. The phase transition with γ = 1/3 critical behaviour occurs at λ = -10/9.

C. Two types of necklaces with independent couplings

There is an example of a tensor model with a planar and a branched-polymer phase which is conceptually even simpler: two or more types of necklace interactions. It has been considered partially in [START_REF] Bonzom | Large N limits in tensor models: Towards more universality classes of colored triangulations in dimension d ≥ 2[END_REF] for the case of quartic necklace interactions B 1i 2 N with identical couplings t 12 2 N = t 13 2 N and so forth. This results in geometries with branchedpolymer behaviour as soon as there are more then one type of necklaces.

Here we show that for independent couplings, e.g. t 12 2 N , t 13 2 N , there is only a very narrow regime of branched-polymer behaviour. The two-point function G 2 is the solution to the equation

G 2 (t 12 2 N , t 13 2 N ) = G 2 (t, s) = P ns (t G 2 2 ) + P ns (s t G 2 2 ) -1 (106)
where we use the relative coupling s = t 13 2 N /t 12 2 N and t = t 12 2 N (one can check that indeed all results are symmetric under s → s -1 ). Then we find that G 2 has the expansion

G 2 (t, s) = G bp 2 -32 2s + 5 -5 √ 2s 1/2 5s + 2 -5 √ 2s 1/2 s + 1 - √ 2s 1/2 t bp -t 1/2 + O t bp -t (107)
typical for the branched-polymer regime only for s 1 < s < s 2 with

s 1 = 5 4 3 - √ 5 ≈ 0.955 , s 2 = 1 5 √ 5 + 3 ≈ 1.047 (108)
and the critical values are

t bp = 1 32 s + 1 - √ 2s , G bp 2 = 4 √ 2s s + 1 - √ 2s -8 . ( 109
)
Again the expansion Eq. ( 107) is consistent with the solutions s 1 , s 2 for the transition points as the radiant in the numerator of the expansion coefficient is positive for s 1 < s < s 2 .

IV. CONCLUSIONS

We analysed the critical behaviour of tensor models for all quartic and sextic interactions in dimension 3 and 4, and computed the unique scalings leading to a defined 1/N expansion, with non-trivial leading order, i.e. with infinitely many contributions at leading order.

Most models present a branched-polymer phase, a 2d pure gravity (planar) phase, and an intermediate phase with positive critical exponent, in which limit infinitely many smooth fractal spheres -baby universes -proliferate to form cacti with infinitely many branching points. While this phase structure has, to our knowledge, always been recovered from models with a matrix interaction and a polymerization perturbation, we provide a model generating cacti of spheres of one or another colour converging towards the smooth fractal Brownian sphere in the continuum, except for a very small region of the phase space for which the interference between the two colours is so strong it breaks down to a branched-polymer phase. At the boundaries of this region, the two colors compete to produce proliferating baby-universes. Furthermore, while in the context of matrix models, the polymerization terms were necessarily nonconnected multi-trace interactions coming from higher order curvature perturbations to the Einstein-Hilbert action, in the present case, the considered models are all built from connected interactions, and can all be interpreted as D-dimensional discrete Einstein-Hilbert theories of gravity without higher order curvature perturbations.

We study explicitly the multi-critical behaviour of a random tensor model with one quartic matrix-like "necklace" interaction and polymerization "melonic" perturbations of orders 4 and 6 and see that instead of the entropy exponent γ = 2/4 we would expect from [START_REF] Alvarez-Gaumé | A Proposal for strings at D > 1[END_REF], we obtain the multi-critical branched-polymer exponent γ = 3/4 at the end of the multi-critical branched-polymer line, the planar coupling playing the role of an higher order polymer interaction (we would not be able to reach an order-3 multi-critical branched-polymer exponent from only two quartic polymer interactions and a sextic one). We stress that all the models introduced can be understood on the level of discrete Ddimensional spaces dual to the Feynman stranded diagrams, so that all the critical and multi-critical phases involved have interpretations in terms of statistical systems coupled to random D-dimensional lattices with very specific adjacency conditions.

The results in this paper rely on a bijection between the Feynman stranded diagrams and embedded diagrams, which generalizes the so-called intermediate-field representation. Throughout this bijection, the link with matrixmodels and two-dimensional discrete surfaces is made clear, and we therefore see that the non-planarity of the order 6 stranded interactions is replaced by a planar coloring rule on some cubic embedded vertices. In light of this, a more appropriate criterium for interactions candidates to generate more interesting continuum limits could be the lack of a planar embedded equivalent intermediate-field vertex from the bijection. are in one-to-one correspondence with the faces of color i of G. The interest of this correspondence is thus that we can use ribbon graphs and matrix models tools to identify the leading orders for random tensor models. From stranded graphs to embedded diagrams. A choice has been made for every interaction, of a pairing of every tensor T with a T , and this defines the white vertices in Tables I and II. To indicate the pairing in the tables, we have used the graphical convention [START_REF] Alvarez-Gaumé | A Proposal for strings at D > 1[END_REF]. The pairs (T, T ) in the stranded graph G correspond to the edges around the white vertices in Ψ(G). By following alternatively the propagators and the pairs (T, T ) in the stranded graph G, we partition the pairs (T, T ) into cyclically ordered groups (the propagators are crossed from the T to the T ). All the edges corresponding to pairs (T, T ) in the same group are connected to a common black vertex in Ψ(G). These edges are cyclically ordered, due to the ordering of the pairs in (T, T ) in the group. This defines the diagram Ψ(G).

If as a convention the tensor T (resp. T ) leaves on the left (resp. right) of the edge, when seen from the white vertex, then the propagators in G correspond to the corners around black vertices in Ψ(G), and they are followed counterclockwise when the propagators are crossed from the T to the T in G. The white vertices in Tables I and II render the internal colored structure of the stranded interaction : the corners around white vertices in M i correspond to strands of color i between tensors inside the interaction in G. Moreover, if a color i is not included in the set of colors of an edge, it means that an edge of color i goes between the two tensors of the corresponding pair (T, T ). Thus the corners around black vertices in M i also correspond to strands of color i between tensors in G. This explains why the faces around M i correspond to faces of color i in G. From embedded diagrams to stranded graphs. Given an embedded diagram M, replace each white vertex by the stranded representation, by matching the pairs (T, T ) with the corresponding edges. Be sure to respect the convention that the T (resp. T ) are on the left (resp. right) of the edge when seen from the white vertex. For every corner around a black vertex in M, add a propagator between the corresponding T and T in the stranded graph.

2 N

 2 (T, T ) = N a,b,c,d,e,f,g,h=1 T abcd Tabef T ghef Tghcd .

  Figure[START_REF] Gurau | The 1/N expansion of colored tensor models in arbitrary dimension[END_REF] is an example of a contribution to the 2-point function of a quartic rank-4 model. The corresponding diagrams Γ contain a certain number F c (Γ) of closed curves for each colour c, which we name faces. We denote

e g2N 4

 4 Tr B 12 2 N (T, T )+g3N 4 Tr B 13 2 N (T, T )+g4N 4 Tr B 14 2 N (T, T )

  All U (N )3 -invariant interactions of order 4 and 6 in rank-3 tensors. The non-local structure of the interaction is shown as stranded diagram, in the coloured-graph representations common in the tensor-model literature and in the intermediate-field representation explained below in section II.

  ,c,d,e,f,g,h,i,j,k,l=1

  An example of a stranded diagram and the corresponding intermediate-field embedded diagram are shown below,

  valent vertices d, e, and f have a coloured internal structure, and all incident edges carry all colours {123}, as shown in[START_REF] Tutte | A census of planar maps[END_REF].

  [START_REF] Aldous | The Continuum Random Tree III[END_REF] , M 13 , M 14 ) = dM 12 dM 13 dM 14 /Z 0 e -N 2

	2 Tr(M 2 12 +M 2 13 +M 2 14

  FIG.1. Left: Two-point function G2(t, s) for two type 2N quartic necklace interactions various coupling ratios s according to Eq. (106).The purple curve shows planar critical points valid for s < s1 and s > s2 while the magenta curve indicates the critical points of the branched-polymer regime for sc,1 < s < sc,2. Right: For comparison, two-point functions G2(t, λ) for type 3K and type 3M interactions in similar ratio of couplings t = t3 k and λ =

	1.7 G2								2.0	F
				s = s1	s = s2
										λ = 8
	1.6								1.8
	1.5	s = 2							1.6
	1.4						s = 0		1.4			λ = -10 / 9
	1.3								1.2
	1.2								1.0
															λ = -7
	0.02	0.03	0.04	0.05	0.06	0.07	0.08	t	0.00	0.02	0.04	0.06	0.08	0.10	0.12	t
	4t 3k +t 3 M
	2t 3k

  10 At these phase-transition values of s we find at critical coupling t pt 1 =

	√	5 + 3/96 and t pt 2 = 5/96
		G 2 (t, s 1,2 ) =	4 3	√	5 -1 -48(c 1,2 ) 1/3 t pt 1,2 -t	2/3 + O t pt 1,2 -t	(110)
	with expansion coefficients c 1 = 3 17	√	5 -38 and c 2 = 6/	√	5 -11/25.

To get rid of the symmetry factor, we can just consider stranded diagrams with one distinguished propagator.

We stress however that such tensor models can be of interest as Sachdev-Ye-Kitaev (SYK) like tensor models when adding a time parameter[START_REF] Kitaev | A simple model of quantum holography[END_REF][START_REF] Witten | An SYK-Like Model Without Disorder[END_REF][START_REF] Gurau | The complete 1/N expansion of a SYK-like tensor model[END_REF][START_REF] Bonzom | Diagrammatics of a colored SYK model and of an SYK-like tensor model, leading and next-to-leading orders[END_REF]. In some cases, it is believed that they could serve as toy models to study the quantum properties of black-holes[START_REF] Maldacena | Remarks on the Sachdev-Ye-Kitaev model[END_REF][START_REF] Polchinski | The spectrum in the Sachdev-Ye-Kitaev model[END_REF].

In many cases, there is a one-to-one correspondence between planar diagrams and families of plane trees which carry labelings[START_REF] Cori | Planar maps are well labeled trees[END_REF][START_REF] Schaeffer | Conjugaison d'arbres et cartes combinatoires aléatoires[END_REF]. This additional information is enough to escape the universality class of branched polymers.

Non-local refers here to the fact that the stranded interaction is non-connected.

In the sense that stranded interactions are connected.

In the quartic case, the intermediate-field vertex would have two incident edges, and it can just be replaced with an edge.

More precisely, the correspondence is similar to that described in the Appendix, but instead of choosing one given pairing of the (T, T ) for all the interactions in the graph, we choose the pairing locally for every subgraph as in Eq. (57), making sure that the paired tensors are those linked by a propagator with a G 2 blob. For technical details see[START_REF] Lionni | Colored discrete spaces: higher dimensional combinatorial maps and quantum gravity[END_REF], Sec. 3.2.6.

In order to stay concise, we do not explain in detail the correspondence for non-leading diagrams, as we only use the correspondence for leading-order diagrams in the following. We stress however that this correspondence can be extended to all orders.

This can be computed from[START_REF] Bender | The Number of Degree-Restricted Rooted Maps on the Sphere[END_REF].

2

For s < (3 -√ 5)/5 ≈ 0, 153 and s > 5(3 -√ 5)/4 ≈ 6.545 the expansion Eq. (107) is defined as well but describes singularities of non-physical (G 2 (0, s) = 1) solutions of Eq. (106).

ACKNOWLEDGEMENTS

We would like to thank S. Dartois, V. Bonzom, V. Rivasseau and J. Ben Geloun for useful discussions. The work leading to this publication was supported by the German Academic Exchange Service (DAAD) with funds from the German Federal Ministry of Education and Research (BMBF) and the People Programme (Marie Curie Actions) of the European Union's Seventh Framework Programme (FP7/2007-2013) under REA grant agreement n • 605728 (P.R.I.M.E. -Postdoctoral Researchers International Mobility Experience).

We typically count diagrams with one distinguished edge, to which an orientation is given. The reason is that this breaks symmetries, as mentioned previously, and therefore leading order diagrams can then be counted more easily. Cutting that distinguished edge, one obtains a contribution to the 2-point function, and conversely. Counting rooted diagrams therefore goes back to counting leading-order contributions to the 2-point function.

Non-separable embedded diagrams [START_REF] Goulden | Combinatorial Enumeration[END_REF] are such that no face has two corners incident to the same vertex. The distinguished edge is either a cut-edge of colour i, in which case there is one 2-point function at each extremity, and the contribution is t i G 2 2 , or it belongs to a leading-order non-separable planar component, in which case all the edges have the same colour set {1i}. In each corner on each vertex is a (possibly trivial) rooted leading-order 2-point function. We denote P ns the generating function of rooted non-separable planar embedded diagrams,

where c n is the number of non-separable embedded diagrams with n edges. It is given by the following parametrization [START_REF] Goulden | Combinatorial Enumeration[END_REF][START_REF] Tutte | A census of planar maps[END_REF],

In an embedded diagram, the number of corners is twice the number of edges. Because we need to add one 2-point function G 2 per corner of the non-separable component containing the root, we need to replace the parameter z counting edges of the non-separable planar component with t 1i 2 N G 2 2 . We therefore obtain the following equation for the leading-order two-point function G 2 , which depends on the coupling constants

The reduction of this equation to the case with a single necklace and melonic interactions has been investigated in [START_REF] Bonzom | Large N limits in tensor models: Towards more universality classes of colored triangulations in dimension d ≥ 2[END_REF], as well as the case where all necklaces are given the same coupling. In Section III, we study in detail the critical regimes when considering quartic necklaces with independent couplings, and show that the same critical regimes are involved.

B. The sextic case in dimension three

For K 3,3 -type interactions the total number of faces of leading-order Feynman diagrams is

as shown in [START_REF] Bonzom | Colored triangulations of arbitrary dimensions are stuffed Walsh maps[END_REF] where we recall that B is the number of interaction bubbles, so that when putting a N 3 scaling in the action, corresponding to s = 1, the amplitude of leading-order Feynman diagrams contributing to the free energy scales as N 3 . Vacuum leading-order diagrams are such that bubbles are locally as in one of the two situations illustrated in Eq. ( 57), as proven in [START_REF] Bonzom | Colored triangulations of arbitrary dimensions are stuffed Walsh maps[END_REF],

in which the G 2 blobs indicate the insertion of a possibly trivial contribution to the leading-order 2-point function.

As explained previously, we count the diagrams contributing to G 2 . It corresponds to counting diagrams with one This is a very specific kind of hard-dimer model on leading-order quartic necklaces (it is similar to a dimer model on random lattices with black and white vertices such that distinguished edges only link black and white vertices, and such that every vertex is incident to a distinguished edge). Coupling a model to hard dimers is in general similar to considering multi-critical behaviour for models with higher order interactions. This suggests that the critical behaviour from interactions of type 3 N could in principle be obtained from the quartic necklaces regimes. Each one of the various models involving interactions of type 3 N and 2 N should correspond to a very specific kind of hard-dimer model on random lattices with coloured vertices and particular restrictions on how the distinguished edges can be, which could exhibit very interesting critical behaviour. We will however leave their study for further work.

III. MULTI-CRITICAL BEHAVIOUR OF THE MODELS

Enhanced coloured tensor models exhibit a rich multi-critical structure at leading order in the 1/N expansion. Critical behaviour of branched-polymers is the standard case in tensor models, but enhancing necklace-type interactions allows for planar regimes as well. When both types of interactions are present in the action, there are regimes in coupling space for each behaviour. The boundary between them obeys the special transition behaviour for which "cacti" dominate which are made of an infinite number of discrete spheres of diverging areas -baby universes-with tree-like incidence relations. These are characterized by a critical exponent γ = 1/3. Moreover, for several couplings the regimes have further internal structure of multi-critical subspaces. In this section we explain this structure for the models under consideration and illustrate it with explicit examples. As a new result we find that there is generically no multi-criticality on the boundary between the regimes induced by trees.

A. Universality transition between planar and branched-polymer regimes

All the models under consideration yield equations for the two-point function G 2 of the general form

where Q 0 , Q 1 , ... are a finite number of polynomials in G 2 and P ns is the generating function of non-separable planar diagrams, Eq. ( 53). (Strictly speaking, sextic necklaces B 3 N and the necklace-type interaction B 3 N have to be described by a modified generating function P

ns for non-separable rooted diagrams with restricted valency, or dually, faces restricted to hexagons. Nevertheless, there are no reasons to expect qualitative differences.)

We conjecture that solutions G 2 always have the same qualitative behaviour of a planar regime, a branched-polymer regime and a baby-universe transition phase. To this end we have analyzed the critical behaviour of the various models considered here. We provide some general arguments and illustrate this statement with a number of explicit examples in the following parts of the section.

Let us consider the general case of one planar coupling g and a finite number of melonic couplings {t j }. The corresponding equations for

where V covers a finite set of melonic interactions,

and k ≥ 2 results from the interaction of order 2k that is relevant for the planarities. Here, k = 2 for quartic necklaces Eq. ( 55), and k = 3 for type 3 K interaction Eq. (85). Heuristically, the generic two-phase structure of the models is obvious from Eq. ( 88): For the planar coupling g much larger than the melonic couplings {t j }, the V part can be neglected and the equation is simply the equation for planar diagrams [START_REF] Bonzom | Large N limits in tensor models: Towards more universality classes of colored triangulations in dimension d ≥ 2[END_REF]. On the contrary, for g much smaller than the {t j }, the equation reduces effectively to 0 = V (G 2 , t j ) for which the implicit-function scheme theorem [START_REF] Flajolet | Analytic Combinatorics[END_REF]Ch. VII.4] applies yielding critical branched-polymer behaviour.

We have used three methods to determine critical loci. Due to the definition Eq. ( 53) of P ns , the equations for G 2 are parametric. The straightforward way to determine critical loci is given by the vanishing of the Jacobian determinant together with the original equations. This is again a system of equations which has to be solved for G 2 , we find γ = -1/2 critical behaviour

The coefficient of the fractional-exponent expansion term already indicates, that this expansion breaks down for t 2 < t pt 2 where

for g > 0. Indeed, this is the solution we find for the line of phase transition between planar and branched-polymer regime. We have checked the expansion around the transition line (g, t pt 2 , t pt 3 ), for example in t 3 ,

yielding critical behaviour with exponent γ = 1/3. In particular, using the Newton-polygon method [56, ch. VII.7] we find that there are no further multi-critical points on this line. Critical surfaces on which G 2 has γ = 1/2 critical behaviour are fourth roots of the equation

whose solutions are hard to handle explicitly. However, from the branching lines of the solutions we can find an explicit expression for the multi-critical line

Appendix A: Correspondence between stranded Feynman graphs and embedded diagrams

In this appendix, we summarize the one-to-one correspondence between stranded Feynman diagrams of random tensor models, and embedded diagrams, in the case of interactions of order 6. This correspondence generalizes the intermediate field representation for the quartic models [START_REF] Bonzom | Enhancing non-melonic triangulations: A tensor model mixing melonic and planar maps[END_REF][START_REF] Gurau | The 1/ N Expansion of Tensor Models Beyond Perturbation Theory[END_REF][START_REF] Delepouve | Borel summability and the non perturbative 1/N expansion of arbitrary quartic tensor models[END_REF]. It was introduced in [START_REF] Bonzom | Colored triangulations of arbitrary dimensions are stuffed Walsh maps[END_REF] in a more combinatorial language for interactions of any order, and further developed in [START_REF] Lionni | Colored discrete spaces: higher dimensional combinatorial maps and quantum gravity[END_REF]. See also [START_REF] Bonzom | Large N limits in tensor models: Towards more universality classes of colored triangulations in dimension d ≥ 2[END_REF] for a pedagogic presentation. Thanks to this correspondence, we can rewrite arbitrary random tensor models as multi-trace matrix models [START_REF] Bonzom | Colored triangulations of arbitrary dimensions are stuffed Walsh maps[END_REF][START_REF] Lionni | Colored discrete spaces: higher dimensional combinatorial maps and quantum gravity[END_REF].

Let us describe the bijection Ψ between stranded diagrams G with sextic interactions B 1 , . . . , B n , and embedded diagrams M defined as follows:

• The diagrams have white vertices as shown in Tables I and II, corresponding to the interactions B 1 , . . . , B n under consideration. As represented in these tables, the edges attached to these white vertices carry specific sets of colors I ⊂ {1, . . . , D}. The ordering of the edges around white vertices is fixed.

• The other extremities of the edges can only be attached to black vertices. There is no restriction on the number or the kind of edges attached to black vertices. The edges around black vertices are cyclically ordered.

If we take such an embedded diagram M = Ψ(G) and delete all the edges whose sets of colors does not include the color i, as well as all the isolated white vertices, we obtain a fatgraph M i equivalent to a ribbon-graph, whose faces