arxiv:hep-th/ v1 3 Feb 1993

Transcription

1 NBI-HE-9-89 PAR LPTHE 9-49 FTUAM 9-44 November 99 Matrix moel calculations beyon the spherical limit arxiv:hep-th/93004v 3 Feb 993 J. Ambjørn The Niels Bohr Institute Blegamsvej 7, DK-00 Copenhagen Ø, Denmark L. Chekhov L.P.T.H.E., Universitè Pierre et Marie Curie 4, pl. Jussieu, 755, Paris Ceex 05, France C. F. Kristjansen The Niels Bohr Institute Blegamsvej 7, DK-00 Copenhagen Ø, Denmark Yu. Makeenko Institute of Theoretical an Experimental Physics B.Cheremushkinskaya 5, 759 Moscow, Russia Abstract We propose an improve iterative scheme for calculating higher genus contributions to the multi-loop (or multi-point) correlators an the partition function of the hermitian one matrix moel. We present explicit results up to genus two. We evelop a version which gives irectly the result in the ouble scaling limit an present explicit results up to genus four. Using the latter version we prove that the hermitian an the complex matrix moel are equivalent in the ouble scaling limit an that in this limit they are both equivalent to the Kontsevich moel. We iscuss how our results away from the ouble scaling limit are relate to the structure of mouli space. Permanent aress: Steklov Mathematical Institute, Vavilov st.4, GSP Moscow, Russia

2 Introuction The hermitian one matrix moel is interesting for many reasons. The matrices involve being N N, the moel can be solve in the limit when N, as was first emonstrate in Ref. []. Furthermore the moel amits a /N -expansion. This expansion is the simplest example of the t Hooft topological (or genus) expansion [] in the given case of a D = 0 imensional quantum fiel theory. If the leaing orer (planar or genus zero approximation) result for some quantity, say the vacuum energy, is normalize to be of orer, the genus g contribution is of orer O(N g ). Even though the hermitian matrix moel is often consiere a toy moel, its perturbative expansion reveals the combinatorial problems one will encounter in the multi-imensional case: The sum of graphs contributing to a given orer of the genus expansion is convergent, while the attempt of summing over all genera leas to a ivergence. This was one of the original motivations for stuying the genus expansion of the hermitian matrix moel. One more application of the hermitian matrix moel is ue to the fact that it escribes [3] iscretize -imensional ranom surfaces an hence D quantum gravity. In this language a iagram of genus g is equivalent to a (piecewise linear) surface of genus g. The hermitian one matrix moel was extensively stuie from this point of view after the observation [4] that it might not only escribe pure D gravity but also some minimal conformal fiel theories couple to D quantum gravity. The relation given in [4] for genus zero was extene to any orer of the genus expansion in the socalle ouble scaling limit, where the coupling constants approach some critical values in a N-epenent way [5]. In the ouble scaling limit many interesting approaches have been evelope in orer to eal with higher genus contributions. The technique of exactly integrable systems (KV an KP hierarchies) was applie to the moel in Ref. [6]. It provies an algorithm for genus by genus calculations in the ouble scaling limit. Another approach propose by Witten [7] is base on the interpretation of the ouble scaling limit of the matrix moel as a topological theory. The problem is then reuce to calculating intersection inices on mouli space. In this approach it was possible to obtain explicit results in the ouble scaling limit up to genus four [8]. A further evelopment of this approach was the representation [9] of the generating function for the intersection inices as a hermitian one matrix moel in an external fiel (the Kontsevich moel). For the hermitian matrix moel away from the ouble scaling limit, the original genus zero solution [] was obtaine by solving the large-n integral sale point equation for the spectral ensity. To exten the results to higher genera, the orthogonal Although the relation conjecture in [4] turne out to be wrong, the paper was nevertheless instrumental for the later evelopments.

3 polynomial technique was applie in Ref. [0]. While this metho works well for the partition function in the case of simple potentials an lower genera, it is ifficult to implement if one wants an explicit calculation of multi-point correlators an if the potential is generic. Another metho of solving the hermitian matrix moel is base on the loop equations [] which were first written own for the hermitian matrix moel in Ref. []. In this approach one solves the genus zero equation an uses the solution in an iterative genus by genus proceure. In the scheme originally propose in Ref. [] an elaborate in Ref. [3], one etermines the higher genera contributions by solving algebraically the higher loop equations an imposing on all correlators the same analyticity structure as that of their genus zero versions. Following this line of action genus one correlators for φ 4 an φ 6 potentials were calculate in ref. [4]. However, this metho involves the entire chain of loop equations an is in practise only applicable when one only wants to calculate a few terms in the genus expansion for a potential with a small number of terms. An alternative iterative proceure for calculating higher genus contributions to the correlators of the the hermitian matrix moel has been propose recently [5]. It is base only on the first in the chain of loop equations, is applicable for an arbitrary potential an can be pursue orer by orer of the genus expansion without any problems. A key point in this scheme is a change of variables to the so calle moments of the potential. The scheme is close in spirit to the iterative solution of the unitary an hermitian (with cubic potential) one matrix moels in an external fiel [6] 3. The genus one partition function an correlators were calculate for an even potential in Ref. [8] an for an arbitrary potential in Ref. [5], where it was also proven that genus g contribution to the n-loop correlator epens on at most (3g + s) lower moments. An analogous algorithm for the complex matrix moel was evelope in Ref. [9] an the genus an 3 correlators were explicitly calculate for an arbitrary potential. In the present paper we propose an improve algorithm for an iterative calculation of higher genus contributions to the multi-loop correlators an to the partition function of the hermitian one matrix moel. We perform calculations up to an incluing genus two. Our algorithm iffers from that of Ref. [5] by a reefinition of the moments, which has the avantage that it allows us to evelop a version which gives irectly the result in the ouble scaling limit where we perform explicit calculations up to an incluing genus four. We compare the analysis to that of the complex matrix moel an prove in particular that the hermitian an complex matrix moels are equivalent in 3 For the Kontsevich moel similar moments has been use in Ref. [7]. 3

