Date

2001-02-14

Description

Comment: 10 pages, 5 figures, Latex, psfig

We consider the model of random surfaces with extrinsic curvature term

embedded into 3d Euclidean lattice $Z^3$. On a 3d Euclidean lattice it has

equivalent representation in terms of transfer matrix $K(Q_{i},Q_{f})$, which

describes the propagation of loops $Q$. We study the spectrum of the transfer

matrix $K(Q_{i},Q_{f})$ on finite dimensional lattices. The renormalisation

group technique is used to investigate phase structure of the model and its

critical behaviour.

embedded into 3d Euclidean lattice $Z^3$. On a 3d Euclidean lattice it has

equivalent representation in terms of transfer matrix $K(Q_{i},Q_{f})$, which

describes the propagation of loops $Q$. We study the spectrum of the transfer

matrix $K(Q_{i},Q_{f})$ on finite dimensional lattices. The renormalisation

group technique is used to investigate phase structure of the model and its

critical behaviour.

Type

Identifier

doi:10.1142/S0217732301004182

Mod.Phys.Lett.A16:1069-1078,2001

Database

Link to record

Show preview

Hide preview

ar
X

iv :c

on d-

m at

/0 10

22 45

v1 [

co nd

-m at.

sta t-m

ec h]

1 4 F

eb 20

01

NRCPS-HE-2001-4

The Spectrum of the Loop Transfer Matrix on Finite Lattice

Georgios Daskalakis and George K.Savvidy

National Research Center Demokritos, Ag. Paraskevi, GR-15310 Athens, Hellenic Republic

Abstract

We consider the model of random surfaces with extrinsic curvature term embedded into 3d Euclidean lattice Z3. On a 3d Euclidean lattice it has equivalent representation in terms of transfer matrix K(Qi, Qf), which describes the propagation of loops Q. We study the spectrum of the transfer matrix K(Qi, Qf) on finite dimensional lattices. The renormalisation group technique is used to investigate phase structure of the model and its critical behaviour.

1 Introduction

Various models of random surfaces built out of triangles embedded into continuous space Rd and surfaces built out of plaquettes embedded into Euclidean lattice Zd have been considered in the literature [1, 2]. These models are based on area action and suffer the problem of non-scaling behaviour of the string tension and the dominance of branched polymers [3]. Several studies have analyzed the physical effects produced by rigidity of the surface introduced by adding dimension-less extrinsic curvature term to the area action [4]. Comprehensive review of the work in this area up to 1997 can be found in [3].

In this article we shall consider a model of random surfaces solely based on the with extrinsic curvature term embedded into 3d Euclidean lattice Z3 [5]. The corresponding transfer matrix has the form [6]

K(Q1, Q2) = exp{−β [k(Q1) + 2l(Q1△Q2) + k(Q2)] }, (1)

where Q1 and Q2 are closed polygon-loops on a two-dimensional lattice, k(Q) is the curvature and l(Q) is the length of the polygon-loop Q 1. This transfer matrix describes the propagation of the initial loop Q1 to the final loop Q2.

The spectrum of the transfer matrix which depends only on symmetric difference of initial and final loops Q1 △Q2

K˜(Q1, Q2) = exp{−2βl(Q1 △Q2) }, (2)

has been evaluated analytically in terms of correlation functions of the 2d Ising model in [8].

Our aim in this article is to consider the full transfer matrix (1) which depends not only on symmetric difference of loops Q1 △ Q2, but also on individual curvature terms k(Q1) + k(Q2). Because analytical solution of this problem is not known yet, we shall consider numerical evaluation of the transfer matrix on finite lattices in order to gain new insight into spectral properties of the transfer matrix in general case and critical properties of the model.

In the next sections we shall numerically compute the spectrum of transfer matrix on finite Euclidean lattices T 2 of the sizes N ⋆M , where N,M = 1, 2, 3, 4. Number of loops on these lattices γ = 2NM will be 2, 4, 16, 64, 512, 4096, 65536 and transfer matrices would be of the size γ × γ = 2NM × 2NM .

2 Basic Formulae

The partition function is defined as [6]

Z(β) = ∑

{Q1,Q2,...,Qτ}

Kβ(Q1, Q2) · · ·Kβ(Qτ , Q1) = trK τ β , (3)

where Kβ(Q1, Q2) is the transfer matrix of size γ × γ, defined above (1). Q1 and Q2 are closed polygon-loops on a two-dimensional lattice T 2 of size N ×M and γ = 2NM is the total number of polygon-loops on a lattice T 2. The transfer matrix (1) can be viewed

1We shall use the word ”loop” for the ”polygon-loop”.

2

as describing the propagation amplitude of the polygon-loop Q1 at time τ to another polygon-loop Q2 at the time τ + 1.

The functional k(Q) is the total curvature of the polygon-loop Q which is equal to the number of corners of the polygon ( the vertices with self-intersection are not counted!) and l(Q) is the length of Q which is equal to the number of its links. The length functional l(Q1 △Q2) is defined as [6]

l(Qi) + l(Qi+1)− 2 · l(Qi ∩Qi+1) = l(Qi △Qi+1), (4)

where the polygon-loop Q1 △ Q2 ≡ Q1 ∪ Q2 \ Q1 ∩ Q2 is a union of links Q1 ∪ Q2 without common links Q1 ∩ Q2. The operation △ maps two polygon-loops Q1 and Q2 into a polygon-loop Q = Q1 △Q2.

