Symmetric configurations highlighted by collective quantum coherence

Recent developments in quantum gravity have shown the Lorentzian treatment to be a fruitful approach towards the emergence of macroscopic spacetimes. In this paper, we discuss another related aspect of the Lorentzian treatment: we argue that collective quantum coherence may provide a simple mechanism for highlighting symmetric configurations over generic non-symmetric ones. After presenting the general framework of the mechanism, we show the phenomenon in some concrete simple examples in the randomly connected tensor network, which is tightly related to a certain model of quantum gravity, i.e., the canonical tensor model. We find large peaks at configurations invariant under Lie-group symmetries as well as a preference for charge quantization, even in the Abelian case. In future study, this simple mechanism may provide a way to analyze the emergence of macroscopic spacetimes with global symmetries as well as various other symmetries existing in nature, which are usually postulated.


Introduction
In recent years there have been some developments which show that the Lorentzian treatment of quantum gravity can be very successful in describing spacetime. One is the success of causal dynamical triangulation (CDT) [1]. CDT takes into account a causality condition, distinguishing Lorentzian spacetimes from the Euclidean ones, and has been shown to generate macroscopic spacetimes similar to the de Sitter spacetime, which is a good approximation of our actual universe. This is in contrast with the Euclidean counterpart, in which the emergence of macroscopic spaces seems to be more non-trivial [2,3]. Another development is the argument by van Raamsdonk that quantum correlations are the essence of connectivity of classical spacetimes, which has been derived in the context of AdS/CFT correspondence [4]. This seems to imply that the quantum correlations of a model for quantum gravity are crucial for the emergence of macroscopic spacetimes. It was also argued that Lorentzian path integral has some advantages over the Euclidean one in describing the wave function of universe in the mini-superspace approximation [5].
In this paper, we discuss another related aspect of Lorentzian treatment in quantum gravity: we argue that quantum coherence may provide a mechanism for the preference of configurations which are invariant under part of the underlying symmetries. A number of models of quantum gravity assume that spaces (spacetimes) are generated from the dynamics of certain 'building-blocks'. The number of building-blocks of a macroscopic space is generally huge, and their dynamics would thus be complicated and chaotic. On the other hand, there exist various symmetries in nature such as gauge symmetries, Lorentz invariance, de Sitter symmetry, and so on. It is a difficult question to see how such symmetries can emerge in general from such probably chaotic dynamics of building-blocks. This would also be tightly related to the dynamics of the emergence of macroscopic spacetimes, since the de Sitter symmetry should be present almost instantly after the Big Bang (and also approximately in the present era), according to cosmological models [6]. Of course, it would be natural to impose that theories of quantum gravity should have underlying symmetries which can interchange building-blocks. However, it would be difficult to devise an explicit scenario of spontaneous breaking of the underlying symmetry to the existing symmetries in nature, because there exists a huge hierarchy between them. In this paper, we consider the problem rather in an opposite direction. Considering the recent emphasis on Lorentzian treatment in quantum gravity, we propose that quantum coherence may provide a mechanism to highlight symmetric configurations over generic non-symmetric ones.
The phenomenon itself is quite common in a broad range of physics. One can find a similar phenomena even in high school text books. One example is X-rays diffracted by a crystal. This system has strong peaks in particular directions, each of these peaks represents particular discrete translational symmetries of the crystal. Though our formalism presented in this paper is described in a different manner suited for our purpose, the aspect is very similar to the text book phenomenon: strong peaks appear in association with part of the underlying symmetries.
The main purpose of this paper is to argue that, due to collective quantum coherence, peak patterns like the ones mentioned above associated to parts of the underlying symmetries can actually occur in a large variety of systems, and we explicitly show the phenomenon in a model tightly related to quantum gravity. After presenting the general framework of the mechanism in a manner suitable for our purpose, we consider the grand partition function of the randomly connected tensor network [7,8,9] with a Lorentzian modification: the exponent of an integrand is adjusted to be purely imaginary. It is known that the grand partition function of the randomly connected tensor network has a similar expression as the exact wave functions [10] of the canonical tensor model [11,12], which is a tensor model [13,14,15] in the Hamilton formalism and studied as a model of quantum gravity. Though our final goal would be to use this mechanism in the context of the canonical tensor model, we restrict ourselves in this paper to general arguments and concrete demonstrations in simpler related situations and leave the application to the canonical tensor model for future study. This paper is organized as follows. In Section 2, we provide the general argument of the mechanism in a setup which is common in a broad range of systems with underlying symmetries. We qualitatively argue that collective quantum coherence highlights configurations invariant under part of underlying symmetries of a system, and discuss some conditions which enhance the phenomenon. In Section 3, we explain the model we consider as an example. It is defined by a Lorentzian modification of the grand partition function of the randomly connected tensor network, which uses a tensor acting on an N -dimensional vectorspace. In Section 4, we consider the N = 3 case of the model, and find strong peaks on SO(2) symmetric configurations. In Section 5, the N = 4 case is considered. In this case, it turns out that there are two possibilities of highlighted symmetries, SO(2) and SO(3). We observe strong peaks on the configurations with those symmetries similarly to the results for N = 3. We also observe charge quantization in the case of SO (2). The final section is devoted to a summary and future prospects. In the appendix, we show some details of the computations in the N = 3 case.