4 the ouble scaling limit. We formulate a limiting proceure which allows us to obtain the Kontsevich from the hermitian matrix moel an use it to prove in a very irect manner the equivalence between the ouble scaling limit of the hermitian matrix moel an the Kontsevich moel. We iscuss a possible interpretation of our results for the hermitian matrix moel as representing a particular iscretization of mouli space which preserves basic geometrical properties such as intersection inices, an which allows for a etaile stuy of the bounary of mouli space [0]. We conjecture that the moment-technique evelope in the present paper might be an efficient tool in the stuy of this fine structure of mouli space. The paper is organize as follows: In Section we introuce the basic ingreients of the iterative proceure, the loop insertion operator, the loop equation an the new variables the moments. In Section 3 we escribe the iterative proceure an present the results for the average of loop operators an the partition function for the hermitian matrix moel up to genus two. In Section 4 we emonstrate in etail how our iterative proceure works in the ouble scaling limit an calculate the correlators an the partition function up to genus four. In Section 5 we show that the Kontsevich moel can be obtaine as a certain limit of the hermitian matrix moel, prove the equivalence of the ouble scaling limit of the hermitian matrix moel an the Kontsevich moel an iscuss the possible connection between our results for the hermitian matrix moel away from the ouble scaling limit an the structure of the iscretize mouli space. Finally Section 6 contains a short iscussion of future perspectives an in the Appenix we evelop a ouble scaling version of the iterative proceure for the complex matrix moel propose in [9], an prove that the ouble scaling limits of the complex an hermitian moels coincie. Main efinitions The iterative proceure we are going to escribe is base on three main ingreients, the loop insertion operator, the first in the chain of loop equations of our moel an a suitable change of variables. It allows us to calculate the genus g contribution to the n-loop correlator for (in principle) any g an any n an (also in practise) for any potential. The possibility of going to an arbitrarily high genus is provie by the loop equation while the possibility of consiering an arbitrarily high number of loops is provie by the loop insertion operator. The change of variables allows us to consier an arbitrary potential. 4

5 . Loop insertion operator The hermitian one matrix moel is efine by the partition function Z[g, N] = e NF = N N where the integration is over hermitian N N matrices an V (φ) = Expectation values (or averages) are efine in the usual way as j= φ exp( N Tr V (φ)) (.) g j j φj (.) Q(φ) = Z φ exp( N Tr V (φ)) Q(φ). (.3) We introuce the generating functional (the -loop average) W(p) = N Tr φ k /p k+ (.4) k=0 an the s loop correlator (s ) W(p,...,p s ) = N s k,...,k s= Tr φ k... Tr φ ks conn. /p k +...p ks+ s (.5) where conn. refers to the connecte part. One can rewrite the last two equations as follows W(p,...,p s ) = N Tr s p φ... Tr. (.6) p s φ conn. As is explaine in Ref. [4], these quantities are associate with the (Laplace images of the) sum over iscretize open surfaces with s bounaries. The correlators can be obtaine from the free energy, F, by application of the loop insertion operator, V (p) : W(p,...,p s ) = V (p s ) V (p s )... F V (p ) (.7) where V (p) j. (.8) p j+ g j j= Notice that Eq. (.7) can be rewritten as W(p,...,p s ) = V (p s ) V (p s )... V (p ) W(p ) (.9) 5

6 which shows that if the -loop correlator is known for an arbitrary potential, all multiloop correlators can be calculate. This is why it is name above as the generating functional. With the normalization chosen above, the genus expansion of the correlators reas Similarly we have W(p,...,p s ) = g=0 N g W g(p,...,p s ) (s ). (.0) F = for the genus expansion of the free energy. g=0 N gf g (.). The loop equation The first in the chain of loop equations for the hermitian matrix moel can conveniently be written as [] C ω πi V (ω) p ω W(ω) = (W(p)) + NW(p, p) (.) where V (ω) = j g j ω j /j an C is a curve which encloses the singularities of W(ω) an not the point ω = p. This contour integration acts as a projector picking up the coefficients of p. Due to Eq. (.9) the secon term on the right han sie of the loop equation (.) is expresse via W(p), so that (.) is a close equation which etermines this quantity. To leaing orer in /N one can ignore the last term in Eq. (.) an assuming that the singularities of W(ω) consist of only one cut on the real axis an that W(p) behaves as /p as p one fins [] W 0 (p) = C ω πi V { } / (ω) (p x)(p y) (.3) p ω (ω x)(ω y) where x an y are etermine by the matrix moel potential in the following way C ω πi V (ω) (ω x)(ω y) = 0, C ω πi ωv (ω) (ω x)(ω y) =. (.4) This genus zero solution will be use below in the iterative proceure of solving the loop equation. Inserting the genus expansion (.0) in Eq. (.) it appears that W g (p), g, obeys the following equation { ˆK W0 (p) } W g (p) = g g = W g (p) W g g (p) + 6 V (p) W g (p) (.5)

7 where ˆK is a linear operator, namely ˆKf(p) C ω πi V (ω) f(ω). (.6) p ω In Eq. (.5) W g (p) is expresse entirely in terms of W gi (p), g i < g. This is what makes it possible to evelop the iterative proceure mentione in the Introuction..3 The new variables To characterize the matrix moel potential we introuce instea of the couplings g j the moments M k an J k efine by M k = J k = C C ω πi ω πi V (ω) (ω x) k+/ (ω y) / k, (.7) V (ω) (ω x) / (ω y) k+/ k. (.8) These moments epen on the coupling constants g j s both explicitly an via x an y which are etermine by Eq. (.4): M k = g k + ( x + (k + )y) g k+ +..., (.9) J k = g k + ( (k + )x + y) g k (.0) Notice that M k an J k epen explicitly only on g j with j k. There are several motivations for introucing these new variables. First, as we shall see below, for each term in the genus expansion of the partition function an the correlators, the epenence on the infinite set of coupling constants will arrange into a simple function of a finte number of the moments. Moreover, these new variables reflect more irectly than the coupling constants the possible critical behaviour of the matrix moel. Let us briefly escribe how this comes about. Performing the contour integral in (.3) by taking resiuals at ω = p an ω = one fins W 0 (p) = { } V (p) M(p) (p x)(p y) (.) where M(p) is a polynomial in p of egree two less than V (p). As alreay iscusse in the Introuction, W 0 (p) can also be etermine by an analysis of the matrix moel in the eigenvalue picture []. Requiring that W 0 (p) is analytic in the complex plane except for a branch cut [y, x] on the real axis an behaves as /p for p correspons to requiring that the eigenvalue ensity, ρ(λ), has support only on the interval [y, x] an is normalize to. The eigenvalue ensity is in this situation given by ρ(λ) = π M(λ) (λ y)(x λ) y λ x. (.) 7

