arxiv:hep-th/ v2 19 May 1992

Transcription

1 RIMS Correlation Functions of the XXZ model for < 1 MichioJimbo a, Kei Miki b, Tetsuji Miwa c and Atsushi Nakayashiki d arxiv:hep-th/ v2 19 May 1992 a Department of Mathematics,Faculty of Science, Kyoto University, Kyoto 606, Japan b Department of Mathematical Science, Faculty of Engineering Science, Osaka University, Toyonaka, Osaka 560, Japan c Research Institute for Mathematical Sciences, Kyoto University, Kyoto 606, Japan d The Graduate School of Science and Technology, Kobe University, Rokkodai, Kobe 657, Japan Dedicated to Professor Chen Ning Yang on the occasion of his seventieth birthday Abstract. A new approach to the correlation functions is presented for the XXZ model in the anti-ferroelectric regime. The method is based on the recent realization of the quantum affine symmetry using vertex operators. With the aid of a boson representation for the latter, an integral formula is found for correlation functions of arbitrary local operators. As a special case it reproduces the spontaneous staggered polarization obtained earlier by Baxter.

2 2 M. Jimbo et al. 0. Introduction In this article we consider the one-dimensional infinite spin-chain H XXZ = 1 σ x 2 k+1 σk x + σy k+1 σy k + σz k+1 k σz k= 0.1 known classically as the XXZ model. We shall limit ourselves strictly to the antiferroelectric regime < 1. Our aim is to find an exact expression for the spin correlation functions σ γ1 i 1 σ γn i n, where σ γ i = σi x, σy i, σz i and = vac vac denotes the ground state average. Here the ground state vac means one of the two ground states in the anti-ferroelectric regime. In Baxter s paper [1] they are denoted by ±. The present paper is based entirely on the framework of the recent work [2], which discribes the realization of the quantum affine symmetry for 0.1 using the q-deformed vertex operators. Let us recall below the contents of [2] which are relevant to the subsequent discussions. The Hamiltonian 0.1 is an operator acting on the infinite tensor product V = V V V of the two dimensional space V = C 2. The space V admits also an action of the quantum affine algebra U = U q ŝl 2 via the iterated coproduct. A naïve computation shows that the algebra U provides an exact symmetry for 0.1. Namely if = q + q 1 /2, then [H XXZ, U ] = 0 where U = U qŝl 2 denotes the subalgebra without the grading operator d [2], while d plays the role of the boost operator. However the actions of 0.1 and U are both defined only formally, and the issue is how to extract the theory free from the difficulties of divergence. The basic idea in [2] is to replace the formal object V by the level 0 U-module F λ,µ = V λ V µ a HomV µ, V λ, 0.2 where V λ λ = Λ 0, Λ 1 denotes the level 1 highest weight U-module, and V λ a signifies its dual the superfix indicates that the U-module structure is given via the antipode a. The choice λ = µ resp. λ µ is responsible for the even resp. odd particle sector. V Λ 0 resp. V Λ 1 means that we are working in the boundary condition σ z 2k = 1, σz 2k+1 = 1 resp. σz 2k = 1, σz 2k+1 = 1 for k >> 0, and V Λ 0 a resp. V Λ 1 a means the same thing for k << 0. The tensor product should be completed in the q-adic sense to allow for infinite sums see [2] for a precise treatment. To make contact with the naïve picture of V, one utilizes the embedding of V λ into the half infinite tensor product V V V [3]. This is supplied by iterating the vertex operators Φ µv λ : V λ V µ V. 0.3 It is conjectured [2] that there is a unique normalization of 0.3 which makes the infinite iteration convergent. Similarly V µ a embeds to the other half infinite tensor product V V V, giving altogether the embedding F λ,µ V V V. 0.4 Eq. 0.4 provides the principle of interpreting the notions defined in the picture V by pulling them back to F λ,µ. The translation operator T is the first such example. The Hamiltonian 0.1 itself is defined in terms of T and the grading operator d; see [2]. In this paper we follow the same principle to formulate the local operators σ γ k as acting on F λ,µ.