The mechanism
The setup of our mechanism is common in physics. We consider a quantity which can schematically be represented by where φ and Q represent sets of multiple variables. Such quantities appear as path integrals, partition functions, wave functions, etc., with corresponding appropriate expressions of φ, Q, S(φ, Q), and C dφ. For instance in the context of quantizing general relativity, a wave function of a universe could be expressed in this form by considering C dφ and Q as integration over metrics on manifolds and metrics on their boundaries respectively. Later in this paper, we consider a discrete model defined by an expression like (1). For our purpose, it is essentially important that the integrand of (1) takes complex values with non-trivial phases in general. Thus, our setup must have Lorentzian characteristics rather than Euclidean: the integrand of (1) should not take just positive real values.
The class of theories we consider are assumed to have group symmetries in the following sense. For convenience, let us discretely label the variables φ, Q as φ a (a = 1, 2, · · · , N ) and Q i (i = 1, 2, · · · , M ) respectively. In general it is also possible to consider continuous labels such as spacetime coordinates. Let us assume that group transformations of φ and Q are given by where R(g) andR(g) are representations of a group G for elements g. For our purpose, it is necessary that neither of the representations, R norR, are trivial. We do not assume R or R to be irreducible, and they are allowed to contain trivial representations in general. We assume that the integration and S are invariant as for any g ∈ G. In particular, the assumptions imply the invariance of Ψ as Now, an interesting physical question would be what the locations of maxima of |Ψ(Q)| 2 are. For the simplicity of the following discussions, let us assume S(φ, Q) to be real. This is common in physics, but the assumption may not be necessary as long as the integrand can take various complex values with non-trivial phases. To illustrate the mechanism, let us consider a crude approximation of (1). This approximation is based on the principle of stationary phase. The principle states that the major contribution of a highly oscillatory integral comes from points where the phase function (in our case S(φ, Q)) is stationary [16]. This approximation would thus be given by only picking up stationary points as where φ = φ σ (σ = 1, 2, · · · , n crit ) are the critical points defined by and A σ are the prefactors obtained by performing Gaussian integrations around each critical point 1 . Here, for simplicity, it is assumed that the critical points are isolated for generic values of Q. If M and N are sufficiently large and S(φ, Q) is not too simple 2 , expression 1 The expression with the sum over critical points can be sophisticated by a method based on the Picard-Lefschetz theory. For instance, see Section 3 of [17] and references therein for more details.
2 For instance, S(φ, Q) should not be a quadratic function of φ and Q.
(5) suggests that Ψ(Q) has the following properties. One is that, for generic values of Q, there would exist a substantial number, n crit , of critical points, and there would be no strong correlations among values of phase S(φ σ , Q) of different critical points. Therefore, Ψ(Q) will be suppressed by mutual cancellations among the contributions of all the critical points. It is also likely that, for generic values of Q, Ψ(Q) behaves almost randomly as a function of Q due to the uncorrelated values of S(φ σ , Q) among the critical points. These properties would make it unlikely that there exist sensible observables in such generic regions of Q. Though such situations would not be generally true for all the possible cases, the kinds of cancellations and randomness mentioned above are expected to naturally occur, when M and N (and n crit as well) are sufficiently large and S(φ, Q) is not too simple. 3 On the other hand, there exist values of Q where |Ψ(Q)| 2 seems to take larger values. Let us consider a value of Q, say Q H , which is invariant under a subgroup H of G: Then, for such a Q H , critical points exist along the trajectories of the group action H: This is an immediate consequence of the invariance in (3). The representation R(H) is generally reducible, and therefore a critical point φ σ may be contained in the trivial part of the representation. In that case, the critical point is usually an isolated single point (except if it happens to coincide with another critical point, which might or might not have some deeper reasons). Otherwise, they form a set of critical points, on which the group action H is nontrivially represented. For the simplicity of our discussions below, let us assume both H and G to be Lie groups. The discussions below can obviously be generalized to finite groups and a similar mechanism will hold if the orders of the groups are large enough. In the case of Lie groups, the set of critical points on which H acts non-trivially form a continuous set of critical points. Taking into account the differences of the two classes of critical points, Ψ(Q H ) can now be represented by where the first term denotes the contributions from the continuous critical points, while the latter contains the contributions from the isolated ones. Here A σ in the first term is determined by the Gaussian integrations over the transverse directions to the group trajectories. One important matter is that the phase function is constant along each group trajectory because The constancy of the phase similarly holds even when the prefactor A σ is included. This is actually an exact property beyond the Gaussian approximation above, because the whole system of φ is invariant under H for Q = Q H .
In (9), there exist two main reasons for which |Ψ(Q H )| 2 may become relatively large due to the continuous critical points. One is that phases are constant along each group trajectory. Therefore, continuous critical points on each trajectory contribute coherently to Ψ(Q H ). The other is that, when Q reaches Q H , some of the isolated critical points are connected by a group orbit generated by H to form a continuous set of critical points. Therefore, the number of critical points is enhanced from a finite number to a continuous infinite. No divergences are generally caused by that though for compact groups, as can be seen in the expression (9). Expression (9) suggests some conditions for large enhancement. One is that the trajectories of continuous critical points should occupy large subspaces in the space of φ. Another is that the number of sets of continuous critical points should not be too large to avoid mutual cancellations among them. These conditions may be rephrased in that there exist a small number (ideally one) of continuous sets of critical points whose orbits R(H)φ σ have large volumes. As for the dimension of the subgroup H, there seems to exist a tension between the following two factors, and it seems hard to say which factor wins. One is that configurations invariant under smaller dimensional subgroups are more common, because it is easier for configurations to satisfy the symmetry condition (7) for smaller dimensional subgroups. On the other hand, enhancement tends to become larger for larger dimensional subgroups, because the dimensions of group orbits in the space of φ are larger.
The preferences discussed above would be consistent with the present form of the actual universe. Comparing to the large supposed symmetry which can interchange 'building-blocks' of spaces, the existing symmetries in nature are extremely small. On the other hand, the representation space of the existing symmetries is the universe itself (for instance, consider the translational symmetry), and therefore the dimension of the representations is extremely large: it should be in the order of the number of 'building-blocks' of the universe.
Lastly, we want to comment on an essential difference between the Lorentzian treatment and the Euclidean one. In the Euclidean treatment, preferred configurations are usually obtained by minimization of S(φ, Q), often referred to as ground states. Here, whether preferred configurations are symmetric or not is just an outcome, but is not relevant in determining them. On the other hand, in the Lorentzian treatment, quantum coherence is the essence rather than values of S(φ, Q), and symmetries play important roles in determining preferred configurations. This seems to imply an important paradigm change from Euclidean to Lorentzian: The collective behavior of the critical points is more important than the value of a single special point.