8 This function vanishes uner normal circumstances as a square root at both ens of its support. Critical behaviour arises when some of the roots of M(λ) approach x or y. For a non-symmetric potential the so calle m th multi-critical point is reache when (m ) extra zeros accumulate at either x or y. Comparing (.3) an (.) it appears that M k (k ) M(λ) λ k, J k (k ) M(λ) λ=x λ k (.3) λ=y so the conition for being at the m th multi-critical point is simply M = M =... = M k = 0 an M k 0 if the extra zeros accumulate at x an J = J =... = J k = 0, J k 0 if the extra zeros accumulate at y. For the symmetric potential (V (φ) = V ( φ)) the eigenvalues will be istribute symmetrically aroun zero an hence confine to an interval of the type [ z, z] in the one arc case. In a situation like this the m th multi-critical point is characterize by the eigenvalue ensity having (m ) extra zeros at both z an z. Reexamining (.7) an (.8) it appears that for the symmetric potential M k = ( ) k+ J k, (V (φ) = V ( φ)). (.4) Hence in this case we have only one set of moments an the conition for being at the m th multi-critical point is the vanishing of the first (m ) of these. This formalism obviously allows for a treatment of a more general situation where m extra zeros accumulate at x an n extra zeros accumulate at y. Such multi-critical moels have been stuie in Ref. []. We will however restrict ourselves to the traitional moels. The superiority of the moments efine in (.7) an (.8) as compare to the coupling constants will become even more clear when we consier the ouble scaling limit in Section 4.. By then it will also become evient why these new moments are more convenient than those originally introuce in Ref. [5]. 3 The iterative proceure Our iterative solution of the loop equation results in a certain representation of the free energy an the correlators in terms of the moments. In this section we escribe the structure of W g (p) an F g an show that the iterative proceure can be conveniently formulate by referring to the eigenvectors of a linear operator, ˆK. We prove that the genus g contribution to the s-loop correlator epens at most on (3g + s) lower moments ( (3g ) for the partition function) for g > 0. We perform explicit calculations up to genus two. 8

9 3. The structure of F g an W g (p) Our iterative proceure of solving the loop equation results in the following representation of the genus g contribution to the free energy F g = α j >, β i > α...α s ; β...β l α, β, γ g M α...m αs J β...j βl M α J β γ g, (3.) where = x y is the istance between the enpoints of the support of the eigenvalue istribution, the brackets enote rational numbers an α, β an γ are non-negative integers. The inices α,...,α s, β,..., β l take values in the interval [, 3g ] an the summation is over sets of inices. In particular F g epens on at most (3g ) moments. Furthermore, since nothing allows us to istinguish between x an y, F g must be invariant uner the interchange of the two. Hence one gets F : α...α s ; β...β l α, β, γ = ( ) γ β...β l ; α...α s β, α, γ. (3.) There are certain restrictions on the integers which enter Eq. (3.). Let us enote by N M an N J the total powers of M s an J s, respectively, i.e. N M = s α, N J = l β. (3.3) Then it hols that N M 0, N J 0 (3.4) an F g : N M + N J = g, (3.5) F g : si= (α i ) + l j= (β j ) + γ = 4g 4. (3.6) The relation (3.5) follows from the fact that the partition function Z = e g N g F g is invariant uner simultaneous rescalings of N an the eigenvalue ensity, ρ(λ); N k N, ρ(λ) ρ(λ). The relation (3.6) follows from the invariance of Z uner rescalings of k the type N k N, g j k j g j. Finally the following inequality must be fulfille: F g : s l (α i ) + (β j ) 3g 3 (3.7) i= j= This requirement becomes more transparent when we consier the ouble scaling limit in Section 4.. In combination with Eq. (3.6) it gives g γ 4g 4. (3.8) 9

10 To explain the structure of W g (p), let us introuce the basis vectors χ (n) (p) an Ψ (n) (p): where χ (n) (p) = M { Ψ (n) (p) = J { } n x (p) χ (k) (p)m n k+ Φ (n) Φ (n) y k= n (p) k= Ψ (k) (p)j n k+ } n, (3.9) n (3.0) Φ (n) x (p) = (p x) n {(p x)(p y)} / n 0, (3.) Φ (n) y (p) = (p y) n {(p x)(p y)} / n 0. (3.) It is easy to show for the operator ˆK efine by Eq. (.6) that { ˆK W0 (p) } χ (n) (p) = { ˆK W0 (p) } Ψ (n) (p) = (p x) n n, (3.3) (p y) n n (3.4) an that the kernel of { ˆK W0 (p) } is spanne by Φ (0) x (p) = Φ(0) y (p). Since W g (p) can be obtaine from F g accoring to Eq. (.7), the representation (3.) implies 4 W g (p) = 3g n= { A (n) g χ (n) (p) + D g (n) Ψ (n) (p) } g. (3.5) We o not a any multiple of Φ (0) x (p) or Φ(0) y (p). Doing so woul contraict the bounary conition W(p) /p for p since this behaviour was alreay obtaine for genus zero. We note that this structure of W g (p) is in agreement with the assumption [3] that W g (p) is analytic in the complex plane except for a branch cut [y, x] on the real axis. The coefficients A g (n) are of the same structure as F g an the relation (3.4) still hols. However in this case the inices α,..., α s, β,...,β l take values in the interval [, 3g n]. Hence W g (p) epens on at most (3g ) moments. The invariance of the partition function uner the rescalings escribe above has the following implications for the structure of A (n) g : A (n) g : N M + N J = g, (3.6) A g (n) : si= (α i ) + l j= (β j ) + γ = 4g n. (3.7) 4 A similar formula is proven in Ref. [5] for a ifferent efinition of the moments an of the basis vectors. 0

