Search

We describe our implementation of fermionic tensor network contraction on arbitrary lattices within both a globally ordered and locally ordered formalism. We provide a pedagogical description of these two conventions as implemented for the quimb library. Using hyperoptimized approximate contraction strategies, we present benchmark fermionic projected entangled pair states simulations of finite Hubbard models defined on the three-dimensional diamond lattice and random regular graphs., Comment: 9 pages, 9 figures

We introduce a change of perspective on tensor network states that is defined by the computational graph of the contraction of an amplitude. The resulting class of states, which we refer to as tensor network functions, inherit the conceptual advantages of tensor network states while removing computational restrictions arising from the need to converge approximate contractions. We use tensor network functions to compute strict variational estimates of the energy on loopy graphs, analyze their expressive power for ground-states, show that we can capture aspects of volume law time evolution, and provide a mapping of general feed-forward neural nets onto efficient tensor network functions. Our work expands the realm of computable tensor networks to ones where accurate contraction methods are not available, and opens up new avenues to use tensor networks., Comment: 4+9 pages, 13 figures

We use the finite projected entangled pair state (PEPS) method to investigate the global phase diagram of the spin-1/2 square-lattice $J_1$-$J_2$-$J_3$ antiferromagnetic (AFM) Heisenberg model. The ground state phase diagram is established with a rich variety of phases: AFM, gapless quantum spin liquid, valence-bond solid (VBS), stripe, and incommensurate spiral phases. The nature of the VBS region is revealed, containing a plaquette VBS and a mixed columnar-plaquette VBS, with the emergence of short-range incommensurate spin correlations in some region. The long-range incommensurate magnetic phase is also explicitly characterized as a planar spiral with incommensurate spatial periodicities. Most interestingly, there exists several multicritical points connecting different phases. These findings elucidate the true nature of the long-standing square-lattice $J_1$-$J_2$-$J_3$ antiferromagnet at zero-temperature. Our results also pave the way to accurately simulate complex two-dimensional quantum systems that may host nonuniform features by means of finite PEPS., Comment: 13 pages; 19 figures

Motivated by the novel phenomena observed in the layered material $\rm SrCu_2(BO_3)_2$, the Shastry-Sutherland model (SSM) has been extensively studied as the minimal model for $\rm SrCu_2(BO_3)_2$. However, the nature of its quantum phase transition from the plaquette valence-bond solid (PVBS) to antiferromagnetic (AFM) phase is under fierce debate, posing a challenge to understand the underlying quantum criticality. Via the state-of-the-art tensor network simulations, we study the ground state of the SSM on large-scale size up to $20 \times 20$ sites. We identify the continuous transition nature accompanied by an emergent O(4) symmetry between the PVBS and AFM phase, which strongly suggests a deconfined quantum critical point (DQCP). Furthermore, we map out the phase diagram of an extended SSM that can be continuously tuned to the SSM, which demonstrates the same DQCP phenomena along a whole critical line. Our results indicate a compelling scenario for understanding the origin of the proposed proximate DQCP in recent experiments of $\rm SrCu_2(BO_3)_2$., Comment: 5+6 pages; 4+5 figures; 3 tables

The emergence of exotic quantum phenomena in frustrated magnets is rapidly driving the development of quantum many-body physics, raising fundamental questions on the nature of quantum phase transitions. Here we unveil the behaviour of emergent symmetry involving two extraordinarily representative phenomena, i.e., the deconfined quantum critical point (DQCP) and the quantum spin liquid (QSL) state. Via large-scale tensor network simulations, we study a spatially anisotropic spin-1/2 square-lattice frustrated antiferromagnetic (AFM) model, namely the $J_{1x}$-$J_{1y}$-$J_2$ model, which contains anisotropic nearest-neighbor couplings $J_{1x}$, $J_{1y}$ and the next nearest neighbor coupling $J_2$. For small $J_{1y}/J_{1x}$, by tuning $J_2$, a direct continuous transition between the AFM and valence bond solid phase is observed.(Of course, the possibility of weakly first order transition can not be fully excluded.) With growing $J_{1y}/J_{1x}$, a gapless QSL phase gradually emerges between the AFM and VBS phases. We observe an emergent O(4) symmetry along the AFM--VBS transition line, which is consistent with the prediction of DQCP theory. Most surprisingly, we find that such an emergent O(4) symmetry holds for the whole QSL--VBS transition line as well. These findings reveal the intrinsic relationship between the QSL and DQCP from categorical symmetry point of view, and strongly constrain the quantum field theory description of the QSL phase. The phase diagram and critical exponents presented in this paper are of direct relevance to future experiments on frustrated magnets and cold atom systems., Comment: 5+7 pages, 4+11 figures

A quantum spin liquid (QSL) is a novel phase of matter with long-range entanglement where localized spins are highly correlated with the vanishing of magnetic order. Such exotic quantum states provide the opportunities to develop new theoretical frameworks for many-body physics and have the potential application in realizing robust quantum computations. Here we show that a gapless QSL can naturally emerge from a deconfined quantum critical point (DQCP), which is originally proposed to describe Landau forbidden continuous phase transition between antiferromagnetic (AFM) and valence-bond solid (VBS) phases. Via large-scale tensor network simulations of a square-lattice spin-1/2 frustrated Heisenberg model, both QSL state and DQCP-type AFM-VBS transition are observed. With tuning coupling constants, the AFM-VBS transition vanishes and instead, a gapless QSL phase gradually develops in between. Remarkably, along the phase boundaries of AFM-QSL and QSL-VBS transitions, we always observe the same correlation length exponents $\nu\approx 1.0$, which is intrinsically different from the one of the DQCP-type transition, indicating new types of universality classes. Our results explicitly demonstrate a new scenario for understanding the emergence of gapless QSL from an underlying DQCP. The discovered QSL phase survives in a large region of tuning parameters and we expect its experimental realizations in solid state materials or quantum simulators.