3 Correlation Functions of the XXZ model for < 1 3 In the language of 0.2 the ground state vacuum vector is given by the identity element of F λ,λ = HomV λ, V λ. The inner product of two vectors f, g F λ,λ is given as the trace f g = tr V λq 2ρ fg tr V λ q 2ρ 0.5 where ρ = Λ 0 +Λ 1 we have changed the normalization from [2], see Sect. 1 below. Thus the correlation functions vac σ γ1 i 1 σ γn i n vac can be expressed as the trace of products of the vertex operators 0.3. To perform the evaluation of these traces, we invoke the bosonization method. The realization of the q-deformed currents on level one modules was done in [4] using ordinary bosons. In the same spirit we derive the formulas for the vertex operators 0.3 in terms of bosons. This leads to an explicit formula for the correlators in terms of certain integrals of meromorphic functions. We verify that in the simplest case of the one-point function σk z this formula reproduces the known result for the spontaneous staggered polarization due to Baxter [1]. Jacques Perk noted a strong similarity between our formula and the formula for the Ising model spin-spin correlation functions given in [5,6]. The bosonization of the vertex operators also enables us to compute the n-point functions discussed in Sect. 6.8 of [2]. They are the matrix elements not the trace of the product of the vertex operators with respect to the highest weight vectors. The conjectural formula 6.39 of [2] is thus proved. In this paper, we do not go into details on this matter. The plan of the paper is as follows. In Sect. 1 we formulate the local operators and their correlators using the scheme mentioned above. In our algebraic formulation the algebra U, intertwiners, etc. are a priori defined over the base field F = Qq q an indeterminate, see the remark at the end of Sect. 2; we shall see however that the resulting formulas are meaningful for complex values of q with q < 1. In Sect. 2 we outline the bosonization of the vertex operators. The basic ingredients are the Drinfeld realization of the algebra U and the result of [4], which we recall briefly. Unlike the classical case q = 1 only one of the components Φ z of 0.3 has the exponential form, while the other one does not but is given as the q-commutator of Φ z with a generator f 1 of U. Sect. 3 is devoted to the formula for the correlators.

