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:2005.06104 · Physical review. B./Physical review. B

Time-dependent variational principle with ancillary Krylov subspace

Mingru Yang, Steven R. White

Published Sep 29, 2020• 136 citations
arXivHTMLPDF

Abstract

We propose an improved scheme to do the time-dependent variational principle (TDVP) in finite matrix product states (MPSs) for two-dimensional systems or one-dimensional systems with long-range interactions. We present a method to represent the time-evolving state in a MPS with its basis enriched by state averaging with global Krylov vectors. We show that the projection error is significantly reduced so precise time evolution can still be obtained even if a larger time step is used. Combined with the one-site TDVP, our approach provides a way to dynamically increase the bond dimension while still preserving unitarity for real-time evolution. Our method can be more accurate and exhibit slower bond dimension growth than the conventional two-site TDVP.

DOI
10.1103/physrevb.102.094315

Related papers

Time-Dependent Variational Principle with Controlled Bond Expansion for Matrix Product States

Jheng-Wei Li, Andreas Gleis, Jan von Delft

We present a controlled bond expansion (CBE) approach to simulate quantum dynamics based on the time-dependent variational principle (TDVP) for matrix product states. Our method alleviates the numerical difficulties of the standard, fixed-rank one-site TDVP integrator by increasing bond dimensions on the fly to reduce the projection error. This is achieved in an economical, local fashion, requiring only minor modifications of standard one-site TDVP implementations. We illustrate the performance and accuracy of CBE-TDVP with several numerical examples on finite quantum lattices, including new results on bipolaron formation in the Peierls-Hubbard model and spin pumping via adiabatic flux insertion in a chiral spin liquid.

arXiv:2208.10972•Physical Review Letters•Jul 9, 2024•26 citations

Stochastic Adaptive Single-Site Time-Dependent Variational Principle

Yihe Xu, Zhaoxuan Xie, Xiaoyu Xie +2 more

In recent years, the time-dependent variational principle (TDVP) method based on the matrix product state (MPS) wave function formulation has shown its great power in performing large-scale quantum dynamics simulations for realistic chemical systems with strong electron-vibration interactions. In this work, we propose a stochastic adaptive single-site TDVP (SA-1TDVP) scheme to evolve the bond-dimension adaptively, which can integrate the traditional advantages of both the high efficiency of the single-site TDVP (1TDVP) variant and the high accuracy of the two-site TDVP (2TDVP) variant. Based on the assumption that the level statistics of entanglement Hamiltonians, which originate from the reduced density matrices of the MPS method, follows a Poisson or Wigner distribution, as generically predicted by random-matrix theory, additional random singular values are generated to expand the bond-dimension automatically. Tests on simulating the vibrationally resolved quantum dynamics and absorption spectra in the pyrazine molecule and perylene bisimide (PBI) J-aggregate trimer as well as a spin-1/2 Heisenberg chain show that it can be automatic and as accurate as 2TDVP but reduce the computational time remarkably.

arXiv:2110.12703•JACS Au•Jan 18, 2022•23 citations

Performance of the time-dependent variational principle for matrix product states in the long-time evolution of a pure state

Shimpei Goto, Ippei Danshita

The projection of time-dependent variational principle (TDVP) for matrix product states enables us to perform long-time simulations of one-dimensional quantum systems with the conservation of the total energy and the norm of wave functions. We compare long-time dynamics after a quantum quench simulated by TDVP with those by the exact diagonalization method in order to evaluate the performance of TDVP. We show that in a nonintegrable model the projection of TDVP clearly improves the long-time behaviors of global observables included in the Hamiltonian, such as the kinetic and interaction energies. In contrast, this projection can lead to larger error for other observables than that caused by the truncation of states.

arXiv:1809.01400•Physical review. B./Physical review. B•Feb 21, 2019•47 citations

Diversifying time evolution of matrix product states using BUGs

Kornél Kapás, Dominik Sulz, Charlotte Verhoeven +3 more