The nature of the zero-temperature phase diagram of the spin-$1/2$ $J_1$-$J_2$ Heisenberg model on a square lattice has been debated in the past three decades, which may hold the key to understand high temperature superconductivity. By using the state-of-the-art tensor network method, specifically, the finite projected entangled pair state (PEPS) algorithm, to simulate the global phase diagram the $J_1$-$J_2$ Heisenberg model up to $24\times 24$ sites, we provide very solid evidences to show that the nature of the intermediate nonmagnetic phase is a gapless quantum spin liquid (QSL), whose spin-spin and dimer-dimer correlations both decay with a power law behavior. There also exists a valence-bond solid (VBS) phase in a very narrow region $0.56\lesssim J_2/J_1\leq0.61$ before the system enters the well known collinear antiferromagnetic phase. The physical nature of the discovered gapless QSL and potential experimental implications are also addressed. We stress that we make the first detailed comparison between the results of PEPS and the well-established density matrix renormalization group (DMRG) method through one-to-one direct benchmark for small system sizes, and thus give rise to a very solid PEPS calculation beyond DMRG. Our numerical evidences explicitly demonstrate the huge power of PEPS for precisely capturing long-range physcis for highly frustrated systems, and also demonstrate the finite PEPS method is a very powerful approach to study strongly corrleated quantum many-body problems.

In a recent comments [arXiv:1909.12788 (2019)], Zhao et al. argue that the definition of dimer orders used in our paper [Phys. Rev. B 98, 241109 (2018)] may not rule out valence bond solid (VBS) orders of $J_1$-$J_2$ model on the open boundary conditions (OBC). In this reply, we show that their argument does not apply to our case.

Based on the scheme of variational Monte Carlo sampling, we develop an accurate and efficient two-dimensional tensor-network algorithm to simulate quantum lattice models. We find that Monte Carlo sampling shows huge advantages in dealing with finite projected entangled pair states, which allows significantly enlarged system size and improves the accuracy of tensor network simulation. We demonstrate our method on the square-lattice antiferromagnetic Heisenberg model up to $32 \times 32$ sites, as well as a highly frustrated $J_1-J_2$ model up to $24\times 24$ sites. The results, including ground state energy and spin correlations, are in excellent agreement with those of the available quantum Monte Carlo or density matrix renormalization group methods. Therefore, our method substantially advances the calculation of 2D tensor networks for finite systems, and potentially opens a new door towards resolving many challenging strongly correlated quantum many-body problems.

Recently, there has been significant progress in solving quantum many-particle problem via machine learning based on the restricted Boltzmann machine. However, it is still highly challenging to solve frustrated models via machine learning, which has not been demonstrated so far. In this work, we design a brand new convolutional neural network (CNN) to solve such quantum many-particle problems. We demonstrate, for the first time, of solving the highly frustrated spin-1/2 J$_1$-J$_2$ antiferromagnetic Heisenberg model on square lattices via CNN. The energy per site achieved by the CNN is even better than previous string-bond-state calculations. Our work therefore opens up a new routine to solve challenging frustrated quantum many-particle problems using machine learning.

The spin-1/2 $J_1$-$J_2$ Heisenberg model on square lattices are investigated via the finite projected entangled pair states (PEPS) method. Using the recently developed gradient optimization method combining with Monte Carlo sampling techniques, we are able to obtain the ground states energies that are competitive to the best results. The calculations show that there is no N\'eel order, dimer order and plaquette order in the region of 0.42 $\lesssim J_2/J_1\lesssim$ 0.6, suggesting a single spin liquid phase in the intermediate region. Furthermore, the calculated staggered spin, dimer and plaquette correlation functions all have power law decay behaviours, which provide strong evidences that the intermediate nonmagnetic phase is a single gapless spin liquid state.

Recently, the tensor network states (TNS) methods have proven to be very powerful tools to investigate the strongly correlated many-particle physics in one and two dimensions. The implementation of TNS methods depends heavily on the operations of tensors, including contraction, permutation, reshaping tensors, SVD and son on. Unfortunately, the most popular computer languages for scientific computation, such as Fortran and C/C++ do not have a standard library for such operations, and therefore make the coding of TNS very tedious. We develop a Fortran2003 package that includes all kinds of basic tensor operations designed for TNS. It is user-friendly and flexible for different forms of TNS, and therefore greatly simplifies the coding work for the TNS methods.

The projected entangled pair states (PEPS) methods have been proved to be powerful tools to solve the strongly correlated quantum many-body problems in two-dimension. However, due to the high computational scaling with the virtual bond dimension $D$, in a practical application PEPS are often limited to rather small bond dimensions, which may not be large enough for some highly entangled systems, for instance, the frustrated systems. The optimization of the ground state using time evolution method with simple update scheme may go to a larger bond dimension. However, the accuracy of the rough approximation to the environment of the local tensors is questionable. Here, we demonstrate that combining the time evolution method with simple update, Monte Carlo sampling techniques and gradient optimization will offer an efficient method to calculate the PEPS ground state. By taking the advantages of massive parallel computing, we can study the quantum systems with larger bond dimensions up to $D$=10 without resorting to any symmetry. Benchmark tests of the method on the $J_1$-$J_2$ model give impressive accuracy compared with exact results.