arxiv: v4 [cond-mat.stat-mech] 7 Aug 2012

Transcription

1 Coarse granng tensor renormalzaton b the hgher-order sngular value decomposton arxv: v4 [cond-mat.stat-mech] 7 Aug 2012 Z. Y. Xe 1, J. Chen 2, M. P. Qn 2, J. W. Zhu 3, L. P. Yang 4 and T. Xang 2,1 1 Insttute of Theoretcal Phscs, Chnese Academ of Scences, P.O. Bo 2735, Bejng , Chna 2 Insttute of Phscs, Chnese Academ of Scences, P.O. Bo 603, Bejng , Chna 3 Insttute of Computatonal Mathematcs and Scentfc/Engneerng Computng, Academ of Mathematcs and Sstems Scence, Chnese Academ of Scences, Bejng , Chna and 4 Bejng Computatonal Scence Research Center, Bejng, , Chna (Dated: August 8, 2012) We propose a novel coarse granng tensor renormalzaton group method based on the hgherorder sngular value decomposton. Ths method provdes an accurate but low computatonal cost technque for studng both classcal and quantum lattce models n two- or three-dmensons. We have demonstrated ths method usng the Isng model on the square and cubc lattces. B keepng up to 16 bond bass states, we obtan b far the most accurate numercal renormalzaton group results for the 3D Isng model. We have also appled the method to stud the ground state as well as fnte temperature propertes for the two-dmensonal quantum transverse Isng model and obtan the results whch are consstent wth publshed data. PACS numbers: Cc, w, Hk I. INTRODUCTION The smulaton of two or hgher dmensonal quantum lattce models remans a great challenge. Ths has stmulated great nterest on the nvestgaton of renormalzaton group (RG) methods for the tensor-network states The use of the tensor-network state as a varatonal wave functon for the classcal lattce model was frst consdered b Nshno and coworkers The, and recentl Garca-Saez et. al 17, proposed a number of RG approaches to stud the thermodnamc propertes of the Isng and other models. However, due to the heav computatonal cost, the mamal truncated tensor dmenson, D, that can be handled wth ther methods s small (between 2 and 5) n 3D, and consequentl the accurac of the results the obtaned s low n comparson wth the Monte Carlo ones. In 2007, Levn and Nave 3 proposed a coarse graned tensor renormalzaton group (TRG) method for studng two dmensonal (2D) classcal models based on the sngular value decomposton (SVD) of matr. Later we proposed a second renormalzaton group (SRG) method 7,8 to globall optmze the truncaton scheme and mprove sgnfcantl the accurac of the TRG. The applcaton of these methods n classcal and quantum lattce models has acheved great success 4,5,7 9. However, t s dffcult to etend these methods to 3D, not just due to the ncrease of the order of local tensors, but also due to the change of lattce topolog n the coarse granng process 18. In ths paper, we ntroduce a novel coarse granng TRG method based on the hgher-order sngular value decomposton (HOSVD) 20 to stud phscal propertes of 2D or 3D lattce models. We wll frst dscuss about a smple TRG method based on the HOSVD (abbrevated as HOTRG hereafter), and then dscuss about a more sophstcated method that ncorporates the second renormalzaton effect of envronment tensors to the HOTRG. The HOSVD takes nto account more accuratel the nterpla between dfferent components of a tensor. It provdes a better scheme to truncate a local tensor than the SVD. Ths paper s arranged as follows. In Sec. II, a detaled ntroducton to the HOTRG and HOSRG for the 2D statstcal lattce models s gven. We have taken the 2D Isng model to show how accurate the HOTRG and HOSRG can be n comparson wth other methods. In Sec. III, we have etended the HOTRG and HOSRG to the 3D statstcal lattce models. For the 3D Isng model, we have obtaned b far the most accurate numercal renormalzaton group results for the 3D Isng model. Our accurac s comparable wth the best Monte Carlo results. In Sec. IV, we have appled HOTRG to 2D quantum lattce models. Our prelmnar results show that the HOTRG provdes a powerful tool for studng the ground state and thermodnamc propertes of 2D quantum lattce models. A summar s gven n Sec. V. II. TWO-DIMENSIONAL SYSTEMS A. HOTRG Let us start b takng the Isng model, H = j σ zσ j z, (1) as an eample to show how the method works n 2D frst. σ z s the paul matr at ste. An etenson of the method to 3D wll be descrbed later. The partton functon of the 2D Isng model can be represented as a

