Skip to main content

Comparison of polynomial approximations to speed up planewave-based quantum Monte Carlo calculations

Cornell Affiliated Author(s)

Author

W.D. Parker
C.J. Umrigar
D. Alfè
F.R. Petruzielo
R.G. Hennig
J.W. Wilkins

Abstract

The computational cost of quantum Monte Carlo (QMC) calculations of realistic periodic systems depends strongly on the method of storing and evaluating the many-particle wave function. Previous work by Williamson et al. (2001) [35] and Alfè and Gillan, (2004) [36] has demonstrated the reduction of the O(N3) cost of evaluating the Slater determinant with planewaves to O(N2) using localized basis functions. We compare four polynomial approximations as basis functions - interpolating Lagrange polynomials, interpolating piecewise-polynomial-form (pp-) splines, and basis-form (B-) splines (interpolating and smoothing). All these basis functions provide a similar speedup relative to the planewave basis. The pp-splines have eight times the memory requirement of the other methods. To test the accuracy of the basis functions, we apply them to the ground state structures of Si, Al, and MgO. The polynomial approximations differ in accuracy most strongly for MgO, and smoothing B-splines most closely reproduce the planewave value for of the variational Monte Carlo energy. Using separate approximations for the Laplacian of the orbitals increases the accuracy sufficiently to justify the increased memory requirement, making smoothing B-splines, with separate approximation for the Laplacian, the preferred choice for approximating planewave-represented orbitals in QMC calculations. © 2015 Elsevier Inc.

Date Published

Journal

Journal of Computational Physics

Volume

287

Number of Pages

77-87,

URL

https://www.scopus.com/inward/record.uri?eid=2-s2.0-84923044488&doi=10.1016%2fj.jcp.2015.01.037&partnerID=40&md5=815f6ff7bcc2f02ebac2a86f682ae35c

DOI

10.1016/j.jcp.2015.01.037

Group (Lab)

Cyrus Umrigar Group

Funding Source

1542776
EP/K038249/1

Download citation