The model
Our final aim for developing this formalism is to apply the highlighting mechanism discussed in the previous section to the canonical tensor model [11,12], a model of quantum gravity, to explore the possibility of emergence of macroscopic spacetimes. On the other hand, the mechanism itself will generally hold in various settings with underlying symmetries. Therefore, to concentrate on the mechanism itself, we restrict ourselves to show the validity of the mechanism in some concrete simple examples connected to the canonical tensor model: we consider the grand partition function of the randomly connected tensor networks [7,8,9] with a Lorentzian modification. In fact, this setting is not so far from the canonical tensor model [11,12], since its known exact wave functions [10] have similar expressions.
The grand partition function [7] of the randomly connected tensor network is defined by where the integration region C is tilted from a real plane to a complex one to make the integration convergent [7], and with short-hand notations, φ 2 = φ a φ a and Qφ 3 = Q abc φ a φ b φ c . Here Q abc is assumed to be real and totally symmetric.
In this paper, we rather discuss a quantity which is similar but has a Lorentzian modification of multiplying the exponent by i, where R N denotes the whole N -dimensional real space, and > 0 is a regularization parameter assuring the convergence of the integration for any real values of Q. In addition to the interest in the exact wave functions of the canonical tensor model, Ψ T (Q) may be regarded as the grand partition function of the randomly connected tensor network with imaginary parameters.
Note that Ψ T (Q) is obviously symmetric under the orthogonal group G = O(N ), and satisfies the conditions required in the previous section, where the representations are given by with R(g) a b being the matrix of the (N -dimensional) vector representation of G = O(N ) for g ∈ G. The regularization parameter will generally be taken to be small positive values, but we will not discuss the details of its vanishing limit → 0 + . This is because Ψ T (Q) may have some singular behaviors in the limit → 0 + , when some of the components of Q vanish. Analyzing the details of such behavior is not in the range of our main purpose of this paper: we are rather interested in showing the enhancement of (13) for symmetric Q, and, to see this, it is enough to take sufficiently small.
4 Simplest non-trivial case: N = 3 We will now discuss a concrete example of the mechanism at work. We will use the model defined by (13) in the case of N = 3. In this case (13) is invariant under transformations of the orthogonal group G = O(3), and the totally symmetric real tensor Q abc has 10 independent parameters. As mentioned above (7) we will now consider a subgroup H ⊂ G. The subgroup H = SO(2) will be used with the following representation of its Lie algebra: The generator of H = SO(2) can always be put into this form by performing an appropriate G = O(3) transformation. Demanding (7) to hold leads to the following requirement for Q: This restricts the tensor a lot, as there are only two independent parameters left: others = 0.
(12) is now written as The critical points can now easily be computed using (6). There are two isolated critical points and a continuous set of critical points which is invariant under the group action of H: So for x > 0, 2x > 3y or x < 0, 2x < 3y there is a continuous set of critical points with topology S 1 . 4 If we make r 0 large by taking x ∼ 0 within these regions with finite y, the continuous set of critical points will occupy a relatively large subspace in the space of φ, and |Ψ(Q)| 2 will become large.
To see this explicitly, we can compute the contributions of the critical points by performing the Gaussian integrations around them, corresponding to (9). For large r 0 we obtain where sgn(x) denotes the sign of x. The explicit calculation of this can be found in the appendix. The last term is the contribution of the continuous set of critical points (21), and clearly dominates over the isolated critical points for large r 0 (or x ∼ 0). The second term appears to have a pole at 2x − 3y = 0, but this is just a consequence of the approximate treatment with the Gaussian integration. Indeed there is no such a singular behavior of Ψ, as can be checked by an exact formula below.
Because of the simplicity of the current model, it is also possible to obtain an exact expression of the important contribution. The calculation is shown in the appendix, and the result is Here f (x, y, φ 3 ) is defined in (43), and λ is an arbitrary real number with sgn(λ) = sgn(y). This is obtained by considering a contour integration shown in Fig.4 to single out the contribution of a pole. This expression does not anymore depend on the regularization parameter , since the first term is a non-singular convergent integration with an exponential damping behavior of the integrand ∼ exp(−3yλφ 2 3 ). The second term is exactly the contribution of the continuous critical points for large r 0 , so the stationary phase approximation above seems to work really well for a large space of critical points.