2 2 (a) M (n) T (n+1) (b) ' ' ' M (n) ' 1 ' 2 T (n+1) ' 1 ' 2 ' FIG. 1: (a) A HOTRG contracton of the tensor network state along the as on the square lattce. (b) Steps of contracton and renormalzaton of two local tensors. The ntal tensor T (0) = T. translaton nvarant tensor network state 8, Z = Tr ' T, (2) where runs over all the lattce stes and Tr s to sum over all bond ndces, and the local tensor T s defned at each lattce ste as shown n Fg. 1(a), T = α W α, W α, W α, W α, (3) where W s a 2 2 matr defned b ( ) cosh(1/t), snh(1/t) W =, (4) cosh(1/t), snh(1/t) and T s the temperature. To coarse gran, we contract the lattce alternatvel along the horzontal (-as) and vertcal (-as) drectons. Ths scheme of coarse granng s smple to mplement. Fg. 1(a), as an eample, shows how the contracton along the -as s done. At each step, two stes are contracted nto a sngle ste n the coarse graned lattce (Fg. 1(b)), and the lattce sze s reduced b a factor of 2. The contracted tensor at each coarse graned lattce ste s defned b M (n) = 1 1 T(n) 2, (5) 2 where = 1 2, = 1 2, and the superscrpt n denotes the n th teraton. The bond dmenson of M (n) along the -as s the square of the correspondng bond dmenson of. To truncate M (n) nto a lower rank tensor, we frst do a HOSVD for ths tensor 20 M (n) = jkls jkl U L UR j UU k UD l, (6) where U s are the untar matrces. S s the core tensor of M (n), whch possesses the followng propertes for an nde, sa nde j: (1) all orthogonalt, S :,j,:,: S :,j,:,: = 0, f j j, (7) where S :,j,:,: S :,j,:,: s the nner-product of these two sub-tensors. (2) pseudo-dagonal, S :,j,:,: S :,j,:,:, f j < j, where S :,j,:,: s the norm of ths sub-tensor whch s the square root of all elements square sum. These norms pla a smlar role as the sngular values of a matr. In M (n), the two vertcal bonds, and, do not need to be renormalzed. Thus n the practcal calculaton, U U and U D are not needed to be determned. Moreover, the rght bond of M (n) s lnked drectl to the left bond of an dentcal tensor on the rght neghborng ste, thus to truncate an one of the horzontal bonds of M (n) wll automatcall truncate the other horzontal bond. The truncaton can be done b comparng the values of and ε 1 = >D ε 2 = j>d S,:,:,: 2 (8) S :,j,:,: 2. (9) If ε 1 < ε 2, we truncate the frst dmenson of S or the second dmenson of U L to D. Otherwse, we truncate the second dmenson of S or the second dmenson of U R to D. Ths knd of truncaton scheme provdes a smple and optmal appromaton to mnmze the truncatonerror 21,22. Ithasbeensuccessfullappledtoman felds such as data compresson, mage processng, pattern recognton, and etc 23. After the truncaton, we can update the local tensor usng the followng formula T (n+1) = j M (n) j U(n+1) j, (10) where = U L (or U R ) f ε 1 s smaller (or larger) than ε 2. The above HOTRG calculaton can be repeated teratvel untl the free energ and other phscal quanttes calculated are converged. The cost of the calculaton scales as D 7 n the computer tme and D 4 n the memor space. Ths s comparable wth the cost of TRG 7,8.

