Correlation function of the two-dimensional Ising model on a finite lattice. II A. I. Bugrij, O. Lisovyy, Bogolyubov Institute for Theoretical Physics Metrolohichna str., 14-b, Kyiv-143, 03143, Ukraine arxiv:0708.3643v1 [nlin.si] 7 Aug 007 Laboratoire de Mathématiues et Physiue Théoriue CNRS/UMR 6083, Université de Tours, Parc de Grandmont, 3700 Tours, France Abstract We calculate the two-point correlation function and magnetic susceptibility in the anisotropic D Ising model on a lattice with one infinite and the other finite dimension, along which periodic boundary conditions are imposed. Using exact expressions for a part of lattice form factors, we propose the formulas for arbitrary spin matrix elements, thus providing a possibility to compute all multipoint correlation functions in the anisotropic Ising model on cylindrical and toroidal lattices. The scaling limit of the corresponding expressions is also analyzed. 1 The model Two-point correlation function and magnetic susceptibility of the isotropic Ising model on a cylinder were calculated in [1, ]. Analogous results can be obtained in the anisotropic case as well. The complications due to the presence of two different coupling constants can be overcome, since the matrix whose determinant yields the correlation function continues to have a Toeplitz form. The computation idea has already been presented in [1], and below we will often use the results of that paper, omitting detailed calculations. The calculation of multipoint correlation functions can be reduced to the problem of finding Ising spin matrix elements the so-called form factors) in the orthonormal basis of transfer matrix eigenstates. Although neither the first nor the second problem has been solved on a finite lattice, the language of matrix elements turns out to be more convenient when constructing the corresponding expressions. A part of form factors can be found from the formulas for the two-point correlation function. Since the structure of these expressions is relatively simple, we generalize them to the case of arbitrary matrix elements. The hamiltonian of the anisotropic Ising model on a rectangular lattice is defined as H[σ] = r σr) J x x + J y y ) σr), where the two-dimensional vector r = r x, r y ) labels the lattice sites: r x = 1,,..., M, r y = 1,,..., N; the spins σr) take on the values ±1. The parameters J x and J y 1
determine the coupling energies of adjacent spins in the horizontal and vertical direction. The operators of shifts by one lattice site, x and y, are given by x σr x, r y ) = σr x + 1, r y ), y σr x, r y ) = σr x, r y + 1), where for periodic boundary conditions one has and for antiperiodic ones R) x NS) x ) M ) = R) N y = 1, ) M ) = NS) N y = 1. If the lattice is periodic in both directions, the partition function of the model at the temperature β 1 Z = e βh[σ] can be written as a sum of four terms Z = 1 [σ] Q NS, NS) + Q NS, R) + Q R, NS) Q R, R)), 1) each of them being proportional to the pfaffian of the operator 0 1 1 + t x x 1 ˆD = 1 0 1 1 + t y y 1 t x x 1 0 1, 1 1 t y y 1 0 Q = cosh K x cosh K y ) MN Pf ˆD, with different boundary conditions for x, y. Here, we have introduced the dimensionless parameters K x = βj x, K y = βj y, t x = tanh K x, t y = tanhk y. One can verify that when the torus degenerates into a cylinder M N), the partition function is determined by the antiperiodic term: Z = Q NS, NS). Let us consider the two-point correlation function in the Ising model on the cylinder σ0, 0)σr x, r y ) = Z 1 [σ] σ0, 0)σr x, r y )e βh[σ] = i P Z def Z = QNS, NS) ip def Q NS, NS), ) where Z def denotes the partition function of the Ising model with a defect: the coupling parameters K x, K y should be replaced by K x iπ/, K y iπ/ along a path that connects the correlating spins see Fig. 1 in [1], where the numbering of lattice sites and the locations of correlating spins were described; the bold line in the figure is the path along which the couplings are modified). The exponent P in the formula ) is eual to the number of steps
along the defect line: P = r x + r y in the case of a shortest path. When the correlating spins are located along a line parallel to the cylinder axis i. e. r y = 0), the ratio of pfaffians in the right hand side of ) can be expressed in terms of the determinant of a Toeplitz matrix σ0, 0)σr x, 0) = det A, 3) A kk = π π dp e ipk k) πn ) NS) t x 1 + t y + t y 1 ) e ip + t xe ip ) 1 + t x ) ) ) 1 + t y tx 1 t y cos p ty 1 t x ) cos, 4) whose size r x r x is determined by the distance between the correlating spins. Here and below the superscripts NS) and R) in sums and products imply that the corresponding operations are performed with respect to Neveu-Schwartz = π N j + 1 ), j = 0, 1,..., N 1) or Ramond = π j, j = 0, 1,..., N 1) values of uasimomenta. N Our goal is to transform 3) 4) into a representation with an explicit dependence on the distance. Ferromagnetic phase In the translationally invariant case the pfaffians Pf ˆD can be easily calculated, Q =,p cosh K x cosh K y sinh K x cos sinh K y cosp) 1/. 5) The product over any of the two uasimomentum components in 5) can be found in an explicit form. For instance, the term in 1), which corresponds to periodic boundary conditions along the x axis and antiperiodic ones along the y axis, may be written as Q R,NS) = sinh K y ) MN/ R) e Nγ)/ 1 + e Nγ) ) = = sinh K x ) MN/ NS) e Mγp)/ 1 e Mγp) ), 6) where the functions γ) and γp) are determined by the relations cosh γ) = 1 + t x )1 + t y ) t x 1 t y) cosh γp) = 1 + t x )1 + t y ) t y 1 t x) and the conditions γ), γp) > 0. We will call the domain of values of K x, K y, where sinh K x sinh K y > 1 t y1 t x ) t x 1 t y) t x1 t y ) t y 1 t x) cos, 7) cosp 8) 3
and hence γ0) = ln t x +K y ) the ferromagnetic region of parameters. Notice that in this case the numerator of the integrand in 4) can be represented in the following factorized form in terms of the function γ): t x 1 + t y ) + t y 1 ) e ip + t xe ip) = 1 t y)e γ0) 1 e γπ) e ip )1 e γ0) e ip ). Then, using the identity 1 N NS) 1 cosh θ cos Nθ tanh = sinh θ, 9) one may compute the sum over the discrete Neveu-Schwartz spectrum in 4). As a result, we can express matrix elements A kk in terms of contour integrals A kk = 1 dz k Az), 10) πi z zk where Az) = z =1 [ ] 1 e γπ) z)1 e γ0) z 1 ) 1 1 e γ0) z)1 e γπ) z 1 ) Tz), 11) Tz) = tanh Nγp), z = e ip. 1) Recall that in order to calculate the determinant A by the Wiener-Hopf method [3] one has to represent the kernel Az) in a factorized form Az) = Pz)Qz 1 ), 13) where the functions Pz) and Qz) are analytic inside the unit circle. Setting Wz) = Pz) Qz 1 ) 14) the determinant can be written as follows: where g l r x ) = A = e hrx) g l r x ), 15) l=0 hr x ) = r x [ln P0) + ln Q0)] + 1 πi 1)l l!l!πi) l z i <1 l dz i z rx i ) l 1 j=i+1 z =1 dz lnqz 1 ) ln Pz), 16) z [z i 1 z j 1 ) z i z j ) ] 1 z i 1 z j ) 4 Wz i 1 ). Wz 1 i ) 17)
We will see that the sum 15) contains a finite number of terms, i. e. there exists l 0 such that g l 0 for all l > l 0. Let us now calculate the integral 17). To do this, we first rewrite the function Tz) defined by 1) in the factorized form: Tz) = tanh Nγp) R) cosh γp) cos) = NS) cosh γp) cos) R) e γ) 1 e γ) z)1 e γ) z 1 ) = NS) e γ) 1 e γ) z)1 e γ) z 1 ) = P T z)q T z 1 ) = e Λ 1 Q T z) Q T z 1 ). Here, we have introduced the notation [ ] Λ 1 = 1 NS) R) γ) γ), Next let us define R) 1 e γ) z) Q T z) = NS) 1 e γ) z) P 0 z) = and write Wz) as follows: Wz) = P 0z) Q 0 z 1 ) 1 [ ] 1 e γπ) z 1, Q 1 e γ0) 0 z) = z P T z) Q T z 1 ) = P 0z) Q T z) Q 0 z 1 ) Tz) = 1 = 1, P T z) = e Λ 1 Q T z). 18) = [ ] 1 e γ0) z 1, 1 e γπ) z = t x1 t y ) Nγp) t y 1 t x) cosh γπ) cos p) coth Q T z) sinh γp) e γ0) γπ) Λ 1, where z = e ip. Taking into account that one can write Wz) as coth Nγp) sinh γp) = 1 N t y 1 t x) t x 1 t y ) R) 1 Wz) = cosh γπ) cosp) Q γ0) γπ) Λ T z) e 1 1 N 5 cosh γ) cosp, 19) R) 1 cosh γ) cos p.
Analogous computation for W 1 z 1 ) gives W 1 z 1 ) = cosh γ0) cosp) Q γπ) γ0) T z) e 1 N R) 1 cosh γ) cosp. Using the last two formulas, we can calculate the integral 17). Let us interchange the order of integration and summation over discrete uasimomenta. Then all the singularities inside the integration contours are exhausted by a finite number of poles at the points z k) = e γ k)), k) = πk N, k = 0, 1,..., N 1, corresponding to the Ramond values of uasimomentum. The integral 17) can then be computed by residues, and the result reads g l r x ) = e l/λ l)!n l R) l e rxγ i) η i ) sinh γ i ) G l [], where The function G l [] is given by e η i) = R) 1 e ) γ) γ i) NS) 1 e γ) γ i) ). 0) G l [] = C l l = Cl l l l 1 j=i+1 [ e γ i 1 ) e γ j 1) ) e γ i ) e γ j) ) ] [ cosh γπ) cosh γ i 1 ) [ ty 1 t x ) ] l t x 1 t y) l e γ i) l 1 j=i+1 [1 + cos i 1 )1 cos i )]. ) )] cosh γ i ) cosh γ0) ) = 1 e γ i 1 ) γ j ) [cos i 1 cos j 1 ) cos i cos j ) ] l sinh γ i)+γ j ) It can be symmetrized with respect to the permutations i j 1. It turns out that the result of such a symmetrization coincides with the even part of the function [ ty 1 t x ) ] l l t x 1 t y) 6 sin i j sinh γ i)+γ j ),
so that g l r x ) may be written as g l r x ) = e l/λ l)!n l [ ty 1 t x) t x 1 t y ) ] l R)e rx γ i) η i ) sinh γ i ) F l [], g 0 = 1, 1) F n [] = n sin i j sinh γ i)+γ j ). ) The formula ) implies that for l > N F l [] 0. By virtue of the relation 15), we then obtain the following representation for the correlation function in the ferromagnetic region: [N/] σ0, 0)σr x, 0) ) = ξξ T e rx /Λ g n r x ), 3) where ξ T = R) p n=0 ξ = 1 sinh Kx sinh K y ) 1 4, 4) R) p R) sinh γ)+γp) 3 Paramagnetic phase NS) sinh γ)+γp) NS) p NS) sinh γ)+γp) 1 4. 5) In the paramagnetic region of parameters, determined by the condition sinh K x sinh K y < 1, we have the ineuality ln t x + K y < 0, and hence γ0) = ln t x K y. Because of this, the kernel of the Toeplitz matrix whose determinant yields the correlation function should be rearranged: Az) = σ0, 0)σr x, 0) +) = A, A kk = 1 dz k 1 Az), 6) πi z zk z =1 [ 1 e γπ) z)1 e γ0) z) ] 1 1 e γ0) z 1 )1 e γπ) z 1 ) Tz). 7) To calculate the determinant of such a matrix, we must modify the method described in the previous section. For the factorized kernel Az) = Pz)Qz 1 ), 7
where the functions ln Pz) and ln Qz) are analytic inside the unit circle z < 1, we have the following general formula [1] f l+1 r x ) = A = 1)rx e hrx+1) Q 0) 1) l l!l + 1)!πi) l+1 l 1 j=i+1 z i <1 z i z j ) l j=0 l+1 f l+1 r x ), 8) l=0 j=i+1 dz i z rx i ) l+1 ) z i Wz i ) ) z i+1 Wzi+1 1 ) i=0 z i 1 z j 1 ). 9) 1 z i z j+1 ) Here, Wz) and hr x ) are defined by the formulas 14) and 16). Using the transformations analogous to those described above, one finds zwz) = z e γ0)+γπ) /Λ Q T z)cosh γ0) cos p)cosh γπ) cos p) 1 N R) 1 cosh γ) cosp, z 1 W 1 z 1 ) = z 1 e γ0)+γπ) Q Tz) 1 N R) 1 cosh γ) cos p. Using the last two relations, the integral 9) can be computed by residues: f l+1 = [ ty 1 t x ) ] l +l γ0)+γπ) l/λ+ e R) l+1 e rxγ i) η i ) t x 1 t y ) l + 1)!N l+1 sinh γ i ) l+1 F[] sinh γ i)+γ j ), where F[] = Cl+1 l ll+1) l 1 j=i+1 cos i cos j ) l+1 j=i+1 cos i 1 cos j 1 ) sin i. We again note that after symmetrization over the permutations i j+1, the function F[] coincides with the even part of the function l+1 sin i j. 8
The correlation function in the paramagnetic region of parameters is then given by g l+1 r x ) = [N 1)/] σ0, 0)σr x, 0) +) = ξξ T e rx /Λ e l+1)/λ l + 1)!N l+1 F n [] = [ ty 1 t x) t x 1 t y ) n ] l+1) / sin i j sinh γ i)+γ j ) n=0 g n+1 r x ), 30) R)e rx γ i) η i ) sinh γ i ) Fl+1 [], 31), F 1 = 1. 3) Comparing the formulas 1) 5) and 30) 3) for the correlation function with the isotropic case, we see that all the difference consists in the redefinition of the function γ) given by 7) and in the appearance of the multipliers t y 1 t x) t x 1 t y ) in the appropriate powers in the function g n r x ). In the thermodynamic limit N, these expressions reduce to the known formulas [4]. 4 Magnetic susceptibility Since we have two independent coupling parameters in the anisotropic Ising model on the cylinder, we can consider two limit cases: K x 0, K y 0 and K x 0, K y 0. The two-dimensional Ising lattice splits into mutually non-interacting closed one-row Ising chains of the length N in the first case, and into N chains of infinite length in the second case. Such limit cases provide primitive tests for the form factor representations obtained above. To perform such checks, it would also be desirable to have an expression for the correlation function σ0)σr) of two spins with arbitrary location on the lattice, i. e. for r x 0 and r y 0. The same problem arises when one tries to write the lattice analog of the Lehmann representation for the two-point Green function or to compute the magnetic susceptibility: in both cases we have to perform the Fourier transformation eir σ0)σr). Meanwhile, the matrix whose determinant gives the correlation function has a non-toeplitz form for r y 0, and the method used above to find form factor expansions becomes inapplicable. The answer can nevertheless be obtained in this case as well if we formulate the problem in terms of spin matrix elements. Recall that the N -dimensional space in which the Ising transfer matrix T acts can be splitted into two invariant subspaces the Neveu-Schwartz sector NS) and the Ramond 9
sector R). The eigenvalues corresponding to the eigenvectors from each subspace are given by λ NS) = sinh K x ) N/ exp 1 { π ) ) 3π ±γ ± γ ±... ± γ π π ) }, 33) N N N λ R) = sinh K x ) N/ exp 1 { ) π ±γ 0) ± γ ±... ± γ π π )}, 34) N N and we have different selection rules in different regions of parameters of the model: the number of minus signs in the NS sector 33) is even in both ferromagnetic and paramagnetic region, while the number of minus signs in the R sector 34) is even in the ferromagnetic region and odd in the paramagnetic one. Note that the expressions 33), 34) for the tranfer matrix eigenvalues can be obtained from the representation 1) for the partition function, if we expand the products 5): Z = Z NS) + Z R), Z NS) = 1 QNS,NS) + Q R,NS) ) = i λ NS i ) M, Z R) = 1 QNS,R) Q R,R) ) = i λ R i ) M. Because of translation invariance, the transfer matrix can be diagonalized simultaneously with the translation operator. This additional reuirement determines the orthonormal basis of eigenvectors up to permutations. It is convenient to interpret the elements of this basis in terms of NS and R multiparticle states with appropriate values of uasimomenta: they are generated from the vacua NS and R by the action of fermionic creation operators. Let us consider for definiteness the ferromagnetic region and show how to express the two-point correlation function in terms of the matrix elements NS 1,..., K ˆσ0, 0) p 1,...,p L R of Ising spin the elements NS NS and R R vanish due to Z -symmetry of the model): σ0, 0)σr x, r y ) = NS ˆσ0, 0)ˆσr x, r y ) NS = = NS ˆσ0, 0)T rxˆσ0, r y )T rx NS = = [N/+1] K=0 e rx/λ R) NS ˆσ0, 0) p 1... p K R K)! p 1...p K K K exp r x γp j ) + ir y p j ). 35) 10
Comparing this expression with the formulas 1) 3), we find spin matrix elements between the NS-vacuum and an arbitrary R eigenstate: [ ] NS ˆσ0, 0) p 1... p L R ty 1 t L = ξξ x) / L T t x 1 t y) e ηp j) Λ 1 N sinh γp j ) L sin p i p j sinh γp i)+γp j ) 36) The formula 35) then implies that form factor expansion of the correlation function σ0, 0)σr x, r y ) is obtained from 1) 3) by the substitution e rx γ i) η i ) sinh γ i ) e rx γ i)+ir y i η i ) sinh γ i ) in the function g n. An analogous result may be found in the paramagnetic case as well. Thus we have. where σ0, 0)σr x, r y ) ) σ0, 0)σr x, r y ) +) g n r x, r y ) = e n/λ n!n n = ξξ T e rx /Λ [ ] ty 1 t n x) / t x 1 t y ) [N/] l=0 [N 1)/] = ξξ T e rx /Λ R) n l=0 g l r x, r y ), 37) g l+1 r x, r y ), 38) e rx γ j)+ir y j η j ) sinh γ j ) Fn []. 39) As an illustration, consider the paramagnetic expansion for K y = 0. Notice that in this limit cosh γ) = coth K x, γ) = t x = const, and, therefore ξ = ξ T = 1, Λ 1 = η) = 0, σ0, 0)σr x, r y ) = 1 b) tx ) rx e iry = t x ) rx δ 0 ry. N As expected, the spins from different lattice rows are uncorrelated. Magnetic susceptibility of the two-dimensional Ising model on a finite M N lattice in zero field may be written as a sum of correlation functions: β 1 χ = M 1 N 1 σ0, 0)σr x, r y ). r x=0 r y=0 In order to find the susceptibility on the cylinder, we can use the expression β 1 χ = r x= r y=0 N 1 σ0, 0)σr x, r y ) 11
with some precautions []. Because of simple structure of the form factor expansions 37) 39) this sum is easily computed, and we obtain the following expressions for the susceptibility in the ferromagnetic and paramagnetic region of parameters: where R) β 1 χ ) = ξξ T n [N/] l=0 χ l, χ n = e n/λ n!n n 1 e η j) 1 sinh γ j ) F n[] coth and δ is the Kronecker symbol. 5 Scaling limit β 1 χ +) = ξξ T [N 1)/] l=0 [ ] ty 1 t n x) / t x 1 t y) n ) n γ j ) + Λ 1 δ χ l+1, ) j, 0 mod π, An important stage in the study of the two-dimensional Ising model is the analysis of its scaling limit [5, 6]. A new effect with respect to the case of the infinite plane is that the cylindrical parameters ξ T, Λ 1, η, which tend to zero for fixed distance from the critical point N, sinh K 1 = const 0), do not vanish in the scaling limit on the cylinder. Although the formulas for correlation functions obviously become more involved, it is possible to generalize several results obtained for the infinite plane to the case of the isotropic Ising model on the cylinder [7]. For instance, it was shown [8] that the Ising model correlation functions on the cylinder satisfy some integrable euations generalizing Painlevé III and Painlevé V euations obtained in the planar case. When considering the scaling limit of the anisotropic model on the cylinder, we do not encounter much difficulty provided we take into account some subtleties in the definition of scaling variables. Since in this case there are two parameters K x, K y, instead of a critical point we have the critical line sinh K x sinh K y 1 = 0. One has certain freedom in the definition of the distance from the critical line, which can be used to ensure that the results obtained in the anisotropic case in the scaling limit coincide with those obtained for the isotropic K x = K y ) model. It turns out that this condition is satisfied with the following choice of scaling variables: r x, r y, N, γ0) 0, γ0)r x = x = const, γ0)r y sinh K x = y = const, γ0)n sinh K x = β = const. 40) Because of the presence of the exponential factors e rx γ)+iry and e Nγ) in the corresponding sums and integrals, we can suppose that 1 and, therefore, ) ) γ) = γ0) 1 +, γ) = γ0) sinhk x 1 +. γ0) sinh K x γ0) 1
In order to find the scaling limit asymptotics of the cylindrical parameters Λ, ξ T and η), it is convenient to use integral representations Λ 1 = 1 π η) = 1 π π 0 π 0 ξ T = N π dp ln coth Nγp) dp π 0 cosp e γ) cosh γ) cos p euivalent to the formulas 18), 0), 5). Denoting, ln coth Nγp) dp d γ p) γ ) sinh Nγp) sinh Nγ) ln sinp + )/ sinp )/, ω) = 1 + and using the above relations, we obtain the answer for the scaled correlation function of the D anisotropic Ising model on the cylinder: σ0, 0)σr x, r y ) ) σ0, 0)σr x, r y ) +) g n x, y, β) = 1 n!β n = ξ ξ T β) e x / Λβ) = ξ ξ T β) e x / Λβ) [l] n Here, the summation is performed over integer l = [l] l 1 =, g n x, y, β), 41) n=0 g n+1 x, y, β), 4) n=0 e x ω j)+iy j η j,β)... ω j ) l n=, F n [l]. 43) the uasimomenta take on the bosonic values j = πl j, and the cylindrical parameters β and form factors are given by Λ 1 β) = η j, β) = ξ T β) = 1 π 1 π 0 0 β π dp ln coth βωp) dp 0 ω j ) βωp) ln coth, ω j ) + p, dp d ω p) ω ) sinh β ωp) sinh βω) ln p + p, 13
F n [l] = 1 n i j ω i ) + ω j ). We see that the spin-spin correlation function on the cylinder in the scaling limit defined by 40) is given by the same functions of renormalized coordinates as in the isotropic model. Therefore, all results of [8] determinant representations of the correlation fuctions, differential euations) apply to the anisotropic case as well. 6 Spin matrix elements In order to calculate multipoint correlation functions on the cylinder and torus, it is necessary and sufficient) to have formulas not only for the form factors NS ˆσ0, 0) p 1...p L R, but also for all other spin matrix elements. For instance, the two-point correlation function in the Ising model on the periodic lattice of size M N in the ferromagnetic region of parameters can be written as follows: σ0, 0)σr x, r y ) = Tr {ˆσ0, 0)T rxˆσ0, r y )T M rx } Tr T M = = { [N/+1] [N/+1] NS) K=0 L=0 1... K p 1...p L e M rx) KP γ j ) r x Λ 1 P + L γp j ) R) NS 1... K ˆσ0, 0) p 1...p L R e iry ) + e rx KP KP LP j ir y p j K)!L)! γ j ) M r x) Λ 1 P ) + L } / ) γp j ) M T Tr, λ 0 where λ 0 denotes the largest transfer matrix eigenvalue, which corresponds to the eigenvector NS. In [9], we have proposed a formula for arbitrary spin matrix elements of the isotropic Ising model on a finite periodic lattice. This formula can be easily generalized to the anisotropic case, using two hints from the above: assume that all the difference due to anisotropy consists in the redefinition of the function γ) in the formula for spin matrix element, and in the appearance of the factor t y1 t x ) to some power see, e. g., the formula 36)); t x 1 t y ) also assume that in the scaling limit 40) all multipoint correlation functions are given by the same formulas [7, 9] as in the isotropic case. Omitting the calculations, we present only the final expression for the spin matrix element, which follows from the above assumptions: NS 1... K ˆσ0, 0) p 1...p L R = ξξ T K e η j)+λ 1 N sinh γ j ) L e ηp j) Λ 1 N sinh γp j ) 14
[ ty 1 t x ) ] K L) t x 1 t y) K sin i j sinh γ i)+γ j ) L sin p i p j sinh γp i)+γp j ) 1 i K 1 j L sinh γ i)+γp j ) sin i p j. 44) As in the isotropic model case, this representation was verified for finite-row Ising chains with N = 1,, 3, 4. The scaling limit of the expression 44) obviously coincides with the classical result [11] in the limit of continuous infinite plane, and with the corresponding results on the continuous infinite cylinder [10]. 7 Discussion Lattice systems with cylindrical geometry recently became important for nanoelectronics. However, the elementary cell in these systems is hexagonal. The star-triangle transformation brings the Ising model defined on such a lattice to the Ising model on a triangular lattice. The simple form of the expressions obtained in the present paper gives a hope that similar results can also be obtained in the latter case. Exact expressions for the correlation functions are also known for several 1D uantum systems, in particular, for the XXZ Heisenberg spin chain [1]. In these models, matrix elements of local operators on the hamiltonian eigenstates were calculated using the Bethe ansatz. One can try to use our results to analyze the behaviour of such integrable onedimensional models in a finite volume and/or at non-zero temperature. We thank S. Z. Pakuliak, V. N. Roubtsov and V. N. Shadura for their interest to this work and useful discussions. This work was supported by the INTAS program under grant No. 03513350. References [1] A. I. Bugrij, Correlation function of the two-dimensional Ising model on the finite lattice. I, Theor. Math. Phys. 17, 001), 58 548; hep-th/0011104. [] A. I. Bugrij, O. Lisovyy, Magnetic susceptibility of the two-dimensional Ising model on a finite lattice, JETP 94, 00), 1140 1148 [Zh. Eksp. Teor. Fiz. 11, 00), 138 1338]; hep-th/010670. [3] B. M. McCoy, T. T. Wu, The Two-Dimensional Ising Model, Harvard University Press, 1973). [4] K. Yamada, On the spin-spin correlation function in the Ising suare lattice and the zero field susceptibility, Prog. Theor. Phys. 71, 1984), 1416 1418. [5] J. Palmer, C. A. Tracy, Two-dimensional Ising correlations: convergence of the scaling limit, Adv. Appl. Math., 1981), 39 388. 15
[6] T. T. Wu, B. M. McCoy, C. A. Tracy, and E. Barouch, Spin-spin correlation functions for the two-dimensional Ising model: Exact theory in the scaling region, Phys. Rev. B13, 1976), 316 374. [7] A. I. Bugrij, Form factor representation of the correlation function of the two dimensional Ising model on a cylinder, in Integrable Structures of Exactly Solvable Two-Dimensional Models of Quantum Field Theory, eds. S. Pakuliak and G. von Gehlen, NATO Sci. Ser. II Math. Phys. Chem. 35, Kluwer Acad. Publ., Dordrecht, 001), 65 93; hep-th/0107117. [8] O. Lisovyy, Nonlinear differential euations for the correlation functions of the D Ising model on the cylinder, Adv. Theor. Math. Phys. 5, 001), 909 9; hep-th/0108015. [9] A. I. Bugrij, O. Lisovyy, Spin matrix elements in D Ising model on the finite lattice, Phys. Letts. A319, 003), 390 394; arxiv: 0708.365. [10] P. Fonseca, A. Zamolodchikov, Ising field theory in a magnetic field: analytic properties of the free energy, J. Stat. Phys. 110, 003), 57 590; hep-th/011167. [11] B. Berg, M. Karowski, P. Weisz, Construction of Green s functions from an exact S matrix, Phys. Rev. D19, 1979), 477 479. [1] N. Kitanine, J. M. Maillet, V. Terras, Form factors of the XXZ Heisenberg spin 1 finite chain, Nucl. Phys. B554, 1999), 647 678; math-ph/980700. 16