Bridging the gap between quantum Monte Carlo and F12-methods

• Published in: Chemical physics 2012-06, Vol.401, p.36-44
• Main Authors: Chinnamsetty, Sambasiva Rao, Luo, Hongjun, Hackbusch, Wolfgang, Flad, Heinz-Jürgen, Uschmajew, André
• Format: Article
• Language: English
• Subjects: F12 method; Jastrow factor; Quantum Monte Carlo; Tensor product approximation; Wavelets
• ISSN: 0301-0104
• DOI: 10.1016/j.chemphys.2011.09.019

Tensor product approximation of pair-correlation functions: τ(x,y)≈∑κ=1κuk(1)(x1,y1)uk(2)(x2,y2)uk(3)(x3,y3) Pair-correlation function τ(x,y)∣∣x·y∣ =∣x∣∣y∣ of the He atom and corresponding tensor product approximation errors. [Display omitted] Tensor product approximation of pair-correlation functions opens a new route from quantum Monte Carlo (QMC) to explicitly correlated F12 methods. Thereby one benefits from stochastic optimization techniques used in QMC to get optimal pair-correlation functions which typically recover more than 85% of the total correlation energy. Our approach incorporates, in particular, core and core-valence correlation which are poorly described by homogeneous and isotropic ansatz functions usually applied in F12 calculations. We demonstrate the performance of the tensor product approximation by applications to atoms and small molecules. It turns out that the canonical tensor format is especially suitable for the efficient computation of two- and three-electron integrals required by explicitly correlated methods. The algorithm uses a decomposition of three-electron integrals, originally introduced by Boys and Handy and further elaborated by Ten-no in his 3d numerical quadrature scheme, which enables efficient computations in the tensor format. Furthermore, our method includes the adaptive wavelet approximation of tensor components where convergence rates are given in the framework of best N-term approximation theory.