3 3 1 2 a k l E (n) j 1 j 2 E (n+1) FIG. 2: (color onlne) Graphcal representaton of Eq. (13) for determnng the envronment tensor E (n) kaj 1 1 from E (n+1) jkl n the backward teraton. ThekestepntheaboveHOTRGteratonstodetermne the four untar matrces on the rght hand sde of Eq. (6). In our calculaton, we determne these matrces b takng the sngular value decomposton of matrces. As an eample, let us consder how to evaluate U L. We frst convert M nto a matr M, M, = M wththefrstndeastherowndeandtherestndces (,, ) as the column nde of ths matr. Then from the theor of HOSVD, we know that U L s equal to the left untar matr of M under the sngular value decomposton. Thus U L can be smpl determned from the canoncal transformaton of the untar matr M M M M = U L Λ L (U L ), (11) where Λ L s the egenvalue of M M. Furthermore, t can be shown that S,:,:,: 2 = Λ L. (12) The cost for evaluatng these U-matrces scales wth D 6. B. HOSRG The HOTRG s a local optmzaton method. It mnmzes the error n the truncaton of a local tensor. However, t gnores the renormalzaton effect of envronment. To develop a global optmzaton method, t s necessar to consder the envronment contrbuton n the renormalzaton of local tensors. In Refs. 7,8, we proposed a SRG appraoch to ncorporate the envronment contrbuton n the optmzaton of local tensors. Ths knd of SRG approach can be also used to mprove the performance of HOTRG, whch leads to a global optmzed HOTRG method, referred as HOSRG below. The HOSRG follows the same coarse granng steps as n the HOTRG. However, at each step, one needs to j calculate a bond denst matr defned on a bond whose bass space wll be truncated. Ths bond denst matr s defned b tracng out all envronment tensors. The SRG ntroduced n Ref. 7,8 s an nfnte lattce algorthm. In the calculaton of the bond denst matr, the sze of the envronment s alwas assumed to be nfnte. At each step of coarse granng, a combned forward and backward teraton s performed. In the forward teraton, the TRG s appled to determne all transformaton matrces and local tensors. A backward teraton s then performed to determne the bond denst matr. Ths scheme can be readl etended to the HOSRG. We have done ths knd ofcalculatonand fnd that t doesprovde much more accurate results than the HOTRG. Smlar as n the DMRG, one can also ntroduce a fnte lattce algorthm to perform the HOSRG calculaton. In ths case, the whole sze of the sstem s fed and the number of tensors n the envronment s reduced at each step of coarse granng. Agan the second renormalzaton effect of the envronment s handled b performng forward-backward teratons. But now ths knd of forward-backward teratons can be repeated for man tmes, smlar as to do a fnte sze sweepng n the DMRG. Ths provdes a self-consstent approach to treat the sstem as well as envronment tensors. It can further mprove the accurac of the HOSRG. Below we gve an ntroducton to ths method. In the frst round of forward teraton, we carr out a standard HOTRG calculaton to determne teratvel all the transformaton matrces U (n) and local tensors. Ths teraton ends when the sstem reaches a desred sze, sa 2 N lattce ponts wth N = One then carres out a backward teraton to calculate the envronment tensor E (n) teratvel, startng from E (N+1) whch s set to be an unt tensor. The teraton formula for determnng E (n) s gven b E (n) kaj 1 1 = jl 2j 2a E (n+1) jkl 2j 2al U(n+1) 1 2, U(n+1) j 1j 2,j. (13) A graphcal representaton of ths equaton s shown n Fg. 2. Ths backward teraton s termnated after E (2) s determned. The above forward-backward teraton determnes all coarse graned sstem and envronment tensors that are needed for carrng out the second renormalzaton calculaton. From these tensors we can perform another forward and backward teraton to mprove ther accurac. But startng from ths round of teraton, the HOTRG s no longer needed. Instead, at each step of coarse granng we evaluate the bond denst matr ρ zw, usng the followng formula ρ zw, (n) = E(n+2) jkl 1 2 j 1j 2j U(n+2) k 1k 2k U(n+2) l 1l 2l 1k 1a T(n) 2al 1 zj 1k 2b T(n) wj 2bl 2. (14) The repeated nde summaton s assumed. A graphcal representaton of ths formula s shown n Fg. 3.

4 4 1 2 l 1 a k U (n+2) k 1 k 2 l z b w j 1 j 2 l 2 E (n+2) U (n+2) j (a) M (n) (b) U (n) V (n) z z' z z' M(n) ' 1 ' 1 ' 2 ' 2 V (n) ' U (n) ' FIG. 3: (color onlne) Graphcal representaton of Eq. (14) for determnng the bond denst matr ρ (n) zw, through E (n+2) jkl. Relatve error of free energ TRG HOTRG SRG HOSRG Temperature FIG. 4: (color onlne) Comparson of the relatve errors of free energ wth respect to the eact results for the 2D Isng model obtaned bvarous methods wth D = 24. The crtcal temperature T c = 2/ln(1+ 2). T (n+1) z T (n+1) z' ' ' FIG.5: (a)ahotrgcoarse granngstepalongthez-ason the cubc lattce. (b) Steps of contracton and renormalzaton of two local tensors. result s alread less than 10 7 even at the crtcal temperature, much more accurate than the TRG result 7,8. The HOSRG also performs better than the SRG. But the dfference n the results obtaned b these two methods s relatvel small around the crtcal pont. The HOTRG s lessaccuratethanthetwosrgmethods, buttscomputatonall economc. The dfference between TRG/SRG and HOTRG/HOSRG les manl n the bass truncaton scheme. The formers basedon the SVD, whlethe latter s based on the HOSVD. The above comparson ndcates that the HOSVD scheme works better. To dagonalze ths bond denst matr 25 ( ρ (n) = Λ ), (15) wecanfndtsegenpar,(λ, ). Sameasnthedenst matr renormalzaton group 24, the egenvalues of ths denst matr determne the probabltes of the correspondng egenvectors n the vrtual bond bass space. B keepng the largest D egenvalues of Λ and the correspondng egenvectors of, one can update the local tensor T (n+1) usng Eq. (10). After fnshng ths forward teraton, we can take a backward teraton to update all envronment tensors wth Eq.(13). Ths forwardbackward teraton s then repeated untl all sstem and envronment tensors are converged. Fg. 4 compares the relatve errors of free energ wth respect to the rgorous soluton 26 for the 2D Isng model obtaned wth four dfferent methods. B keepng just 24 states, we fnd that the relatve error of the HOTRG III. THREE-DIMENSIONAL SYSTEMS TheaboveHOTRGand HOSRGmethods canbe readl etended to 3D. Ths s an advantage of the coarse granng scheme proposed here. On the cubc lattce, a full ccle of lattce contracton needs to be done n three steps, along the -as, -as, and z-as, respectvel. At each step, two neghborng tensors wll be combned to form a sngle coarse graned tensor and the lattce sze s reduced b a factor of 2. As an eample, Fg. 5 shows how the tensors are contracted along the z-as. The HOSVD of the coarse graned local tensor (Fg. 5(b)) can be smlarl done as for the 2D case. But the local tensor now has s bond ndces and a HOSVD for a hgher order tensor should be done. Moreover, the bass spaces for both the -as and -as bonds need to be renormalzed. Thus we should determne from the core tensor and the untar matrces of M (n) not onl the transformaton matr for the -

