Electrostatically embedded symmetry-adapted perturbation theory.
Autor: | Glick CS; School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA.; Center for Computational Molecular Science and Technology, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA., Alenaizan A; School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA.; Center for Computational Molecular Science and Technology, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA., Cheney DL; Molecular Structure and Design, Bristol-Myers Squibb Company, P.O. Box 5400, Princeton, New Jersey 08543, USA., Cavender CE; Department of Biochemistry and Biophysics, University of Rochester School of Medicine and Dentistry, Rochester, New York 14642, USA., Sherrill CD; School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA.; Center for Computational Molecular Science and Technology, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA.; School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0765, USA. |
---|---|
Jazyk: | angličtina |
Zdroj: | The Journal of chemical physics [J Chem Phys] 2024 Oct 07; Vol. 161 (13). |
DOI: | 10.1063/5.0221974 |
Abstrakt: | Symmetry-adapted perturbation theory (SAPT) is an ab initio approach that directly computes noncovalent interaction energies in terms of electrostatics, exchange repulsion, induction/polarization, and London dispersion components. Due to its high computational scaling, routine applications of even the lowest order of SAPT are typically limited to a few hundred atoms. To address this limitation, we report here the addition of electrostatic embedding to the SAPT (EE-SAPT) and ISAPT (EE-ISAPT) methods. We illustrate the embedding scheme using water trimer as a prototype example. Then, we show that EE-SAPT/EE-ISAPT can be applied for efficiently and accurately computing noncovalent interactions in large systems, including solvated dimers and protein-ligand systems. In the latter application, particular care must be taken to properly handle the quantum mechanics/molecular mechanics boundary when it cuts covalent bonds. We investigate various schemes for handling charges near this boundary and demonstrate which are most effective in the context of charge-embedded SAPT. (© 2024 Author(s). Published under an exclusive license by AIP Publishing.) |
Databáze: | MEDLINE |
Externí odkaz: |