No Arabic abstract
An infinite projected entangled-pair state (iPEPS) is a variational tensor network ansatz for 2D wave functions in the thermodynamic limit where the accuracy can be systematically controlled by the bond dimension $D$. We show that for the doped Hubbard model in the strongly correlated regime ($U/t=8$, $n=0.875$) iPEPS yields lower variational energies than state-of-the-art variational methods in the large 2D limit, which demonstrates the competitiveness of the method. In order to obtain an accurate estimate of the energy in the exact infinite $D$ limit we introduce and test an extrapolation technique based on a truncation error computed in the iPEPS imaginary time evolution algorithm. The extrapolated energies are compared with accurate quantum Monte Carlo results at half filling and with various other methods in the doped, strongly correlated regime.
The excitation ansatz for tensor networks is a powerful tool for simulating the low-lying quasiparticle excitations above ground states of strongly correlated quantum many-body systems. Recently, the two-dimensional tensor network class of infinite entangled pair states gained new ground state optimization methods based on automatic differentiation, which are at the same time highly accurate and simple to implement. Naturally, the question arises whether these new ideas can also be used to optimize the excitation ansatz, which has recently been implemented in two dimensions as well. In this paper, we describe a straightforward way to reimplement the framework for excitations using automatic differentiation, and demonstrate its performance for the Hubbard model at half filling.
We present a scheme to perform an iterative variational optimization with infinite projected entangled-pair states (iPEPS), a tensor network ansatz for a two-dimensional wave function in the thermodynamic limit, to compute the ground state of a local Hamiltonian. The method is based on a systematic summation of Hamiltonian contributions using the corner transfer-matrix method. Benchmark results for challenging problems are presented, including the 2D Heisenberg model, the Shastry-Sutherland model, and the t-J model, which show that the variational scheme yields considerably more accurate results than the previously best imaginary time evolution algorithm, with a similar computational cost and with a faster convergence towards the ground state.
We analyze a criterion which guarantees that the ground states of certain many body systems are stable under perturbations. Specifically, we consider PEPS, which are believed to provide an efficient description, based on local tensors, for the low energy physics arising from local interactions. In order to assess stability in the framework of PEPS, one thus needs to understand how physically allowed perturbations of the local tensor affect the properties of the global state. In this paper, we show that a restricted version of the Local Topological Quantum Order (LTQO) condition provides a checkable criterion which allows to assess the stability of local properties of PEPS under physical perturbations. We moreover show that LTQO itself is stable under perturbations which preserve the spectral gap, leading to nontrivial examples of PEPS which possess LTQO and are thus stable under arbitrary perturbations.
Variational Monte Carlo studies employing projected entangled-pair states (PEPS) have recently shown that they can provide answers on long-standing questions such as the nature of the phases in the two-dimensional $J_1 - J_2$ model. The sampling in these Monte Carlo algorithms is typically performed with Markov Chain Monte Carlo algorithms employing local update rules, which often suffer from long autocorrelation times and interdependent samples. We propose a sampling algorithm that generates independent samples from a PEPS, bypassing all problems related to finite autocorrelation times. This algorithm is a generalization of an existing direct sampling algorithm for unitary tensor networks. We introduce an auxiliary probability distribution from which independent samples can be drawn, and combine it with importance sampling in order to evaluate expectation values accurately. We benchmark our algorithm on the classical Ising model and on variational optimization of two-dimensional quantum spin models.
We present an extension of a framework for simulating single quasiparticle or collective excitations on top of strongly correlated quantum many-body ground states using infinite projected entangled pair states, a tensor network ansatz for two-dimensional wave functions in the thermodynamic limit. Our approach performs a systematic summation of locally perturbed states in order to obtain excited eigenstates localized in momentum space, using the corner transfer matrix method, and generalizes the framework to arbitrary unit cell sizes, the implementation of global Abelian symmetries and fermionic systems. Results for several test cases are presented, including the transverse Ising model, the spin-$frac{1}{2}$ Heisenberg model and a free fermionic model, to demonstrate the capability of the method to accurately capture dispersions. We also provide insight into the nature of excitations at the $k=(pi,0)$ point of the Heisenberg model.