skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Multipole moments in the effective fragment potential method

Abstract

In the effective fragment potential (EFP) method the Coulomb potential is represented using a set of multipole moments generated by the distributed multipole analysis (DMA) method. Misquitta, Stone, and Fazeli recently developed a basis space-iterated stockholder atom (BS-ISA) method to generate multipole moments. This study assesses the accuracy of the EFP interaction energies using sets of multipole moments generated from the BS-ISA method, and from several versions of the DMA method (such as analytic and numeric grid-based), with varying basis sets. Both methods lead to reasonable results, although using certain implementations of the DMA method can result in large errors. With respect to the CCSD(T)/CBS interaction energies, the mean unsigned error (MUE) of the EFP method for the S22 data set using BS-ISA–generated multipole moments and DMA-generated multipole moments (using a small basis set and the analytic DMA procedure) is 0.78 and 0.72 kcal/mol, respectively. Here, the MUE accuracy is on the same order as MP2 and SCS-MP2. The MUEs are lower than in a previous study benchmarking the EFP method without the EFP charge transfer term, demonstrating that the charge transfer term increases the accuracy of the EFP method. Regardless of the multipole moment method used, it is likelymore » that much of the error is due to an insufficient short-range electrostatic term (i.e., charge penetration term), as shown by comparisons with symmetry-adapted perturbation theory.« less

Authors:
 [1];  [2];  [3]; ORCiD logo [1]
  1. Ames Lab. and Iowa State Univ., Ames, IA (United States)
  2. Purdue Univ., West Lafayette, IN (United States)
  3. Univ. of London, London (England)
Publication Date:
Research Org.:
Ames Laboratory (AMES), Ames, IA (United States)
Sponsoring Org.:
USDOE Office of Science (SC), Basic Energy Sciences (BES) (SC-22)
OSTI Identifier:
1355900
Report Number(s):
IS-J-9325
Journal ID: ISSN 1089-5639
Grant/Contract Number:
CHE-1465154; AC02-07CH11358
Resource Type:
Journal Article: Accepted Manuscript
Journal Name:
Journal of Physical Chemistry. A, Molecules, Spectroscopy, Kinetics, Environment, and General Theory
Additional Journal Information:
Journal Volume: 121; Journal Issue: 9; Journal ID: ISSN 1089-5639
Publisher:
American Chemical Society
Country of Publication:
United States
Language:
English
Subject:
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CHEMISTRY

Citation Formats

