Ground-State Properties of Quantum Many-Body Systems: Entangled-Plaquette States and Variational Monte Carlo
Abstract
We propose a new ansatz for the ground-state wave function of quantum many-body systems on a lattice. The key idea is to cover the lattice with plaquettes and obtain a state whose configurational weights can be optimized by means of a Variational Monte Carlo algorithm. Such a scheme applies to any dimension, without any “sign” instability. We show results for various two dimensional spin models (including frustrated ones). A detailed comparison with available exact results, as well as with variational methods based on different ansatzs is offered. In particular, our numerical estimates are in quite good agreement with exact ones for unfrustrated systems, and compare favorably to other methods for frustrated ones.
pacs:
02.70.Ss, 05.50+qI Introduction
The study of the ground-state (GS) properties of quantum many-body systems is one of the most challenging tasks of theoretical physics. Exact results can be obtained numerically only for systems of relatively small size (i.e., few particles); this limitation is particularly severe, e.g., when studying phase transitions, wherein the emergence of long-range order can only be established by carrying out an extrapolation of the estimates to the thermodynamic limit. Two main numerical techniques, namely Density Matrix Renormalization Group (DMRG)dmrg () and Quantum Monte Carlo (QMC)qmc () have successfully been employed to investigate quantum spin models on a lattice. These computational approaches, however, find their optimal applicability under different specific constraints. DMRG yields extremely accurate results in one dimension (1D) even for very large systems, but fails in describing the properties of the quantum GS in higher dimension due to the unfavorable scaling with the system size of the computational resources needed.liang () QMC, on the other hand, is the method of choice for quantum systems obeying Bose statistics in any dimension, but suffers from the notorious “sign problem” for Fermi systems.
Generalizations of the variational family of Matrix-Product States (MPS) underlying the DMRG have recently been investigated to go beyond the discussed limitations of DMRG itself and QMC. Specifically the most natural extension of MPS is given by Projected-Entangled Pair States (PEPS)peps () which efficiently approximate ground states of local Hamiltonianshas () and have been used to simulate, otherwise intractable, 2D quantum systems.peps1 (); orusetal () Other variational families of states have also been proposed and tested in 2D.vidalmera (); Levin (); Weng (); LevinWen (); cincio (); vidal2 (); vidl () Despite the promising results obtained so far, it seems very hard to use those methods in 3D [or even for systems with periodic boundary conditions (PBC)] due to the unfavorable scaling of the required computer resources.
A new possibility has recently emerged to combine the main advantages of DMRG and Monte Carlo in order to build new algorithms.SBS (); SV () Some of themSBS (); Snew () can be used in 2D and it seems that the one based on String-Bond StatesSBS () may be used even for some 3D systems.Alessandro ()
In this paper we introduce a new numerical technique that combines the strengths of QMC with an extension of PEPS to simulate lattices models, overcoming some of their limitations. We also test this technique with non-trivial models, and compare it with other techniques, including PEPS. Specifically, we propose a new class of states called Entangled-Plaquette States (EPS) which allow an accurate characterization of the GS of quantum spin systems by means of a simple Variational Monte Carlo (VMC) algorithm (see Ref. VMC, for a general review).
The basic idea underlying EPS can be schematically described as follows: Assume to cut a lattice in several sub-blocks (e.g. small PEPS) and extract for each sub-block the GS wave function. The wave function of the original system, expressed as the product of the sub-block wave functions yields reasonable (up to corrections which scale with the sub-block boundaries) estimates of GS energy and short-range (of the order of the sub-block size) correlations. These estimates dramatically improve if the sub-block size is increased and, more importantly, if overlapping sub-blocks (i.e., entangled plaquettes) are employed. The latter is the crucial point which allows, accounting for its correlated nature, a description of the quantum GS much more accurate than that obtainable with a simple non-overlapping-plaquette product state (i.e., in a mean-field fashion).o1 () A GS wave function whose weights are the product of variational parameters in one to one correspondence to the spin configuration of each entangled plaquette, appears, consequently, the natural choice for a variational ansatz. Therefore, our numerical approach is based on a family of states, namely EPS, which share many analogies with PEPS, and takes advantage of Monte Carlo sampling to estimate physical observables of interest. It can be applied to systems of any spatial dimensionality and, as a pure variational method (i.e., not involving imaginary time projection), is sign problem free.
We test our numerical protocol on a variety of quantum spin models on a square lattice comprising as many as sites. Our energy estimates are in excellent agreement with those (“exact” in practice) obtained by QMC for a system of lattice hard core bosons. In the presence of nearest-neighbor repulsion at the Heisenberg point, we find an extrapolated (to infinite lattice size) value of the energy per site which differs from the QMC resultsand () by less than 10, and is more accurate than VMC estimates obtained with a Jastrow wave function.cep90 () In the case of a frustrated antiferromagnet (i.e., the so-called model), for which a sign problem exists in QMC, our energy estimates (whose error relative to the exact ones is less than 10) compare favorably with those obtained with PEPS, or fixed-node Green Function Monte Carlo (GFMC).massimo ()
Ii Methodology
Consider a collection of spins arranged on a square lattice with (without loss of generality, here and in the following we will refer to this specific case). Provided a trial state where and , the energy expectation value on the given state is:
(1) |
where
(2) |
and (real weights are assumed for simplicity). According to the variational principle, is an upper bound of the GS energy which can be evaluated by minimizing Eq. 1 with respect to the weights. At this point we have to make an ansatz for the wave function: imagine to cover the lattice with plaquettes (say one of dimension for each site) and assign a coefficient to all the possible spin configurations of any single plaquette. Given a global spin configuration , its weight can be expressed as follows:
(3) |
where depends only on the spin state of the sites belonging to the plaquette. With this choice the analytic expression of the derivative of Eq. 1 with respect to is:
(4) |
The multidimensional summations in Eq. 1 and 4 can exactly be evaluated only for small . For large systems () , one has to employ the Monte Carlo method. Specifically, the energy, as well as its derivatives, can be estimated from the same sample. Moreover, the only quantity depending on the plaquette coefficient with respect to which the derivative is taken is . By using Eq. 3 it turns out that is simply equal to (easy and fast to compute). Similarly to other works,SBS (); SV (), the steps of the basic variational algorithm adopted to calculate the GS energy are: i) Start from a randomly chosen initial configuration; ii) generate a large set of new configurations by flipping one or more spins via the Metropolis algorithm;metro () iii) evaluate the energy and its gradient vector; iv) update all the of a small step against the gradient direction; v) iterate from ii) until convergence of the energy is reached. It is worth mentioning that expectation values of physical observables other than the energy can be evaluated according to Eq. 1 and 2 when is replaced by a generic operator .
For a single-spin flip the acceptance probability is given by:
(5) |
where the flip is proposed for the spin and the products run over all the plaquettes which include such a spin.
In a typical calculation, we start with plaquettes; once the energy has converged, the size of the plaquettes is increased to improve the estimate. The number of coefficients which need be stored in memory for a spin-1/2 system is , and can be reduced taking into account problem-dependent symmetries. For each optimization step, a few thousands updates are necessary to get a rough estimate of the energy and efficiently move the coefficients along the gradient. At the later stages of the simulation, to reach the optimal energy value, an important role is played by the gradient step which has to be carefully tuned. An example of how the error of the GS energy relative to the GFMC result decreases as a function of the plaquettes size is illustrated in Fig. 1.
The relative error is already small (less than ) for plaquettes and reaches a value when plaquettes are used. Numerical data refer to a system of lattice bosons (at half filling) which interact via an infinite on-site (hard core) repulsion on a lattice with PBC. Since the total magnetization along is a good quantum number, we performed the calculation in the canonical ensemble (i.e., in the sector). Consequently we chose to update the configuration by flipping pairs of spins and for which . The expression of the acceptance probability for the pair update is a straightforward generalization of Eq. 5.
Iii Results
Estimates of the GS energy per site of a system of lattice hard core bosons, for three different lattice sizes, are presented and compared to GFMC results in Tab. 1. Although the GFMC method is not purely variational (i.e., the GS wave function is projected out from a trial state via imaginary time evolution) our numerical data are in excellent agreement with GFMC ones, even for the largest system considered in this work ().
This work | GFMC | |
---|---|---|
8 | -1.1004(4) | -1.1008(3) |
10 | -1.0988(4) | -1.0998(1) |
20 | -1.0952(7) | -1.0967(4) |
As a second case study we investigate the quantum spin Hamiltonian:
(6) |
where the first (second) summation runs over nearest- (next-nearest-) neighbor sites. For and the above Hamiltonian describes the antiferromagnetic Heisenberg model, which is isomorphic to a system of hard core bosons with nearest-neighbor repulsion of strength =. This model has been extensively studied numerically in the past (see for instance Refs. sand, and cep90, ). GS energies per site for various lattice sizes are presented in Tab 2 and compared to exact calculationsexact () and Stocastic Series Expansion (SSE) data.sand ()
This work | SSE^{1}^{1}1Ref. sand, | Exact^{2}^{2}2Ref. exact, | |
---|---|---|---|
4 | -0.7016(1) | -0.701777(7) | -0.7018 |
6 | -0.6785(2) | -0.678873(4) | -0.6789 |
8 | -0.6724(3) | -0.673487(4) | - |
10 | -0.6699(3) | -0.671549(4) | - |
The energy per site is a monotonic increasing function of the lattice size. The simple formulaman ()
(7) |
where ,beta () gives an extrapolated energy per site , which is in good agreement with the most accurate QMC result:sand () and lower than other extrapolated values obtained by purely variational methods. For example, on the basis of a Jastrow wave function, Trivedi and Ceperley found .cep90 ()
Calculations carried on with open boundary conditions (OBC) yield GS energies lower than those obtained with the general PEPS, or SBS method. For example, we get for while the PEPS result is and the SBS one .
This work | SSE^{3}^{3}3Ref. sand, | |
---|---|---|
4 | 0.1807(4) | 0.17962(1) |
6 | 0.1550(4) | 0.152568(9) |
8 | 0.1432(4) | 0.13760(1) |
10 | 0.1352(5) | 0.12855(2) |
Estimates of the spin-spin correlation function computed at the maximum distance on the lattice according to the formula where are shown in Tab 3. We find an extrapolated value of the staggered magnetization defined by of 0.324(1) which is in reasonable agreement with reported in Ref sand, and with other estimates.man () It has to be mentioned that the discrepancy between our result and the SSE one might be due to a violation of the so-called “area law”.area ()
This work | Exact^{4}^{4}4Ref. exact, | |
---|---|---|
0.0 | -0.6785(2) | -0.6789 |
0.1 | -0.6377(1) | -0.6381 |
0.2 | -0.5985(1) | -0.5990 |
0.3 | -0.5616(1) | -0.5625 |
0.4 | -0.5277(1) | -0.5297 |
0.5 | -0.4985(2) | -0.5038 |
0.6 | -0.4860(2) | -0.4932 |
0.7 | -0.5255(1) | -0.5300 |
0.8 | -0.5843(1) | -0.5865 |
0.9 | -0.6453(1) | -0.6491 |
1.0 | -0.7091(1) | -0.7144 |
This work | PEPS^{5}^{5}5data provided by V. Murg | |
---|---|---|
0.0 | -0.61567(8) | -0.61506 |
0.1 | -0.5865(1) | -0.5845 |
0.2 | -0.5571(1) | -0.5555 |
0.3 | -0.5290(2) | -0.5274 |
0.4 | -0.5028(1) | -0.5016 |
0.5 | -0.4781(1) | -0.4779 |
0.6 | -0.4523(2) | -0.4508 |
0.7 | -0.4553(1) | -0.4541 |
0.8 | -0.4928(2) | -0.4906 |
0.9 | -0.5371(1) | -0.5344 |
1.0 | -0.5800(1) | -0.5792 |
By including the next-nearest-neighbor interaction () in the Hamiltonian the system becomes frustrated and is believed to undergo a phase transition for . This model cannot be simulated by QMC due to the sign problem (arising in turn from the underlying Fermi statistics). Our method, instead, is purely variational and can be applied without the occurrence of any sign instability. With the aim of comparing our estimates with exact results (available only for ), we computed the GS energy, as a function of , of the model on a square lattice. Numerical values are shown in Tab. 4; the error of our estimates relative to the exact result is shown in Fig. 2. This quantity never exceeds 10, moreover it has to be mentioned that the GS energies computed by means of EPS compare favorably to those obtained by GFMC with the fixed-node approximation massimo () and, except in a narrow region of , where the the agreement is however remarkable, to variational results obtained with a BCS-type ansatz.sor1 (); unp ()
Also in this case, for OBC, we test our scheme against the PEPS one. The estimated energies are lower than those computed with PEPS even for a lattices (see Tab. 5), where the PEPS approach performs at its best.
Iv Conclusions and outlook
In this work we have presented Entangled-Plaquette States: an ansatz for the GS wave function of quantum many-body systems which allows the accurate estimate of physical observables by means of Variational Monte Carlo. Our approach not only gives accurate results for unfrustrated systems (where other methods are in principle “exact” ) but, most importantly, applied to systems for which QMC simulations suffer from a sign problem, yields estimates whose accuracy seems to be not obtainable with different techniques. Therefore our method appears as a very promising avenue to investigate the effects of Fermi statistics (e.g., frustration). The extension of our computational approach to 3D can be easily implemented (i.e., taking cubic plaquettes) and several possible improvements (e.g., cover the lattice with plaquettes of different shapes and sizes), are being currently investigated.
Acknowledgments
We acknowledge discussions with F. Verstraete and M. Wolf, and thank V. Murg for providing us with the PEPS data. This work has been supported by the DFG (FOR 635) and the EU project QUEVADIS, and the Canadian NSERC through the grant G121210893.
Note added.— After this paper was accepted a related approach has been proposed by J. Hitesh et al., arxiv:0907.4646.
References
- (1) S. R. White, Phys. Rev. Lett. 69, 2863 (1992).
- (2) D. M. Ceperley and B. J. Alder, Phys. Rev. Lett. 45, 566 (1980).
- (3) S Liang and H. Pang, Phys. Rev. B 49, 9214 (1994).
- (4) F. Verstraete and J. I. Cirac, cond-mat/0406426.
- (5) M. B. Hastings, Phys. Rev. B 76, 035114 (2007).
- (6) V. Murg, F. Verstraete and J. I. Cirac, Phys. Rev. A 75, 033605 (2007); arxiv: 0901.2019.
- (7) J. Jordan, R. Orus, G. Vidal, F. Verstraete and J. I. Cirac, Phys. Rev. Lett. 101, 250602 (2008).
- (8) G. Vidal, Phys. Rev. Lett. 99, 220405 (2007), ibidem 101, 110501 (2008).
- (9) M. Levin and C. P. Nave, Phys. Rev. Lett. 99, 120601 (2007).
- (10) H. C. Jiang, Z. Y. Weng and T. Xiang, Phys. Rev. Lett. 101, 090603 (2008).
- (11) Z.-C. Gu, M. Levin and X.-G. Wen, Phys. Rev. B 78, 205116 (2008).
- (12) L. Cinicio, J. Dziarmaga and M. M. Rams, Phys. Rev. Lett. 100, 240603 (2008).
- (13) G Evenbly and G. Vidal, Phys. Rev. B 79, 144108 (2009); Phys. Rev. Lett. 102, 180406 (2009).
- (14) G Evenbly and G. Vidal, arxiv: 0904.3383.
- (15) N. Schuch, M. M. Wolf, F. Verstraete and J. I. Cirac, Phys. Rev. Lett. 100, 040501 (2008).
- (16) A. W. Sandvik and G. Vidal, Phys. Rev. Lett. 99, 220602 (2007).
- (17) L. Wang, Y.-J. Kao and A. W. Sandvik, arxiv:0901.0214v1.
- (18) A. Sfondrini et al., in preparation.
- (19) W. M. C. Foulkes, L. Mitas, R. J. Needs and G. Rajagopal, Rev. Mod. Phys. 73, 33 (2001).
- (20) L. Isaev, G. Ortiz and J. Dukelsky Phys. Rev. B 79, 024409 (2009); arxiv:0903.1630.
- (21) A. W. Sandvik, Phys Rev B 56, 11678 (1996).
- (22) N. Trivedi and D. M. Ceperley, Phys. Rev. B 41, 4552 (1990).
- (23) M. Boninsegni, Phys. Lett. A 216, 313 (1996).
- (24) N. Metropolis, A. Rosenbluth, M. N. Rosenbluth, A. H. Teller and E. Teller, J. Chem. Phys 21, 1087 (1953).
- (25) H. J. Schulz, T. A. L. Ziman and D. Poilblanc, J. Physique I 6, 75 (1996).
- (26) E. Manousakis Rev. Mod. Phys. 63, 1 (1991).
- (27) P. Hasenfratz and F. Niedermayer, Z. Phys. B 92, 91 (1993).
- (28) R. W. Chhajlany, P. Tomczak and A. Wòjcik, Phys. Rev Lett 99, 167204 (2007).
- (29) L. Capriotti, F. Becca, A. Parola and S. Sorella, Phys. Rev Lett 87, 097201 (2001).
- (30) F. Becca, L. Capriotti, A. Parola and S. Sorella, arxiv:0905.4854