Infinite projected entangled pair states (iPEPS), the tensor network ansatz for two-dimensional systems in the thermodynamic limit, already provide excellent results on ground-state quantities using either imaginary-time evolution or variational optimisation. Here, we show (i) the feasibility of real-time evolution in iPEPS to simulate the dynamics of an infinite system after a global quench and (ii) the application of disorder-averaging to obtain translationally invariant systems in the presence of disorder. To illustrate the approach, we study the short-time dynamics of the square lattice Heisenberg model in the presence of a bi-valued disorder field.
Infinite projected entangled pair states (iPEPS), the tensor network ansatz for two-dimensional systems in the thermodynamic limit, already provide excellent results on ground-state quantities using either imaginary-time evolution or variational optimisation. Here, we show (i) the feasibility of real-time evolution in iPEPS to simulate the dynamics of an infinite system after a global quench and (ii) the application of disorder-averaging to obtain translationally invariant systems in the presence of disorder. To illustrate the approach, we study the short-time dynamics of the square lattice Heisenberg model in the presence of a bi-valued disorder field.
A typical quantum state obeying the area law for entanglement on an infinite two-dimensional (2D) lattice can be represented by a tensor network Ansatz, known as an infinite projected entangled pair state (iPEPS), with a finite bond dimension D. Its real-imaginary time evolution can be split into small time steps. An application of a time step generates a new iPEPS with a bond dimension k times the original one. The new iPEPS does not make optimal use of its enlarged bond dimension kD; hence, in principle, it can be represented accurately by a more compact Ansatz, preferably with the original D. In this work we show how the more compact iPEPS can be optimized variationally to maximize its overlap with the new iPEPS. To compute the overlap we use the corner-transfer-matrix renormalization group. By simulating sudden quench of the transverse field in the 2D quantum Ising model with the proposed algorithm, we provide a proof of principle that real-time evolution can be simulated with iPEPS. A similar proof is provided with the same model for imaginary-time evolution of purification of its thermal states.
An infinite projected entangled-pair state (iPEPS) is a powerful variational tensor network ansatz for two-dimensional ground states in the thermodynamic limit, and can be seen as a natural generalization of a matrix-product state to two dimensions. One of the main challenges in iPEPS simulations is the optimization of the tensors, i.e., finding the optimal variational parameters, in order to have the best representation of the ground state of a given Hamiltonian. The author presents a variational optimization scheme, in which the energy is minimized in an iterative way by sweeping over all the tensors in the ansatz, in a similar spirit as done in the density-matrix renormalization group method. Benchmark results for challenging problems are presented that 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.
Infinite projected entangled pair states (iPEPS) provide a convenient variational description of infinite, translationally-invariant two-dimensional quantum states. However, the simulation of local excitations is not directly possible due to the translationally-invariant ansatz. Furthermore, as iPEPS are either identical or orthogonal, expectation values between different states as required during the evaluation of non-equal-time correlators are ill-defined. Here, we show that by introducing auxiliary states on each site, it becomes possible to simulate both local excitations and evaluate non-equal-time correlators in an iPEPS setting under real-time evolution. We showcase the method by simulating the t-J t − J model after a single hole has been placed in the half-filled antiferromagnetic background and evaluating both return probabilities and spin correlation functions, as accessible in quantum gas microscopes.
An infinite projected entangled pair state (iPEPS) is a tensor network ansatz to represent a quantum state on an infinite 2D lattice whose accuracy is controlled by the bond dimension D. Its real, Lindbladian, or imaginary time evolution can be split into small time steps. Every time step generates a new iPEPS with an enlarged bond dimension ${D}^{\ensuremath{'}}>D$, which is approximated by an iPEPS with the original D. In P. Czarnik and J. Dziarmaga, Phys. Rev. B 98, 045110 (2018), an algorithm was introduced to optimize the approximate iPEPS by maximizing directly its fidelity to the one with the enlarged bond dimension ${D}^{\ensuremath{'}}$. In this paper, we implement a more efficient optimization employing a local estimator of the fidelity. For imaginary time evolution of a thermal state's purification, we also consider using unitary disentangling gates acting on ancillas to reduce the required D. We test the algorithm simulating Lindbladian evolution and unitary evolution after a sudden quench of transverse field hx in the 2D quantum Ising model. Furthermore, we simulate thermal states of this model and estimate the critical temperature with good accuracy: 0.1 for hx=2.5 and 0.5 for the more challenging case of hx=2.9 close to the quantum critical point at hx=3.04438(2).
The infinite projected entangled pair states (iPEPS) algorithm [J. Jordan et al., Phys. Rev. Lett. 101, 250602 (2008)] has become a useful tool in the calculation of ground-state properties of two-dimensional quantum lattice systems in the thermodynamic limit. Despite its many successful implementations, the method has some limitations in its present formulation which hinder its application to some highly entangled systems. The purpose of this paper is to unravel some of these issues, in turn enhancing the stability and efficiency of iPEPS methods. For this, we first introduce the fast full update scheme, where effective environment and iPEPS tensors are both simultaneously updated (or evolved) throughout time. As we shall show, this implies two crucial advantages: (i) dramatic computational savings and (ii) improved overall stability. In addition, we extend the application of the local gauge fixing, successfully implemented for finite-size PEPS [M. Lubasch et al., Phys. Rev. B 90, 064425 (2014)], to the iPEPS algorithm. We see that the gauge fixing not only further improves the stability of the method but also accelerates the convergence of the alternating least-squares sweeping in the (either ``full'' or ``fast full'') tensor update scheme. The improvement in terms of computational cost and stability of the resulting ``improved'' iPEPS algorithm is benchmarked by studying the ground-state properties of the quantum Heisenberg and transverse-field Ising models on an infinite square lattice.
We investigate the imaginary-time relaxation critical dynamics of the two-dimensional transverse-field Ising model using infinite projected entangled pair states (iPEPS) with the full-update strategy. Simulating directly in the thermodynamic limit, we explore the relaxation process near the critical point with two types of initial states: a fully polarized state and a product state with a small magnetization. For the fully polarized state, the magnetization shows a power law scaling M∝τ−β/(νz) in the imaginary-time evolution, from which both the critical point and critical exponent can be determined with high accuracy. For the nearly paramagnetic state, the relaxation process exhibits a behavior of M∝τθ with θ=0.1958 being the critical initial-slip exponent, which is in good agreement with that obtained from the dynamic scaling of the self-correlation in quantum Monte Carlo method. These universal features emerge well before the system converges to the ground state, demonstrating the efficiency of imaginary-time evolution for probing quantum criticality. Our results demonstrate that iPEPS can serve as a robust and scalable method for studying dynamical critical phenomena in two-dimensional quantum many-body systems.