The eigenvalues of the transfer matrix Kβ(Q1, Q2) define all statistical properties of the system and can be found as a solution of the following integral equation in the loop space Π [8] ∑

{Q2}

Kβ(Q1, Q2) Ψ(Q2) = Λ(β) Ψ(Q1), (5)

where Ψ(Q) is a function on the loop space. The eigenvalues define the partition function (3)

Z3dNM(β) = Λ τ 0(β,NM) + ...+ Λ

τ γ−1(β,NM), (6)

and in the thermodynamical limit the free energy is equal to

− β f3d(β) = lim τNM→∞

1

τNM ln Z3d(β) = lim

N,M→∞

1

NM ln Λ0(β,MN), (7)

where Λ0(β,NM) is the largest eigenvalue defined at temperature β for the lattice of the size N ×M . The correlation length is defined by the ratio of eigenvalues Λ1(β,NM) and Λ0(β,NM) as

ξNM(β) = −1/ln[Λ1(β,NM)/Λ0(β,NM)],

and grows if the eigenvalue Λ1 approaches to the eigenvalue Λ0 at some critical tempera- ture βc. By the Frobenius-Perron theorem Λ0(β) is simple and Λ0(β) > Λ1(β) ≥ Λ2(β) ≥ ...

Because we can analyze numerically only finite-dimensional matrices we shall use renormalisation group approach to critical phenomena [9] in order to extract critical in- dexes. The renormalisation group equation

β¯ = β¯(β) (8)

is defined through the correlation function

ξNM(β)

N +M =

ξN¯M¯(β¯)

N¯ + M¯ , N +M > N¯ + M¯. (9)

The critical point βc is a fixed point β ∗ of the renormalisation group equation (8)

β∗ = β¯(β∗). (10)

The renormalisation group equation (8) can be expanded near fixed point β∗ into the series of the following form

β¯(β)− β∗ = g1(β − β ∗) + g2(β − β

∗)2 + ...

3

20 40 60 80 100

1

2

3

4

Figure 1: The typical behaviour of the eigenvalues. These are eigenvalues of the transfer matrix K as a function of ω = e−4β for the size 4× 4.

and the critical index ν is equal to

ν = ln[(N +M)/(N¯ + M¯)]

lng1 . (11)

We shall use these formulas to define fixed point and the critical indexes. In the next sections we shall explicitly construct transfer matrices on finite Euclidean

lattices T 2 of the sizes N ⋆M , where N,M = 1, 2, 3, 4. Number of loops on these lattices γ = 2NM will be 2, 4, 16, 64, 512, 4096, 65536 and transfer matrices would be of the size γ × γ = 2NM × 2NM .

3 General properties of the spectrum

We shall fix the boundary spins on a lattice T 2, let us say, in up direction. The convention is that we shall not count the boundary spins, thus if the lattice has the size 2× 3 then it has six internal spins 2. In this setting when all spins are in up direction then in the dual picture we shall have just an empty loop. When one spin inside lattice is turned down, then in the dual picture we shall have simply one loop consisting of four bonds, a box-loop, and so on.

Let us begin from the simplest lattice of the size 1× 1 and then move to more com- plicated ones in order to understand general properties of the spectrum. For this lattice which has just one spin we have empty loop ∅ and one box-loop 2. The transfer matri- ces K and K˜ describing transitions between these two loops can be computed by using definition (1) and (2)

K =

( 1 e−12β

e−12β e−8β

) , K˜ =

( 1 e−8β

e−8β 1

) , (12)

and then it is easy to find its eigenvalues:

2Λ0,1 = 1 + e −8β ±

√ 1 + 4e−24β + e−16β − 2e−8β, Λ˜0,1 = 1± e

−8β . (13)

For the lattice of the size 1 × 2 we have four loops: the empty loop, two one-box loops

2And ten boundary spins which we shall not count.

4

Figure 2: The spectrum of the transfer matrices K of the size 16 × 16 as a function of ω = e−2β . The qualitative character of the spectrum is the same as in Figure 1, but now we can see that at high temperatures some of the eigenvalues become negative.

and one two-box loop. The transfer matrices are therefore four by four:

K =

1 e−12β e−12β e−16β

e−12β e−8β e−20β e−16β

e−12β e−20β e−8β e−16β

e−16β e−16β e−16β e−8β

, K˜ =

1 e−8β e−8β e−12β

e−8β 1 e−12β e−8β

e−8β e−12β 1 e−8β

e−12β e−8β e−8β 1

, (14)

The eigenvalues of the matrix K˜ have been found in [8]

Λ˜0 = 1 + 2e −8β + e−12β , Λ˜1,2 = 1− e

−12β , Λ˜3 = 1− 2e −8β + e−12β , (15)

and for the matrix K we have found eigenvalues on a computer (see Figure 1). In both cases all eigenvalues are positive. As the last example of explicitly depicted transfer matrix

5

we shall present it for the lattice 2× 2. It has sixteen loops and is of the size 16× 16

K =

