Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:1701.07035 · Physical review. B./Physical review. B

Variational optimization algorithms for uniform matrix product states

Valentin Zauner-Stauber, Laurens Vanderstraeten, Matthew Fishman, Frank Verstraete, Jutho Haegeman

Published Jan 25, 2018• 267 citations
arXivHTMLPDF

Abstract

We combine the density matrix renormalization group (DMRG) with matrix product state tangent space concepts to construct a variational algorithm for finding ground states of one-dimensional quantum lattices in the thermodynamic limit. A careful comparison of this variational uniform matrix product state algorithm (VUMPS) with infinite density matrix renormalization group (IDMRG) and with infinite time evolving block decimation (ITEBD) reveals substantial gains in convergence speed and precision. We also demonstrate that VUMPS works very efficiently for Hamiltonians with long-range interactions and also for the simulation of two-dimensional models on infinite cylinders. The new algorithm can be conveniently implemented as an extension of an already existing DMRG implementation.

DOI
10.1103/physrevb.97.045145

Related papers

Unifying time evolution and optimization with matrix product states

Jutho Haegeman, Christian Lubich, Ivan Valer'evich Oseledets +2 more

We show that the time-dependent variational principle provides a unifying framework for time-evolution methods and optimization methods in the context of matrix product states. In particular, we introduce a new integration scheme for studying time evolution, which can cope with arbitrary Hamiltonians, including those with long-range interactions. Rather than a Suzuki-Trotter splitting of the Hamiltonian, which is the idea behind the adaptive time-dependent density matrix renormalization group method or time-evolving block decimation, our method is based on splitting the projector onto the matrix product state tangent space as it appears in the Dirac-Frenkel time-dependent variational principle. We discuss how the resulting algorithm resembles the density matrix renormalization group (DMRG) algorithm for finding ground states so closely that it can be implemented by changing just a few lines of code and it inherits the same stability and efficiency. In particular, our method is compatible with any Hamiltonian for which ground-state DMRG can be implemented efficiently. In fact, DMRG is obtained as a special case of our scheme for imaginary time evolution with infinite time step.

arXiv:1408.5056•Physical review. B./Physical review. B•Oct 10, 2016•800 citations

Parallel time-dependent variational principle algorithm for matrix product states

Paul Secular, Nikita Gourianov, Michael Lubasch +3 more

Combining the time-dependent variational principle (TDVP) algorithm with the parallelization scheme introduced by Stoudenmire and White for the density matrix renormalization group (DMRG), we present the first parallel matrix product state (MPS) algorithm capable of time evolving one-dimensional (1D) quantum lattice systems with long-range interactions. We benchmark the accuracy and performance of the algorithm by simulating quenches in the long-range Ising and XY models. We show that our code scales well up to 32 processes, with parallel efficiencies as high as $86%$. Finally, we calculate the dynamical correlation function of a 201-site Heisenberg XXX spin chain with $1/{r}^{2}$ interactions, which is challenging to compute sequentially. These results pave the way for the application of tensor networks to increasingly complex many-body systems.

arXiv:1912.06127•Physical review. B./Physical review. B•Jun 5, 2020•40 citations

Symmetry-protected local minima in infinite DMRG

Robert N. C. Pfeifer

The infinite density matrix renormalization group (iDMRG) algorithm is a highly successful numerical algorithm for the study of low-dimensional quantum systems, and is also frequently used to initialize the more popular finite DMRG algorithm. Implementations of both finite and infinite DMRG frequently incorporate support for the protection and exploitation of symmetries of the Hamiltonian. In common with other variational tensor network algorithms, convergence of iDMRG to the ground state is not guaranteed, with the risk that the algorithm may become stuck in a local minimum. In this paper, I demonstrate the existence of a particularly harmful class of physically irrelevant local minima affecting both iDMRG and to a lesser extent also infinite time-evolving block decimation (iTEBD), for which the ground state is compatible with the protected symmetries of the Hamiltonian but cannot be reached using the conventional iDMRG or iTEBD algorithms. I describe a modified iDMRG algorithm which evades these local minima, and which also admits a natural interpretation on topologically ordered systems with a boundary.

arXiv:1505.06266•Physical Review B•Nov 23, 2015•1 citation

Infinite size density matrix renormalization group, revisited

Ian P. McCulloch

