A Fast and Efficient Algorithm for Slater Determinant Updates in Quantum Monte Carlo Simulations
arXiv:0906.4354 · doi:10.1063/1.3142703
Abstract
We present an efficient low-rank updating algorithm for updating the trial wavefunctions used in Quantum Monte Carlo (QMC) simulations. The algorithm is based on low-rank updating of the Slater determinants. In particular, the computational complexity of the algorithm is O(kN) during the k-th step compared with traditional algorithms that require O(N^2) computations, where N is the system size. For single determinant trial wavefunctions the new algorithm is faster than the traditional O(N^2) Sherman-Morrison algorithm for up to O(N) updates. For multideterminant configuration-interaction type trial wavefunctions of M+1 determinants, the new algorithm is significantly more efficient, saving both O(MN^2) work and O(MN^2) storage. The algorithm enables more accurate and significantly more efficient QMC calculations using configuration interaction type wavefunctions.
References in corpus (5)
- Full optimization of Jastrow-Slater wave functions with application to the first-row atoms and homonuclear diatomic molecules
- Neutral and charged excitations in carbon fullerenes from first-principles many-body theories
- Energies of the first row atoms from quantum Monte Carlo
- Spectroscopic data for the LiH molecule from pseudopotential quantum Monte Carlo calculations
- Density-density functionals and effective potentials in many-body electronic structure calculations