5 5 l V (n+1) b n u E (n) k d m j f E (n+1) FIG. 6: (Color onlne) Graphcal representaton for the determnaton of the envronment tensor E (n) mnjuk from E(n+1) lrfbud n 3D. drectonbonds, U (n), but alsothetransformatonmatr for the -drecton bonds, V (n). After that the dmensons for both -as and -as bonds are truncated and thelocaltensorsupdatedusngu (n) andv (n). Thecontracton and renormalzaton of tensors along other two drectons can be smlarl done. Ths three-step teraton can then be repeated untl the results are converged. After the above HOTRG teraton, one can also do a backward teraton to evaluate the envronment tensors and carr out the HOSRG calculaton n 3D. A graphcal representaton for teratvel determnng the envronment tensor n ths backward teraton s shown n Fg. 6. A seres of forward-backward teratons s then performed to take nto account the second renormalzaton effect of the envronment to the coarse graned tensors. In the subsequent forward teratons, we evaluate and dagonalze the bond denst matr (see Fg. 7) and update the coarse graned tensors. The envronment tensors are evaluated agan n the backward teraton. In the 3D calculaton, the computatonal tme scales wth D 11 and the memor scales wth D 6. Ths cost n the computatonal resource s sgnfcantl smaller than n other 3D numercal RG methods 11 17,19. We have studed the 3D Isng model usng the HOTRG for D up to 16. The temperature dependence of the nternal energ U and the specfc heat C for the 3D Isng model obtaned b the HOTRG wth D = 14 s shown n Fg. 8 and compared wth the Monte Carlo result 27. Our result of the specfc heat agreeswth the MonteCarloone. At the crtcal temperature, T c = , the nternal energ s found to be U c = for D = 14. Ths value of U c, as shown n Table I, also agrees well wth other publshed data. From the temperature dependence of the specfc heat around the crtcal pont, one can estmate the crtcal eponent of the specfc heat wth the formula C t α (16) r where t = 1 T/T c. However, as the specfc heat data are obtaned smpl from the numercal dervatve of the nternal energ, the accurac of the specfc heat data s much less than that of the nternal energ, especall around the crtcal pont. Ths causes a bg error n the determnaton of the eponent α wth the above formula. Ths problem can be solved b drectl evaluatng ths eponent from the temperature dependence of the nternal energ. From the temperature ntegraton of the specfc heat, t s smple to show that the nternal energ should ehbt the followng crtcal behavor U = U c +at+bt 1 α, (17) where a and b are unknown parameters whch can be determned b fttng. Fg. 9 shows the fttng curves for the nternal energ around the crtcal pont obtaned wth Eq. (17). The crtcal eponent s found to be α = and for the temperature hgher and lower than the crtcal value, respectvel. These values of the crtcal eponent are consstent wth the result obtaned from the seres epanson 28, 0.104, and the Monte Carlo calculaton 29, Fg. 10 shows the temperature dependence of the spontaneous magnetzaton M obtaned b the HOTRG wth D=14. Our data agree well wth the Monte Carlo results 35. From the sngular behavor of M, we fnd that the crtcal temperature T c = for D = 14. Furthermore, b fttng the data of M n the crtcal regme wth the formula M t γ, (18) wefnd thatthe eponentγ = ,consstentwth the Monte Carlo 29 (0.3262) and seres epanson 36 (0.3265) results. l V (n+1) b V (n+2) U (n+2) d u z w E (n+2) FIG. 7: (Color onlne) Graphcal representaton for the determnaton of the bond denst matr ρ (n) zw, from the envronment tensor E (n+2) lrfbud n 3D. f r

