Time-dependent variational principle with ancillary Krylov subspace

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.