No Arabic abstract
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.
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 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.
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.
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.
The recently developed stochastic gradient method combined with Monte Carlo sampling techniques [PRB {bf 95}, 195154 (2017)] offers a low scaling and accurate method to optimize the projected entangled pair states (PEPS). We extended this method to the fermionic PEPS (fPEPS). To simplify the implementation, we introduce a fermi arrow notation to specify the order of the fermion operators in the virtual entangled EPR pairs. By defining some local operation rules associated with the fermi arrows, one can implement fPEPS algorithms very similar to that of standard PEPS. We benchmark the method for the interacting spinless fermion models, and the t-J models. The numerical calculations show that the gradient optimization greatly improves the results of simple update method. Furthermore, much larger virtual bond dimensions ($D$) and truncation dimensions ($D_c$) than those of boson and spin systems are necessary to converge the results. The method therefore offer a powerful tool to simulate fermion systems because it has much lower scaling than the direct contraction methods.