0 6 6 6 6 8 8 8 8 11 11 11 11 11 11 10 6 4 10 12 10 8 12 8 12 9 13 11 13 11 13 12 6 10 4 10 12 8 12 12 8 13 9 13 11 13 11 12 6 12 10 4 10 12 8 12 8 9 13 11 13 11 13 12 6 10 12 10 4 12 8 8 12 13 9 13 11 13 11 12 8 8 8 12 12 4 12 12 12 11 11 9 13 13 9 10 8 12 12 8 8 12 4 12 12 11 11 13 9 9 13 10 8 8 12 12 8 12 12 4 12 9 11 13 13 9 9 10 8 12 8 8 12 12 12 12 4 11 11 9 9 13 13 10 11 9 13 9 13 11 11 9 11 6 14 10 14 10 14 9 11 13 9 13 9 11 11 11 11 14 6 14 10 14 10 9 11 11 13 11 13 9 13 13 9 10 14 6 12 14 12 9 11 13 11 13 11 13 9 13 9 14 10 12 6 12 14 9 11 11 13 11 13 13 9 9 13 10 14 14 12 6 12 9 11 13 11 13 11 9 13 9 13 14 10 12 14 12 6 9 10 12 12 12 12 10 10 10 10 9 9 9 9 9 9 4

, (16)

where we show only the exponents. The eigenvalues of the matrix K˜ have been found in [8]

Λ˜0 = 1 + 4e −8β + 6e−12β + 5e−16β , Λ˜1,2,3,4 = 1 + 2e

−8β − 3e−16β , Λ˜5,6,7,8 = 1− 2e

−12β + 2e−16β , Λ˜9,10 = 1− 2e −12β + e−16β ,

Λ˜11,12,13,14 = 1− 2e −8β + e−16β , Λ˜15 = 1− 4e

−8β + 6e−12β − 3e−16β .

(17)

The eigenvalues of matrix K are shown on Figure 2. The general properties of the spec- trum at low and high temperatures can be easily understood. At β → 0 the eigenvalues are :

Λ0 = 2 NM , Λ1 = ..... = Λ2NM−1 = 0

and when β →∞ they are :

Λ0 = 1, Λ1 = ..... = Λ2NM−1 = 0.

These properties of the eigenvalues at the ends of the spectrum can be seen on Figures 1 and 2. These Figures also allow to see the behaviour of the spectrum at the intermediary temperatures, the largest eigenvalue Λ0 and the second one Λ1 approach to each other. The rest of the eigenvalues follow the behaviour of the Λ1. As we will see minimal distance between eigenvalue Λ0 and Λ1 decreases when the size of the lattice increases signaling the possibility of a phase transition.

4 Numerical calculation of eigenvalues Λ0 and Λ1

To generate the transfer matrices of the larger size we have developed a special fast algorithm to create polygon-loops Q and then to compute the elements of the transfer matrixK. This part of the program is working fast enough to generate matrices of the size 65536 and higher. To compute eigenvalues Λ0 and Λ1 we have used standard technique, that is acting by the matrix K on an arbitrary vector Ψ(Q) many, m, times. In the limit

6

Figure 3: The ratio Λ1/Λ0 as a function of β for the matrices of sizes 16 × 16, 64 × 64, 256× 256, 512× 512, 4096× 4096.

m → ∞ it converges to Ψ0(Q). In the same manner, acting on the orthogonal vector, one can get Ψ1(Q) and therefore corresponding eigenvalues. The real problem to handle higher size matrices is mostly connected with the memory of the computer.

The ratio of the eigenvalues Λ0/Λ1 is shown on Figure 3. This ratio increases with the size of the matrix K and reaches the value 0.6 for the matrix of the size 4096× 4096. To find fixed point β∗ of the renormalization group equation we have constructed the curves ξNM(β)/(N +M) for different lattices N ∗M . As one can see on Figure 4 there are two fixed points at β∗1 ≃ 0.17 and at β

∗ 2 ≃ 0.21. At the stable fixed point β

∗ 2 ≃ 0.21

we have computed the value g1 ≃ 1.99 and the critical index ν ≃ 0.59. This should be compared with the results obtained in [13] by low temperature expansion and by Monte- Carlo simulation of the corresponding system in [14, 15]. Our value is of the same order of magnitude. To confirm these results one should go to high-dimensional matrices.

It is also interesting to see the behaviour of eigenfunctions. The eigenfunctions Ψ(Q) for the matrix K˜(β) are already known and are independent of the temperature [8]. This is because transfer matrices K˜(β) commute with each other at different temperatures and the system is integrable [8]. This is a general property of integrable systems [7, 10, 11, 12]. The situation with the full transfer matrices K(β) is not known, therefore we should check if the eigenfunctions are temperature dependent. For that reason we shall compute the scalar product

( ΨΛ(Q)(β ′) ⋆ΨΛ(Q)(β) )

between eigenfunction ΨΛ(Q)(β) at different temperatures β ′ and β in order to see if

the scalar product depends on temperature. For the first eigenvalue Λ0 and β ′ = 0 the

result of the calculations is shown on Figure 5. It is clearly seen that the scalar product is temperature dependent and that it drastically changes near the critical points. This result drives us to the conclusion that full transfer matrix represents a non-integrable system.

7

Figure 4: The ratio ξNM(β)/(N +M) as a function of β for the matrices of sizes 16× 16 and 512× 512. There are two fixed points at β∗1 ≃ 0.17 and at β

∗ 2 ≃ 0.21.

Figure 5: Temperature dependents of the scalar product {Ψ0(Q)(0) ⋆Ψ0(Q)(β)} for the matrices of the sizes 16× 16, 64× 16, 512× 512, 4096× 4096.