Time evolution of tensor-network states is a central tool for studying nonequilibrium dynamics in strongly correlated quantum systems. In the past decade, using the time-dependent variational principle (TDVP) discretized by projector-splitting has become the gold standard for evolving matrix product states (MPS) with long-range Hamiltonians. More recently, a conceptually different class of discretizations of the TDVP, the class of Basis Update and Galerkin (BUG) integrators, has been developed within the framework of dynamical low-rank approximation, offering an alternative, more flexible approach. In this work, we formulate two BUG integrators explicitly for MPS. The BUG integrators result in "single-site" MPS integrators that are robust with respect to small singular values, avoid backward-in-time substeps, and naturally allow for the adaptive choice of bond dimensions. We consider the augmented BUG integrator of first order and the midpoint BUG integrator of second order and clarify their relation to the projector-splitting discretization of the TDVP. Using representative models of strongly correlated quantum systems, like the XY spin model with long-range interaction and two-dimensional lattice models of fermions and spins of various lengths, we perform a systematic numerical comparison between BUG and projector-splitting time discretizations of the continuous-time TDVP, focusing on scalability, numerical complexity, and accuracy. We find that BUG integrators, due to their reduced numerical cost and the flexibility of the concept, provide a competitive and, in certain regimes, advantageous alternative to the well-established time-stepping by projector-splitting in TDVP, especially in models with large local dimensions, and in cases where observables are extracted with high temporal resolution.

cond-mat.str-elCross-listedarXiv:2609.12848•Sep 14, 2026•0 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

Improving matrix product state methods for long-range interactions

Mingru Yang

We study one dimensional models of diatomic molecules where both the electrons and nuclei are treated as quantum particles, going beyond the usual Born-Oppenheimer approximation. The continuous system is approximated by a grid which computationally resembles a ladder, with the electrons living on one leg and the nuclei on the other. To simulate DMRG efficiently with this system, a three-site algorithm has been implemented. We also use a compression method to treat the long-range interactions between charged particles. We find that 1D diatomic molecules with spin-1/2 nuclei in the spin-triplet state will unbind when the mass of the nuclei reduces to only a few times larger than the electron mass, while the molecule with nuclei in the singlet state always binds, given the two electrons in their singlet state in both cases. \nWe propose an improved scheme to do the time dependent variational principle (TDVP) in finite matrix product states (MPS) for two-dimensional systems or one-dimensional systems with long range interactions. We present a method to represent the time-evolving state in a MPS with its basis enriched by state-averaging with global Krylov vectors. We show that the projection error is significantly reduced so that precise time evolution can still be obtained even if a larger time step is used. Combined with the one-site TDVP, our approach provides a way to dynamically increase the bond dimension while still preserving unitarity for real time evolution. Our method can be more accurate and exhibit slower bond dimension growth than the conventional two-site TDVP.\n \nWe apply our improved TDVP method to investigate the spin squeezing dynamics of the XXZ model of $1/r^\\alpha$ interaction in two dimension. Comparing with the spin squeezing parameter and other observables obtained from discrete truncated Wigner approximation (DTWA), we verify the validity of DTWA and unveil the potential for this method to study dynamics of large-scale spin systems. Our results confirm the existence of large collective regime when $\\alpha > 2$, which can be a guide for future experimental realizations. Combined with the purification method, our improved TDVP method is proved to be also useful to study the thermalization of the long-range interacting system.