I revisit the infinite-size variant of the Density Matrix Renormalization Group (iDMRG) algorithm for obtaining a fixed-point translationally invariant matrix product wavefunction in the context of one-dimensional quantum systems. A crucial ingredient of this algorithm is an efficient transformation for obtaining the matrix elements of the wavefunction as the lattice size is increased, and I introduce here a versatile transformation that is demonstrated to be much more effective than previous versions. The resulting algorithm has a surprisingly close relationship to Vidal's Time Evolving Block Decimation for infinite systems, but allows much faster convergence. Access to a translationally invariant matrix product state allows the calculation of correlation functions based on the transfer matrix, which directly gives the spectrum of all correlation lengths. I also show some advantages of the Matrix Product Operator (MPO) technique for constructing expectation values of higher moments, such as the exact variance $$.

arXiv:0804.2509•Apr 16, 2008•97 citations

Tangent-space methods for uniform matrix product states

Laurens Vanderstraeten, Jutho Haegeman, Frank Verstraete

In these lecture notes we give a technical overview of tangent-space methods for matrix product states in the thermodynamic limit. We introduce the manifold of uniform matrix product states, show how to compute different types of observables, and discuss the concept of a tangent space. We explain how to variationally optimize ground-state approximations, implement real-time evolution and describe elementary excitations for a given model Hamiltonian. Also, we explain how matrix product states approximate fixed points of one-dimensional transfer matrices. We show how all these methods can be translated to the language of continuous matrix product states for one-dimensional field theories. We conclude with some extensions of the tangent-space formalism and with an outlook to new applications.

arXiv:1810.07006•SciPost Physics Lecture Notes•Jan 15, 2019•230 citations

Finding the ground states of symmetric infinite-dimensional Hamiltonians: explicit constrained optimizations of tensor networks

S. N. Saadatmand

Abstract Understanding extreme non-locality in many-body quantum systems can help resolve questions in thermostatistics and laser physics. The existence of symmetry selection rules for Hamiltonians with non-decaying terms on infinite-size lattices can lead to finite energies per site, which deserves attention. Here, we present a tensor network approach to construct the ground states of nontrivial symmetric infinite-dimensional spin Hamiltonians based on constrained optimizations of their infinite matrix product states description, which contains no truncation step, offers a very simple mathematical structure, and other minor advantages at the cost of slightly higher polynomial complexity in comparison to an existing method. More precisely speaking, our proposed algorithm is in part equivalent to the more generic and well-established solvers of infinite density-matrix renormalization-group and variational uniform matrix product states, which are, in principle, capable of accurately representing the ground states of such infinite-range-interacting many-body systems. However, we employ some mathematical simplifications that would allow for efficient brute-force optimizations of tensor-network matrices for the specific cases of highly-symmetric infinite-size infinite-range models. As a toy-model example, we showcase the effectiveness and explain some features of our method by finding the ground state of the U(1)-symmetric infinite-dimensional antiferromagnetic XX Heisenberg model.

