TY - JOUR

T1 - Electrostatically embedded many-body expansion for large systems, with applications to water clusters

AU - Dahlke, Erin E.

AU - Truhlar, Donald G.

PY - 2007/1

Y1 - 2007/1

N2 - The use of background molecular charge to incorporate environmental effects on a molecule or active site is widely employed in quantum chemistry. In the present article we employ this practice in conjunction with many-body expansions. In particular, we present electrostatically embedded two-body and three-body expansions for calculating the energies of molecular clusters. The system is divided into fragments, and dimers or trimers of fragments are calculated in a field of point charges representing the electrostatic potential of the other fragments. We find that including environmental point charges can lower the errors in the electrostatically embedded pairwise additive (EE-PA) energies for a series of water clusters by as much as a factor of 10 when compared to the traditional pairwise additive approximation and that for the electrostatically embedded three-body (EE-3B) method the average mean unsigned error over nine different levels of theory for a set of six tetramers and one pentamer is only 0.05 kcal/mol, which is only 0.4% of the mean unsigned net interaction energy. We also test the accuracy of the EE-PA and EE-3B methods for a cluster of 21 water molecules and find that the errors relative to a full MP2/aug'-cc-pVTZ calculation to be only 2.97 and 0.38 kcal/mol, respectively, which are only 1.5% and 0.2%, respectively, of the net interaction energy. This method offers the advantage over some other fragment-based methods in that it does not use an iterative method to determine the charges and thus provides substantial savings for large clusters. The method is convenient to adapt to a variety of electronic structure methods and program packages, it has N 2 or N3 computational scaling for large systems (where N is the number of fragments), it is easily converted to an O(N) method, and its linearity allows for convenient analytic gradients.

AB - The use of background molecular charge to incorporate environmental effects on a molecule or active site is widely employed in quantum chemistry. In the present article we employ this practice in conjunction with many-body expansions. In particular, we present electrostatically embedded two-body and three-body expansions for calculating the energies of molecular clusters. The system is divided into fragments, and dimers or trimers of fragments are calculated in a field of point charges representing the electrostatic potential of the other fragments. We find that including environmental point charges can lower the errors in the electrostatically embedded pairwise additive (EE-PA) energies for a series of water clusters by as much as a factor of 10 when compared to the traditional pairwise additive approximation and that for the electrostatically embedded three-body (EE-3B) method the average mean unsigned error over nine different levels of theory for a set of six tetramers and one pentamer is only 0.05 kcal/mol, which is only 0.4% of the mean unsigned net interaction energy. We also test the accuracy of the EE-PA and EE-3B methods for a cluster of 21 water molecules and find that the errors relative to a full MP2/aug'-cc-pVTZ calculation to be only 2.97 and 0.38 kcal/mol, respectively, which are only 1.5% and 0.2%, respectively, of the net interaction energy. This method offers the advantage over some other fragment-based methods in that it does not use an iterative method to determine the charges and thus provides substantial savings for large clusters. The method is convenient to adapt to a variety of electronic structure methods and program packages, it has N 2 or N3 computational scaling for large systems (where N is the number of fragments), it is easily converted to an O(N) method, and its linearity allows for convenient analytic gradients.

UR - http://www.scopus.com/inward/record.url?scp=34248152064&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=34248152064&partnerID=8YFLogxK

U2 - 10.1021/ct600253j

DO - 10.1021/ct600253j

M3 - Article

AN - SCOPUS:34248152064

VL - 3

SP - 46

EP - 53

JO - Journal of Chemical Theory and Computation

JF - Journal of Chemical Theory and Computation

SN - 1549-9618

IS - 1

ER -