We present an efficient implementation of dynamically pruned quantum dynamics, both in coordinate space and in phase space. We combine the ideas behind the biorthogonal von Neumann basis (PvB) with the orthogonalized momentum-symmetrized Gaussians (Weylets) to create a new basis, projected Weylets, that takes the best from both methods. We benchmark pruned time-dependent dynamics using phase-space-localized PvB, projected Weylets, and coordinate-space-localized DVR bases, with real-world examples in up to six dimensions. For the examples studied, coordinate-space localization is the most important factor for efficient pruning and the pruned dynamics is much faster than the unpruned, exact dynamics. Phase-space localization is useful for more demanding dynamics where many basis functions are required. There, projected Weylets offer a more compact representation than pruned DVR bases.
REFERENCES
We use std::unordered_set of the C++ programming language standard library as implemented in the GNU compiler collection.112
Brown and Carrington noted problems with the condition number of G, namely numbers after an adaption of the width of each von Neumann function for a grid of size 48.88 We do not experience such problems.