6 Internal Energ U Temperature Specfc Heat Internal Energ U HOTRG Fttng curve ( =0.1023, T > T c ) Fttng curve ( =0.1137, T < T c ) Temperature FIG. 8: (Color onlne) The nternal energ and the specfc heat for the 3D Isng model obtaned b the HOTRG wth D = 14. The Monte Carlo result (black curve) obtaned from an emprcal ft formula gven n Ref. [27] s shown for comparson. Fg. 11 shows the crtcal temperature T c determned from the sngular ponts of the nternal energ as well as the magnetzaton for D up to 16. The values of T c obtaned from these two quanttes agree wth each other. ForD = 16,Tcobtanedfromthenternalenergandthe magnetzaton are and , respectvel. The relatve dfference s less than But T c does not var monotoncall wth D. It becomes converged onl when D 13, ndcatng the mportance of keepng a large D n the 3D TRG calculaton. The error n T c, estmated fromthe dfference between the valuesoft c for D = 15 and D = 16, s also less than Our results agree wth the Monte Carlo data The above dscusson ndcates that the HOTRG works ver well n 3D. The accurac of the results can be further mproved b applng the HOSRG. However, the HOSRG calculaton costs much more CPU tme. A thorough stud wth the HOSRG on the 3D Isng model s stll n progress and the results wll be publshed sepa- TABLE I: Comparson of the nternal energ at the crtcal temperature U c for the 3D Isng model obtaned b dfferent methods. method U c HOTRG(D = 16) (3) Seres epanson (1) Seres epanson (1) Seres epanson (15) Monte Carlo (4) Monte Carlo (8) Monte Carlo (4) FIG. 9: (Color onlne) The nternal energ (D = 14) and ts fttng curves wth Eq. (17) around the crtcal pont for the 3D Isng model. α s the crtcal eponent for the specfc heat. ratel. IV. GROUND STATE AND THERMODYNAMICS OF 2D QUANTUM LATTICE MODELS A d-dmensonal quantum lattce model s equvalent to a (d+1)-dmensonal classcal model, the HOTRG and HOSRG methods above ntroduced can be also etended to stud the ground state and thermodnamc propertes of d-dmensonal quantum lattce models. For one- Magnetzaton M HOTRG Monte Carlo 10 0 ft (T c -T)/T c Temperature FIG. 10: (color onlne) Temperature dependence of the magnetzaton for the 3D Isng model (D = 14). The Monte Carlo result s from Ref. [35]. Inset: Logarthmc plot the magnetzaton around the crtcal pont. The slope of the fttng curve gves the crtcal eponent of the magnetzaton, γ =

7 7 T c from U from M D Ground State Energ E E 0 M z M h Magnetzaton FIG. 11: (Color onlne) The crtcal temperature T c as a functon of the bond dmenson D for the 3D Isng model obtaned from the nternal energ (U) and magnetzaton (M), respectvel. dmensonal quantum lattce models, there are alread man mature methods for studng the ground state as well as the thermodnamc propertes. For eample, the ground state can be studed b the DMRG 24 and the thermodnamcs can be studed b the quantum transfer matr renormalzaton group (TMRG) 44,45. Here we wll onl dscuss how to appl the HOTRG/HOSRG to a 2D quantum lattce model. As an eample, we wll take the 2D quantum Isng model wth a transverse feld to show how these methods work. The Hamltonan of ths model s defned b H = j σ zσ j z h σ. (19) TABLE II: Comparson of the crtcal pont T c for the 3D Isng model obtaned b dfferent methods. method HOTRG(D=16, from U) HOTRG(D=16, from M) Monte Carlo Monte Carlo Monte Carlo Monte Carlo Seres epanson CTMRG TPVA CTMRG TPVA Algebrac varaton T c FIG. 12: (Color onlne) The ground state energ E 0, the magnetzaton M σ and M z σ z versus the appled feld h for the 2D quantum Isng model obtaned b the HOTRG wth D = 14. We start b representng the partton functon of ths model as a tensor-network model n the 2+1 dmensons. B usng the Trotter-Suzuk decomposton formula, we can epress the partton functon as 19 where Z = Tre βh Tr [ e τhz e τh] L +O(τ 2 ) (20) H z = j σ zσ j z, (21) H = h σ. (22) β = Lτ sthenversetemperatureandτ sasmalltrotter parameter. Ths partton functon can be also epressed as a product of evoluton matr V where Z = TrV L, (23) V = e τh/2 e τhz e τh/2 (24) TABLE III: The crtcal feld h c for the 2D quantum Isng model wth transverse feld obtaned b dfferent methods. method h c HOTRG(D=14) Monte Carlo Seres Epanson PEPS VDMA TERG PEPS CTM