8

We would like to thank Thordur Jonsson for stimulating discussions and A.Peppas for his collaboration in early stages of this work and calculation of eigenvalues for low dimensional matrices. This work was supported in part by the EEC Grant no. HPRN- CT-1999-00161 and the Grant from Ministry of Research and Technology: ”Bilateral Cooperation between Hellenic Republic and Georgia”.

References

[1] D.Weingarten. Nucl.Phys.B210 (1982) 229 A.Maritan and C.Omero. Phys.Lett. B109 (1982) 51 T.Sterling and J.Greensite. Phys.Lett. B121 (1983) 345 B.Durhuus,J.Fro¨hlich and T.Jonsson. Nucl.Phys.B225 (1983) 183 J.Ambjørn,B.Durhuus,J.Fro¨hlich and T.Jonsson. Nucl.Phys.B290 (1987) 480 T.Hofsa¨ss and H.Kleinert. Phys.Lett. A102 (1984) 420 M.Karowski and H.J.Thun. Phys.Rev.Lett. 54 (1985) 2556 F.David. Europhys.Lett. 9 (1989) 575

[2] D.Gross. Phys.Lett. B138 (1984) 185 V.A.Kazakov. Phys.Lett. B150 (1985) 282 F.David. Nucl.Phys. B257 (1985) 45 J.Ambjørn, B.Durhuus and J.Fro¨hlich. Nucl.Phys. B257 (1985) 433 V.A.Kazakov, I.K.Kostov and A.A.Migdal. Phys.Lett. B157 (1985) 295

[3] J.Ambjorn, B.Durhuus and T.Jonsson. Quantum geometry. Cambridge Monographs on Mathematical Physics. Cambridge 1998; J.Phys.A 21 (1988) 981

[4] W.Helfrich. Z.Naturforsch. C28 (1973) 693; J.Phys.(Paris) 46 (1985) 1263 L.Peliti and S.Leibler. Phys.Rev.Lett. 54 (1985) 1690 A.Polykov. Nucl.Phys.B268 (1986) 406 D.Forster. Phys.Lett. 114A (1986) 115 H.Kleinert. Phys.Lett. 174B (1986) 335 T.L.Curtright and et.al. Phys.Rev.Lett. 57 (1986)799; Phys.Rev. D34 (1986) 3811 F.David. Europhys.Lett. 2 (1986) 577 P.O.Mazur and V.P.Nair. Nucl.Phys. B284 (1987) 146 E.Braaten and C.K.Zachos. Pys.Rev. D35 (1987) 1512 E.Braaten, R.D.Pisarski and S.M.Tye. Phys.Rev.Lett. 58 (1987) 93 P.Olesen and S.K.Yang. Nucl.Phys. B283 (1987) 73 R.D.Pisarski. Phys.Rev.Lett. 58 (1987) 1300

[5] R.V. Ambartzumian, G.K. Savvidy , K.G. Savvidy and G.S. Sukiasian. Phys. Lett. B275 (1992) 99 G.K. Savvidy and K.G. Savvidy. Mod.Phys.Lett. A8 (1993) 2963 G.K. Savvidy and K.G. Savvidy. Int. J. Mod. Phys. A8 (1993) 3993 G.K.Savvidy and F.J.Wegner. Nucl.Phys.B413(1994)605 G.K. Savvidy and K.G. Savvidy. Phys.Lett. B324 (1994) 72

9

[6] G.K. Savvidy and K.G. Savvidy. Phys.Lett. B337 (1994) 333; Mod.Phys.Lett. A11 (1996) 1379. G.K. Savvidy, K.G. Savvidy and P.K.Savvidy Phys.Lett. A221 (1996) 233

[7] H.A.Kramers and G.H.Wannier. Phys.Rev. 60 (1941) 252 L.Onsager. Phys.Rev. 65 (1944) 117 M.Kac and J.C.Ward. Phys.Rev. 88 (1952) 1332 C.A.Hurst and H.S.Green. J.Chem.Phys. 33 (1960) 1059

[8] T.Jonsson and G.K.Savvidy. Phys.Lett.B449 (1999) 254 T.Jonsson and G.K.Savvidy. Nucl.Phys. B575 (2000) 661 G.K.Savvidy. JHEP 0009 (2000) 044

[9] M.P.Nightingale. Physica 83A (1976) 561

[10] R.J.Baxter, Exactly solved models in statistical mechanics. Academic Press, London 1982

[11] A.B. Zamolodchikov. Commun.Math.Phys.79 (1981) 489 V.V.Bazhanov and R.J.Baxter. J.Stat.Phys. 69(1992) 453

[12] A.Polyakov. Gauge fields and String. (Harwood Academic Publishers, 1987) E.Fradkin, M.Srednicky and L.Susskind. Phys.Rev. D21 (1980) 2885 C.Itzykson. Nucl.Phys. B210 (1982) 477 A.Casher, D.Foerster and P.Windey. Nucl.Phys. B251 (1985) 29

[13] R. Pietig and F.J. Wegner. Nucl.Phys. B466 (1996) 513; Nucl.Phys.B525 (1998) 549