11 We also have an analogue of (3.7) for A (n) g. It reas A (n) g : k l (α i ) + (β j ) 3g n. (3.8) i= j= Again we refer to Section 4. for further explanations. However, we note that combining (3.7) an (3.8) one gets again the boun (3.8) on γ. As was the case for F g, W g (p) must be invariant uner the interchange of x an y. This means that D (n) g must appear from A (n) g by the replacements, J M. s.) Furthermore, W g (p) (We note that we o not have a relation like (3.) for the A (n) g must be an o function of p for a symmetric potential. This implies that A (n) g = ( ) n D (n) g. That the structure of F g an W g (p) actually is as escribe in this section can be proven by inuction. We will not go through the proof here. Instea we refer to Refs. [5, 9]. In the latter reference a formula somewhat similar to (3.5) was proven for the complex matrix moel. However the strategy of the proof will be evient from Sections 3. an 3.3 where we escribe the iterative proceures which allow us to calculate W g (p) an F g for any g starting from W 0 (p). 3. The iterative proceure for etermining W g (p) Accoring to Eq. (.), we nee to calculate W 0 (p, p) in orer to start the iterative proceure. To o this we write the loop insertion operator as V (p) = V (p) + x V (p) x + y V (p) y (3.9) where V (p) = j. (3.0) p j+ g j The erivatives x/v (p) an y/v (p) can be obtaine from (.4) an rea j= x V (p) = M Φ () x, y V (p) = J Φ () y. (3.) Using the relation one fins [3] V (p) V (ω) = p W 0 (p, p) = p ω, (3.) (x y) 6(p x) (p y). (3.3)

12 This enables us to fin W (p) an we see that it is of the form (3.5) with A () = 8, D () = 8, A() = 6, (3.4) D() = 6. (3.5) Carrying on the iteration process is straightforwar. In each step one must calculate the right han sie of the loop equation (.). Decomposing the result obtaine into fractions of the type (p x) n, (p y) n allows one to ientify immeiately the coefficients A g (n) an D g (n). To calculate W g (p, p) it is convenient to write the loop insertion operator as V (p) = n M n + V (p) M n j J j + x V (p) J j V (p) x + y V (p) y (3.6) where M n V (p) = (p x) n / (p y) 3/ (n + /)Φ (n+) x + ( n ( ) n i M i x y i= (p) ) n i+ ( ) n + ( ) n J x y y V (p) x +(n + /)M n+ V (p). (3.7) The erivatives x/v (p) an y/v (p) are given by (3.) an the function Φ (n) x was efine in (3.). Of course J n /V (p) just appears from M n /V (p) by the replacements y x an J M. We note that there is no simplification of the algorithm in the case of the symmetric potential. We can only put x = y = z at the en of the calculation. This complication stems of course from the fact that we have to keep the o coupling constants in the loop insertion operator until all ifferentiations have been performe. Only hereafter they can be put equal to zero. The same is not true in the ouble scaling limit however. We will come back to this later. By taking a closer look at the loop insertion operator (3.6) an bearing in min the results (3.3), (3.4) an (3.5), it is easy to convince oneself that A (n) g an D g (n) epen only on x an y via (x y) an have the structure shown in equation (3.). Furthermore it appears that W g (p,...,p s ) epens on at most (3g + s) parameters. The

13 results for g = obtaine with the ai of Mathematica rea A () = A () J 67 J J 5 J J + 49 J J 57 J J M 8 4 J + 49 M M M M + M 8 4 J M J M 64 3 J M + 67 M M 5 M M, 57 8 J = J M 8 3 J 49 M M 56 4 M M 3 M 8 3 J M + J M 8 J M 67 M 8 3 M M M 3, A (3) = 49 M 4 56 M 5 M M + 67 M 8 3 M M J 8 3 J M 8 J, M A (4) = 49 M 3 8 M M, A (5) = M ; D () = A () (M J, ), D () = A () (M J, ), D (3) = A (3) (M J, ), D (4) = A (4) (M J, ), D (5) = A (5) (M J, ). (3.8) The genus two contribution to W(p) is now etermine by Eq. (3.5). 3.3 The iterative proceure for F g In this section we present an algorithm which allows one to etermine F g, as soon as the result for W g (p) is known. The strategy consists in writing the basis vectors χ (n) (p) an Ψ (n) (p) as erivatives with respect to V (p). It is easy to verify that the following relations hol χ () (p) = x V (p), (3.9) 3

14 Ψ () (p) = χ () (p) = Ψ () (p) = y V (p), (3.30) { V (p) 3 ln M } 3 ln, (3.3) { V (p) 3 ln J } 3 ln. (3.3) Combining this with the results (3.4) an (3.5) one immeiately fins which coincies with the expression of Ref. [5]. F = 4 ln M 4 ln J ln (3.33) 6 In the general case things are not quite as simple. The basis vectors can not be written as total erivatives. This is of course in accorance with the fact that the A an D coefficients now have a more complicate epenence on the potential (cf. Section 3. ). However, a rewriting of the basis vectors allows one to ientify relatively simply W g (p) as a total erivative. In the case of χ (n) (p) this rewriting reas where Φ (n) x χ (n) (p) = M { n ( ) n i n M n n shoul be written as } ( ) n i Φ (i) y x M i i= V (p) x y } n V (p) χ (k) M n k+, n (3.34) k= Φ x (n) = { n ( ) n i n i= +M n x V (p) n } ( ) n i y x M i V (p) x y Φ (i) M n V (p) n, (3.35) Φ () x = M x V (p). (3.36) The basis vector χ () (p) shoul of course still be written as in (3.9). The rewriting of the Ψ (n) s is analogous to that of the χ (n) s. It can be obtaine by performing the replacements J M an x y in the formulas above. 4

15 By means of these rewritings we have been able to etermine F. The result reas F 9 = 7680 J M + 8 J J 3 8 M M J + 64 J M 3 M 3 64 J M J 3 40 J 4 M 40 M M 3 9 M J 3 9 J 3 + J M 64 J M 7 8 J M 4 + J 3 60 J 5 9 J J 3 8 J J J 3 M 3 60 M M M 3 8 M 4 35 M M 3. (3.37) It is obvious from the formulas above that F g epens for a non-symmetric potential on at most (3g ) moments. Furthermore, for a symmetric potential, F g is a sum of two ientical terms an epens on only at most 3g moments. A consequence of this oubling for the ouble scaling limit will be iscusse in Section The ouble scaling limit It is easy to etermine which terms in the explicit expressions for F g an W g (p) etermine in the previous section that contribute in the ouble scaling limit. However it is rather time consuming to etermine F g an W g (p) away from the ouble scaling limit. In this section we evelop an algorithm which gives us irectly the result in the ouble scaling limit. Using this algorithm we calculate the correlators an the partition function explicitly up to genus four an escribe their general structure. 4. Multi-critical points Let us consier first the case of the symmetric potential. We hence have x = y = z an J k = ( ) k+ M k for all values of k. As is mentione earlier, the m th multi-critical point is characterize by the eigenvalue ensity having (m ) extra zeros accumulating at both z an z, an the conition for being at this point is the vanishing of the first (m ) moments. For the m th multi-critical moel the ouble scaling limit of the correlators is obtaine by fixing the ratio of any given coupling an, say g, to its critical value an setting The moments then scale as p = z c + aπ, (4.) z = z c aλ /m. (4.) J k M k a m k, k [, m ]. (4.3) 5