Asymmetric configuration
We showed that for small x the behavior of Ψ(x, y) is dominated by the continuous critical points. We will now show that the symmetric configuration is the preferred configuration.
Let us first note that the situation of (21) is rather special. If Q does not have a certain symmetry, then such a continuous set of critical points will in general not exist. In the case of N = 3 it can actually be shown explicitly that any other additional non-symmetric terms in (18) will remove the continuous set of critical points in favor of a number of isolated points. If an asymmetric term is parametrized by z, taking the limit z → 0 moves these isolated points towards the circle defined by (21), and for z = 0 they are linked by a constant trajectory of S, which forms the circle of continuous critical points.
From this and earlier discussions we expect that adding an asymmetric term disrupts the symmetry and Ψ will be suppressed. The simplest way to disrupt the symmetry would be to allow Q 113 and Q 223 to take different values: others = 0.
The integral is then given by After performing the Gaussian integration over φ 1 and φ 2 one finds Here we have introduced x = x 1 +x 2 2 and z = x 1 −x 2 2 . This parametrization is mainly chosen because one can now identify the symmetric configuration by simply setting z = 0. For z = 0 this integral reduces to (43). Although it proves quite difficult to solve this integral analytically, it can be evaluated numerically quite well.
As can be seen in Fig.1, |Ψ| 2 clearly peaks at z = 0, which indicates that the system prefers symmetric configurations. We also studied the region x < 0 with y = −1 numerically, but have not observed any interesting behavior. These seem to confirm our expectations, as there are either no real solutions to (21) in that case, or the solution has very small r 0 .
The above numerical study is just one of the possible asymmetric terms, one more will now be investigated before making some more general (but local) statements. Instead of (24)  (22) or (23).
we now take All other terms (of course with the exception of permutation of indices) are again put to zero. The following integral now needs to be solved: The φ 2 integration can be done readily, which gives This double integration can be done numerically and the result is shown in Fig.2. The result is similar to Fig.1, supporting the generality of the phenomenon.