[14] G.K.Bathas et.al. Mod.Phys.Lett. A10 (1995) 2695 D.Johnson and R.K.P.C.Malmini, Phys.Lett. B378 (1996) 87 M.Baig, D.Espriu, D.Johnson and R.K.P.C.Malmini, J.Phys. A30 (1997) 407; J.Phys. A30 (1997) 7695 G.Koutsoumbas et.al. Phys.Lett. B410 (1997) 241

[15] A.Cappi, P.Colangelo, G.Gonella and A.Maritan, Nucl. Phys. B370 (1992) 659 G.Gonnella, S.Lise and A.Maritan, Europhys. Lett. 32 (1995) 735 E.N.M.Cirillo and G.Gonella. J.Phys.A: Math.Gen.28 (1995) 867

10

iv :c

on d-

m at

/0 10

22 45

v1 [

co nd

-m at.

sta t-m

ec h]

1 4 F

eb 20

01

NRCPS-HE-2001-4

The Spectrum of the Loop Transfer Matrix on Finite Lattice

Georgios Daskalakis and George K.Savvidy

National Research Center Demokritos, Ag. Paraskevi, GR-15310 Athens, Hellenic Republic

Abstract

We consider the model of random surfaces with extrinsic curvature term embedded into 3d Euclidean lattice Z3. On a 3d Euclidean lattice it has equivalent representation in terms of transfer matrix K(Qi, Qf), which describes the propagation of loops Q. We study the spectrum of the transfer matrix K(Qi, Qf) on finite dimensional lattices. The renormalisation group technique is used to investigate phase structure of the model and its critical behaviour.

1 Introduction

Various models of random surfaces built out of triangles embedded into continuous space Rd and surfaces built out of plaquettes embedded into Euclidean lattice Zd have been considered in the literature [1, 2]. These models are based on area action and suffer the problem of non-scaling behaviour of the string tension and the dominance of branched polymers [3]. Several studies have analyzed the physical effects produced by rigidity of the surface introduced by adding dimension-less extrinsic curvature term to the area action [4]. Comprehensive review of the work in this area up to 1997 can be found in [3].

In this article we shall consider a model of random surfaces solely based on the with extrinsic curvature term embedded into 3d Euclidean lattice Z3 [5]. The corresponding transfer matrix has the form [6]

K(Q1, Q2) = exp{−β [k(Q1) + 2l(Q1△Q2) + k(Q2)] }, (1)

where Q1 and Q2 are closed polygon-loops on a two-dimensional lattice, k(Q) is the curvature and l(Q) is the length of the polygon-loop Q 1. This transfer matrix describes the propagation of the initial loop Q1 to the final loop Q2.

The spectrum of the transfer matrix which depends only on symmetric difference of initial and final loops Q1 △Q2

K˜(Q1, Q2) = exp{−2βl(Q1 △Q2) }, (2)

has been evaluated analytically in terms of correlation functions of the 2d Ising model in [8].

Our aim in this article is to consider the full transfer matrix (1) which depends not only on symmetric difference of loops Q1 △ Q2, but also on individual curvature terms k(Q1) + k(Q2). Because analytical solution of this problem is not known yet, we shall consider numerical evaluation of the transfer matrix on finite lattices in order to gain new insight into spectral properties of the transfer matrix in general case and critical properties of the model.

In the next sections we shall numerically compute the spectrum of transfer matrix on finite Euclidean lattices T 2 of the sizes N ⋆M , where N,M = 1, 2, 3, 4. Number of loops on these lattices γ = 2NM will be 2, 4, 16, 64, 512, 4096, 65536 and transfer matrices would be of the size γ × γ = 2NM × 2NM .

2 Basic Formulae

The partition function is defined as [6]

Z(β) = ∑

{Q1,Q2,...,Qτ}

Kβ(Q1, Q2) · · ·Kβ(Qτ , Q1) = trK τ β , (3)

where Kβ(Q1, Q2) is the transfer matrix of size γ × γ, defined above (1). Q1 and Q2 are closed polygon-loops on a two-dimensional lattice T 2 of size N ×M and γ = 2NM is the total number of polygon-loops on a lattice T 2. The transfer matrix (1) can be viewed

1We shall use the word ”loop” for the ”polygon-loop”.

2

as describing the propagation amplitude of the polygon-loop Q1 at time τ to another polygon-loop Q2 at the time τ + 1.

The functional k(Q) is the total curvature of the polygon-loop Q which is equal to the number of corners of the polygon ( the vertices with self-intersection are not counted!) and l(Q) is the length of Q which is equal to the number of its links. The length functional l(Q1 △Q2) is defined as [6]

l(Qi) + l(Qi+1)− 2 · l(Qi ∩Qi+1) = l(Qi △Qi+1), (4)

where the polygon-loop Q1 △ Q2 ≡ Q1 ∪ Q2 \ Q1 ∩ Q2 is a union of links Q1 ∪ Q2 without common links Q1 ∩ Q2. The operation △ maps two polygon-loops Q1 and Q2 into a polygon-loop Q = Q1 △Q2.

The eigenvalues of the transfer matrix Kβ(Q1, Q2) define all statistical properties of the system and can be found as a solution of the following integral equation in the loop space Π [8] ∑

{Q2}

Kβ(Q1, Q2) Ψ(Q2) = Λ(β) Ψ(Q1), (5)

where Ψ(Q) is a function on the loop space. The eigenvalues define the partition function (3)

Z3dNM(β) = Λ τ 0(β,NM) + ...+ Λ

τ γ−1(β,NM), (6)