16 Furthermore, it is well known that the genus g (g ) contribution to the free energy has the following scaling behaviour F g a ( g)(m+/). (4.4) Bearing in min that the structure of F g is as shown in Eq. (3.), one fins that the following relation must hol s l (m α j ) + (m β i ) (α + β)(m ) m( g) g. (4.5) j= i= Since the free energy shoul look the same for all multi-critical moels, we have N J + N M g, (4.6) s l (α j ) + (β i ) 3g 3. (4.7) j= i= We alreay know from the analysis in Section 3. that the equality sign must hol in (4.6). Only terms for which max(s, l) m an for which the equality sign hols also in (4.7) will contribute in the ouble scaling limit. From Eq. (3.6) it follows that these terms will have γ = g. (4.8) A similar analysis can be carrie out for the generating functional. Here it is known that the genus g contribution to W(p) has the following scaling behaviour W g (π, Λ) a ( g)(m+/) (4.9) with the exception of W 0 (π, Λ) which also contains a non-scaling part. From (3.9) an (3.0) it appears that the basis vectors scale as χ (n) Ψ (n) a m n+/. (4.0) Remembering that the structure of the A an D coefficients is as shown in (3.), one fins that the following relation must hol for both A (n) g an D (n) g s l (m α j ) + (m β i ) (α + β)(m ) m( g) + n g. (4.) j= i= again provie max(s, l) m. Since also the generating functional shoul look the same for all multi-critical moels, it follows that A (n) g conitions an D (n) g must satisfy the following N J + N M g, (4.) 6

17 s l (α j ) + (β i ) 3g n +. (4.3) j= i= Here we recognize the homogeneity requirement (3.8). We know from Section 3. that the equality sign always hols in (4.). Only terms for which max(s, l) m an for which the equality sign hols also in (4.3) will contribute in the ouble scaling limit. We note this means that γ = g for these terms (cf. Eq. (3.7)). Let us turn now to the case of the non-symmetric potential, an let us consier the m th multi-critical point assuming that the extra zeros accumulate at x. To obtain the ouble scaling limit of the correlators, we first fix as before the ratio between any given coupling an the first one to its critical value. Then we scale p an x in the following way keeping however y fixe: p = x c + aπ, (4.4) x = x c aλ /m. (4.5) Uner these circumstances the J-moments o not scale but M k behaves as M k a m k k [, m ] (4.6) Using again (4.4) we fin that the only terms which survive the ouble scaling limit for F g are those for which N M = g, (4.7) s (α j ) = 3g 3. (4.8) j= Comparing Eq. (4.7) with Eq. (3.5), we see that all the J-epenence isappears when the prescription for taking the ouble scaling limit is as in (4.4) an (4.5). Furthermore, comparing (4.8) with (3.6) we fin that all remaining terms have γ = g as in the symmetric case. The relations (4.7) an (4.8) will be exploite in Section 5.. Let us carry out the analysis of W g (p) also for the non-symmetric case since this will give us an aitional information about the structure of the D an A coefficients. For the basis vectors we have χ (n) a m n+/, Ψ (n) a /. (4.9) By exploiting again the known scaling behaviour (4.9) of W g (p), one can just as in the symmetric case erive certain homogeneity conitions that the A an D coefficients must fulfill. One fins A (n) g : N M g, s (α j ) 3g n ; (4.0) j= 7

18 D (n) g : N M g, s (α j ) 3g. (4.) j= The conitions are ifferent for the A an D coefficients ue to the ifferent scaling behaviour of the basis vectors. Terms that are important in the ouble scaling limit are only those for which s m an for which the equality sign hols in both relations in (4.0) or both relations in (4.). Comparing the secon relation in (4.) with (4.3) we see that all D-terms vanish in the ouble scaling limit. Furthermore, comparing (4.0) with (3.6) we see that all J epenent A terms isappear. Hence in the ouble scaling limit everything is expresse only in terms of the M s an, an all terms have γ = g. As is mentione earlier, D g (n) can be obtaine from A (n) g by performing the replacements x y an J M Hence we also have that the following homogeneity requirements shoul be fulfille A (n) g : N J g, l (β j ) 3g ; (4.) j= D (n) g : N J g, l (β j ) 3g n. (4.3) j= These relations coul of course be erive by analyzing the scaling behaviour assuming that the extra zeros accumulate at y. This implies fixing x an replacing (4.4) an (4.5) with p = y c aπ, (4.4) y = y c + aλ /m. (4.5) In this situation only terms for which l m an for which the equality sign hols in both relations in (4.) or both relations in (4.3) will contribute in the ouble scaling limit. For symmetry reasons, now all A terms plus M-epenent D terms isappear. Also the free energy is expresse uner these circumstances entirely in terms of the J s an. We note that inepenently of the etails of the prescription for taking the ouble scaling limit, in this limit the multi-loop correlator W g (p,...,p s ) epens on at most 3g +s moments an the free energy F g epens on at most 3g moments. 4. Determining W g (p) in the.s.l. The homogeneity conitions erive in the previous sections allow us to etermine whether a given term contributes to the ouble scaling limit or not. As an example, let us consier the case of a non-symmetric potential where the critical behaviour is associate with the enpoint x. Using (4.0) we see that only A (5), the first term 8