Perturbations
By adding certain asymmetric terms, we have seen that the symmetric configurations seem to be preferred. To be more general one can look at general small perturbations around the symmetric configurations to make local statements. Let the tensor given in (17) be denoted by Q 0 and its corresponding action, (18), by S 0 . Ψ(Q) can be written, at least formally, for small perturbations Q = Q 0 + δQ as Considering the first order term in δQ we get an integral like Note that any term which is odd in φ 1 or φ 2 will vanish, since the whole integrand is then odd in φ 1 or φ 2 . There is actually just one term (which is not present in the unperturbed function) which is not odd but also vanishes: This means that ∂Ψ/∂Q abc | Q=Q 0 δQ abc = 0 for any perturbations δQ transverse to the symmetric configurations. Therefore, Ψ(Q 0 ) is a local extremum with respect to the transverse directions. In general the second derivative will not be negative, as is what one would expect for a local maximum. However, in the region where the continuous critical points are important, Q = Q 0 becomes a local maximum of |Ψ(Q)| 2 with respect to the transverse directions. This can be qualitatively shown by assuming that the contributions of the continuous set of critical points φ σ h dominate over the other isolated contributions in the summation (9): where we used the short-hand notation δQφ 3 ≡ δQ abc φ a φ b φ c , and · · · H denotes taking an average over the group trajectory of H. With ∂Ψ/∂Q abc | Q=Q 0 δQ abc = 0, (32) proves |Ψ(Q)| 2 takes a local maximum at Q = Q 0 with respect to the transverse directions to the symmetric configurations. This qualitative proof of the relation between the classicality of φ and the maxima of |Ψ| 2 is generally valid for any symmetric configurations and N . Here, the existence of the imaginary unit in the exponent is essentially important for the appearance of the minus sign in (32), and therefore the relation should be considered to be a characteristic of the Lorentzian treatment. The stationary points of the action are often regarded as the 'classical path' of a system, and the above argument shows that in the 'classical limit' the system generally prefers symmetric configurations.