and in the thermodynamical limit the free energy is equal to

− β f3d(β) = lim τNM→∞

1

τNM ln Z3d(β) = lim

N,M→∞

1

NM ln Λ0(β,MN), (7)

where Λ0(β,NM) is the largest eigenvalue defined at temperature β for the lattice of the size N ×M . The correlation length is defined by the ratio of eigenvalues Λ1(β,NM) and Λ0(β,NM) as

ξNM(β) = −1/ln[Λ1(β,NM)/Λ0(β,NM)],

and grows if the eigenvalue Λ1 approaches to the eigenvalue Λ0 at some critical tempera- ture βc. By the Frobenius-Perron theorem Λ0(β) is simple and Λ0(β) > Λ1(β) ≥ Λ2(β) ≥ ...

Because we can analyze numerically only finite-dimensional matrices we shall use renormalisation group approach to critical phenomena [9] in order to extract critical in- dexes. The renormalisation group equation

β¯ = β¯(β) (8)

is defined through the correlation function

ξNM(β)

N +M =

ξN¯M¯(β¯)

N¯ + M¯ , N +M > N¯ + M¯. (9)

The critical point βc is a fixed point β ∗ of the renormalisation group equation (8)

β∗ = β¯(β∗). (10)

The renormalisation group equation (8) can be expanded near fixed point β∗ into the series of the following form

β¯(β)− β∗ = g1(β − β ∗) + g2(β − β

∗)2 + ...

3

20 40 60 80 100

1

2

3

4

Figure 1: The typical behaviour of the eigenvalues. These are eigenvalues of the transfer matrix K as a function of ω = e−4β for the size 4× 4.

and the critical index ν is equal to

ν = ln[(N +M)/(N¯ + M¯)]

lng1 . (11)

We shall use these formulas to define fixed point and the critical indexes. In the next sections we shall explicitly construct transfer matrices on finite Euclidean

lattices T 2 of the sizes N ⋆M , where N,M = 1, 2, 3, 4. Number of loops on these lattices γ = 2NM will be 2, 4, 16, 64, 512, 4096, 65536 and transfer matrices would be of the size γ × γ = 2NM × 2NM .

3 General properties of the spectrum

We shall fix the boundary spins on a lattice T 2, let us say, in up direction. The convention is that we shall not count the boundary spins, thus if the lattice has the size 2× 3 then it has six internal spins 2. In this setting when all spins are in up direction then in the dual picture we shall have just an empty loop. When one spin inside lattice is turned down, then in the dual picture we shall have simply one loop consisting of four bonds, a box-loop, and so on.

Let us begin from the simplest lattice of the size 1× 1 and then move to more com- plicated ones in order to understand general properties of the spectrum. For this lattice which has just one spin we have empty loop ∅ and one box-loop 2. The transfer matri- ces K and K˜ describing transitions between these two loops can be computed by using definition (1) and (2)

K =

( 1 e−12β

e−12β e−8β

) , K˜ =

( 1 e−8β

e−8β 1

) , (12)

and then it is easy to find its eigenvalues:

2Λ0,1 = 1 + e −8β ±

√ 1 + 4e−24β + e−16β − 2e−8β, Λ˜0,1 = 1± e

−8β . (13)

For the lattice of the size 1 × 2 we have four loops: the empty loop, two one-box loops

2And ten boundary spins which we shall not count.

4

Figure 2: The spectrum of the transfer matrices K of the size 16 × 16 as a function of ω = e−2β . The qualitative character of the spectrum is the same as in Figure 1, but now we can see that at high temperatures some of the eigenvalues become negative.

and one two-box loop. The transfer matrices are therefore four by four:

K =

1 e−12β e−12β e−16β

e−12β e−8β e−20β e−16β

e−12β e−20β e−8β e−16β

e−16β e−16β e−16β e−8β

, K˜ =

1 e−8β e−8β e−12β

e−8β 1 e−12β e−8β

e−8β e−12β 1 e−8β

e−12β e−8β e−8β 1

, (14)

The eigenvalues of the matrix K˜ have been found in [8]

Λ˜0 = 1 + 2e −8β + e−12β , Λ˜1,2 = 1− e

−12β , Λ˜3 = 1− 2e −8β + e−12β , (15)

and for the matrix K we have found eigenvalues on a computer (see Figure 1). In both cases all eigenvalues are positive. As the last example of explicitly depicted transfer matrix

5

we shall present it for the lattice 2× 2. It has sixteen loops and is of the size 16× 16

K =

0 6 6 6 6 8 8 8 8 11 11 11 11 11 11 10 6 4 10 12 10 8 12 8 12 9 13 11 13 11 13 12 6 10 4 10 12 8 12 12 8 13 9 13 11 13 11 12 6 12 10 4 10 12 8 12 8 9 13 11 13 11 13 12 6 10 12 10 4 12 8 8 12 13 9 13 11 13 11 12 8 8 8 12 12 4 12 12 12 11 11 9 13 13 9 10 8 12 12 8 8 12 4 12 12 11 11 13 9 9 13 10 8 8 12 12 8 12 12 4 12 9 11 13 13 9 9 10 8 12 8 8 12 12 12 12 4 11 11 9 9 13 13 10 11 9 13 9 13 11 11 9 11 6 14 10 14 10 14 9 11 13 9 13 9 11 11 11 11 14 6 14 10 14 10 9 11 11 13 11 13 9 13 13 9 10 14 6 12 14 12 9 11 13 11 13 11 13 9 13 9 14 10 12 6 12 14 9 11 11 13 11 13 13 9 9 13 10 14 14 12 6 12 9 11 13 11 13 11 9 13 9 13 14 10 12 14 12 6 9 10 12 12 12 12 10 10 10 10 9 9 9 9 9 9 4

