1
|
Abstract
Vibrational energies are partitioned into the contributions of molecular parts called segments, for instance, residues in proteins. The fragment molecular orbital method is used to facilitate vibrational calculations of large systems at the DFTB and HF-3c levels. The vibrational analysis is combined with the partitioning of the electronic energy, yielding free-energy contributions of segments to the binding energy, pinpointing hot spots for drug discovery and other studies. The analysis is illustrated on two protein-ligand complexes in solution.
Collapse
Affiliation(s)
- Dmitri G Fedorov
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science and Technology (AIST), Central 2, Umezono 1-1-1, Tsukuba 305-8568, Japan
| |
Collapse
|
2
|
Tripathy V, Saha A, Raghavachari K. Electrostatically embedded molecules-in-molecules approach and its application to molecular clusters. J Comput Chem 2021; 42:719-734. [PMID: 33586802 DOI: 10.1002/jcc.26492] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/17/2020] [Revised: 11/25/2020] [Accepted: 01/16/2021] [Indexed: 11/06/2022]
Abstract
We report the application of our fragment-based quantum chemistry model MIM (Molecules-In-Molecules) with electrostatic embedding. The method is termed "EE-MIM (Electrostatically Embedded Molecules-In-Molecules)" and accounts for the missing electrostatic interactions in the subsystems resulting from fragmentation. Point charges placed at the atomic positions are used to represent the interaction of each subsystem with the rest of the molecule with minimal increase in the computational cost. We have carefully calibrated this model on a range of different sizes of clusters containing up to 57 water molecules. The fragmentation methods have been applied with the goal of reproducing the unfragmented total energy at the MP2/6-311G(d,p) level. Comparative analysis has been carried out between MIM and EE-MIM to gauge the impact of electrostatic embedding. Performance of several different parameters such as the type of charge and levels of fragmentation are analyzed for the prediction of absolute energies. The use of background charges in subsystem calculations improves the performance of both one- and two-layer MIM while it is noticeably important in the case of one-layer MIM. Embedded charges for two-layer MIM are obtained from a full system calculation at the low-level. For one-layer MIM, in the absence of a full system calculation, two different types of embedded charges, namely, Geometry dependent (GD) and geometry independent (GI) charges, are used. A self-consistent procedure is employed to obtain GD charges. We have further tested our method on challenging charged systems with stronger intermolecular interactions, namely, protonated ammonia clusters (containing up to 30 ammonia molecules). The observations are similar to water clusters with improved performance using embedded charges. Overall, the performance of NPA charges as embedded charges is found to be the best.
Collapse
Affiliation(s)
- Vikrant Tripathy
- Department of Chemistry, Indiana University, Bloomington, Indiana, USA
| | - Arjun Saha
- Department of Chemistry, University of Southern California, Los Angeles, California, USA
| | | |
Collapse
|
3
|
Nishimoto Y, Fedorov DG. The fragment molecular orbital method combined with density-functional tight-binding and periodic boundary conditions. J Chem Phys 2021; 154:111102. [PMID: 33752370 DOI: 10.1063/5.0039520] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Abstract
The density-functional tight-binding (DFTB) formulation of the fragment molecular orbital method is combined with periodic boundary conditions. Long-range electrostatics and dispersion are evaluated with the Ewald summation technique. The first analytic derivatives of the energy with respect to atomic coordinates and lattice parameters are formulated. The accuracy of the method is established in comparison to numerical gradients and DFTB without fragmentation. The largest elementary cell in this work has 1631 atoms. The method is applied to elucidate the polarization, charge transfer, and interactions in the solution.
Collapse
Affiliation(s)
- Yoshio Nishimoto
- Graduate School of Science, Kyoto University, Kitashirakawa Oiwakecho, Sakyoku, Kyoto 606-8502, Japan
| | - Dmitri G Fedorov
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science and Technology (AIST), 1-1-1 Umezono, Tsukuba, Ibaraki 305-8568, Japan
| |
Collapse
|
4
|
Abstract
High-order charge transfer is incorporated into the fragment molecular orbital (FMO) method using a charge transfer state with fractional charges. This state is used for a partition analysis of properties based on segments that may be different from fragments in FMO. The partition analysis is also formulated for calculations without fragmentation. All development in this work is limited to density-functional tight-binding. The analysis is applied to a water cluster, crambin (PDB: 1CBN), and two complexes of Trp-cage (1L2Y) with ligands. The contributions of functional groups in ligands are obtained, providing useful information for drug discovery.
Collapse
Affiliation(s)
- Dmitri G Fedorov
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science and Technology (AIST), Central 2, Umezono 1-1-1, Tsukuba 305-8568, Japan
| |
Collapse
|
5
|
Affiliation(s)
- Dmitri G. Fedorov
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science and Technology (AIST), Central 2, Umezono 1-1-1, Tsukuba 305-8568, Japan
| |
Collapse
|
6
|
Kaliakin DS, Fedorov DG, Alexeev Y, Varganov SA. Locating Minimum Energy Crossings of Different Spin States Using the Fragment Molecular Orbital Method. J Chem Theory Comput 2019; 15:6074-6084. [DOI: 10.1021/acs.jctc.9b00641] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Affiliation(s)
- Danil S. Kaliakin
- Department of Chemistry, University of Nevada, Reno, 1664 N. Virginia Street, Reno, Nevada 89557-0216, United States
| | - Dmitri G. Fedorov
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science and Technology (AIST), Central 2, Umezono 1-1-1, Tsukuba 305-8568, Japan
| | - Yuri Alexeev
- Computational Science Division and Argonne Leadership Computing Facility, Argonne National Laboratory, Argonne, Illinois 60439, United States
| | - Sergey A. Varganov
- Department of Chemistry, University of Nevada, Reno, 1664 N. Virginia Street, Reno, Nevada 89557-0216, United States
| |
Collapse
|
7
|
Abstract
Based on induced solvent charges, a new model of solvent screening is developed in the framework of the fragment molecular orbital combined with the polarizable continuum model. The developed model is applied to analyze interactions in a prototypical zwitterionic system, sodium chloride in water, and it is shown that the large underestimation of the interaction in the original solvent screening based on local charges is successfully corrected. The model is also applied to a complex of the Trp-cage (PDB: 1L2Y ) miniprotein with an anionic ligand, and the physical factors determined protein-ligand binding in solution are unraveled.
Collapse
Affiliation(s)
- Dmitri G Fedorov
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat) , National Institute of Advanced Industrial Science and Technology (AIST) , Central 2, Umezono 1-1-1 , Tsukuba 305-8568 , Japan
| |
Collapse
|
8
|
Fedorov DG, Brekhov A, Mironov V, Alexeev Y. Molecular Electrostatic Potential and Electron Density of Large Systems in Solution Computed with the Fragment Molecular Orbital Method. J Phys Chem A 2019; 123:6281-6290. [DOI: 10.1021/acs.jpca.9b04936] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Dmitri G. Fedorov
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science and Technology (AIST), Central 2, Umezono 1-1-1, Tsukuba, 305-8568, Japan
| | - Anton Brekhov
- Department of Chemistry, Lomonosov Moscow State University, Moscow, 119991, Russian Federation
| | - Vladimir Mironov
- Department of Chemistry, Lomonosov Moscow State University, Moscow, 119991, Russian Federation
| | - Yuri Alexeev
- Argonne Leadership Computing Facility and Computational Science Division, Argonne National Laboratory, Argonne, Illinois, 60439, United States
| |
Collapse
|
9
|
Vuong VQ, Nishimoto Y, Fedorov DG, Sumpter BG, Niehaus TA, Irle S. The Fragment Molecular Orbital Method Based on Long-Range Corrected Density-Functional Tight-Binding. J Chem Theory Comput 2019; 15:3008-3020. [PMID: 30998360 DOI: 10.1021/acs.jctc.9b00108] [Citation(s) in RCA: 27] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
The presently available linear scaling approaches to density-functional tight-binding (DFTB) based on the fragment molecular orbital (FMO) method are severely impacted by the problem of artificial charge transfer due to the self-interaction error (SIE), which hampers the simulation of zwitterionic systems such as biopolymers or ionic liquids. Here we report an extension of FMO-DFTB where we included a long-range corrected (LC) functional designed to mitigate the DFTB SIE, called the FMO-LC-DFTB method, resulting in a robust method which succeeds in simulating zwitterionic systems. Both energy and analytic gradient are developed for the gas phase and the polarizable continuum model of solvation. The scaling of FMO-LC-DFTB with system size N is shown to be almost linear, O( N1.13-1.28), and its numerical accuracy is established for a variety of representative systems including neutral and charged polypeptides. It is shown that pair interaction energies between fragments for two mini-proteins are in excellent agreement with results from long-range corrected density functional theory. The new method was employed in long time scale (1 ns) molecular dynamics simulations of the tryptophan cage protein (PDB: 1L2Y ) in the gas phase for four different protonation states and in stochastic global minimum structure searches for 1-ethyl-3-methylimidazolium nitrate ionic liquid clusters containing up to 2300 atoms.
Collapse
Affiliation(s)
- Van Quan Vuong
- Bredesen Center for Interdisciplinary Research and Graduate Education , University of Tennessee , Knoxville , Tennessee 37996 , United States
| | - Yoshio Nishimoto
- Fukui Institute for Fundamental Chemistry , Kyoto University , Kyoto 606-8501 , Japan
| | - Dmitri G Fedorov
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat) , National Institute of Advanced Industrial Science and Technology (AIST) , Tsukuba 305-8568 , Japan
| | - Bobby G Sumpter
- Center for Nanophase Materials Sciences and Computational Sciences and Engineering Division , Oak Ridge National Laboratory , Oak Ridge , Tennessee 37831 , United States
| | - Thomas A Niehaus
- Univ Lyon, Université Claude Bernard Lyon 1 , CNRS, Institut Lumière Matière , F-69622 Villeurbanne , France
| | - Stephan Irle
- Bredesen Center for Interdisciplinary Research and Graduate Education , University of Tennessee , Knoxville , Tennessee 37996 , United States.,Center for Nanophase Materials Sciences and Computational Sciences and Engineering Division , Oak Ridge National Laboratory , Oak Ridge , Tennessee 37831 , United States.,Chemical Sciences Division , Oak Ridge National Laboratory , Oak Ridge , Tennessee 37831 , United States
| |
Collapse
|
10
|
Nakata H, Bai S. Development of a new parameter optimization scheme for a reactive force field based on a machine learning approach. J Comput Chem 2019. [DOI: 10.1002/jcc.25841] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2018] [Revised: 02/17/2019] [Accepted: 03/27/2019] [Indexed: 12/14/2022]
|
11
|
Malloum A, Fifen JJ, Dhaouadi Z, Nana Engo SG, Conradie J. Structures, relative stability and binding energies of neutral water clusters, (H2O)2–30. NEW J CHEM 2019. [DOI: 10.1039/c9nj01659g] [Citation(s) in RCA: 34] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
We have revised the structures of neutral water clusters, (H2O)n=2–30, with the affordable M06-2X functional, presenting up to 25 isomers for each cluster size.
Collapse
Affiliation(s)
- Alhadji Malloum
- Department of Physics
- Faculty of Science
- University of Ngaoundere
- Ngaoundere
- Cameroon
| | - Jean Jules Fifen
- Department of Physics
- Faculty of Science
- University of Ngaoundere
- Ngaoundere
- Cameroon
| | - Zoubeida Dhaouadi
- Laboratoire de Spectroscopie Atomique Moléculaire et Applications
- Faculté des Sciences de Tunis
- Université de Tunis El Manar
- Tunis
- Tunisia
| | - Serge Guy Nana Engo
- Department of Physics
- Faculty of Science
- University of Ngaoundere
- Ngaoundere
- Cameroon
| | - Jeanet Conradie
- Department of Chemistry
- University of the Free State
- Bloemfontein
- South Africa
| |
Collapse
|
12
|
Nakata H, Fedorov DG. Simulations of infrared and Raman spectra in solution using the fragment molecular orbital method. Phys Chem Chem Phys 2019; 21:13641-13652. [DOI: 10.1039/c9cp00940j] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023]
Abstract
Calculation of IR and Raman spectra in solution for large molecular systems made possible with analytic FMO/PCM Hessians.
Collapse
Affiliation(s)
| | - Dmitri G. Fedorov
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat)
- National Institute of Advanced Industrial Science and Technology (AIST)
- Tsukuba
- Japan
| |
Collapse
|
13
|
Nakata H, Fedorov DG. Analytic second derivatives for the efficient electrostatic embedding in the fragment molecular orbital method. J Comput Chem 2018; 39:2039-2050. [DOI: 10.1002/jcc.25360] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/20/2018] [Revised: 04/27/2018] [Accepted: 04/29/2018] [Indexed: 01/09/2023]
Affiliation(s)
- Hiroya Nakata
- Department of Fundamental Technology Research; Research and Development Center Kagoshima, Kyocera, 1-4 Kokubu Yamashita-cho; Kirishima-shi Kagoshima, 899-4312 Japan
| | - Dmitri G. Fedorov
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science and Technology (AIST), 1-1-1 Umezono; Tsukuba Ibaraki, 305-8568 Japan
| |
Collapse
|
14
|
Nishimoto Y, Fedorov DG. Adaptive frozen orbital treatment for the fragment molecular orbital method combined with density-functional tight-binding. J Chem Phys 2018; 148:064115. [DOI: 10.1063/1.5012935] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022] Open
Affiliation(s)
- Yoshio Nishimoto
- Fukui Institute for Fundamental Chemistry, Kyoto University, 34-4 Takano Nishihiraki-cho, Sakyo-ku, Kyoto 606-8103, Japan
| | - Dmitri G. Fedorov
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science and Technology (AIST), 1-1-1 Umezono, Tsukuba, Ibaraki 305-8568, Japan
| |
Collapse
|
15
|
Fedorov DG, Kitaura K. Pair Interaction Energy Decomposition Analysis for Density Functional Theory and Density-Functional Tight-Binding with an Evaluation of Energy Fluctuations in Molecular Dynamics. J Phys Chem A 2018; 122:1781-1795. [DOI: 10.1021/acs.jpca.7b12000] [Citation(s) in RCA: 37] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Dmitri G. Fedorov
- Research
Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science and Technology (AIST), Central 2, Umezono 1-1-1, Tsukuba 305-8568, Japan
| | - Kazuo Kitaura
- Advanced
Institute for Computational Science (AICS), RIKEN, 7-1-26 Minatojima-Minami-Machi,
Chuo-ku, Kobe, Hyogo 650-0047, Japan
- Fukui
Institute for Fundamental Chemistry, Kyoto University, Takano-Nishihiraki-cho
34-4, Sakyou-ku, Kyoto 606-8103, Japan
| |
Collapse
|
16
|
Dohn AO, Jónsson EÖ, Levi G, Mortensen JJ, Lopez-Acevedo O, Thygesen KS, Jacobsen KW, Ulstrup J, Henriksen NE, Møller KB, Jónsson H. Grid-Based Projector Augmented Wave (GPAW) Implementation of Quantum Mechanics/Molecular Mechanics (QM/MM) Electrostatic Embedding and Application to a Solvated Diplatinum Complex. J Chem Theory Comput 2017; 13:6010-6022. [PMID: 29083921 DOI: 10.1021/acs.jctc.7b00621] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
A multiscale density functional theory-quantum mechanics/molecular mechanics (DFT-QM/MM) scheme is presented, based on an efficient electrostatic coupling between the electronic density obtained from a grid-based projector augmented wave (GPAW) implementation of density functional theory and a classical potential energy function. The scheme is implemented in a general fashion and can be used with various choices for the descriptions of the QM or MM regions. Tests on H2O clusters, ranging from dimer to decamer show that no systematic energy errors are introduced by the coupling that exceeds the differences in the QM and MM descriptions. Over 1 ns of liquid water, Born-Oppenheimer QM/MM molecular dynamics (MD) are sampled combining 10 parallel simulations, showing consistent liquid water structure over the QM/MM border. The method is applied in extensive parallel MD simulations of an aqueous solution of the diplatinum [Pt2(P2O5H2)4]4- complex (PtPOP), spanning a total time period of roughly half a nanosecond. An average Pt-Pt distance deviating only 0.01 Å from experimental results, and a ground-state Pt-Pt oscillation frequency deviating by <2% from experimental results were obtained. The simulations highlight a remarkable harmonicity of the Pt-Pt oscillation, while also showing clear signs of Pt-H hydrogen bonding and directional coordination of water molecules along the Pt-Pt axis of the complex.
Collapse
Affiliation(s)
- A O Dohn
- Faculty of Physical Sciences and Science Institute, University of Iceland , 107 Reykjavı́k, Iceland
| | - E Ö Jónsson
- Faculty of Physical Sciences and Science Institute, University of Iceland , 107 Reykjavı́k, Iceland
| | - G Levi
- Department of Chemistry, Technical University of Denmark , 2800 Kgs. Lyngby, Denmark
| | - J J Mortensen
- CAMD, Department of Physics, Technical University of Denmark , 2800 Kgs. Lyngby, Denmark
| | - O Lopez-Acevedo
- Department of Applied Physics, Aalto University , 02150 Espoo, Finland
| | - K S Thygesen
- CAMD, Department of Physics, Technical University of Denmark , 2800 Kgs. Lyngby, Denmark
| | - K W Jacobsen
- CAMD, Department of Physics, Technical University of Denmark , 2800 Kgs. Lyngby, Denmark
| | - J Ulstrup
- Department of Chemistry, Technical University of Denmark , 2800 Kgs. Lyngby, Denmark
| | - N E Henriksen
- Department of Chemistry, Technical University of Denmark , 2800 Kgs. Lyngby, Denmark
| | - K B Møller
- Department of Chemistry, Technical University of Denmark , 2800 Kgs. Lyngby, Denmark
| | - H Jónsson
- Faculty of Physical Sciences and Science Institute, University of Iceland , 107 Reykjavı́k, Iceland.,Department of Applied Physics, Aalto University , 02150 Espoo, Finland
| |
Collapse
|
17
|
Liu KY, Herbert JM. Understanding the many-body expansion for large systems. III. Critical role of four-body terms, counterpoise corrections, and cutoffs. J Chem Phys 2017; 147:161729. [DOI: 10.1063/1.4986110] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Affiliation(s)
- Kuan-Yu Liu
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - John M. Herbert
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| |
Collapse
|
18
|
Sylvetsky N, Kesharwani MK, Martin JML. The aug-cc-pVnZ-F12 basis set family: Correlation consistent basis sets for explicitly correlated benchmark calculations on anions and noncovalent complexes. J Chem Phys 2017; 147:134106. [DOI: 10.1063/1.4998332] [Citation(s) in RCA: 41] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Nitai Sylvetsky
- Department of Organic Chemistry, Weizmann Institute of Science, 76100 Reḥovot, Israel
| | - Manoj K. Kesharwani
- Department of Organic Chemistry, Weizmann Institute of Science, 76100 Reḥovot, Israel
| | - Jan M. L. Martin
- Department of Organic Chemistry, Weizmann Institute of Science, 76100 Reḥovot, Israel
| |
Collapse
|
19
|
|
20
|
Yuan D, Li Y, Ni Z, Pulay P, Li W, Li S. Benchmark Relative Energies for Large Water Clusters with the Generalized Energy-Based Fragmentation Method. J Chem Theory Comput 2017; 13:2696-2704. [DOI: 10.1021/acs.jctc.7b00284] [Citation(s) in RCA: 29] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Affiliation(s)
- Dandan Yuan
- School
of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic
Chemistry of Ministry of Education, Institute of Theoretical and Computational
Chemistry, Nanjing University, Nanjing 210023, P. R. China
| | - Yunzhi Li
- School
of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic
Chemistry of Ministry of Education, Institute of Theoretical and Computational
Chemistry, Nanjing University, Nanjing 210023, P. R. China
| | - Zhigang Ni
- School
of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic
Chemistry of Ministry of Education, Institute of Theoretical and Computational
Chemistry, Nanjing University, Nanjing 210023, P. R. China
| | - Peter Pulay
- Department
of Chemistry and Biochemistry, University of Arkansas, Fayetteville, Arkansas 72701, United States
| | - Wei Li
- School
of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic
Chemistry of Ministry of Education, Institute of Theoretical and Computational
Chemistry, Nanjing University, Nanjing 210023, P. R. China
| | - Shuhua Li
- School
of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic
Chemistry of Ministry of Education, Institute of Theoretical and Computational
Chemistry, Nanjing University, Nanjing 210023, P. R. China
| |
Collapse
|
21
|
Chen GD, Weng J, Song G, Li ZH. Generalized Switch Functions in the Multilevel Many-Body Expansion Method and Its Application to Water Clusters. J Chem Theory Comput 2017; 13:2010-2020. [DOI: 10.1021/acs.jctc.7b00144] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Affiliation(s)
- Guo Dong Chen
- Collaborative Innovation Center of Chemistry for Energy Material, Shanghai Key Laboratory of Molecular Catalysis & Innovative Materials, Department of Chemistry, Fudan University, Shanghai 200433, China
| | - Jingwei Weng
- Collaborative Innovation Center of Chemistry for Energy Material, Shanghai Key Laboratory of Molecular Catalysis & Innovative Materials, Department of Chemistry, Fudan University, Shanghai 200433, China
| | - Guoliang Song
- Collaborative Innovation Center of Chemistry for Energy Material, Shanghai Key Laboratory of Molecular Catalysis & Innovative Materials, Department of Chemistry, Fudan University, Shanghai 200433, China
| | - Zhen Hua Li
- Collaborative Innovation Center of Chemistry for Energy Material, Shanghai Key Laboratory of Molecular Catalysis & Innovative Materials, Department of Chemistry, Fudan University, Shanghai 200433, China
| |
Collapse
|
22
|
Affiliation(s)
- Jinfeng Liu
- State
Key Laboratory of Natural Medicines, Department of Basic Medicine
and Clinical Pharmacy, China Pharmaceutical University, Nanjing, 210009, China
| | - Lian-Wen Qi
- State
Key Laboratory of Natural Medicines, Department of Basic Medicine
and Clinical Pharmacy, China Pharmaceutical University, Nanjing, 210009, China
| | - John Z. H. Zhang
- School
of Chemistry and Molecular Engineering, East China Normal University, Shanghai, 200062, China
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai, 200062, China
- Department
of Chemistry, New York University, New York, New York 10003, United States
| | - Xiao He
- School
of Chemistry and Molecular Engineering, East China Normal University, Shanghai, 200062, China
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai, 200062, China
| |
Collapse
|
23
|
Abstract
To complement our study of the role of finite precision in electronic structure calculations based on a truncated many-body expansion (MBE, or "n-body expansion"), we examine the accuracy of such methods in the present work. Accuracy may be defined either with respect to a supersystem calculation computed at the same level of theory as the n-body calculations, or alternatively with respect to high-quality benchmarks. Both metrics are considered here. In applications to a sequence of water clusters, (H2O)N=6-55 described at the B3LYP/cc-pVDZ level, we obtain mean absolute errors (MAEs) per H2O monomer of ∼1.0 kcal/mol for two-body expansions, where the benchmark is a B3LYP/cc-pVDZ calculation on the entire cluster. Three- and four-body expansions exhibit MAEs of 0.5 and 0.1 kcal/mol/monomer, respectively, without resort to charge embedding. A generalized many-body expansion truncated at two-body terms [GMBE(2)], using 3-4 H2O molecules per fragment, outperforms all of these methods and affords a MAE of ∼0.02 kcal/mol/monomer, also without charge embedding. GMBE(2) requires significantly fewer (although somewhat larger) subsystem calculations as compared to MBE(4), reducing problems associated with floating-point roundoff errors. When compared to high-quality benchmarks, we find that error cancellation often plays a critical role in the success of MBE(n) calculations, even at the four-body level, as basis-set superposition error can compensate for higher-order polarization interactions. A many-body counterpoise correction is introduced for the GMBE, and its two-body truncation [GMBCP(2)] is found to afford good results without error cancellation. Together with a method such as ωB97X-V/aug-cc-pVTZ that can describe both covalent and non-covalent interactions, the GMBE(2)+GMBCP(2) approach provides an accurate, stable, and tractable approach for large systems.
Collapse
Affiliation(s)
- Ka Un Lao
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Kuan-Yu Liu
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Ryan M Richard
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - John M Herbert
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| |
Collapse
|
24
|
Nakata H, Fedorov DG. Efficient Geometry Optimization of Large Molecular Systems in Solution Using the Fragment Molecular Orbital Method. J Phys Chem A 2016; 120:9794-9804. [DOI: 10.1021/acs.jpca.6b09743] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Hiroya Nakata
- Department
of Fundamental Technology Research, R and D Center Kagoshima, Kyocera, 1-4 Kokubu Yamashita-cho, Kirishima-shi, Kagoshima 899-4312, Japan
| | - Dmitri G. Fedorov
- Research
Center for Computational Design of Advanced Functional Materials, National Institute of Advanced Industrial Science and Technology, 1-1-1
Umezono, Tsukuba, Ibaraki 305-8568, Japan
| |
Collapse
|
25
|
Nakata H, Nishimoto Y, Fedorov DG. Analytic second derivative of the energy for density-functional tight-binding combined with the fragment molecular orbital method. J Chem Phys 2016; 145:044113. [DOI: 10.1063/1.4959231] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022] Open
Affiliation(s)
- Hiroya Nakata
- Department of Fundamental Technology Research, R and D Center Kagoshima, Kyocera, 1-4 Kokubu Yamashita-cho, Kirishima-shi, Kagoshima 899-4312, Japan
| | - Yoshio Nishimoto
- Fukui Institute for Fundamental Chemistry, Kyoto University, 34-4 Takano Nishihiraki-cho, Sakyo-ku, Kyoto 606-8103, Japan
| | - Dmitri G. Fedorov
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science and Technology (AIST), 1-1-1 Umezono, Tsukuba, Ibaraki 305-8568, Japan
| |
Collapse
|
26
|
Cisneros G, Wikfeldt KT, Ojamäe L, Lu J, Xu Y, Torabifard H, Bartók AP, Csányi G, Molinero V, Paesani F. Modeling Molecular Interactions in Water: From Pairwise to Many-Body Potential Energy Functions. Chem Rev 2016; 116:7501-28. [PMID: 27186804 PMCID: PMC5450669 DOI: 10.1021/acs.chemrev.5b00644] [Citation(s) in RCA: 258] [Impact Index Per Article: 32.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2015] [Indexed: 12/17/2022]
Abstract
Almost 50 years have passed from the first computer simulations of water, and a large number of molecular models have been proposed since then to elucidate the unique behavior of water across different phases. In this article, we review the recent progress in the development of analytical potential energy functions that aim at correctly representing many-body effects. Starting from the many-body expansion of the interaction energy, specific focus is on different classes of potential energy functions built upon a hierarchy of approximations and on their ability to accurately reproduce reference data obtained from state-of-the-art electronic structure calculations and experimental measurements. We show that most recent potential energy functions, which include explicit short-range representations of two-body and three-body effects along with a physically correct description of many-body effects at all distances, predict the properties of water from the gas to the condensed phase with unprecedented accuracy, thus opening the door to the long-sought "universal model" capable of describing the behavior of water under different conditions and in different environments.
Collapse
Affiliation(s)
| | - Kjartan Thor Wikfeldt
- Science
Institute, University of Iceland, VR-III, 107, Reykjavik, Iceland
- Department
of Physics, Albanova, Stockholm University, S-106 91 Stockholm, Sweden
| | - Lars Ojamäe
- Department
of Chemistry, Linköping University, SE-581 83 Linköping, Sweden
| | - Jibao Lu
- Department
of Chemistry, The University of Utah, Salt Lake City, Utah 84112-0850, United States
| | - Yao Xu
- Lehrstuhl
Physikalische Chemie II, Ruhr-Universität
Bochum, 44801 Bochum, Germany
| | - Hedieh Torabifard
- Department
of Chemistry, Wayne State University, Detroit, Michigan 48202, United States
| | - Albert P. Bartók
- Engineering
Laboratory, University of Cambridge, Trumpington Street, Cambridge CB21PZ, United Kingdom
| | - Gábor Csányi
- Engineering
Laboratory, University of Cambridge, Trumpington Street, Cambridge CB21PZ, United Kingdom
| | - Valeria Molinero
- Department
of Chemistry, The University of Utah, Salt Lake City, Utah 84112-0850, United States
| | - Francesco Paesani
- Department
of Chemistry and Biochemistry, University
of California San Diego, La Jolla, California 92093, United States
| |
Collapse
|
27
|
Fiedler B, Coriani S, Friedrich J. Molecular Dipole Moments within the Incremental Scheme Using the Domain-Specific Basis-Set Approach. J Chem Theory Comput 2016; 12:3040-52. [PMID: 27300371 DOI: 10.1021/acs.jctc.6b00076] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
We present the first implementation of the fully automated incremental scheme for CCSD unrelaxed dipole moments using the domain-specific basis-set approach. Truncation parameters are varied, and the accuracy of the method is statistically analyzed for a test set of 20 molecules. The local approximations introduce small errors at second order and negligible ones at third order. For a third-order incremental CCSD expansion with a CC2 error correction, a cc-pVDZ/SV domain-specific basis set (tmain = 3.5 Bohr), and the truncation parameter f = 30 Bohr, we obtain a mean error of 0.00 mau (-0.20 mau) and a standard deviation of 1.95 mau (2.17 mau) for the total dipole moments (Cartesian components of the dipole vectors). By analyzing incremental CCSD energies, we demonstrate that the MP2 and CC2 error correction schemes are an exclusive correction for the domain-specific basis-set error. Our implementation of the incremental scheme provides fully automated computations of highly accurate dipole moments at reduced computational cost and is fully parallelized in terms of the calculation of the increments. Therefore, one can utilize the incremental scheme, on the same hardware, to extend the basis set in comparison to standard CCSD and thus obtain a better total accuracy.
Collapse
Affiliation(s)
- Benjamin Fiedler
- Institute for Chemistry, Technische Universität Chemnitz , Straße der Nationen 62, D-09111 Chemnitz, Germany
| | - Sonia Coriani
- Dipartimento di Scienze Chimiche e Farmaceutiche, Università degli Studi di Trieste , Via L. Giorgieri 1, I-34127 Trieste, Italy.,Aarhus Institute of Advanced Studies, Aarhus University , Høegh-Guldbergs Gade 6B, DK-8000 Aarhus C, Denmark
| | - Joachim Friedrich
- Institute for Chemistry, Technische Universität Chemnitz , Straße der Nationen 62, D-09111 Chemnitz, Germany
| |
Collapse
|
28
|
|
29
|
Frey JA, Holzer C, Klopper W, Leutwyler S. Experimental and Theoretical Determination of Dissociation Energies of Dispersion-Dominated Aromatic Molecular Complexes. Chem Rev 2016; 116:5614-41. [DOI: 10.1021/acs.chemrev.5b00652] [Citation(s) in RCA: 53] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Jann A. Frey
- Departement
für Chemie und Biochemie, Universität Bern, Freiestrasse 3, CH-3012 Bern, Switzerland
| | - Christof Holzer
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), Fritz-Haber-Weg 2, D-76131 Karlsruhe, Germany
| | - Wim Klopper
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), Fritz-Haber-Weg 2, D-76131 Karlsruhe, Germany
| | - Samuel Leutwyler
- Departement
für Chemie und Biochemie, Universität Bern, Freiestrasse 3, CH-3012 Bern, Switzerland
| |
Collapse
|
30
|
Fedorov DG, Kitaura K. Subsystem Analysis for the Fragment Molecular Orbital Method and Its Application to Protein-Ligand Binding in Solution. J Phys Chem A 2016; 120:2218-31. [PMID: 26949816 DOI: 10.1021/acs.jpca.6b00163] [Citation(s) in RCA: 47] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
A subsystem analysis is derived incorporating interfragment interactions into the fragment properties, such as energies or charges. The relative stabilities of three alanine isomers, the α-helix, the β-turn, and the extended form are studied and the differences in fragment properties are elucidated. The analysis is further elaborated for studies of binding energies. The binding of the Trp-cage protein (PDB: 1L2Y ) to two ligands is studied in detail. Binding energies defined for each fragment can be used as a convenient descriptor for analyzing contributions to binding in solution.
Collapse
Affiliation(s)
- Dmitri G Fedorov
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science and Technology (AIST) , Central 2, Umezono 1-1-1, Tsukuba, 305-8568, Japan
| | - Kazuo Kitaura
- Graduate School of System Informatics, Kobe University , 1-1 Rokkodai-cho, Nada-ku, Kobe 657-8501, Japan
| |
Collapse
|
31
|
Pruitt SR, Nakata H, Nagata T, Mayes M, Alexeev Y, Fletcher G, Fedorov DG, Kitaura K, Gordon MS. Importance of Three-Body Interactions in Molecular Dynamics Simulations of Water Demonstrated with the Fragment Molecular Orbital Method. J Chem Theory Comput 2016; 12:1423-35. [DOI: 10.1021/acs.jctc.5b01208] [Citation(s) in RCA: 33] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Spencer R. Pruitt
- Argonne
Leadership Computing Facility, Argonne National Laboratory, 9700 S. Cass
Avenue, Lemont, Illinois 60439, United States
| | - Hiroya Nakata
- Department of Fundamental Technology Research, R&D Center Kagoshima, Kyocera Corporation, 1-4 Kokubu Yamashita-cho, Kirishima-shi, Kagoshima 899-4312, Japan
| | - Takeshi Nagata
- Nanosystem Research
Institute, National Institute of Advanced Industrial Science and Technology, 1-1-1 Umenzono, Tsukuba, Ibaraki 305-8568, Japan
| | - Maricris Mayes
- Department
of Chemistry and Biochemistry, University of Massachusetts Dartmouth, 285 Old Westport Road, Dartmouth, Massachusetts 02747-2300, United States
| | - Yuri Alexeev
- Argonne
Leadership Computing Facility, Argonne National Laboratory, 9700 S. Cass
Avenue, Lemont, Illinois 60439, United States
| | - Graham Fletcher
- Argonne
Leadership Computing Facility, Argonne National Laboratory, 9700 S. Cass
Avenue, Lemont, Illinois 60439, United States
| | - Dmitri G. Fedorov
- Nanosystem Research
Institute, National Institute of Advanced Industrial Science and Technology, 1-1-1 Umenzono, Tsukuba, Ibaraki 305-8568, Japan
| | - Kazuo Kitaura
- Graduate
School of System Informatics, Kobe University, 1-1 Rokkodai-cho, Nada-ku, Kobe 657-8501, Japan
| | - Mark S. Gordon
- Department
of Chemistry and Ames Laboratory, Iowa State University, 201 Spedding
Hall, Ames, Iowa 50011, United States
| |
Collapse
|
32
|
Abstract
The GEBF method is demonstrated to be more accurate than the EE-MB method for medium-sized water clusters.
Collapse
Affiliation(s)
- Dandan Yuan
- School of Chemistry and Chemical Engineering
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education
- Institute of Theoretical and Computational Chemistry
- Nanjing University
- Nanjing
| | - Xiaoling Shen
- School of Chemistry and Chemical Engineering
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education
- Institute of Theoretical and Computational Chemistry
- Nanjing University
- Nanjing
| | - Wei Li
- School of Chemistry and Chemical Engineering
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education
- Institute of Theoretical and Computational Chemistry
- Nanjing University
- Nanjing
| | - Shuhua Li
- School of Chemistry and Chemical Engineering
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education
- Institute of Theoretical and Computational Chemistry
- Nanjing University
- Nanjing
| |
Collapse
|
33
|
Nakata H, Fedorov DG, Nagata T, Kitaura K, Nakamura S. Simulations of Chemical Reactions with the Frozen Domain Formulation of the Fragment Molecular Orbital Method. J Chem Theory Comput 2015; 11:3053-64. [DOI: 10.1021/acs.jctc.5b00277] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Hiroya Nakata
- Department
of Biomolecular Engineering, Tokyo Institute of Technology, 4259 Nagatsuta-cho, Midori-ku, Yokohama, Kanagawa 226-8501, Japan
- Research
Cluster for Innovation, RIKEN, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
- Japan Society for the Promotion of Science, Kojimachi Business Center Building, 5-3-1 Kojimachi, Chiyoda-ku, Tokyo 102-0083, Japan
| | - Dmitri G. Fedorov
- Nanosystem
Research Institute, National Institute of Advanced Industrial Science and Technology, 1-1-1 Umezono, Tsukuba, Ibaraki 305-8568, Japan
| | - Takeshi Nagata
- Nanosystem
Research Institute, National Institute of Advanced Industrial Science and Technology, 1-1-1 Umezono, Tsukuba, Ibaraki 305-8568, Japan
- Graduate
School of System Informatics, Kobe University, 1-1 Rokkodai-cho, Nada-ku, Kobe 657-8501, Japan
| | - Kazuo Kitaura
- Graduate
School of System Informatics, Kobe University, 1-1 Rokkodai-cho, Nada-ku, Kobe 657-8501, Japan
| | - Shinichiro Nakamura
- Research
Cluster for Innovation, RIKEN, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
| |
Collapse
|
34
|
Miliordos E, Xantheas SS. An accurate and efficient computational protocol for obtaining the complete basis set limits of the binding energies of water clusters at the MP2 and CCSD(T) levels of theory: Application to (H2O)m, m = 2-6, 8, 11, 16, and 17. J Chem Phys 2015; 142:234303. [DOI: 10.1063/1.4922262] [Citation(s) in RCA: 48] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Evangelos Miliordos
- Physical Sciences Division, Pacific Northwest National Laboratory, 902 Battelle Boulevard, P.O. Box 999, MS K1-83, Richland, Washington 99352, USA
| | - Sotiris S. Xantheas
- Physical Sciences Division, Pacific Northwest National Laboratory, 902 Battelle Boulevard, P.O. Box 999, MS K1-83, Richland, Washington 99352, USA
| |
Collapse
|
35
|
Xu X, Zhang W, Tang M, Truhlar DG. Do Practical Standard Coupled Cluster Calculations Agree Better than Kohn–Sham Calculations with Currently Available Functionals When Compared to the Best Available Experimental Data for Dissociation Energies of Bonds to 3d Transition Metals? J Chem Theory Comput 2015; 11:2036-52. [DOI: 10.1021/acs.jctc.5b00081] [Citation(s) in RCA: 94] [Impact Index Per Article: 10.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Xuefei Xu
- Department
of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Wenjing Zhang
- Department
of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
- College
of Chemistry and Molecular Engineering, Zhengzhou University, Zhengzhou, Henan Province 450001, People’s Republic of China
| | - Mingsheng Tang
- College
of Chemistry and Molecular Engineering, Zhengzhou University, Zhengzhou, Henan Province 450001, People’s Republic of China
| | - Donald G. Truhlar
- Department
of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|
36
|
Affiliation(s)
- Krishnan Raghavachari
- Department of Chemistry, Indiana University, 800 E. Kirkwood Avenue, Bloomington, Indiana 47405, United States
| | - Arjun Saha
- Department of Chemistry, Indiana University, 800 E. Kirkwood Avenue, Bloomington, Indiana 47405, United States
| |
Collapse
|
37
|
Analytic second derivative of the energy for density functional theory based on the three-body fragment molecular orbital method. J Chem Phys 2015; 142:124101. [DOI: 10.1063/1.4915068] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
38
|
Vogiatzis KD, Klopper W, Friedrich J. Non-covalent Interactions of CO2 with Functional Groups of Metal–Organic Frameworks from a CCSD(T) Scheme Applicable to Large Systems. J Chem Theory Comput 2015; 11:1574-84. [DOI: 10.1021/ct5011888] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Affiliation(s)
- Konstantinos D. Vogiatzis
- Department
of Chemistry, University of Minnesota, 207 Pleasant Street Southeast, Minneapolis, Minnesota 55455-0431, United States
- Institute
of Physical Chemistry, Karlsruhe Institute of Technology, Fritz-Haber-Weg
2, D-76131 Karlsruhe, Germany
| | - Wim Klopper
- Institute
of Physical Chemistry, Karlsruhe Institute of Technology, Fritz-Haber-Weg
2, D-76131 Karlsruhe, Germany
| | - Joachim Friedrich
- Institute
of Chemistry, Chemnitz University of Technology, Strasse der Nationen 62, D-09111 Chemnitz, Germany
| |
Collapse
|
39
|
Alipour M. Relative energies of water nanoclusters (H2O)20: comparison of empirical and nonempirical double-hybrids with generalized energy-based fragmentation approach. NEW J CHEM 2015. [DOI: 10.1039/c5nj00817d] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The applicability of recently developed parameterized and parameter-free double-hybrids for predicting the relative energies of water nanoclusters has been examined.
Collapse
Affiliation(s)
- Mojtaba Alipour
- Department of Chemistry
- College of Sciences
- Shiraz University
- Shiraz
- Iran
| |
Collapse
|
40
|
Wang B, Yang KR, Xu X, Isegawa M, Leverentz HR, Truhlar DG. Quantum mechanical fragment methods based on partitioning atoms or partitioning coordinates. Acc Chem Res 2014; 47:2731-8. [PMID: 24841937 DOI: 10.1021/ar500068a] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
Conspectus The development of more efficient and more accurate ways to represent reactive potential energy surfaces is a requirement for extending the simulation of large systems to more complex systems, longer-time dynamical processes, and more complete statistical mechanical sampling. One way to treat large systems is by direct dynamics fragment methods. Another way is by fitting system-specific analytic potential energy functions with methods adapted to large systems. Here we consider both approaches. First we consider three fragment methods that allow a given monomer to appear in more than one fragment. The first two approaches are the electrostatically embedded many-body (EE-MB) expansion and the electrostatically embedded many-body expansion of the correlation energy (EE-MB-CE), which we have shown to yield quite accurate results even when one restricts the calculations to include only electrostatically embedded dimers. The third fragment method is the electrostatically embedded molecular tailoring approach (EE-MTA), which is more flexible than EE-MB and EE-MB-CE. We show that electrostatic embedding greatly improves the accuracy of these approaches compared with the original unembedded approaches. Quantum mechanical fragment methods share with combined quantum mechanical/molecular mechanical (QM/MM) methods the need to treat a quantum mechanical fragment in the presence of the rest of the system, which is especially challenging for those parts of the rest of the system that are close to the boundary of the quantum mechanical fragment. This is a delicate matter even for fragments that are not covalently bonded to the rest of the system, but it becomes even more difficult when the boundary of the quantum mechanical fragment cuts a bond. We have developed a suite of methods for more realistically treating interactions across such boundaries. These methods include redistributing and balancing the external partial atomic charges and the use of tuned fluorine atoms for capping dangling bonds, and we have shown that they can greatly improve the accuracy. Finally we present a new approach that goes beyond QM/MM by combining the convenience of molecular mechanics with the accuracy of fitting a potential function to electronic structure calculations on a specific system. To make the latter practical for systems with a large number of degrees of freedom, we developed a method to interpolate between local internal-coordinate fits to the potential energy. A key issue for the application to large systems is that rather than assigning the atoms or monomers to fragments, we assign the internal coordinates to reaction, secondary, and tertiary sets. Thus, we make a partition in coordinate space rather than atom space. Fits to the local dependence of the potential energy on tertiary coordinates are arrayed along a preselected reaction coordinate at a sequence of geometries called anchor points; the potential energy function is called an anchor points reactive potential. Electrostatically embedded fragment methods and the anchor points reactive potential, because they are based on treating an entire system by quantum mechanical electronic structure methods but are affordable for large and complex systems, have the potential to open new areas for accurate simulations where combined QM/MM methods are inadequate.
Collapse
Affiliation(s)
- Bo Wang
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, 207 Pleasant St. SE, Minneapolis, Minnesota 55455-0431, United States
| | - Ke R. Yang
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, 207 Pleasant St. SE, Minneapolis, Minnesota 55455-0431, United States
| | - Xuefei Xu
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, 207 Pleasant St. SE, Minneapolis, Minnesota 55455-0431, United States
| | - Miho Isegawa
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, 207 Pleasant St. SE, Minneapolis, Minnesota 55455-0431, United States
| | - Hannah R. Leverentz
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, 207 Pleasant St. SE, Minneapolis, Minnesota 55455-0431, United States
| | - Donald G. Truhlar
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, 207 Pleasant St. SE, Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|
41
|
Affiliation(s)
- Ryan M. Richard
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Ka Un Lao
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - John M. Herbert
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| |
Collapse
|
42
|
Fiedler L, Leverentz HR, Nachimuthu S, Friedrich J, Truhlar DG. Nitrogen and Sulfur Compounds in Atmospheric Aerosols: A New Parametrization of Polarized Molecular Orbital Model Chemistry and Its Validation against Converged CCSD(T) Calculations for Large Clusters. J Chem Theory Comput 2014; 10:3129-39. [DOI: 10.1021/ct5003169] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Luke Fiedler
- Department
of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, United States
| | - Hannah R. Leverentz
- Department
of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, United States
| | - Santhanamoorthi Nachimuthu
- Department
of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, United States
| | - Joachim Friedrich
- Institute
for Chemistry, Chemnitz University of Technology, Straße der Nationen 62, 09111, Chemnitz, Germany
| | - Donald G. Truhlar
- Department
of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, United States
| |
Collapse
|
43
|
Marenich AV, Ho J, Coote ML, Cramer CJ, Truhlar DG. Computational electrochemistry: prediction of liquid-phase reduction potentials. Phys Chem Chem Phys 2014; 16:15068-106. [PMID: 24958074 DOI: 10.1039/c4cp01572j] [Citation(s) in RCA: 310] [Impact Index Per Article: 31.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
This article reviews recent developments and applications in the area of computational electrochemistry. Our focus is on predicting the reduction potentials of electron transfer and other electrochemical reactions and half-reactions in both aqueous and nonaqueous solutions. Topics covered include various computational protocols that combine quantum mechanical electronic structure methods (such as density functional theory) with implicit-solvent models, explicit-solvent protocols that employ Monte Carlo or molecular dynamics simulations (for example, Car-Parrinello molecular dynamics using the grand canonical ensemble formalism), and the Marcus theory of electronic charge transfer. We also review computational approaches based on empirical relationships between molecular and electronic structure and electron transfer reactivity. The scope of the implicit-solvent protocols is emphasized, and the present status of the theory and future directions are outlined.
Collapse
Affiliation(s)
- Aleksandr V Marenich
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, 207 Pleasant Street S.E., Minneapolis, MN 55455-0431, USA.
| | | | | | | | | |
Collapse
|
44
|
Wang K, Li W, Li S. Generalized Energy-Based Fragmentation CCSD(T)-F12a Method and Application to the Relative Energies of Water Clusters (H2O)20. J Chem Theory Comput 2014; 10:1546-53. [DOI: 10.1021/ct401060m] [Citation(s) in RCA: 54] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Kedong Wang
- School
of Physics and Electronical Engineering, Henan Normal University, Xinxiang, 453007, People’s Repubic of China
| | | | | |
Collapse
|