19 of A (4) an the first two terms of A (3) of the long list of complicate expressions for the D an A coefficients of W (p) given in Section 3. survive in the ouble scaling limit. Obviously it woul be convenient to have an algorithm which gives as output only terms which are relevant for the continuum limit. An algorithm which gives us all potentially relevant terms can be obtaine by a slight moification of the one escribe in Section 3.. By potentially relevant we mean relevant for m multi-critical moels with m sufficiently large. It is easy to convince oneself that terms which are relevant for the ouble scaling limit of W g (p) can only appear from terms which are relevant for the ouble scaling limit of the right han sie of the loop equation. Let us assume that we know the scaling relevant versions of A (n) G an D(n) G for G =,...,g an let us assume that we want to calculate the ouble scaling limit of the right han sie of the loop equation for W g+. First of all we note that all terms in the basis vectors χ (n) (p) an Ψ (n) (p) show the same scaling behaviour so none of them can be ignore. However we can replace all occurrences of (p y) with (x c y) = c in this limit. It is known from earlier analysis [4] that for a m multi-critical moel W g (p, p) scales as W g (π, π) a g(m+/). (4.6) Using equation (4.9) it is easy to show that all the proucts in the sum in (.) have the same scaling behaviour. Hence, no superfluous terms appear from the sum if we start from the ouble scale versions of the W G, G =,...g. However, if we calculate W g (p, p) by letting the loop insertion operator as written in (3.6) act on W g (p), irrelevant terms will appear even if we start from ouble scale version of W g (p). By comparing (4.9) an (4.6), we see that only operators in /V (p) which lower the power of a by m + 3/ give rise to the relevant terms. Therefore, we shoul iscar all other operators. By examining carefully each term in (3.6), one fins that only the following part of the loop insertion operator remains in the ouble scaling limit: where M n V (p) x V (p) V (p) x = n M n + x V (p) M n V (p) x { = (n + /) Φ (n+) x (p) M } n+ Φ () x M (p) (.s.l.) (4.7) (.s.l.), (4.8) = Φ x (p) (.s.l.) (4.9) M an Φ (n) x (p) = (p x) n { c (p x)} / (.s.l.). (4.30) 9

20 We are now in a position to calculate the ouble scaling limit of the right han sie of the loop equation for W g+ provie we know the ouble scale versions of W (p),..., W g (p). We see that, since all the y-epenence has isappeare, we o not have to perform any ecomposition of the result. This is in accorance with the outcome of the analysis of the previous section that all D terms isappear in the ouble scaling limit. There it was also foun that all J-epenent A terms woul vanish. This also appears from the formulas above. No J terms ever appear if we o not start out with any an we o not. The starting point of the iteration proceure is of course as before W 0 (p, p) but now we shoul take only the part of it that contributes in the ouble scaling limit. From equation (3.3) we fin 5 W (NS) 0 (p, p) = 6 (.s.l.). (4.3) (p x) There is an interesting feature of the A coefficients. For all genera it hols that A () A () g g = = 0. This is because the smallest negative power of (p x) appearing in any basis vector is 3/. By multiplication of two basis vectors or by application of the loop insertion operator this power will be lowere by at least 3/. Hence terms of the type (p x) an (p x) never turn up on the right han sie of the loop equation. We have calculate the A coefficients as they look in the ouble scaling limit for g =, 3 an 4 following the iterative proceure outline above. The results for g = an g = 3 rea A () = 0, A () = 0, A (3) = 49 M 56 c M 4 5 M 3 3 c M 3, A (4) = 49 M 8 c M 3, A (5) = A () 3 = 0, A () 3 = 0, c M, 3 M 3 A (3) 3 = 5355 M c M M M M 3 c M c M + (680 M 3 M M M 5 ) 048 c M 6 55 M 6 5 c M 5, M M 4 5 Here an in the following we use the notation that the superscript (S) refers to the case of the symmetric potential an the superscript (NS) to the case of the non-symmetric potential (where the critical behaviour is associate with the enpoint x). 0

21 A (4) 3 = 5355 M 4 56 c M M M 3 5 c M M M M c M M c M 5, A (5) 3 = 737 M 3 56 c M M M c M M c M 5, A (6) 3 = 6495 M 048 c M M c M 5, A (7) 3 = 303 M 5 c M 5, A (8) 3 = c M 4. (4.3) We remin the reaer that the terms liste above are only potentially relevant. Whether or not they are actually relevant epens on which multi-critical moel one wants to consier. For a m th multi-critical moel all terms involving M k, k > m vanish in the ouble scaling limit. This is true whether one uses a minimal potential or not. We also remin the reaer that we assume that we ha a non-symmetric potential an that the critical behaviour was associate with the enpoint x. In the case where the critical behaviour is associate with the enpoint y all the formulas in this section still hol provie c is replace with c, x with y an M with J. Let us turn now to the case of the symmetric potential. In view of the scaling relations (4.) an (4.) it might seem unnatural to work with terms like (p x) an (p y) an we will see below that there actually exists a way of avoiing this. However, for the moment we will analyze the scaling behaviour in the x, y formalism by (formally) scaling p to z c whenever it occurs in a term like (p y) an to + z c whenever it occurs in a term like (p x). Now let us assume as in the non-symmetric case that we know the scaling relevant versions of A (n) G an D(n) G for G =,..., g an that we want to calculate the ouble scaling limit of the right han sie of the loop equation for W g+ (p). Of course we still have that all terms in the basis vectors show the same scaling behaviour so we still have to keep all of them. However it is easy to convince oneself that in the ouble scaling limit all occurrences of (p y) in χ (n) (p) can be replace by c an all occurrences of (p x) in Ψ (n) can be replace by c. Furthermore, it appears that when we calculate the sum of proucts on the right han sie of the loop equation all proucts that mix χ s an Ψ s will be subominant in the ouble scaling limit. These mixe proucts shoul hence be iscare. Finally to ensure the survival of only ouble scaling relevant terms in W g (p, p), the loop insertion operator shoul be written as V (p) = V (p) x + V (p) y (4.33)