N = 4
In the case of N = 3 there was only one possible Lie subgroup which could function as a partial symmetry. In the case of N = 4, there are three possibilities: SO(2), SO(2) × SO (2) and SO(3). However, it will turn out that the SO(2) × SO(2) reduces to be either trivial or an SO (2). In what follows we will first consider SO(2) × SO(2) and SO(2).
By O(N ) rotations, the SO(2) × SO(2) generators for N = 4 can be put in the following representation: Here t 1 and t 2 are parameters with t 1 > 0 and t 2 ≥ 0, so (33) represents two independent generators in general. However, it turns out that in order to solve (16) there are only two possible values of t 2 : t 2 = 0 or t 2 = 2t 1 5 . Other values of t 2 lead to the trivial configuration (Q abc = 0). In each case, t 2 = 0 and t 2 = 2t 1 , (33) represents only one SO(2) generator, and therefore the symmetry is SO(2), not SO(2) × SO (2). This preference of 'charge quantization' is generally true for any N : one can find more solutions of symmetric configurations of Q by such 'charge quantization' restricting the freedom of Lie group generators.
For t 2 = 0 the following symmetric tensor solves (16): All other terms (of course with the exception of permutation of indices) are again put to zero. The same methods as before can now be used to find the critical points. The specific expressions are nontrivial and not very relevant for the discussion, but the critical points again admit a solution of the form φ 2 1 + φ 2 2 = r 2 0 , φ 3 = constant, φ 4 = constant. This is also what we expect, since the critical points need to admit an SO(2) symmetry represented as the rotations of φ 1,2 .
For t 2 = 2t 1 there is a more restricted set of parameters: All other terms (of course with the exception of permutation of indices) are again put to zero. This gives, besides some isolated critical points, also a set of continuous critical points satisfying φ 2 1 + φ 2 2 = 2(φ 2 3 + φ 2 4 ) = r 2 0 . We have performed some numerical investigations under some perturbations around the symmetric configurations for the above cases, too. We have obtained similar results as the N = 3 case: strong peaks on the symmetric configurations.
The other possible symmetry is SO (3). It turns out that the solution is really similar to SO(2) for N = 3. SO(3) has three independent generators: (36) (16) now gives the following nonzero components (up to permutations of the indices): with all the other terms being zero. The critical points now have the topology of a 2-sphere S 2 , φ 2 1 + φ 2 2 + φ 2 3 = r 2 0 , obtained by a straightforward extension of (21). The quantity (13) can also be reduced to a single integral: where the only relevant difference from the N = 3 case, (43), is the exponent of the denominator. One can also evaluate certain asymmetric configurations by generalizing (37) to Q 114 = Q 224 = Q 334 = Q 114 parameterized by x 1,2,3 , respectively. This will again give just a single integral which is reasonably computable numerically: Again, setting x 1 = x 2 = x 3 in this parametrization gives the SO(3) symmetric configurations. Numerical calculations show that this is indeed the preferred configuration, as given in Fig.3.
In the figure, one can also observe lower peaks emanating in three directions, z 1 = 0, z 2 = 0, and z 1 = z 2 , respectively, which are the SO(2) symmetric configurations.
The N = 4 case demonstrates an interesting aspect of the mechanism explained in section 2. Though the critical points of SO(3) are expected to have a larger individual contribution to Ψ (in other words; a higher peak) because the space of the critical points is larger, the configuration space with SO(2) symmetry is always larger. This interplay between the size of the critical points space and the size of the configuration space with a certain symmetry is expected to be important in the large scale behavior of, for instance, the canonical tensor model.

