Header

UZH-Logo

Maintenance Infos

Second-order moller-plesset perturbation theory in the condensed phase: an efficient and massively parallel gaussian and plane waves approach


Del Ben, Mauro; Hutter, Juerg; VandeVondele, Joost (2012). Second-order moller-plesset perturbation theory in the condensed phase: an efficient and massively parallel gaussian and plane waves approach. Journal of Chemical Theory and Computation, 8(11):4177-4188.

Abstract

A novel algorithm, based on a hybrid Gaussian and plane waves (GPW) approach, is developed for the canonical second order Moller-Plesset perturbation energy (MP2) of finite and extended systems. The key aspect of the method is that the electron repulsion integrals (ia vertical bar lambda sigma) are computed by direct integration between the products of Gaussian basis functions lambda sigma and the electrostatic potential arising from a given occupied virtual pair density ia. The electrostatic potential is obtained in a plane waves basis set after solving the Poisson equation in Fourier space. In particular, for condensed phase systems, this scheme is highly efficient. Furthermore, our implementation has low memory requirements and displays excellent parallel scalability up to 100 000 processes. In this way, canonical MP2 calculations for condensed phase systems containing hundreds of atoms or more than 5000 basis functions can be performed within minutes, while systems up to 1000 atoms and 10 000 basis functions remain feasible. Solid LiH has been employed as a benchmark to study basis set and system size convergence. Lattice constants and cohesive energies of various molecular crystals have been studied with MP2 and double-hybrid functionals.

Abstract

A novel algorithm, based on a hybrid Gaussian and plane waves (GPW) approach, is developed for the canonical second order Moller-Plesset perturbation energy (MP2) of finite and extended systems. The key aspect of the method is that the electron repulsion integrals (ia vertical bar lambda sigma) are computed by direct integration between the products of Gaussian basis functions lambda sigma and the electrostatic potential arising from a given occupied virtual pair density ia. The electrostatic potential is obtained in a plane waves basis set after solving the Poisson equation in Fourier space. In particular, for condensed phase systems, this scheme is highly efficient. Furthermore, our implementation has low memory requirements and displays excellent parallel scalability up to 100 000 processes. In this way, canonical MP2 calculations for condensed phase systems containing hundreds of atoms or more than 5000 basis functions can be performed within minutes, while systems up to 1000 atoms and 10 000 basis functions remain feasible. Solid LiH has been employed as a benchmark to study basis set and system size convergence. Lattice constants and cohesive energies of various molecular crystals have been studied with MP2 and double-hybrid functionals.

Statistics

Citations

47 citations in Web of Science®
59 citations in Scopus®
Google Scholar™

Altmetrics

Downloads

221 downloads since deposited on 13 Mar 2013
22 downloads since 12 months
Detailed statistics

Additional indexing

Item Type:Journal Article, refereed, original work
Communities & Collections:07 Faculty of Science > Department of Chemistry
Dewey Decimal Classification:540 Chemistry
Language:English
Date:2012
Deposited On:13 Mar 2013 14:45
Last Modified:07 Dec 2017 20:35
Publisher:American Chemical Society
ISSN:1549-9618
Publisher DOI:https://doi.org/10.1021/ct300531w

Download

Download PDF  'Second-order moller-plesset perturbation theory in the condensed phase: an efficient and massively parallel gaussian and plane waves approach'.
Preview
Content: Accepted Version
Filetype: PDF
Size: 2MB
View at publisher