, (16)

where we show only the exponents. The eigenvalues of the matrix K˜ have been found in [8]

Λ˜0 = 1 + 4e −8β + 6e−12β + 5e−16β , Λ˜1,2,3,4 = 1 + 2e

−8β − 3e−16β , Λ˜5,6,7,8 = 1− 2e

−12β + 2e−16β , Λ˜9,10 = 1− 2e −12β + e−16β ,

Λ˜11,12,13,14 = 1− 2e −8β + e−16β , Λ˜15 = 1− 4e

−8β + 6e−12β − 3e−16β .

(17)

The eigenvalues of matrix K are shown on Figure 2. The general properties of the spec- trum at low and high temperatures can be easily understood. At β → 0 the eigenvalues are :

Λ0 = 2 NM , Λ1 = ..... = Λ2NM−1 = 0

and when β →∞ they are :

Λ0 = 1, Λ1 = ..... = Λ2NM−1 = 0.

These properties of the eigenvalues at the ends of the spectrum can be seen on Figures 1 and 2. These Figures also allow to see the behaviour of the spectrum at the intermediary temperatures, the largest eigenvalue Λ0 and the second one Λ1 approach to each other. The rest of the eigenvalues follow the behaviour of the Λ1. As we will see minimal distance between eigenvalue Λ0 and Λ1 decreases when the size of the lattice increases signaling the possibility of a phase transition.

4 Numerical calculation of eigenvalues Λ0 and Λ1

To generate the transfer matrices of the larger size we have developed a special fast algorithm to create polygon-loops Q and then to compute the elements of the transfer matrixK. This part of the program is working fast enough to generate matrices of the size 65536 and higher. To compute eigenvalues Λ0 and Λ1 we have used standard technique, that is acting by the matrix K on an arbitrary vector Ψ(Q) many, m, times. In the limit

6

Figure 3: The ratio Λ1/Λ0 as a function of β for the matrices of sizes 16 × 16, 64 × 64, 256× 256, 512× 512, 4096× 4096.

m → ∞ it converges to Ψ0(Q). In the same manner, acting on the orthogonal vector, one can get Ψ1(Q) and therefore corresponding eigenvalues. The real problem to handle higher size matrices is mostly connected with the memory of the computer.

The ratio of the eigenvalues Λ0/Λ1 is shown on Figure 3. This ratio increases with the size of the matrix K and reaches the value 0.6 for the matrix of the size 4096× 4096. To find fixed point β∗ of the renormalization group equation we have constructed the curves ξNM(β)/(N +M) for different lattices N ∗M . As one can see on Figure 4 there are two fixed points at β∗1 ≃ 0.17 and at β

∗ 2 ≃ 0.21. At the stable fixed point β

∗ 2 ≃ 0.21

we have computed the value g1 ≃ 1.99 and the critical index ν ≃ 0.59. This should be compared with the results obtained in [13] by low temperature expansion and by Monte- Carlo simulation of the corresponding system in [14, 15]. Our value is of the same order of magnitude. To confirm these results one should go to high-dimensional matrices.

It is also interesting to see the behaviour of eigenfunctions. The eigenfunctions Ψ(Q) for the matrix K˜(β) are already known and are independent of the temperature [8]. This is because transfer matrices K˜(β) commute with each other at different temperatures and the system is integrable [8]. This is a general property of integrable systems [7, 10, 11, 12]. The situation with the full transfer matrices K(β) is not known, therefore we should check if the eigenfunctions are temperature dependent. For that reason we shall compute the scalar product

( ΨΛ(Q)(β ′) ⋆ΨΛ(Q)(β) )

between eigenfunction ΨΛ(Q)(β) at different temperatures β ′ and β in order to see if

the scalar product depends on temperature. For the first eigenvalue Λ0 and β ′ = 0 the

result of the calculations is shown on Figure 5. It is clearly seen that the scalar product is temperature dependent and that it drastically changes near the critical points. This result drives us to the conclusion that full transfer matrix represents a non-integrable system.

7

Figure 4: The ratio ξNM(β)/(N +M) as a function of β for the matrices of sizes 16× 16 and 512× 512. There are two fixed points at β∗1 ≃ 0.17 and at β

∗ 2 ≃ 0.21.

Figure 5: Temperature dependents of the scalar product {Ψ0(Q)(0) ⋆Ψ0(Q)(β)} for the matrices of the sizes 16× 16, 64× 16, 512× 512, 4096× 4096.

8

We would like to thank Thordur Jonsson for stimulating discussions and A.Peppas for his collaboration in early stages of this work and calculation of eigenvalues for low dimensional matrices. This work was supported in part by the EEC Grant no. HPRN- CT-1999-00161 and the Grant from Ministry of Research and Technology: ”Bilateral Cooperation between Hellenic Republic and Georgia”.

References