4 4 M. Jimbo et al. 1. Vacuum expectation values of local operators Let us recall some notations of [2]. We denote by V = Fv + Fv the 2-dimensional vector space on which the Pauli matrices σ x, σ y, σ z act. We consider V as the 2-dimensional U qŝl2 -module. The XXZ-Hamiltonian formally acts on the infinite tensor product V of V. We label the components of the tensor product by integers k Z from right k to left + k. In our mathematical scheme, we replace the semi-infinite ŝl2 tensor product of the components k 1 by the irreducible highest weight U q - module V Λ i i = 0, 1 Let us define the action of the local operators on V Λ i. Let L EndV V. The operator L naturally acts on the tensor product of the com- }{{} n ponents n k 1 of V. We want to interprete this action as one in End V Λ i. For this purpose, we use the vertex operators [2]: they are the intertwiners Φ µv λ z : V λ V µ V z, ΦµV λ Φ µ λv z : V λ V z V µ, µv zv = Φ λ+ zv v µv + + Φ λ zv v 1.1 Φµ λv zv v ± = Φ µ λv ±zv, where λ = Λ 0, µ = Λ 1 or λ = Λ 1, µ = Λ 0, and V z = V F[z, z 1 ] signifies the U q ŝl2 - module associated to V with the spectral parameter z see [2], Sect. 6. Precisely speaking, we consider the vertex operators as the generating functions of their Fourier components in terms of the spectral parameter z. They are normalized as Φ Λ1V Λ 0 z u Λ0 = u Λ1 v +, ΦΛ 0V Λ 1 z u Λ1 = u Λ0 v + +, 1.2 Φ Λ1 Λ1V Φ Λ 0 z/q2, ΦΛ 1 Λ 0V z = q 1 ΦΛ 1V Λ 0+ z/q2, Λ 0V +z = Φ Λ0 Λ 1V +z = q Φ Λ0V Λ 1 z/q2, ΦΛ 0 q 2 ; q 4 q 4 ; q 4 Φλ µv z Φ µv λ z = id V λ. Λ 1V z = Φ Λ0V Λ 1+ z/q2, Here u Λi denotes the highest weight vector of V Λ i. We have used the standard notation z; p = j=0 1 zpj. Set Λ n = Λ n 2 for simplicity of notation. Given an L as above, we define the operator L = ρ i z n,...,z 1 L End V Λ i by L = q2 ; q 4 n q 4 ; q 4 n id V Λi+n L Φ Λi Λ z Λi+n 1 i+1v 1 Φ Λ z i+nv n Λi+nV Λi+1V Φ Λ i+n 1 z n Φ Λ i z 1. In this paper, we use the following convention different from that of [2] for the invariant bilinear form on P = ZΛ 0 ZΛ 1 Zδ: Λ 0, Λ 0 = 0, Λ 0, α 1 = 0, Λ 0, δ = 1, α 1, α 1 = 2, α 1, δ = 0, δ, δ = 0. Note that Λ 1 = Λ 0 + α 1 /2, δ = α 0 + α 1. We set ρ = Λ 0 + Λ 1 as usual, and also α = α 1 for simplicity. We identify P = Zh 0 Zh 1 Zd as a subset of P via,. We have ρ = 2d + α 1 /2.

5 Correlation Functions of the XXZ model for < 1 5 The vacuum expectation value L i z n,...,z 1 of the local operator L is given by L i z n,...,z 1 = tr V Λ i q 2ρ ρ z i n,...,z 1 L. 1.3 q 2ρ tr V Λi This is a consequence of the formula for the invariant inner product in the space V Λ i V Λ i a [2]. For the XXZ model correlator we specialize the spectral parameters to z 1 = = z n. We expect that the formula unspecialized would give the equal-row vertical arrow correlator for the inhomogeneous six-vertex model with z k being the trigonometric spectral parameter of the k-th vertical line. Here we use the usual language of vertex models on the 2-dimensional square lattice, in which the fluctuation variables are described as arrows sitting on vertical or horizontal edges. Note the following selection rule, which is specific to the six vertex model: L z i n,...,z 1 = 0 for L = E ε n ε n E ε 1 ε 1 E ij is a matrix unit such that ε ε n ε ε n. 2. Bosonization Let us first recall Drinfeld s realization of the quantum affine algebra U = U qŝl2 [7]. It is an associative algebra generated by the letters {x ± k k Z}, {a l l Z 0 }, γ ±1/2 and K, satisfying the following defining relations. γ ±1/2 the center of the algebra, [a k, a l ] = δ k+l,0 1 k [2k]γk γ k q q 1, [a k, K] = 0, Kx ± k K 1 = q ±2 x ± k, [a k, x ± l ] = ±1 k [2k]γ k /2 x ± k+l, x ± k+1 x± l q ±2 x ± l x± k+1 = q±2 x ± k x± l+1 x± l+1 x± k, [x + k, x l ] = 1 q q 1 γk l/2 ψ k+l γ l k/2 ϕ k+l, where [n] = q n q n /q q 1 and {ψ r, ϕ s r Z 0, s Z 0 } are related to {a l l Z 0 } by k=0 k=0 ψ k z k = K exp q q 1 a k z k, k=1 ϕ k z k = K 1 exp q q 1 a k z k. The standard Chevalley generators {e i, f i, t i } are given by the identification We use the coproduct k=1 t 0 = γk 1, t 1 = K, e 1 = x + 0, f 1 = x 0, e 0t 1 = x 1, t 1 1 f 0 = x + 1. e i = e i 1 + t i e i, f i = f i t 1 i + 1 f i, t i = t i t i.

6 6 M. Jimbo et al. Next recall the construction of the level one irreducible integrable highest weight representations of U in terms of bosons [4]. Let P = Z α 2, Q = Zα be the weight/root lattice of sl2, and let F[P], F[Q] be their group algebras. The basis elements of F[P] are written multiplicatively as e nα n 1 2 Z. As an F[Q]-module, F[P] = F[P] 0 F[P] 1 where F[P] i = F[Q]e iα/2. Introduce a U -module structure on the space W = F[a 1, a 2, ] F[P] in the following way. First let {a k }, e β and α act on W as We let also a k = the left multiplication by a k 1 for k < 0, = [a k, ] 1 for k > 0, e β1 f e β2 = f e β1+β2, α f e β = α, βf e β. K = 1 q α, γ = q id. The actions of the generators {x ± n } are given through the generating functions X ± z = n Z x± n z n 1 as follows. X + z = exp X z = exp a n [n] q n/2 z n exp a n [n] q n/2 z n e α z α, a n [n] qn/2 z n exp a n [n] qn/2 z n e α z α. With these actions W becomes a U -module. The submodules F[a 1, a 2, ] F[P] i are isomorphic to the irreducible highest weight modules V Λ i with the highest weight vectors u Λ0 = 1 1 and u Λ1 = 1 e α/2. The grading operator d is introduced by da n1 i 1 a nr i r e β = r n j i j + j=1 β, β 2 Λ i, Λ i a n 1 i 2 1 a nr i r e β. We remark that the trace of p d X + zx z on V Λ 0 is calculated in [8]. We shall describe the vertex operators in terms of the representation constructed above, and then compute a similar trace for them see Sect. 3. Let the components Φ ± z = that it intertwines the action of x 0 Φ µv λ± z be defined as in 1.1. From the condition we find Φ + z = [ Φ z, x 0 ] q, 2.1 where [X, Y ] q = XY qy X. The intertwining relations with the Chevalley generators determine Φ ± z uniquely, but it is not easy to get the expression of the vertex operators from them only. We better use the relations with Drinfeld s generators a n. We lack the formulas for the coproduct of them, in general. However we have the following partial information which will suffice for our purpose [9].

7 Correlation Functions of the XXZ model for < 1 7 Proposition 2.1. For k 0 and l > 0 we have k 1 x + k = x+ k γk + γ 2k K x + k + γ k+3i/2 ψ k i γ k i x + i mod N N+, 2 i=0 l 1 x + l = x+ l γ l + K 1 x + l + γ l i/2 ϕ l+i γ l+i x + i mod N N+, 2 i=1 a l = a l γ l/2 + γ 3l/2 a l mod N N +, a l = a l γ 3l/2 + γ l/2 a l mod N N +. Here N ± and N± 2 are left F[γ±, ψ r, ϕ s r, s Z 0 ]-modules generated by {x ± m k Z} and {x ± m x± n m, n Z} respectively. By using Proposition 2.1 and noting that N ± v ± = 0, N + v F[z, z 1 ]v +, we get the exact relations [a k, Φ 7k/2 [k] z] = q k zk Φ z k > 0, [a k, Φ z] = q [ Φ z, X + w] = 0, 5k/2 [k] k z k Φ z k > 0, t 1 Φ zt 1 1 = q Φ z. These conditions along with the normalization 1.2 determine the form of Φ z completely. Explicitly we have, on V Λ i i = 0, 1, Φ z = exp a n [2n] q7n/2 z n exp a n [2n] q 5n/2 z n e α/2 q 3 z α+i/ It can be checked directly that the operators given by 2.1, 2.2 enjoy the correct intertwining properties. Let us make a comment on the base field F = Qq. In the above expression of the bosonization of the vertex operators, the square root of q appears. But, it can be absorbed in the boson oscillators if we change the definition of their commutation relations. Note also that the combination α + i/2 in the power of q 3 z produces only an integer power acting on V Λ i. In any event, the appearance of q 1/2 is only superficial, and our theory is free from the ambiguity in the choice of the square root.

8 8 M. Jimbo et al. 3. Calculation of correlators As explained in Sect.1, for the evaluation of the correlators it is necessary to calculate the trace of the product of the vertex operators. Using the boson realization described in Sect.2, we shall treat the trace of the form trx d y α O where x, y are complex parameters with x < 1. Since q 2ρ = q 4d α, the choice x = q 4, y = q 1 will be relevant to the correlators. The calculation is simplified by the technique of Clavelli and Shapiro [10], Appendix C. Their prescription is as follows. Introduce a copy of bosons {b n } satisfying [a m, b n ] = 0 and the same commutation relations as the a n. Let ã n = a n 1 x n + b n n > 0, a n + b n x n n < 0. 1 For a linear operator O = O{a n } on the Fock space F a = F[a 1, a 2, ], let Õ = O{ã n } be the operator on F a F b F b = F[b 1, b 2, ] obtained by substituting ã n for a n. We have then n 2 tr Fa x do = 0 Õ 0 1 xn, where d = n>0 [n][2n] a na n and 0 Õ 0 denotes the usual expectation value with respect to the Fock vacuum 0 = 1 1 F a F b, 0 0 = 1. In order to apply their method to our case, we need the following. Set fz = zx; x q 2 zq 2 x n ; q 4 x/z; x, gz = zq 4 x n ; q 4. Let us denote the operators X z and Φ ± z with ã n substituted for a n by Jz and φ ± z. Then the explicit expressions of Jz and φ z after the normal ordering are Jz = f1exp q n/2 [n] zn a n z n b n exp q n/2 [n] φ z = g1exp exp They satisfy the following relations, z n a n xz n b n 1 x n e α z α, q 7n/2 z n a n q 5n/2 z n b n [2n] q 5n/2 z n a n q 7n/2 xz n b n e α/2 [2n]1 x n q 3 z α+i/2 on V Λ i. Jξ 1 Jξ 2 =: Jξ 1 Jξ 2 : φ z 1 φ z 2 =: φ z 1 φ z 2 : g Jξφ z =: Jξφ z : φ zjξ =: Jξφ z : φ + z = q q 1 q 2 < w <1 1 ξ 2 ξ 1 z1 g 1 q2 ξ 2 z2 ξ 1, z 2 z 1 x 1 1 q 2 w 1 fw, q 1 w 1 wfw, dξ 2πi : Jξφ z : f ξ2 ξ 1 f ξ1 ξ 2, w 1 w1 q 2 w 1 fw.

9 Correlation Functions of the XXZ model for < 1 9 Here w = ξ/q 2 z and : : denotes the normal ordering with respect to the bosons a n and b n We do not normal order the operators e nα and α. In the above equations, all factors of the form 1/1 z should be understood as n=0 zn. This fact specifies the contour for the expression φ +. Using the formulas above, 1.3 is evaluated as follows. Set P εn,,ε1 ε z n,,z 1 x, y i = q2 ; q 4 n n,,ε 1 q 4 ; q 4 n tr V Λ i x d y α ΦΛ i Λi+n 1 Λ i+1v ε z 1 Φ Λ 1 i+nv ε z n Φ Λi+nV Λi+1V n Λ i+n 1ε n z n Φ Λ iε 1 z 1 tr V Λix d y α. 3.1 We have L i z n,...,z 1 = P εn,,ε1 z n,, z 1 q 4, q 1 i. ε n,,ε 1 for L = E ε n ε n E ε 1 ε 1 see the remark at the end of Sect. 1. Introduce the following notations A = {a 1,,a s } = {j ε j = 1}, B = {b 1,,b t } = {j ε j = +1}, s + t = n, a i < a j, b i < b j for i < j, hz = q 2 z; x xz; x q 2 z 1 ; x xz 1 ; x. We prepare the integration variables ξ a a A, ζ b b B and set η j = ξ aj 1 j s, = ζ bn+1 j s < j n, η = j η j and z = j z j. Then we have P εn,,ε1 ε n,,ε 1 z n,, z 1 x, y i = 1 t q a A a<j n a A a+ b B b nn+1/2 a A z j q 2 ξ a z j ξ a b B b<j n h1n j<k hz j/z k hη j /η k j,k hη j/z k C a ζ b q 2 z j ζ b z j dξ a dζ b 2πiξ a z a b B C b 2πiζ b z b j<k η k η j η k q 2 η j m Z+i/2 z/ η2m y 2m x m2 i/4 x; x tr V Λix d y α. 3.2 Note that the last factor of the above equation can be rewritten into z η i y z/ η 2 x 1+i ; x 2 η/y z 2 x 1 i ; x 2 y 2 x 1+i ; x 2 y 2 x 1 i ; x 2. The contours of integration should be chosen as follows. Both C a and C b are anticlockwise, and the z i 1 i n lie inside of C a and outside of C b. Other relevant poles with small multipliers q, x e.g., q 2 z i are inside, and those with large multipliers are outside of the contours C a and C b. By specializing 3.2 to the case x = q 4 and y = q 1, we obtain the integral formula for the correlators. The case n = 1 gives the spontaneous staggered polarization P 0 = P P 1 1 [1]. Note that Pε ε 1 i = P ε ε 1 1 i. By using 3.2 for P and P we get P 0 = q 2 ; q 2 2 / q2 ; q 2 2.

10 10 M. Jimbo et al. Acknowledgement. We wish to thank N. H. Jing for sending his preprint [8], and H. Ooguri and J. H. H. Perk for discussions. We are grateful to Brian Davies and Omar Foda for their collaborations in our previous works, on which is based the present paper. References [1] Baxter R J, Spontaneous staggered polarization of the F model, J. Stat. Phys., 9, [2] Davies B, Foda O, Jimbo M, Miwa T and Nakayashiki A, Diagonalization of the XXZ Hamiltonian by vertex operators, RIMS preprint, 873, [3] Foda O and Miwa T, Corner transfer matrices and quantum affine algebras, Int. J. Mod. Phys. A, 7, supplement 1A [4] Frenkel I B and Jing N H, Vertex representations of quantum affine algebras, Proc. Nat l. Acad. Sci. USA, 85, [5] Baxter R J, Solvable eight-vertex model on an arbitrary planar lattice, Phil. Trans. Roy. Soc., A289, [6] Au-Yang H and Perk J H H, Critical Correlations in a Z-invariant Inhomogeneous Ising Model, Physica, 144A, [7] Drinfeld V G, A new realization of Yangians and quantum affine algebras, Soviet Math. Doklady, 36, [8] Jing N H, On a trace of q analog vertex operators, in Quantum Groups, Spring Workshop on Quantum Groups Argonne National Laboratory 16 April 11 May 1990, eds. T. Curtright, D. Fairlie and C. Zachos. World Scientific, Singapore 1992 [9] Chari V and Pressley A, Quantum affine algebras, Commun. Math. Phys., 142, [10] Clavelli L and Shapiro J A, Pomeron Factorization in General Dual Models, Nucl. Phys. B, 57,