In this work, a simple and fundamental numeric scheme dubbed as ab initio optimization principle (AOP) is proposed for the ground states of translational invariant strongly correlated quantum lattice models. The idea is to transform a nondeterministic-polynomial-hard ground-state simulation with infinite degrees of freedom into a single optimization problem of a local function with finite number of physical and ancillary degrees of freedom. This work contributes mainly in the following aspects: (1) AOP provides a simple and efficient scheme to simulate the ground state by solving a local optimization problem. Its solution contains two kinds of boundary states, one of which play the role of the entanglement bath that mimics the interactions between a supercell and the infinite environment, and the other gives the ground state in a tensor network (TN) form. (2) In the sense of TN, a novel decomposition named as tensor ring decomposition (TRD) is proposed to implement AOP. Instead of following the contraction-truncation scheme used by many existing TN-based algorithms, TRD solves the contraction of a uniform TN in an opposite way by encoding the contraction in a set of self-consistent equations that automatically reconstruct the whole TN, making the simulation simple and unified; (3) AOP inherits and develops the ideas of different well-established methods, including the density matrix renormalization group (DMRG), infinite time-evolving block decimation (iTEBD), network contractor dynamics, density matrix embedding theory, etc., providing a unified perspective that is previously missing in this fields. (4) AOP as well as TRD give novel implications to existing TN-based algorithms: A modified iTEBD is suggested and the two-dimensional (2D) AOP is argued to be an intrinsic 2D extension of DMRG that is based on infinite projected entangled pair state. This paper is focused on one-dimensional quantum models to present AOP. The benchmark is given on a transverse Ising chain and 2D classical Ising model, showing the remarkable efficiency and accuracy of the AOP.
In this work, a simple and fundamental numeric scheme dubbed as ab initio optimization principle (AOP) is proposed for the ground states of translational invariant strongly correlated quantum lattice models. The idea is to transform a nondeterministic-polynomial-hard ground-state simulation with infinite degrees of freedom into a single optimization problem of a local function with finite number of physical and ancillary degrees of freedom. This work contributes mainly in the following aspects: (1) AOP provides a simple and efficient scheme to simulate the ground state by solving a local optimization problem. Its solution contains two kinds of boundary states, one of which play the role of the entanglement bath that mimics the interactions between a supercell and the infinite environment, and the other gives the ground state in a tensor network (TN) form. (2) In the sense of TN, a novel decomposition named as tensor ring decomposition (TRD) is proposed to implement AOP. Instead of following the contraction-truncation scheme used by many existing TN-based algorithms, TRD solves the contraction of a uniform TN in an opposite way by encoding the contraction in a set of self-consistent equations that automatically reconstruct the whole TN, making the simulation simple and unified; (3) AOP inherits and develops the ideas of different well-established methods, including the density matrix renormalization group (DMRG), infinite time-evolving block decimation (iTEBD), network contractor dynamics, density matrix embedding theory, etc., providing a unified perspective that is previously missing in this fields. (4) AOP as well as TRD give novel implications to existing TN-based algorithms: A modified iTEBD is suggested and the two-dimensional (2D) AOP is argued to be an intrinsic 2D extension of DMRG that is based on infinite projected entangled pair state. This paper is focused on one-dimensional quantum models to present AOP. The benchmark is given on a transverse Ising chain and 2D classical Ising model, showing the remarkable efficiency and accuracy of the AOP.
We develop a strategy for tensor network algorithms that allows to deal very efficiently with lattices of high connectivity. The basic idea is to fine grain the physical degrees of freedom, i.e., decompose them into more fundamental units which, after a suitable coarse graining, provide the original ones. Thanks to this procedure, the original lattice with high connectivity is transformed by an isometry into a simpler structure, which is easier to simulate via usual tensor network methods. In particular this enables the use of standard schemes to contract infinite 2D tensor networks-such as corner transfer matrix renormalization schemes-which are more involved on complex lattice structures. We prove the validity of our approach by numerically computing the ground-state properties of the ferromagnetic spin-1 transverse-field Ising model on the 2D triangular and 3D stacked triangular lattice, as well as of the hardcore and softcore Bose-Hubbard models on the triangular lattice. Our results are benchmarked against those obtained with other techniques, such as perturbative continuous unitary transformations and graph projected entangled pair states, showing excellent agreement and also improved performance in several regimes.
The study of interacting quantum many-body systems is one of the central challenges in modern physics. Interestingly, by understanding the physical properties of such systems from a microscopic description, one can explain many complex quantum phenomena such as the high-temperature superconductivity, fractional quantum Hall effect, etc. However, due to high complexity of these systems, theoretical investigations such as analytical and perturbation theories are limited to studying only a small subset of problems. In the past few decades, numerical simulation has emerged as an indispensable method for studying many phenomena, and has led to immense success and progress in recent years. The main obstacle in addressing quantum many-body systems is that the number of parameters required to represent a state of the system scales exponentially with the system size. However, one hardly needs the full Hilbert space but only a small subspace to describe physical properties of the system. From the point of view in quantum information theory, quantum entanglement is one of the main factors that is responsible for the complexity of studying quantum many-body systems. In particular, the amount of quantum entanglement in such system quantifies the correlation between constituents and therefore it determines how big the subspace is required to represent the system. In this thesis, we are interested in using tensor network states, a new emerging numerical technique to study dynamical and ground state properties of the system. The key ingredient is that a tensor network state which consists of a fairly small number of parameters can represent very well the complete wave function of a system. In one spatial dimension, Matrix Product State (MPS) has become a standard tool, combined with related algorithms such as Time Evolving Block Decimation algorithm (TEBD) to study dynamical properties of a locally-perturbed quantum lattice system. In this thesis, we introduce a technique called ``infinite boundary conditions'' that allows us to use a finite-size MPS to represent the state of an infinite-size system. As a result, one can simulate dynamics of a system for a longer time with a significant computational cost reduction compared to traditional methods. In two spatial dimensions, the tensor network known as Projected Entangled Pair States (PEPS) is emerging as a useful approach. However algorithms for optimizing the matrix elements of these tensors are complex and not very efficient. In particular, the efficiency of the optimization depends critically on the quality of the so-called ``environment'', which is expensive to compute. We have introduced a new method to accelerate the convergence of calculating the PEPS ground state by recycling the environment tensors for many optimization steps, vastly reducing the number of times that the full environment needs to be calculated.
A bang-bang (BB) algorithm prepares the ground state of a lattice quantum many-body Hamiltonian H=H1+H2 by evolving an initial product state alternating between H1 and H2. We optimize the algorithm with tensor networks in one and two dimensions. The optimization has two stages. In stage one, a shallow translationally-invariant circuit is optimized in an infinite lattice. In stage two, the infinite-lattice gate sequence is used as a starting point for a finite lattice where it remains optimal in the bulk. The prepared state requires optimization only at its boundary, within a healing length from lattice edges, and the gate sequence needs to be modified only within the causal cone of the boundary. We test the procedure in the 1D and 2D quantum Ising model near its quantum critical point employing, respectively, the matrix product state (MPS) and the pair-entangled projected state (PEPS). At the boundary already the infinite-lattice sequence turns out to provide a more accurate state than in the bulk.
Tensor networks capture large classes of ground states of phases of quantum matter faithfully and efficiently. Their manipulation and contraction has remained a challenge over the years, however. For most of the history, ground state simulations of two-dimensional quantum lattice systems using (infinite) projected entangled pair states have relied on what is called a time-evolving block decimation. In recent years, multiple proposals for the variational optimization of the quantum state have been put forward, overcoming accuracy and convergence problems of previously known methods. The incorporation of automatic differentiation in tensor networks algorithms has ultimately enabled a new, flexible way for variational simulation of ground states and excited states. In this work we review the state-of-the-art of the variational iPEPS framework, providing a detailed introduction to automatic differentiation, a description of a general foundation into which various two-dimensional lattices can be conveniently incorporated, and demonstrative benchmarking results.
In this paper we explore the practical use of the corner transfer matrix and its higher-dimensional generalization, the corner tensor, to develop tensor network algorithms for the classical simulation of quantum lattice systems of infinite size. This exploration is done mainly in one and two spatial dimensions (1D and 2D). We describe a number of numerical algorithms based on corner matrices and tensors to approximate different ground-state properties of these systems. The proposed methods also make use of matrix product operators and projected entangled pair operators and naturally preserve spatial symmetries of the system such as translation invariance. In order to assess the validity of our algorithms, we provide preliminary benchmarking calculations for the spin-1/2 quantum Ising model in a transverse field in both 1D and 2D. Our methods are a plausible alternative to other well-established tensor network approaches such as iDMRG and iTEBD in 1D, and iPEPS and TERG in 2D. The computational complexity of the proposed algorithms is also considered and, in 2D, important differences are found depending on the chosen simulation scheme. We also discuss further possibilities, such as 3D quantum lattice systems, periodic boundary conditions, and real-time evolution. This discussion leads us to reinterpret the standard iTEBD and iPEPS algorithms in terms of corner transfer matrices and corner tensors. Our paper also offers a perspective on many properties of the corner transfer matrix and its higher-dimensional generalizations in the light of novel tensor network methods.
The research work in this thesis is based on strongly interacting quantum lattice systems. The biggest part of research was conducted using state-of-the-art tensor network simulations. Tensor networks provide efficient and highly accurate representations of quantum states when any simply connected patch of the quantum state is slightly entangled. Matrix Product States (MPS) is a tensor network representation for quantum states which is quasi-exact for one-dimensional systems when entanglement entropy of any bipartition of the state follows ``area-law". Projected Entangled Pair States (PEPS) is the extension of MPS for higher dimensional systems, where entanglement entropy area-law is non-trivial. These formalisms and their relevant ground state optimization techniques, Density Matrix Renormalization Group (DMRG) for MPS and Simple Update (SU) with Tensor Renormalization Group (TRG) for PEPS are thoroughly analysed. Entanglement entropy area-law is usually obeyed by the ground state of local Hamiltonians, while generic highly excited states follow a volume-law (they span a finite part of the Hilbert space). Recently, a class of interacting system was shown to undergo a dynamical phase transition (Many Body Localization) where the entropy of every eigenstate follows an area-law. This transition is achieved when the system is highly disordered and the quantum many body state becomes localized similarly to the free particle Anderson localization. The area-law property makes highly excited eigenstates of Many Body Localized (MBL) systems efficiently represented by the MPS ansatz. We develop a highly optimized algorithm (eDMRG) which goes beyond ground state optimization and successfully target eigenstates in any part of the spectrum. This algorithm is used together with analytical calculations, based on local integrability of MBL systems, to identify the universal behaviour of the entanglement spectrum of highly excited eigenstates in MBL systems. In the second part we study interacting bosons in 2D optical lattices. Strongly interacting bosons are simulated using the Bose-Hubbard (BH) model when the interactions are strictly local and the Extended Bose-Hubbard (EBH) model, when additional dipolar interactions are present. The ground states of BH/EBH Hamiltonians in a hexagonal lattice are studied. The phase spectrum includes various insulating and critical phases which are studied in detail using infinite-PEPS ansatz. Additional results on entanglement entropy scaling with respect to the filling of the lattice are presented. Finally, a strongly interacting Harper-Hofstadter Hamiltonian is realized by combining synthetic fields with the BH model. The homogeneity of the system is then broken by a parabolic trapping potential, similar to the ones used in cold atom experiments. Using time-dependent Gutzwiller ansatz (GA) the expansion dynamics of the cloud in large square lattices are studied. In contrast to the expansion dynamics of the BH model, it is found that the synthetic fields generate a self-trapping effect. Using phenomenology and simulations the dynamics are studied in the hard-core and soft-core limit.