[1] D.Weingarten. Nucl.Phys.B210 (1982) 229 A.Maritan and C.Omero. Phys.Lett. B109 (1982) 51 T.Sterling and J.Greensite. Phys.Lett. B121 (1983) 345 B.Durhuus,J.Fro¨hlich and T.Jonsson. Nucl.Phys.B225 (1983) 183 J.Ambjørn,B.Durhuus,J.Fro¨hlich and T.Jonsson. Nucl.Phys.B290 (1987) 480 T.Hofsa¨ss and H.Kleinert. Phys.Lett. A102 (1984) 420 M.Karowski and H.J.Thun. Phys.Rev.Lett. 54 (1985) 2556 F.David. Europhys.Lett. 9 (1989) 575

[2] D.Gross. Phys.Lett. B138 (1984) 185 V.A.Kazakov. Phys.Lett. B150 (1985) 282 F.David. Nucl.Phys. B257 (1985) 45 J.Ambjørn, B.Durhuus and J.Fro¨hlich. Nucl.Phys. B257 (1985) 433 V.A.Kazakov, I.K.Kostov and A.A.Migdal. Phys.Lett. B157 (1985) 295

[3] J.Ambjorn, B.Durhuus and T.Jonsson. Quantum geometry. Cambridge Monographs on Mathematical Physics. Cambridge 1998; J.Phys.A 21 (1988) 981

[4] W.Helfrich. Z.Naturforsch. C28 (1973) 693; J.Phys.(Paris) 46 (1985) 1263 L.Peliti and S.Leibler. Phys.Rev.Lett. 54 (1985) 1690 A.Polykov. Nucl.Phys.B268 (1986) 406 D.Forster. Phys.Lett. 114A (1986) 115 H.Kleinert. Phys.Lett. 174B (1986) 335 T.L.Curtright and et.al. Phys.Rev.Lett. 57 (1986)799; Phys.Rev. D34 (1986) 3811 F.David. Europhys.Lett. 2 (1986) 577 P.O.Mazur and V.P.Nair. Nucl.Phys. B284 (1987) 146 E.Braaten and C.K.Zachos. Pys.Rev. D35 (1987) 1512 E.Braaten, R.D.Pisarski and S.M.Tye. Phys.Rev.Lett. 58 (1987) 93 P.Olesen and S.K.Yang. Nucl.Phys. B283 (1987) 73 R.D.Pisarski. Phys.Rev.Lett. 58 (1987) 1300

[5] R.V. Ambartzumian, G.K. Savvidy , K.G. Savvidy and G.S. Sukiasian. Phys. Lett. B275 (1992) 99 G.K. Savvidy and K.G. Savvidy. Mod.Phys.Lett. A8 (1993) 2963 G.K. Savvidy and K.G. Savvidy. Int. J. Mod. Phys. A8 (1993) 3993 G.K.Savvidy and F.J.Wegner. Nucl.Phys.B413(1994)605 G.K. Savvidy and K.G. Savvidy. Phys.Lett. B324 (1994) 72

9

[6] G.K. Savvidy and K.G. Savvidy. Phys.Lett. B337 (1994) 333; Mod.Phys.Lett. A11 (1996) 1379. G.K. Savvidy, K.G. Savvidy and P.K.Savvidy Phys.Lett. A221 (1996) 233

[7] H.A.Kramers and G.H.Wannier. Phys.Rev. 60 (1941) 252 L.Onsager. Phys.Rev. 65 (1944) 117 M.Kac and J.C.Ward. Phys.Rev. 88 (1952) 1332 C.A.Hurst and H.S.Green. J.Chem.Phys. 33 (1960) 1059

[8] T.Jonsson and G.K.Savvidy. Phys.Lett.B449 (1999) 254 T.Jonsson and G.K.Savvidy. Nucl.Phys. B575 (2000) 661 G.K.Savvidy. JHEP 0009 (2000) 044

[9] M.P.Nightingale. Physica 83A (1976) 561

[10] R.J.Baxter, Exactly solved models in statistical mechanics. Academic Press, London 1982

[11] A.B. Zamolodchikov. Commun.Math.Phys.79 (1981) 489 V.V.Bazhanov and R.J.Baxter. J.Stat.Phys. 69(1992) 453

[12] A.Polyakov. Gauge fields and String. (Harwood Academic Publishers, 1987) E.Fradkin, M.Srednicky and L.Susskind. Phys.Rev. D21 (1980) 2885 C.Itzykson. Nucl.Phys. B210 (1982) 477 A.Casher, D.Foerster and P.Windey. Nucl.Phys. B251 (1985) 29

[13] R. Pietig and F.J. Wegner. Nucl.Phys. B466 (1996) 513; Nucl.Phys.B525 (1998) 549

[14] G.K.Bathas et.al. Mod.Phys.Lett. A10 (1995) 2695 D.Johnson and R.K.P.C.Malmini, Phys.Lett. B378 (1996) 87 M.Baig, D.Espriu, D.Johnson and R.K.P.C.Malmini, J.Phys. A30 (1997) 407; J.Phys. A30 (1997) 7695 G.Koutsoumbas et.al. Phys.Lett. B410 (1997) 241

[15] A.Cappi, P.Colangelo, G.Gonella and A.Maritan, Nucl. Phys. B370 (1992) 659 G.Gonnella, S.Lise and A.Maritan, Europhys. Lett. 32 (1995) 735 E.N.M.Cirillo and G.Gonella. J.Phys.A: Math.Gen.28 (1995) 867

10

Comments