•Jan 1, 2020•0 citations
We propose an improved scheme to do the time-dependent variational principle (TDVP) in finite matrix product states (MPSs) for two-dimensional systems or one-dimensional systems with long-range interactions. We present a method to represent the time-evolving state in a MPS with its basis enriched by state averaging with global Krylov vectors. We show that the projection error is significantly reduced so precise time evolution can still be obtained even if a larger time step is used. Combined with the one-site TDVP, our approach provides a way to dynamically increase the bond dimension while still preserving unitarity for real-time evolution. Our method can be more accurate and exhibit slower bond dimension growth than the conventional two-site TDVP.
We present a controlled bond expansion (CBE) approach to simulate quantum dynamics based on the time-dependent variational principle (TDVP) for matrix product states. Our method alleviates the numerical difficulties of the standard, fixed-rank one-site TDVP integrator by increasing bond dimensions on the fly to reduce the projection error. This is achieved in an economical, local fashion, requiring only minor modifications of standard one-site TDVP implementations. We illustrate the performance and accuracy of CBE-TDVP with several numerical examples on finite quantum lattices, including new results on bipolaron formation in the Peierls-Hubbard model and spin pumping via adiabatic flux insertion in a chiral spin liquid.
In recent years, the time-dependent variational principle (TDVP) method based on the matrix product state (MPS) wave function formulation has shown its great power in performing large-scale quantum dynamics simulations for realistic chemical systems with strong electron-vibration interactions. In this work, we propose a stochastic adaptive single-site TDVP (SA-1TDVP) scheme to evolve the bond-dimension adaptively, which can integrate the traditional advantages of both the high efficiency of the single-site TDVP (1TDVP) variant and the high accuracy of the two-site TDVP (2TDVP) variant. Based on the assumption that the level statistics of entanglement Hamiltonians, which originate from the reduced density matrices of the MPS method, follows a Poisson or Wigner distribution, as generically predicted by random-matrix theory, additional random singular values are generated to expand the bond-dimension automatically. Tests on simulating the vibrationally resolved quantum dynamics and absorption spectra in the pyrazine molecule and perylene bisimide (PBI) J-aggregate trimer as well as a spin-1/2 Heisenberg chain show that it can be automatic and as accurate as 2TDVP but reduce the computational time remarkably.
The projection of time-dependent variational principle (TDVP) for matrix product states enables us to perform long-time simulations of one-dimensional quantum systems with the conservation of the total energy and the norm of wave functions. We compare long-time dynamics after a quantum quench simulated by TDVP with those by the exact diagonalization method in order to evaluate the performance of TDVP. We show that in a nonintegrable model the projection of TDVP clearly improves the long-time behaviors of global observables included in the Hamiltonian, such as the kinetic and interaction energies. In contrast, this projection can lead to larger error for other observables than that caused by the truncation of states.
Time evolution of tensor-network states is a central tool for studying nonequilibrium dynamics in strongly correlated quantum systems. In the past decade, using the time-dependent variational principle (TDVP) discretized by projector-splitting has become the gold standard for evolving matrix product states (MPS) with long-range Hamiltonians. More recently, a conceptually different class of discretizations of the TDVP, the class of Basis Update and Galerkin (BUG) integrators, has been developed within the framework of dynamical low-rank approximation, offering an alternative, more flexible approach. In this work, we formulate two BUG integrators explicitly for MPS. The BUG integrators result in "single-site" MPS integrators that are robust with respect to small singular values, avoid backward-in-time substeps, and naturally allow for the adaptive choice of bond dimensions. We consider the augmented BUG integrator of first order and the midpoint BUG integrator of second order and clarify their relation to the projector-splitting discretization of the TDVP. Using representative models of strongly correlated quantum systems, like the XY spin model with long-range interaction and two-dimensional lattice models of fermions and spins of various lengths, we perform a systematic numerical comparison between BUG and projector-splitting time discretizations of the continuous-time TDVP, focusing on scalability, numerical complexity, and accuracy. We find that BUG integrators, due to their reduced numerical cost and the flexibility of the concept, provide a competitive and, in certain regimes, advantageous alternative to the well-established time-stepping by projector-splitting in TDVP, especially in models with large local dimensions, and in cases where observables are extracted with high temporal resolution.
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.
We study one dimensional models of diatomic molecules where both the electrons and nuclei are treated as quantum particles, going beyond the usual Born-Oppenheimer approximation. The continuous system is approximated by a grid which computationally resembles a ladder, with the electrons living on one leg and the nuclei on the other. To simulate DMRG efficiently with this system, a three-site algorithm has been implemented. We also use a compression method to treat the long-range interactions between charged particles. We find that 1D diatomic molecules with spin-1/2 nuclei in the spin-triplet state will unbind when the mass of the nuclei reduces to only a few times larger than the electron mass, while the molecule with nuclei in the singlet state always binds, given the two electrons in their singlet state in both cases. \nWe propose an improved scheme to do the time dependent variational principle (TDVP) in finite matrix product states (MPS) for two-dimensional systems or one-dimensional systems with long range interactions. We present a method to represent the time-evolving state in a MPS with its basis enriched by state-averaging with global Krylov vectors. We show that the projection error is significantly reduced so that precise time evolution can still be obtained even if a larger time step is used. Combined with the one-site TDVP, our approach provides a way to dynamically increase the bond dimension while still preserving unitarity for real time evolution. Our method can be more accurate and exhibit slower bond dimension growth than the conventional two-site TDVP.\n \nWe apply our improved TDVP method to investigate the spin squeezing dynamics of the XXZ model of $1/r^\\alpha$ interaction in two dimension. Comparing with the spin squeezing parameter and other observables obtained from discrete truncated Wigner approximation (DTWA), we verify the validity of DTWA and unveil the potential for this method to study dynamics of large-scale spin systems. Our results confirm the existence of large collective regime when alpha>2\\alpha > 2alpha>2, which can be a guide for future experimental realizations. Combined with the purification method, our improved TDVP method is proved to be also useful to study the thermalization of the long-range interacting system.