22 where /V (p) y can be obtaine from /V (p) x given in equation (4.7) (4.30) by replacing x by y, M by J an c by c. The starting point for the iterative proceure is of course the ouble scale version of W 0 (p, p) which here reas W (S) 0 (p, p) = 6 (p x) + (.s.l.). (4.34) 6 (p y) We see that the loop equation actually ecouples into two inepenent equations. Each of these is a ouble scale version of the loop equation for the non symmetric potential. One correspons to the case where the critical behaviour is associate with the enpoint x, the other to the case where the critical behaviour is associate with the enpoint y. The A g (n) s are hence equal to those obtaine in the non-symmetric case where the critical behaviour is associate with the enpoint x an D g (n) = ( ) n A (n) g. This means that the right han sie of the loop equation of the loop equation can be written as the r.h.s. = 3g n= n ca (n) g (4.35) (p z) n which seems to be the natural way of expressing it bearing in min the scaling relations (4.) an (4.). as As a consequence, it becomes also more natural to express the generating functional where Ã(n) g fulfill 3g W g (S) (p) = Ã (n) g χ (n) (p) (4.36) n= = A g (n) n c an where we now have only one set of basis vectors χ(n) (p) which { ˆK W0 (p) } χ (n) (p) = It is easy to show that χ (n) (p) shoul be chosen as where χ (n) (p) = { } Φ (n) n z (p) χ (k) (p)m n k+ c k n M k= Φ (n) z (p) = (p z) n. (4.37) (4.38) (p z) n+/. (4.39) Furthermore, one can easily convince oneself that the right han sie of the loop equation can be obtaine irectly in the form (4.35) if one replaces (p x) by c (p z) in the formulas (4.7) (4.30), i.e. if one carries out the iteration process starting from W S 0 (p, p) = c 6(p z) (.s.l.) (4.40)

23 an using the following expression for the loop insertion operator. V s (p) = n M n + z V (p) M n V (p) z (.s.l.) (4.4) where M n V (p) z V (p) { = (n + /) Φ (n+) z (p) n c M } n+ Φ () z M (p) c (.s.l.), (4.4) = c Φ () z (p) (.s.l.). (4.43) M This observation allows us t show explicitly that the hermitian an the complex matrix moel are equivalent in the ouble scaling limit. 4.3 Determining F g in the ouble scaling limit The starting point of the calculation is the ouble scale version of W g (p) obtaine as escribe in the previous section an the strategy is again to rewrite the basis vectors in a form which allows one to ientify W g (p) as a total erivative. However this time the rewriting in the case of χ (n) is mae with the ai of (4.8) instea of (3.7) an reas χ (n) (p) = { } n M n M n V (p) χ (k) (p)m n k+ k= n. (4.44) We o not nee any expression for χ () (p) since the sum in (4.44) starts at k = an since we know that in the ouble scaling limit A () g = 0, neither we o nee any expression for Φ x (p). The relevant rewriting of Ψ (n) (p) appears from (4.44) when M is replace with J. For the non-symmetric potential where the critical behaviour is associate with the enpoint x, we immeiately fin for genus one F (NS) = 4 ln M (.s.l.). (4.45) This expression can of course alternatively be obtaine by taking the ouble scaling limit of Eq. (3.33). Before presenting the explicit results for F g (NS) in genus, 3 an 4, let us escribe the general structure of F g (NS) as it appears from our iterative solution of the loop equation. We have F (NS) g = α...α s α, γ herm g α j > M α... M αs M α γ (.s.l.) g (4.46) 3

24 where the rational numbers enote by the brackets are ientical to those of Eq. (3.) having the same inices. The α j s obey Eq. (4.8), γ is given by Eq. (4.8) an as it follows from Eq. (4.7). α = g + s (4.47) Let us consier the g = case as an example of what structures are encoe by Eq. (4.46). For g = we have α j [, 3g ] = [, 4] an the relation (4.8) gives s (α j ) = 3g 3 = 3 g =, (4.48) j= so that s [, 3g 3] = [, 3]. The restriction (4.48) amits s = 3 an α = α = α 3 =, s = an α, α =, 3 an s = an α = 4. The proper powers of M in the enominator for such terms are given by α = g + s = + s. The outcome of the iteration process is F (NS) = M 3 60 c M M M 3 8 c M 4 35 M c M 3. (4.49) We note that all allowe values of s an α i s actually appear. By taking a closer look at the loop insertion operator an the basis vectors it is easy to convince oneself that this will be the case for all genera. The expression (4.49) coul of course also have been obtaine by taking the ouble scaling limit of (3.37) following the recipe given in (4.7) an (4.8). We see that the expression for F simplifies rastically in the ouble scaling limit. an The results for g = 3 (where s 6) an g = 4 (where s 9) rea F (NS) 3 = 05 M 6 56 c M M 4 M 3 56 c M M M 3 5 c M M 56 c M 8 3 M M M 3 M M 3 3 c M c M 33 M M M 3 M 5 c M 5 6 c M + 45 M 4 96 c M M M 6 04 c M M c M 5 (4.50) F (NS) 4 = M c M M 7 M c M M c M 3 5 M M 6 M M 3 3 M 3 c M M 4 M 3 M 4 c M 04 3 c M M 5 M M 4 M 3 c M M M 3 M 4 c M c M 4

25 M 3 M c M M 3 M 3 M c M 5467 M 4 M c M M 3 3 M c M M M 3 M c M M M 3 M c M M M 4 M c M M M 3 M c M M c M M M c M M c M M 3 M 4 M c M M 3 M c M M M 4 M c M M M 3 M c M M M M M 9 c M M 5 M 6 c M c M 3 M M 4 M c M M 3 M c M M c M 7. (4.5) Neeless to say that the results for the case where the critical behaviour is associate with the enpoint y can be obtaine from these by performing the replacements M J, c c. It is obvious from the iscussion in Section 4. an formula (4.44) that inepenently of the etails of the prescription for taking the ouble scaling limit, F g epens in this limit on at most 3g moments. Again we stress that the terms liste above are only potentially relevant. For a m th multi-critical moel all terms involving M k, k > m vanish too. This is true whether one uses a minimal potential or not. It is interesting to note that for the symmetric potential we get a sum of two ientical terms an hence F (S) g = F (NS) g (.s.l.) (4.5) which is a well known property of the ouble scaling limit for the hermitian matrix moel [4]. It can be trace back to the fact that the loop equation for a symmetric potential ecouples into two inepenent equations. Equation (4.5) is the reason why the partition function of D quantum gravity is efine as the square root of the ouble scaling limit of the partition function (.) in the case of a symmetric potential. There is no such square root in the non-symmetric case an the proper continuum partition function can be obtaine just as the ouble scaling limit of (.). Of course we have etermine up till now only the coefficients, F g, of the genus expansion of the free energy (Cf. Eq (.)). For a m th multi-critical moel the relevant expansion parameter in the ouble scaling limit is Na m+/. If we efine moments µ k by (Cf. Eq (4.3)) M k = a m k µ k, k [, m] (4.53) 5

