title: Asymptotic Expansion for Electrostatic Embedding Integrals in QM/MM Calculations.
abstract: In QM/MM studies with large MM regions, the calculation of electrostatic embedding integrals can become a computational bottleneck. To overcome this problem, an asymptotic expansion for nuclear attraction-type integrals is developed. As a result, the long-range interactions between the QM and MM atoms reduce to atom-centered multipole moment-like expansions. The algorithm uses a natural spatial division of the molecular structure. To further improve the computational performance, a cutoff radius for the multipole moment-like expansion is introduced. The new code was validated and benchmarked with deMon2k/CHARMM QM/MM calculations on an RNA polymerase II model with almost 350 000 atoms. It is shown that the computational time for the calculation of the embedding integrals in this system can be reduced below 200 s on a small parallel architecture (eight cores) without a loss of accuracy.
authors.affiliation: Departamento de Química, CINVESTAV. Avenida Instituto Politécnico Nacional 2508, A.P. 14-740 México, D.F. 07000, México.
authors.forename: Aurelio
authors.initials: A
authors.lastname: Alvarez-Ibarra
authors.affiliation: Departamento de Química, CINVESTAV. Avenida Instituto Politécnico Nacional 2508, A.P. 14-740 México, D.F. 07000, México.
authors.forename: Andreas M
authors.initials: AM
authors.lastname: Köster
authors.affiliation: Department of Chemistry and Institute for Biocomplexity and Informatics, University of Calgary, 2500 University Drive NW, Calgary, Alberta, Canada T2N 1N4.
authors.forename: Rui
authors.initials: R
authors.lastname: Zhang
authors.affiliation: Department of Chemistry and Institute for Biocomplexity and Informatics, University of Calgary, 2500 University Drive NW, Calgary, Alberta, Canada T2N 1N4.
authors.forename: Dennis R
authors.initials: DR
authors.lastname: Salahub