Bertoni, Colleen, Slipchenko, Lyudmila V., Misquitta, Alston J., and Gordon, Mark S. Multipole moments in the effective fragment potential method. United States: N. p., 2017. Web. doi:10.1021/acs.jpca.7b00682.
Bertoni, Colleen, Slipchenko, Lyudmila V., Misquitta, Alston J., & Gordon, Mark S. Multipole moments in the effective fragment potential method. United States. doi:10.1021/acs.jpca.7b00682.
Bertoni, Colleen, Slipchenko, Lyudmila V., Misquitta, Alston J., and Gordon, Mark S. Fri . "Multipole moments in the effective fragment potential method". United States. doi:10.1021/acs.jpca.7b00682. https://www.osti.gov/servlets/purl/1355900.
@article{osti_1355900,
title = {Multipole moments in the effective fragment potential method},
author = {Bertoni, Colleen and Slipchenko, Lyudmila V. and Misquitta, Alston J. and Gordon, Mark S.},
abstractNote = {In the effective fragment potential (EFP) method the Coulomb potential is represented using a set of multipole moments generated by the distributed multipole analysis (DMA) method. Misquitta, Stone, and Fazeli recently developed a basis space-iterated stockholder atom (BS-ISA) method to generate multipole moments. This study assesses the accuracy of the EFP interaction energies using sets of multipole moments generated from the BS-ISA method, and from several versions of the DMA method (such as analytic and numeric grid-based), with varying basis sets. Both methods lead to reasonable results, although using certain implementations of the DMA method can result in large errors. With respect to the CCSD(T)/CBS interaction energies, the mean unsigned error (MUE) of the EFP method for the S22 data set using BS-ISA–generated multipole moments and DMA-generated multipole moments (using a small basis set and the analytic DMA procedure) is 0.78 and 0.72 kcal/mol, respectively. Here, the MUE accuracy is on the same order as MP2 and SCS-MP2. The MUEs are lower than in a previous study benchmarking the EFP method without the EFP charge transfer term, demonstrating that the charge transfer term increases the accuracy of the EFP method. Regardless of the multipole moment method used, it is likely that much of the error is due to an insufficient short-range electrostatic term (i.e., charge penetration term), as shown by comparisons with symmetry-adapted perturbation theory.},
doi = {10.1021/acs.jpca.7b00682},
journal = {Journal of Physical Chemistry. A, Molecules, Spectroscopy, Kinetics, Environment, and General Theory},
number = 9,
volume = 121,
place = {United States},
year = {Fri Feb 17 00:00:00 EST 2017},
month = {Fri Feb 17 00:00:00 EST 2017}
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record

Save / Share:
  • No abstract prepared.
  • This study compares the results of the general effective fragment potential (EFP2) method to the results of a previous combined coupled cluster with single, double, and perturbative triple excitations [CCSD(T)] and symmetry-adapted perturbation theory (SAPT) study [Sinnokrot and Sherrill, J. Am. Chem. Soc., 2004, 126, 7690] on substituent effects in {pi}-{pi} interactions. EFP2 is found to accurately model the binding energies of the benzene-benzene, benzene-phenol, benzene-toluene, benzene-fluorobenzene, and benzene-benzonitrile dimers, as compared with high-level methods [Sinnokrot and Sherrill, J. Am. Chem. Soc., 2004, 126, 7690], but at a fraction of the computational cost of CCSD(T). In addition, an EFP-based Montemore » Carlo/simulated annealing study was undertaken to examine the potential energy surface of the substituted dimers.« less
  • The design and development of scientific software components to provide an interface to the effective fragment potential (EFP) methods are reported. Multiscale modeling of physical and chemical phenomena demands the merging of software packages developed by research groups in significantly different fields. Componentization offers an efficient way to realize new high performance scientific methods by combining the best models available in different software packages without a need for package readaptation after the initial componentization is complete. The EFP method is an efficient electronic structure theory based model potential that is suitable for predictive modeling of intermolecular interactions in large molecularmore » systems, such as liquids, proteins, atmospheric aerosols, and nanoparticles, with an accuracy that is comparable to that of correlated ab initio methods. The developed components make the EFP functionality accessible for any scientific component-aware software package. The performance of the component is demonstrated on a protein interaction model, and its accuracy is compared with results obtained with coupled cluster methods.« less
  • The charge transfer (CT) interaction, the most time-consuming term in the general effective fragment potential method, is made much more computationally efficient. This is accomplished by the projection of the quasiatomic minimal-basis-set orbitals (QUAMBOs) as the atomic basis onto the self-consistent field virtual molecular orbital (MO) space to select a subspace of the full virtual space called the valence virtual space. The diagonalization of the Fock matrix in terms of QUAMBOs recovers the canonical occupied orbitals and, more importantly, gives rise to the valence virtual orbitals (VVOs). The CT energies obtained using VVOs are generally as accurate as those obtainedmore » with the full virtual space canonical MOs because the QUAMBOs span the valence part of the virtual space, which can generally be regarded as “chemically important.” The number of QUAMBOs is the same as the number of minimal-basis MOs of a molecule. Therefore, the number of VVOs is significantly smaller than the number of canonical virtual MOs, especially for large atomic basis sets. This leads to a dramatic decrease in the computational cost.« less
  • This work presents an extension of the linear response TDDFT/EFP method to the nonlinear-response regime together with the implementation of nonlinear-response TDDFT/EFP in the quantum-chemistry computer package GAMESS. Included in the new method is the ability to calculate the two-photon absorption cross section and to incorporate solvent effects via the EFP method. The nonlinear-response TDDFT/EFP method is able to make correct qualitative predictions for both gas phase values and aqueous solvent shifts of several important nonlinear properties.