Summary and future prospects
An important unresolved question in quantum gravity is how a macroscopic spacetime like the universe can emerge from a, probably chaotic, microscopic theory. A few recent developments suggest the importance of Lorentz treatment in such emergence [1,4,5]. In view of this, we have proposed a simple mechanism which highlights partial symmetries based on collective quantum coherence. Starting from a widely used setup with underlying symmetries, we have shown that configurations symmetric under part of the underlying symmetries are highlighted by collective quantum coherence over generic non-symmetric configurations. Qualitative arguments show that configurations invariant under partial symmetries with large dimensional representations are preferred. This seems to be in accord with the present form of the universe: relatively small dimensional symmetries like the Poincare/de Sitter symmetries represented on the whole universe. To demonstrate this mechanism, we have shown the occurrence of the phenomenon in some simple computable examples in the randomly connected tensor network [7,8,9], which is closely related to the canonical tensor model [11,12], a model of quantum gravity. We have observed strikingly sharp highlighting of configurations invariant under either Abelian or non-Abelian symmetries and also observed the phenomenon of charge quantization even in the Abelian case. Whether macroscopic spacetimes can be generated through the mechanism in the canonical tensor model will be an immediate future study.
We have put our mechanism in a familiar form, treating it in a rather general way. However, this does not imply that it is directly applicable for just any theory. Formally this could for instance be placed in a naive quantum gravity context by considering the metric at a certain time-slice as the configuration (called Q in Section 2) and the histories of the universe as integration variables (called φ in Section 2), where the diffeomorphism invariance is the underlying symmetry. However, the classical equation of motion corresponding to the critical equation (6) cannot be solved in a fully consistent manner for the Minkowski signature, because solutions will generally encounter singularities at certain spacetime locations. This makes the stationary phase analysis unreliable, and one would have to treat the system in a welldefined quantum mechanical manner from the beginning. In fact, there have been some recent developments in the asymptotic safety program [19,20], which studies gravity as a nonperturbative quantum field theory. The program is usually studied in a Euclidean context, but can often be Wick rotated to Lorentzian signature [21]. Considering technical difficulties immediate results would not be expected, but there would be possibilities for our effect to appear in this context.
Applying the mechanism to CDT also seems to be a non-trivial exercise, as the symmetries of the Regge action are usually not that simple, and after performing a Wick rotation the weights of the partition function are real instead of imaginary. The causality condition, restricting the triangulations to be compatible with a space-time decomposition, makes this Wick rotation rigorous and is considered to be tightly related to the success of CDT. The necessity of such restriction may be compared to what one would find after some analytic continuation of the integral (1) to make the weights real. 6 According to the mechanism, one would expect that symmetric configurations would obtain large weights, while non-symmetric configurations would be suppressed. It would also be nice to look into the Regge calculus from a Lorentzian viewpoint: the symmetries of the classical solutions 7 are enhanced for flat geometries or geometries with a small positive cosmological constant [22]. The mechanism would suggest the preference of such geometries with small curvatures.
Though the mechanism might not be always directly applicable in its current form, the mechanism does seem to be more reliable if the configuration and integration variables (φ and Q) are more independent, so that the form of (1) is directly satisfied. This could for instance be a statistical system with a sufficient amount of variables or a quantum mechanical wave function with an integral representation which has a similar form as (1). This is why the mechanism works well in the case of the random tensor network discussed in Section 3 , and we expect the mechanism to work well for the canonical tensor model.
An important characteristic of our mechanism is that peaks can appear only when there exist collective coherent phenomena along orbits of critical points. Therefore, a large peak inevitably contains a large number of degrees of freedom contributing to form it. Due to this collective nature, statistical tools could be useful to further analyse this mechanism, and characterise these preferred configurations. This might give a connection to the thermodynamic description of spacetime discussed in the literatures [23,24,25,26].

