Delayed Slater determinant update algorithms for high efficiency quantum Monte Carlo

Abstract

Within ab initio Quantum Monte Carlo simulations, the leading numerical cost for large systems is the computation of the values of the Slater determinants in the trial wavefunction. Each Monte Carlo step requires finding the determinant of a dense matrix. This is most commonly iteratively evaluated using a rank-1 Sherman-Morrison updating scheme to avoid repeated explicit calculation of the inverse. The overall computational cost is therefore formally cubic in the number of electrons or matrix size. To improve the numerical efficiency of this procedure, we propose a novel multiple rank delayed update scheme. This strategy enables probability evaluation with application of accepted moves to the matrices delayed until after a predetermined number of moves, K. The accepted events are then applied to the matrices en bloc with enhanced arithmetic intensity and computational efficiency via matrix-matrix operations instead of matrix-vector operations. This procedure does not change the underlying Monte Carlo sampling or its statistical efficiency. For calculations on large systems and algorithms such as diffusion Monte Carlo where the acceptance ratio is high, order of magnitude improvements in the update time can be obtained on both multi-core CPUs and GPUs.

References in corpus (15)

• Ab initio quantum Monte Carlo calculations of spin superexchange in cuprates: the benchmarking case of Ca${}_2$CuO${}_3$
• Computing the energy of a water molecule using MultiDeterminants: A simple, efficient algorithm
• Phase Stability of TiO${}_2$ Polymorphs from Diffusion Quantum Monte Carlo
• The effect of electron correlation on the electronic structure and spin-lattice coupling of the high-T${}_c$ cuprates: quantum Monte Carlo calculations
• Discovering correlated fermions using quantum Monte Carlo
• Toward an improved control of the fixed-node error in quantum Monte Carlo: The case of the water molecule
• Simple formalism for efficient derivatives and multi-determinant expansions in quantum Monte Carlo
• Electronic Origin of the Volume Collapse in Cerium
• Quantum Monte Carlo study of the energetics of the rutile, anatase, brookite, and columbite TiO${}_2$ polymorphs
• Diffusion quantum Monte Carlo study of the equation of state and point defects in aluminum
• Phase Stability and Properties of Manganese Oxide Polymorphs: Assessment and Insights from Diffusion Monte Carlo
• Cohesive energy and structural parameters of binary oxides of groups IIA and IIIB from diffusion quantum Monte Carlo
• Towards a systematic assessment of errors in diffusion Monte Carlo calculations of semiconductors: case study of zinc selenide and zinc oxide
• Quantum Monte Carlo analysis of a charge ordered insulating antiferromagnet: the Ti${}_4$O${}_7$ Magnéli phase
• A Fast and Efficient Algorithm for Slater Determinant Updates in Quantum Monte Carlo Simulations

Cited by in corpus (11)

• QMCPACK : An open source ab initio Quantum Monte Carlo package for the electronic structure of atoms, molecules, and solids
• Phaseless Auxiliary-Field Quantum Monte Carlo on Graphical Processing Units
• Excited-state diffusion Monte Carlo calculations: a simple and efficient two-determinant ansatz
• Magnetic, thermodynamic, and dynamical properties of the three-dimensional fermionic Hubbard model: A comprehensive Monte Carlo study
• Extended Metal-Insulator Crossover with Strong Antiferromagnetic Spin Correlation in Half-Filled 3D Hubbard Model
• A Non-stochastic Optimization Algorithm for Neural-network Quantum States
• Nature of the mixed-parity pairing of attractive fermions with spin-orbit coupling in optical lattice
• Ising phase transitions and thermodynamics of correlated fermions in a two-dimensional spin-dependent lattice potential
• QMCkl: A Kernel Library for Quantum Monte Carlo Applications
• Spin-triplet pairing instability in a two-dimensional repulsive Hubbard model
• Reducing the Cost of Energy Differences in Variational Monte Carlo with Spotlight Sampling