编辑推荐

相关资源

上传者其他资源

其它热门资源

本周本月全部

TIDesigns相关资料

**资源大小：**324.95KB**上 传 者：**dingcx2008 （他上传的所有资源）**上传日期：**2015-08-20**资源类型：**应用文档**资源积分：**1分**评 论：**0条**下载次数：**0**参与讨论：**去论坛

文档简介

张量重整化文献，最近几年最新发展起来的

文档预览

PRL 101, 090603 (2008) PHYSICAL REVIEW LETTERS week ending 29 AUGUST 2008 Accurate Determination of Tensor Network State of Quantum Lattice Models in Two Dimensions H. C. Jiang,1 Z. Y. Weng,1 and T. Xiang2,3 1Center for Advanced Study, Tsinghua University, Beijing, 100084, China 2Institute of Physics, Chinese Academy of Sciences, P.O. Box 603, Beijing 100190, China 3Institute of Theoretical Physics, Chinese Academy of Sciences, P.O. Box 2735, Beijing 100190, China (Received 30 June 2008; published 29 August 2008) We have proposed a novel numerical method to calculate accurately physical quantities of the ground state using the tensor network wave function in two dimensions. The tensor network wave function is determined by an iterative projection approach which uses the Trotter-Suzuki decomposition formula of quantum operators and the singular value decomposition of matrix. The norm of the wave function and the expectation value of a physical observable are evaluated by a coarse-grain tensor renormalization group approach. Our method allows a tensor network wave function with a high bond degree of freedom (such as D ¼ 8) to be handled accurately and efﬁciently in the thermodynamic limit. For the Heisenberg model on a honeycomb lattice, our results for the ground state energy and the staggered magnetization agree well with those obtained by the quantum Monte Carlo and other approaches. DOI: 10.1103/PhysRevLett.101.090603 PACS numbers: 05.10.Cc, 71.10.Àw, 75.10.Jm The application of the density matrix renormalization group (DMRG) proposed by White [1] has achieved great success in one dimension [2,3]. However, in two dimensions, the application of the DMRG in both real and momentum space [4,5] has been limited only to small lattices. The error resulting from the DMRG truncation increases extremely fast with increasing size of lattice. To resolve this problem, the tensor network state, which is an extension of the matrix product in one dimension [6], was proposed [7,8]. In a tensor network state, the local tensors interact with each other on all directions of the lattice. This leads to two problems in the treatment of the tensor network state. First, it is difﬁcult to determine accurately all elements of local tensors by any variational approach since the total degree of freedom of a local tensor increases exponentially with the dimension of the tensor. Second, it is difﬁcult to calculate the expectation value of any physical observable even if we know the expression of the tensor network wave function, since the number of summations over the basis conﬁgurations increases exponentially with the lattice size. In this Letter, we propose a novel method to handle the tensor network wave function in two dimensions. We will show that the tensor network wave function jÉi of the ground state can be accurately determined by applying an iterative projection approach. This approach is similar to the time-evolving block decimation method that was used to determine the matrix product wave function of the ground state in one dimension [9,10]. Then we will generalize the classic coarse-grain renormalization group approach proposed by Levin and Nave [11] to the quantum system, and use it to calculate the norm of the wave function and the expectation value of any physical observable. This provides an accurate and efﬁcient tool to determine the expectation values of physical quantities from the tensor network wave function of the ground state. Below we will take the S ¼ 1=2 Heisenberg model on a honeycomb lattice as an example to show how the method works. The Hamiltonian is deﬁned by X H ¼ Hij; (1) hiji Hij ¼ JSiSj À 12h½ðÀÞiSi;z þ ðÀ1ÞjSj;z; (2) where hiji stands for summation over nearest neighboring sites and h is the magnitude of a staggered magnetic ﬁeld. It is straightforward to extend the method to other quantum lattice models with short range interactions in two dimensions. As the Hamiltonian is translational invariant, we assume the tensor network state to have the following form: Y jÉi ¼ Tr xi yi zi Axiyizi ½miBxjyjzi ½mjjmimji: i2b;j2w (3) A schematic representation of this tensor network state is shown in Fig. 1. In Eq. (3) b (w) stands for the black (white) sublattice. mi is the eigenvalue of Siz. Axiyizi ½mi and Bxjyjzj ½mj are the two three-indexed tensors deﬁned on the black and white sublattices, respectively. i ( ¼ x; y; z) are positive diagonal matrices (or vectors) of di- mension D deﬁned on the bond emitted from site i along the direction. The subscripts xi, yi, and zi are the integer bond indices of dimension D (i.e., each running from 1 to D). A bond links two sites. The two bond indices deﬁned from the two end points take the same values. For example, if the bond connecting i and j along the x direction, then xi ¼ xj. The trace is to sum over all spin conﬁgurations f. . . ; mi; mj; . . .g and over all bond indices. The ground state wave function can be determined by applying the projection operator expðÀHÞ to an arbitrary 0031-9007= 08=101(9)=090603(4) 090603-1 Ó 2008 The American Physical Society PRL 101, 090603 (2008) PHYSICAL REVIEW LETTERS week ending 29 AUGUST 2008 initial state jÉi. In the limit ! 1, expðÀHÞjÉi will converge to the ground state of H. However, this projection cannot be done in a single step since terms in H deﬁned by Eq. (2) do not commute with each other. In real calculation, we will take a small and apply this projection operator to jÉi iteratively for many times. Let us start by dividing the Hamiltonian into three parts X H ¼ Hx þ Hy þ Hz; H ¼ Hi;iþ ð ¼ x; y; zÞ: i2black H ( ¼ x; y; z) contains all the interaction terms along the direction only. These terms commute with each other. By applying the Trotter-Suzuki formula, we can express the projection operator as eÀH % eÀHz eÀHy eÀHx þ oð2Þ: (4) This means that each iteration of projection can be done using expðÀHÞ ( ¼ x; y; z) in three separate steps. In the ﬁrst step, the projection is done with Hx. As only two neighboring spins connected by horizontal bonds have interactions in Hx, the resulting projected wave function can be expressed as FIG. 1. Schematic representation of a tensor network state on a honeycomb lattice. The lattice is divided into two sublattices, represented by the black and white dots, respectively. Each vertex, on which a spin state is inhibited, is connected with three neighboring vertices along three directions, labeled by x, y, and z. On each bond, there is a diagonal matrix (or a vector), , where the subscript ¼ x, y, or z is a bond index of dimension D. At each vertex, a tensor representation of the spin state m, Ax;y;z½m for the black sublattice or Bx;y;z½m for the white sublattice, is deﬁned. A tensor network state is a product of all these bond vectors and vertex tensors. YX eÀHx jÉi ¼ Tr hm0im0jjeÀHijjmimjixi yi zi Axiyizi ½miBxjyjzj ½mjjm0im0ji: (5) i2b;j¼iþx mimj From this, a ðD2dÞ Â ðD2dÞ matrix can be deﬁned by XX S ¼ yizim0i;yjzjm0j hm0im0jjeÀHijjmimjiyi zi Axyizi ½mixBxyjzj ½mjyj zj ; (6) mimj x where d ¼ 2 is the total number of states of a S ¼ 1=2 spin. Taking the singular value decomposition for this matrix, one can further express this S matrix as X S ¼ yizimi;yjzjmj Uyizimi;x~xVxT;yjzjmj ; (7) x where U and V are two unitary matrices and ~x is a positive diagonal matrix of dimension D2d. Next we truncate the basis space by keeping only D largest singular values of ~x. Then we set the left ~x as the new x (x ¼ 1; . . . ; D) and update the tensors A and B by the following formula: Axyizi ½mi ¼ Àyi 1Àzi 1Uyizimi;x; (8) Bxyjzj ½mj ¼ Àyj1Àzj 1Vyjzjm0j;x: (9) A ﬂow chart of the above one-step renormalization of the wave function is shown in Fig. 2. The next two steps of projections can be similarly done with Hy and Hz, respectively. This completes one iteration of the projection. By repeating this iteration procedure many times, an accurate ground state wave function can then be projected out. This iteration process is very efﬁcient. The converging speed depends on the truncation error. In our calculation, we take ¼ 10À3 initially and then gradually reduce it to $10À5 to ensure the convergence of the wave function. The number FIG. 2. Flow chart of the one-step renormalization of the wave function. (a) To use expðÀHi;iþxÞ to act on the tensor network state. (b) To evaluate the S matrix deﬁned by Eq. (6). (c) To perform the singular value decomposition for S. (d) To truncate the basis space of ~x and to ﬁnd A~ and B~ with Eqs. (8) and (9), respectively. 090603-2 PRL 101, 090603 (2008) PHYSICAL REVIEW LETTERS week ending 29 AUGUST 2008 TABLE I. The ground state energy per site E and the staggered magnetization M in the zero ﬁeld limit as a function of D. D E M 3 À0:5365 0.249 4 À0:5456 0.228 5 À0:5488 0.220 6 À0:5513 0.206 7 À0:5490 0.216 8 À0:5506 0.212 of iterations used in our calculation is generally around 105–106. Given jÉi, the expectation value of a measurable quan- tity O is deﬁned by hO^ i ¼ hÉjO^ jÉi hÉjÉi : (10) We note that both hÉjÉi and hÉjO^ jÉi are tensor network functions. For example, Y hÉjÉi ¼ Tr T T ; a b i2b;j2w xix0i;yiy0i;ziz0i xjx0j;yjy0j;zjz0j (11) where the trace is to sum over all bond indices. Both Ta and Tb are D2 Â D2 Â D2 tensors. Ta is deﬁned by X Txax0;yy0;zz0 ¼ ðxyzÞ1=2Axyz½mAx0y0z0 ½mð0x0y0zÞ1=2: m (12) Tb is similarly deﬁned. Thus we can apply the tensor renormalization group method proposed by Levin and Nave [11] to evaluate hÉjÉi and hÉjO^ jÉi. To perform the tensor renormalization, we ﬁrst take two Ta and Tb on the two ends of a bond and deﬁne the following D4 Â D4 matrix: X Mll0;kk0 ¼ Tnal0kTnbk0l: n (13) By taking the singular value decomposition, one can also express this matrix as X Mll0;kk0 ¼ Ull0;nÃnVkk0;n; (14) n¼1;...;D4 where U and V are unitary matrices, Ãn is a positive deﬁned diagonal matrix of dimension D4. Again we will truncate the basis space and keep only basis states corre- sponding to the largest D2 singular values of Ã. Then the M matrix can be approximately expressed as X Mll0;kk0 % Sanll0 Sbnkk0 ; n¼1;...;D2 (15) where pﬃﬃﬃﬃﬃﬃ Sanll0 ¼ ÃnUll0;n; (16) FIG. 3. Steps of coarse graining. (a) Forming the M matrix by tracing out the common bond indices of tensors Ta and Tb deﬁned on the two neighboring sites with Eq. (14), and then perform the singular value decomposition and ﬁnd two new tensors Sa and Sb deﬁned by Eqs. (16) and (17). (b) Tracing out all common bond indices of Sa tensors (similarly for the Sb tensors) on a triangle formed by the three closed vertices to form a coarse-grained tensor T~a deﬁned by Eq. (18). pﬃﬃﬃﬃﬃﬃ Sbnkk0 ¼ ÃnVkk0;n (17) are the two vertex tensors deﬁned in the new lattice shown in Fig. 3(a). After the above transformation, the lattice structure is changed [Fig. 3(b)]. Now we replace each smallest triangle by a single lattice point. This introduces a coarse-grained honeycomb lattice with two coarse-grained tensors T~a and T~b deﬁned by X T~ a xyz ¼ SaxikSayjiSazkj ijk (18) X T~ b xyz ¼ SbxikSbyjiSbzkj: ijk (19) This coarse-grain transformation reduces the lattice by a factor of 3 at each iteration. Iterating this procedure, at the end the honeycomb lattice will eventually be reduced to having only 6 sites (Fig. 4). One can then trace out all bond indices to ﬁnd the norm of the wave function. The above coarse-grain tensor renormalization group transformation can be straightforwardly extended to evaluate hÉjO^ jÉi. The difference is that Ta and Tb now may become site dependent and their deﬁnitions are changed. FIG. 4. Tensor renormalization transformation on the honeycomb lattice. (a) The honeycomb lattice. (b) The deformed lattice in the tensor renormalization deﬁned by Eqs. (13)–(15). (c) The squeezed honeycomb lattice after the coarse-grain decimation. 090603-3 PRL 101, 090603 (2008) PHYSICAL REVIEW LETTERS week ending 29 AUGUST 2008 TABLE II. Comparison of our results with those obtained by other approaches for the ground state energy per site E and the staggered magnetization M of the Heisenberg model with h ¼ 0. Method Spin wave [12] Series expansion [13] Monte Carlo [14] Ours D ¼ 8 E À0:5489 À0:5443 À0:5450 À0:5506 M 0.24 0.27 0.22 0:21 Æ 0:01 FIG. 5 (color online). The staggered magnetization MðhÞ as a function of the staggered magnetic ﬁeld, at different D. We have applied the above approach to the spin- 1 2 anti- ferromagnetic Heisenberg model (2). Both ground state energy and the staggered magnetization M deﬁned by MðhÞ ¼ EðhÞ À Eð0Þ (20) h are calculated. In Eq. (20), EðhÞ is the ground state energy in a ﬁnite staggered magnetic ﬁeld h. The lattice size is N ¼ 6 Â 310. The ﬁnite size effect is negligible compared with the truncation error in the tensor renormalization. Table I shows the ground state energy and the staggered magnetization as a function of D for the Heisenberg model with h ¼ 0. The zero ﬁeld staggered magnetization is obtained by extrapolating MðhÞ obtained at ﬁnite h (Fig. 5) to the limit h ! 0. With D ¼ 8, we ﬁnd that the ground state energy E ¼ À0:5506 and the staggered magnetization M ¼ 0:21 Æ 0:01 in the zero ﬁeld limit. They agree well with results obtained by other approaches (see Table II). In conclusion, we have proposed a novel method to treat the tensor network wave function of quantum lattice models in two dimensions. It allows us to determine the tensor network wave function of the ground state accurately and efﬁciently. The ground state energy and the staggered magnetization of the S ¼ 1=2 Heisenberg model on the honeycomb lattice obtained with this method are consistent with those obtained by other methods. By fully considering the symmetry of the Hamiltonian, we believe that a larger local tensor with D $ 20 can be accessed with our method. This can further improve the accuracy of results and make this method even more powerful. We would like to thank Z. C. Gu and X. G. Wen for helpful discussion on the tensor renormalization group method before their work (Ref. [15]) was published. This work was supported by the NSF-China and the National Program for Basic Research of MOST, China. [1] S. R. White, Phys. Rev. Lett. 69, 2863 (1992). [2] R. J. Bursill, T. Xiang, and G. A. Gehring, J. Phys. Condens. Matter 8, L583 (1996); X. Wang and T. Xiang, Phys. Rev. B 56, 5061 (1997). [3] U. Schollwock, Rev. Mod. Phys. 77, 259 (2005). [4] T. Xiang, J. Lou, and Z. B. Su, Phys. Rev. B 64, 104414 (2001). [5] T. Xiang, Phys. Rev. B 53, R10 445 (1996). [6] S. O¨ stlund and S. Rommer, Phys. Rev. Lett. 75, 3537 (1995). [7] H. Niggemann, A. Klu¨mper, and J. Zittartz, Z. Phys. B 104, 103 (1997); Eur. Phys. J. B 13, 15 (2000). [8] F. Verstraete and J. Cirac, arXiv:cond-mat/0407066. [9] G. Vidal, Phys. Rev. Lett. 91, 147902 (2003); 93, 040502 (2004); S. R. White and A. E. Feiguin, Phys. Rev. Lett. 93, 076401 (2004); A. J. Daley et al., J. Stat. Mech. (2004) P04005. [10] G. Vidal, Phys. Rev. Lett. 98, 070201 (2007). [11] M. Levin and C. P. Nave, Phys. Rev. Lett. 99, 120601 (2007). [12] W. Zheng, J. Oitmaa, and C. J. Hamer, Phys. Rev. B 44, 11 869 (1991). [13] J. Oitmaa, C. J. Hamer, and W. Zheng, Phys. Rev. B 45, 9834 (1992). [14] J. D. Reger, J. A. Riera, and A. P. Young, J. Phys. Condens. Matter 1, 1855 (1989). [15] Z. C. Gu and X. G. Wen, arXiv:0806.3509. 090603-4

相关帖子

回到顶部

EEWORLD下载中心所有资源均来自网友分享，如有侵权，请发送举报邮件到客服邮箱bbs_service@eeworld.com.cn 或通过站内短信息或QQ：273568022联系管理员 高进，我们会尽快处理。