Acknowledgements
The work of N.S. is supported in part by JSPS KAKENHI Grant Number 15K05050. The work of D.O. is supported in part by the Hendrik Mullerfonds. N.S. would like to thank Y. Sato and B. Fraser for some initial discussions during his stay in Chulalongkorn University. We would also like to thank S. Hirano and H. Nielsen for some intriguing discussions in the occasions of their visits to Yukawa Institute for Theoretical Physics.

A Detailed calculations for N = 3
The contribution of the isolated critical points given in (19) and (20) can be obtained by perturbing S(φ σ + δφ) up to second order and performing Gaussian integration: Here S σ ab = ∂ 2 S ∂φ a ∂φ b φ=φ σ , and α σ l are the eigenvalues of the matrix S σ ab . For the evaluation of the square root we choose the principal square root, taking the branch cut along the whole negative real axis.
Next we evaluate the contribution from the continuous set of critical points (21). Since these critical points are degenerate, it is convenient to integrate over the degeneracy first. For this, polar coordinates, where φ 1 = r cos(θ), φ 2 = r sin(θ), are most convenient since the critical points are degenerate along the angular direction. Plugging these coordinates into (18) gives a trivial integration over θ, removing the degeneracy. The remaining Gaussian integrals can be evaluated to be: where δφ 3 = φ 3 − φ σ 3 , δr = r − r 0 , and sgn(x) denotes the sign of x. Here the standard formulas, iab 2 ) have been used. In the last step the r 0 → ∞ limit was taken in which we expect the continuous critical point to become most relevant. This limit can be realized by taking x → 0 with sgn(x) = −sgn(y), as can be seen in (21).
Because of the simplicity of the current model, it is also possible to do part of the calculation exact. The integral which needs to be solved is: where the − φ 2 term was added to make sure the integral converges as mentioned in (13). The φ 1 and φ 2 integrals can be done readily since they are just Gaussian integrals, and using ∞ −∞ dke iak 2 = i π a one finds Interestingly, the integrand has a simple pole at φ 3 = − 1+i x . This value of φ 3 indeed agrees with that of the continuous critical point found in (21). This suggests that the contribution of the pole is related to the continuous critical point, as will be shown more explicitly below.
The contribution of the pole of the integration of (43) can be picked out by taking the contour of Fig.4 in the complex φ 3 plane. Here, λ is an arbitrary real number with the same sign as y, the reason of which will become apparent shortly. The residue theorem then gives: So, by taking L → ∞ on account of sgn(λ)=sgn(y), the integral (43) can be rewritten: Ψ(x, y) = iπ Here note that there is no in the expression: we have been able to safely take the → 0 + limit, because the integrand of the first term damps exponentially ∼ exp(−3yλφ 2 3 ) and no regularization is necessary. By comparing with (41), it can be seen that the continuous critical point is the source of the contribution of the pole for r 0 → ∞ (or x ∼ 0).  (45), where λ must have the same sign as y for the limit L → ∞ and → 0 + . There is a pole at φ 3 = − 1+i x shown by * , and only for sgn(x) = −sgn(y) the pole will lie inside the contour. The orientation of the contour depends on the sign of λ, generating the factor of sgn(λ) in (44).