8 8 Internal Energ U h = 2 h = 3 h = T Specfc Heat C T h = 2 h = 3 h = 4 FIG. 13: (Color onlne) Temperature dependence of the nternal energ n three dfferent felds. The sold dots and open crcles are obtaned wth the 3D coarse granng HOTRG wth D = 14 and the magnar tme evoluton approach wth D = 24, respectvel. s the evoluton operator between two neghborng Trotter laers. To nsert the complete bass set between an two of the eponental terms on the rght hand sde of Eq. (20), t s straghtforward to show that V can be epressed as a product of local tensors. From ths, we can epress the partton functon as a 3D tensor-network model Z Tr T lr f b µ ν +O(τ 2 ), (25) where the s-ndeed local tensor s defned b wth and T lrfbµν = σ W = P = 1 2 ( W σlw σrw σfw σbp σµ P σν, (26) ( coshτ, ) snhτ, (27) coshτ, snhτ ) e τh/2, e τh/2 e τh/2, e τh/2. (28) The operator V governs the bass state evoluton along the Trotter drecton. The matr element of V between two sets of bass states, {ν } and {µ }, s defned b tracng out all spatal ndces of local tensors T n a gven Trotter laer {ν } V {µ } Tr T lr f b µ ν, (29) where Tr s to trace out all spatal ndces onl. The transverse feld h appears onl n the P-matr. When h = 0, t s smple to show that PP = I. In FIG. 14: (Color onlne) The specfc heat versus temperature obtaned b the magnar tme evoluton approach wth D = 24 for the 2D quantum Isng model. ths case, the Hamltonan returns to the 2D classcal Isng model and Eq. (20) s reduced to a product of L 2D tensor-network model for the Isng model wth an nverse temperature τ. At zero temperature, L, the partton functon s a product of nfnte man tensors along all three drectons. The 3D HOTRG or HOSRG method ntroduced prevousl can be drectl appled to stud phscal propertes of the ground state. In ths case, there s no need to evaluate the ground state wavefuncton. Ths s an advantage of ths approach n comparson wth other methods whch are based on the tensor-network representaton of ground state wavefuncton. Fg. 12 shows the feld dependence of the ground state energ E 0, the transverse magnetzaton M and the longtudnal magnetzaton M z for the transverse Isng model obtaned b the HOTRG wth D = 14 and τ = In agreement wth other calculatons, we fnd that ths model ehbts a phase transton at a fnte feld. The crtcal feld s found to be h c = , consstent wth other publshed data, as shown n Table III. At fnte temperature, the lattce dmenson of the tensor network model s fnte along the magnar tme drecton. Now two approaches can be used to evaluate the partton functon. The frst s to follow the steps of 3D HOTRG to contract the tensors alternatvel along the three drectons. The contracton along the Trotter drecton s termnated f all bond varables along that drecton are contracted. The teraton s then carred out purel along the two spatal drectons as for a pure 2D classcal model. Ths s an accurate and effcent approach for evaluatng phscal quanttes. However, the number of temperature ponts that can be studed wth ths approachsqute lmted foragvenτ, snce the temperature s reduced b a factor of 2 at each contracton along the Trotter drecton. The second approach s to do the magnar tme evo-