26 we get by replacing M k with µ k for k [, m] an setting M k equal to zero for k > m in the formulas above exactly the coefficient of the expansion in the ouble scaling parameter. A similar statement hols for the results of the previous section. 4.4 A remark on the complex matrix moel The complex matrix moel is efine by the partition function Z C = e N F C = φ φ exp( NV C (φ φ)) (4.54) N N where the integration is over complex N N matrices an V C (φ φ) = Its generating functional is efine by an the s-loop correlator (s ) by W C (p,..., p s ) = N s j= g j j Tr (φ φ) j. (4.55) W C (p) = Tr (φ φ) k /p k+ (4.56) N k=0 k,...,k s= Tr (φ φ) k... Tr (φ φ) ks ) conn. /p k +...ps ks+ (4.57) so that W C (p,...,p s ) = N Tr s p p φ φ... Tr p s. (4.58) p s φ φ conn. In Ref. [9] an iterative proceure which enables one to calculate Wg C (p,...,p s ) for any genus g an any s starting from W0 C (p) was evelope. This proceure is much the same as for the hermitian matrix moel. It is base on the loop equation, the loop insertion operator an a suitable change of variables. In Ref. [9] the stuies were carrie out only away from the ouble scaling limit. In the Appenix we show how the iterative scheme for the complex matrix moel can be moifie to give irectly the result in the ouble scaling limit. The moifie scheme provies us with an exact proof that the hermitian an the complex matrix moel are equivalent in the ouble scaling limit 6. We fin that W C g (p,...,p s ) = F C g = 4 g+s W(S) g (p,...,p s ) (.s.l.), (4.59) 4 g F(S) g (.s.l.) (4.60) 6 The fact that the hermitian an complex matrix moels belong to the same universality class an that the correlators therefore coincie in the ouble scaling limit oes not contraict Ref. [5] where it was shown that the iscrete Virasoro operators for the complex matrix moel o not have continuum limit. 6

27 inepenently of the type of critical behaviour consiere. In particular we see that the partition function of D quantum gravity can also be obtaine as the square root of the ouble scaling limit of the partition function (4.54) provie the integration is over N/ N/ matrices in stea of over N N matrices. This is of course in accorance with the fact that the number of inepenent components of a N/ N/ complex matrix is the same as that of a N N hermitian matrix. 5 Relation to mouli space The equivalence between the ouble scaling limit of the hermitian matrix moel an the Kontsevich moel has been proven in various ways in Refs. [6, 7]. We present in this section another proof irectly in terms of the moments introuce above. For this purpose we escribe a limiting proceure which allows us to obtain the Kontsevich moel irectly from the hermitian matrix moel. We iscuss also the connection between the hermitian matrix moel an a certain iscretization of mouli space. 5. The limiting proceure As we have seen in Section 4.3, the ouble scaling limit of the free energy is given by the expansion (4.46) which encoes in a remarkable way all multi-critical behaviour of the hermitian matrix moel. Since most of the terms originally present in the expression for F g vanish in the ouble scaling limit the question arises whether it is possible to formulate a matrix moel which gives solely the coefficients which are relevant for the ouble scaling limitṫhe answer to this question is yes an the proper moel is known as the Kontsevich moel. The Kontsevich moel is efine by the integral [9] n n Z K [M, n] = X exp { n tr ( )} MX + X3 6 n n X exp { n tr ( )}. (5.) MX there M=iag(m,...,m n ) (at first) is a positive efinite matrix an the integration is over hermitian n n matrices. The genus expansion of the free energy is given by ln Z K [M, n] = g=0 n g F Kont g. (5.) The goal of this section is to show that the Kontsevich moel can be obtaine from the hermitian matrix moel by a certain limiting proceure which looks conceptually ifferent from the ouble scaling limit. We shall see however in Section 5. that this limiting proceure leas to a result which is ientical to the one obtaine in the ouble scaling limit. 7

28 In orer to formulate the limiting proceure let us recall some ientities. Consier the hermitian matrix moel as efine in Eq. (.) with N = ξn. 7 following ientity [8] 8 We have the Z[g, N(ξ)] = e n tr Λ / Z P [Λ, n], N(ξ) = ξn (5.3) where the partition function Z P [Λ, n] is associate with the Penner moel [9], in the external fiel Λ =iag(λ,...,λ n ): Z P [Λ, n] = n n [ X exp n tr ( ΛX )] X ξ log X. (5.4) The ientification is vali provie the set of coupling constants {g} = {g 0, g,...} an the matrix Λ are relate by the Miwa transformation g k = n tr Λ k δ k k, g 0 = n tr log Λ. (5.5) In orer to show how the partition function (5.) can be obtaine from (5.4), we note that Z P [Λ, n] can be reformulate as the so-calle Kontsevich Penner moel [30]: Z KP [η, n] = n n X exp { ξn tr [ ηxηx + (log( X) + X)]} n n X exp { ξn tr [ ηxηx ]} (5.6) X provie Λ = ξ(η + η ). (5.7) We note in aition that the moel (5.6) possesses a remarkable connection to a iscretization of mouli space which will be iscusse in Section 5.4. Here we will be intereste in the connection to the usual continuum mouli space, i.e. to the Kontsevich moel which precisely has the interpretation as the generating functional for intersection inices on mouli space (Cf. to Section 5.3). In orer to highlight this connection we choose a specific parametrization of η an ξ: η = e εm, ξ = ε 3 (5.8) where m is to be ientifie with the matrix, M which appears in the efinition (5.) of the Kontsevich moel, an where ε is presently just an expansion parameter which 7 We shoul point out a subtlety relate to the representations (5.4) an (5.6). We have use a constant ξ in these equations which has the opposite sign of that conventionally use. The analytic continuation of ξ is neee in orer to make contact with the Kontsevich moel. We will not to try to justify it, since the iea is to provie some heuristic explanation of the observe ientity of the expansion coefficients of the hermitian matrix moel an the Kontsevich moel. 8 See also Ref. [8] where the proof was extene to finite n. 8