arXiv:1910.10370•Journal of Physics Condensed Matter•Apr 14, 2020•1 citation
We combine the density matrix renormalization group (DMRG) with matrix product state tangent space concepts to construct a variational algorithm for finding ground states of one-dimensional quantum lattices in the thermodynamic limit. A careful comparison of this variational uniform matrix product state algorithm (VUMPS) with infinite density matrix renormalization group (IDMRG) and with infinite time evolving block decimation (ITEBD) reveals substantial gains in convergence speed and precision. We also demonstrate that VUMPS works very efficiently for Hamiltonians with long-range interactions and also for the simulation of two-dimensional models on infinite cylinders. The new algorithm can be conveniently implemented as an extension of an already existing DMRG implementation.
We show that the time-dependent variational principle provides a unifying framework for time-evolution methods and optimization methods in the context of matrix product states. In particular, we introduce a new integration scheme for studying time evolution, which can cope with arbitrary Hamiltonians, including those with long-range interactions. Rather than a Suzuki-Trotter splitting of the Hamiltonian, which is the idea behind the adaptive time-dependent density matrix renormalization group method or time-evolving block decimation, our method is based on splitting the projector onto the matrix product state tangent space as it appears in the Dirac-Frenkel time-dependent variational principle. We discuss how the resulting algorithm resembles the density matrix renormalization group (DMRG) algorithm for finding ground states so closely that it can be implemented by changing just a few lines of code and it inherits the same stability and efficiency. In particular, our method is compatible with any Hamiltonian for which ground-state DMRG can be implemented efficiently. In fact, DMRG is obtained as a special case of our scheme for imaginary time evolution with infinite time step.
Combining the time-dependent variational principle (TDVP) algorithm with the parallelization scheme introduced by Stoudenmire and White for the density matrix renormalization group (DMRG), we present the first parallel matrix product state (MPS) algorithm capable of time evolving one-dimensional (1D) quantum lattice systems with long-range interactions. We benchmark the accuracy and performance of the algorithm by simulating quenches in the long-range Ising and XY models. We show that our code scales well up to 32 processes, with parallel efficiencies as high as 8686%86. Finally, we calculate the dynamical correlation function of a 201-site Heisenberg XXX spin chain with 1/r21/{r}^{2}1/r2 interactions, which is challenging to compute sequentially. These results pave the way for the application of tensor networks to increasingly complex many-body systems.
The infinite density matrix renormalization group (iDMRG) algorithm is a highly successful numerical algorithm for the study of low-dimensional quantum systems, and is also frequently used to initialize the more popular finite DMRG algorithm. Implementations of both finite and infinite DMRG frequently incorporate support for the protection and exploitation of symmetries of the Hamiltonian. In common with other variational tensor network algorithms, convergence of iDMRG to the ground state is not guaranteed, with the risk that the algorithm may become stuck in a local minimum. In this paper, I demonstrate the existence of a particularly harmful class of physically irrelevant local minima affecting both iDMRG and to a lesser extent also infinite time-evolving block decimation (iTEBD), for which the ground state is compatible with the protected symmetries of the Hamiltonian but cannot be reached using the conventional iDMRG or iTEBD algorithms. I describe a modified iDMRG algorithm which evades these local minima, and which also admits a natural interpretation on topologically ordered systems with a boundary.
I revisit the infinite-size variant of the Density Matrix Renormalization Group (iDMRG) algorithm for obtaining a fixed-point translationally invariant matrix product wavefunction in the context of one-dimensional quantum systems. A crucial ingredient of this algorithm is an efficient transformation for obtaining the matrix elements of the wavefunction as the lattice size is increased, and I introduce here a versatile transformation that is demonstrated to be much more effective than previous versions. The resulting algorithm has a surprisingly close relationship to Vidal's Time Evolving Block Decimation for infinite systems, but allows much faster convergence. Access to a translationally invariant matrix product state allows the calculation of correlation functions based on the transfer matrix, which directly gives the spectrum of all correlation lengths. I also show some advantages of the Matrix Product Operator (MPO) technique for constructing expectation values of higher moments, such as the exact variance $$.
In these lecture notes we give a technical overview of tangent-space methods for matrix product states in the thermodynamic limit. We introduce the manifold of uniform matrix product states, show how to compute different types of observables, and discuss the concept of a tangent space. We explain how to variationally optimize ground-state approximations, implement real-time evolution and describe elementary excitations for a given model Hamiltonian. Also, we explain how matrix product states approximate fixed points of one-dimensional transfer matrices. We show how all these methods can be translated to the language of continuous matrix product states for one-dimensional field theories. We conclude with some extensions of the tangent-space formalism and with an outlook to new applications.
Abstract Understanding extreme non-locality in many-body quantum systems can help resolve questions in thermostatistics and laser physics. The existence of symmetry selection rules for Hamiltonians with non-decaying terms on infinite-size lattices can lead to finite energies per site, which deserves attention. Here, we present a tensor network approach to construct the ground states of nontrivial symmetric infinite-dimensional spin Hamiltonians based on constrained optimizations of their infinite matrix product states description, which contains no truncation step, offers a very simple mathematical structure, and other minor advantages at the cost of slightly higher polynomial complexity in comparison to an existing method. More precisely speaking, our proposed algorithm is in part equivalent to the more generic and well-established solvers of infinite density-matrix renormalization-group and variational uniform matrix product states, which are, in principle, capable of accurately representing the ground states of such infinite-range-interacting many-body systems. However, we employ some mathematical simplifications that would allow for efficient brute-force optimizations of tensor-network matrices for the specific cases of highly-symmetric infinite-size infinite-range models. As a toy-model example, we showcase the effectiveness and explain some features of our method by finding the ground state of the U(1)-symmetric infinite-dimensional antiferromagnetic XX Heisenberg model.