9 9 M h = 4 h = 3 M z h = 2 h = h = FIG. 15: (Color onlne) Transverse magnetzaton as a functon of temperature for the 2D quantum Isng model. The sold dots and open crcles are obtaned wth the 3D coarse granng HOTRG wth D = 14 and the magnar tme evoluton approach wth D = 24, respectvel. T 0.2 h = FIG. 16: (Color onlne) Longtudnal magnetzaton as a functon of temperature for the 2D quantum Isng model. The sold dots and open crcles are obtaned wth the 3D coarse granng HOTRG wth D = 14 and the magnar tme evoluton approach wth D = 24, respectvel. T luton. In ths case, one starts from a one-trotter laer tensor product sstem, whose tensor operator s defned b V. At each step of evoluton, one more Trotter laer (.e. onev operator)saddedtothesstemandthetruncaton for the spatal bond dmenson s done usng the HOTRG. B performng ths magnar tme evoluton up to certan temperature, phscal quanttes are then evaluated at that temperature b tracng out all bond varables wth the HOTRG. The nverse temperature ncreases lnearl wth the evoluton number. Ths allows us to collect more temperature ponts. However, as the truncaton error s accumulated wth the evoluton, the results such obtaned would become less and less accurate wth decreasng temperature. Thus ths approach should be used onl for evaluatng thermodnamc quanttes n hgh temperatures. A smlar magnar tme evoluton approach was recentl proposed for evaluatng thermodnamc quanttes b Ran et al. 46 based on the bond entanglement mean feld approach proposed n Ref. 4. In our calculaton, both approaches have been used. The frst approach (.e. the 3D coarse granng HOTRG) s more accurate than the second one (.e. the magnar tme evoluton approach) especall n low temperatures. It s used for collectng the low temperature data. The second approach allows more temperature ponts to be evaluated and s appled to evaluate thermodnamc quanttes n hgh temperatures. Fgs show the nternal energ, the specfc heat, the transverse and longtudnal magnetzaton as a functon of temperature for the 2D quantum Isng model n three dfferent appled felds, respectvel. The sold and open crcles are results obtaned b the 3D HOTRG (D = 14) and the magnar tme evoluton (D = 24) approaches, respectvel. The results obtaned wth these two approaches agree wth each other n the ntermedated temperature range. In hgher temperature, the truncaton error of the magnar tme evoluton s relatvel small and the results obtaned wth ths method s more accurate. In low temperature the results obtaned b the 3D HOTRG are much more accurate than those obtaned b the magnar tme evoluton, snce the error s accumulated at ever step of coarse granng or tme evoluton and the 3D HOTRG can reach low temperatures n a much fewer coarse granng steps. In all calculatons, the Trotter step s set to τ = When h < h c, as the ground state s spontaneousl smmetr broken wth a fnte magnetc order, a fnte temperature phase transton s epected. Such phase transton s confrmed b our calculaton, whch can be clearl seen from the temperature dependence of the longtudnal magnetzaton M z (Fg. 16). The crtcal temperaturedecreaseswth ncreasngfeldhforh<h c. The specfc heat s obtaned from the numercal dervatve of the nternal energ shown n Fg. 13. It shows a sngular behavor around the crtcal pont for h = 2. The above dscusson ndcates that the HOTRG (ncludng the magnar tme evoluton approach) provdes a smple and powerful method for studng the ground state and fnte temperature propertes of 2D quantum lattce models. B applng the HOSRG to the above calculaton, we beleve that the accurac of results can be further mproved. But ths takes a longer tme to do the calculaton. V. SUMMARY In summar, we have proposed the HOTRG and HOSRG methods for studng phscal propertes of classcalor quantum lattce models n 2Dor3D. B comparsonwththeeactsolutonofthe2disngmodel,wehave shown that the smple HOTRG calculaton can alread

10 10 produce ver accurate numercal results. The HOSRG takes nto account the second renormalzaton effect of the envronment tensors, and t can sgnfcantl mprove the accurac of the HOTRG. These method allow us to retan an unprecedentedl hgh bond dmenson n the bass truncaton and eld b far the most accurate numercal RG results for the 3D Isng model. We have also appled the HOTRG to stud the ground state and thermodnamc propertes of the 2D quantum Isng model wth a transverse magnetc feld. Our results agree well wth all publshed data. Smmetr or good quantum number of the tensor network state for these models can be used to reduce the computatonal and storage cost. Ths wll allow us to retan more bass states to reduce the truncaton error and further mprove the accurac of results. In ths work, we have taken two translaton nvarant tensor network models, namel the Isng model and the transverse Isng model, as eamples to show how the methods work. However, t should be emphaszed that our methods work more generall. The can be etended to a sstem whch s translaton nvarant b shftng two or more lattce stes, or even to a random sstem, such as a spn glass model. B combnng wth the local update method of quantum tensor product wavefuncton ntroduced n Ref. 4, one can also use ths method to stud ground state propertes of 3D quantum lattce models. We wsh to thank H.W.J. Blote for sendng us the Monte Carlo result shown n Fg. 8. We are ndebted to Prof. T. Nshno for helpful dscusson. Ths work was supported b NSFC (Nos and ) and MOST 973 Project (No. 2011CB309703). Electronc address: tang@aph.ph.ac.cn 1 H. Nggemann, A. Klumper, and J. Zttartz, Z. Phs. B 104, 103 (1997). 2 F. Verstraete and J. Crac, arxv: M. Levn and C. P. Nave, Phs. Rev. Lett. 99, (2007). 4 H. C. Jang, Z. Y. Weng, and T. Xang, Phs. Rev. Lett. 101, (2008). 5 Z. C. Gu, M. Levn, X. G. Wen, Phs. Rev. B 78, (2008). 6 J. Jordan, R. Orus, G. Vdal, F. Verstraete, and J. I. Crac, Phs. Rev. Lett. 101, (2008). 7 Z. Y. Xe, H. C. Jang, Q. N. Chen, Z. Y. Weng, and T. Xang, Phs. Rev. Lett. 103, (2009). 8 H. H. Zhao, Z. Y. Xe, Q. N. Chen, Z. C. We, J. W. Ca, and T. Xang, Phs. Rev. B 81, (2010). 9 Q. N. Chen, M. P. Qn, J. Chen, Z. C. We, H. H. Zhao, B. Normand, and T. Xang, Phs. Rev. Lett. 107, (2011) 10 I. Pzorn, L. Wang, F. Verstraete, Phs. Rev. A 83, (2011). 11 T. Nshno and K. Okunsh, J. Phs. Soc. Jpn. 67, 3066 (1998). 12 K. Okunsh, T. Nshno, Prog. Theor. Phs, 103, 541 (2000). 13 T. Nshno, K. Okunsh, Y. Heda, N. Maeshma, Y. Akutsu, Nucl. Phs. B 575, 504 (2000). 14 T. Nshno, Y. Heda, K. Okunsh, N. Maeshma, Y. Akutsu, and A. Gendar, Prog. Theo. Phs, 105, 409 (2001). 15 N. Maeshma, Y. Heda, Y. Akutsu, T. Nshno, and K. Okunsh, Phs. Rev. E 64, (2001). 16 A. Gendar, and T. Nshno, Phs. Rev. B 71, (2005). 17 A. Garca-Saez, and J. I. Latorre, arxv: Z. C. Gu, M. Levn, and X. G. Wen, unpublshed. 19 R. Orus, Phs. Rev. B 85, (2012). 20 L. de Lathauwer, B. de Moor, and J. Vandewalle, SIAM J. Matr Anal. Appl, 21, 1253 (2000). 21 L. de Lathauwer, B. de Moor, and J. Vandewalle, SIAM J. Matr Anal. Appl, 21, 1324 (2000). 22 V. de. Slva, and L. H. Lm, SIAM. J. Matr Anal. Appl. 30, 1084 (2008). 23 D. J. Luo, C. Dng, and H. Huang, arxv: ; G. Bergqvst, E. G. Larsson, IEEE Sgnal Proc. Mag. 27, 151 (2010). 24 S. R. Whte, Phs. Rev. Lett. 69, 2863 (1992). 25 In general, the bond denst matr ma not be hermtan. In ths case one should dagonalze ths matr usng the canoncal transformaton, rather than the untar transformaton. 26 L. Onsager, Phs. Rev. 65, 117 (1944). 27 X. M. Feng, and H. W. J. Blote, Phs. Rev. E 81, (2010). 28 H. Arsue, and T. Fujwara, Phs. Rev. E 67, (2003). 29 M. Hasenbusch, Int. J. of Mod. Phs. C 12, 911 (2001). 30 P. Butera, and M. Com, Phs. Rev. B 60, 6749 (1999). 31 A. J. Lu, and M. E. Fsher, Phsca A 156, 35 (1989). 32 M. F. Skes, D. L. Hunter, D. S. McKenze, and B. R. Heap, J. Phs. A 5, 667 (1972). 33 M. Hasenbusch, and A. Pnn, J. Phs. A 31, 6157 (1998). 34 S. J. K. Jensen, and O. G. Mourtsen, J. Phs. A 15, 2631 (1982). 35 A. L. Talapov, and H. W. J. Blote, J. Phs. A: Math. Gen. 29, 5727 (1996). 36 M. Campostrn, A. Pelssetto, P. Ross, and E. Vcar, Phs. Rev. E 65, (2002). 37 M. Hasenbusch, Phs. Rev. B 82, (2010). 38 Y. J. Deng, and H. W. J. Blote, Phs. Rev. E 68, (2003). 39 R. Gupta, and P. Tamao, Int. J. Mod. Phs. C 7 305(1996). 40 P. Butera, and M. Com, Phs. Rev. B 62, (2000). 41 S. G. Chung, Phs. Lett. A 359, 707 (2006). 42 H.W. J. Blote and Y. Deng, Phs. Rev. E 66, (2002). 43 H. X. He, C. J. Hamer, and J. Otmaa, J. Phs. A 23, 1775 (1990); J. Otmaa, C. J. Hamer, and W. H. Zheng, J. Phs. A 24, 2863 (1991). 44 R. J. Bursll, T. Xang, G. A. Gehrng, J. Phs.: Condensed Matter 8, L583 (1996). 45 X. Wang, T. Xang, Phs. Rev. B 56, 5061 (1997). cond-mat/

11 11 46 S.J. Ran, W. L, X. X, Z. Zhang, G. Su, arxv: , unpublshed. 47 R. Orus, G. Vdal, Phs. Rev. B (2009).