51
|
Generalized Energy-Based Fragmentation Approach for the Accurate Binding Energies and Raman Spectra of Methane Hydrate Clusters. CHINESE J CHEM PHYS 2022. [DOI: 10.1063/1674-0068/cjcp2111256] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022]
|
52
|
Bandyopadhyay P, Priya P, Sadhukhan M. A simple fragment-based method for van der Waals corrections over density functional theory. Phys Chem Chem Phys 2022; 24:8508-8518. [DOI: 10.1039/d2cp00744d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Modeling intermolecular noncovalent interactions between large molecules remains a challenge for the electron structure theory community due to the high cost. Fragment-based methods usually fare well in reducing the cost...
Collapse
|
53
|
Harshan AK, Bronson MJ, Jensen L. Local-Field Effects in Linear Response Properties within a Polarizable Frozen Density Embedding Method. J Chem Theory Comput 2021; 18:380-393. [PMID: 34905917 DOI: 10.1021/acs.jctc.1c00816] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
In this work, we present a polarizable frozen density embedding (FDE) method for calculating polarizabilities of coupled subsystems. The method (FDE-pol) combines a FDE method with an explicit polarization model such that the expensive freeze/thaw cycles can be bypassed, and approximate nonadditive kinetic potentials are avoided by enforcing external orthogonality between the subsystems. To describe the polarization of the frozen environment, we introduce a Hirshfeld partition-based density-dependent method for calculating the atomic polarizabilities of atoms in molecules, which alleviates the need to fit the atomic parameters to a specific system of interest or to a larger general set of molecules. We show that the Hirshfeld partition-based method predicts molecular polarizabilities close to the basis set limit, and thus, a single basis set-dependent scaling parameter can be introduced to improve the agreement against the reference polarizability data. To test the model, we characterized the uncoupled and coupled response of small interacting molecular complexes. Here, the coupled response properties include the perturbation of the frozen system due to the external perturbation which is ignored in the uncoupled response. We show that FDE-pol can accurately reproduce both the exact uncoupled polarizability and the coupled polarizabilities of the supermolecular systems. Using damped response theory, we also demonstrate that the coupled frequency-dependent polarizability can be described by including local field effects. The results emphasize the necessity of including local-field effects for describing the response properties of coupled subsystems, as well as the importance of accurate atomic polarizability models.
Collapse
Affiliation(s)
- Aparna K Harshan
- Department of Chemistry, The Pennsylvania State University, 104 Chemistry Building, University Park 16802, United States
| | - Mark J Bronson
- Department of Chemistry, The Pennsylvania State University, 104 Chemistry Building, University Park 16802, United States
| | - Lasse Jensen
- Department of Chemistry, The Pennsylvania State University, 104 Chemistry Building, University Park 16802, United States
| |
Collapse
|
54
|
Li W, Ma H, Li S, Ma J. Computational and data driven molecular material design assisted by low scaling quantum mechanics calculations and machine learning. Chem Sci 2021; 12:14987-15006. [PMID: 34909141 PMCID: PMC8612375 DOI: 10.1039/d1sc02574k] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/10/2021] [Accepted: 10/12/2021] [Indexed: 12/11/2022] Open
Abstract
Electronic structure methods based on quantum mechanics (QM) are widely employed in the computational predictions of the molecular properties and optoelectronic properties of molecular materials. The computational costs of these QM methods, ranging from density functional theory (DFT) or time-dependent DFT (TDDFT) to wave-function theory (WFT), usually increase sharply with the system size, causing the curse of dimensionality and hindering the QM calculations for large sized systems such as long polymer oligomers and complex molecular aggregates. In such cases, in recent years low scaling QM methods and machine learning (ML) techniques have been adopted to reduce the computational costs and thus assist computational and data driven molecular material design. In this review, we illustrated low scaling ground-state and excited-state QM approaches and their applications to long oligomers, self-assembled supramolecular complexes, stimuli-responsive materials, mechanically interlocked molecules, and excited state processes in molecular aggregates. Variable electrostatic parameters were also introduced in the modified force fields with the polarization model. On the basis of QM computational or experimental datasets, several ML algorithms, including explainable models, deep learning, and on-line learning methods, have been employed to predict the molecular energies, forces, electronic structure properties, and optical or electrical properties of materials. It can be conceived that low scaling algorithms with periodic boundary conditions are expected to be further applicable to functional materials, perhaps in combination with machine learning to fast predict the lattice energy, crystal structures, and spectroscopic properties of periodic functional materials.
Collapse
Affiliation(s)
- Wei Li
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute of Theoretical and Computational Chemistry, School of Chemistry and Chemical Engineering, Nanjing University Nanjing 210023 China
| | - Haibo Ma
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute of Theoretical and Computational Chemistry, School of Chemistry and Chemical Engineering, Nanjing University Nanjing 210023 China
- Jiangsu Key Laboratory of Advanced Organic Materials, Jiangsu Key Laboratory of Vehicle Emissions Control, Nanjing University Nanjing 210023 China
| | - Shuhua Li
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute of Theoretical and Computational Chemistry, School of Chemistry and Chemical Engineering, Nanjing University Nanjing 210023 China
| | - Jing Ma
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute of Theoretical and Computational Chemistry, School of Chemistry and Chemical Engineering, Nanjing University Nanjing 210023 China
- Jiangsu Key Laboratory of Advanced Organic Materials, Jiangsu Key Laboratory of Vehicle Emissions Control, Nanjing University Nanjing 210023 China
| |
Collapse
|
55
|
Sandler I, Sharma S, Chan B, Ho J. Accurate Quantum Chemical Prediction of Gas-Phase Anion Binding Affinities and Their Structure-Binding Relationships. J Phys Chem A 2021; 125:9838-9851. [PMID: 34739245 DOI: 10.1021/acs.jpca.1c06648] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
This paper systematically examines the performance of contemporary wavefunction and density functional theory methods to identify robust and cost-efficient methods for predicting gas-phase anion binding energies. This includes the local coupled cluster LNO-CCSD(T) and DLPNO-CCSD(T), as well as double-hybrid DSD-PBEP86-D3(BJ) and various hybrid functionals M06-2X, B3LYP-D3(BJ), ωB97M-V, and ωB97X-V. The focus is on dual-hydrogen-bonding anion receptors that are commonly found in supramolecular chemistry and organocatalysis, namely, (thio)ureas, deltamides, (thio)squaramides, and croconamides as well as the yet-to-be-explored rhodizonamides. Of the methods examined, M06-2X emerged as the overall best performing method as the other functionals including DSD-PBEP86-D3(BJ) and the local coupled cluster DLPNO-CCSD(T) method displayed systematic errors that increase with the degree of carbonylation of the receptors. Hybrid ONIOM models that employed semiempirical methods (PM7, GFN1-xTB, and GFN2-xTB) and "threefold"-corrected small-basis set potentials (HF-3c, B97-3c, and PBEh-3c) were explored, and the best models resulted in 50- to 500-fold reduction in CPU time compared to W1-local. These calculations provide important insight into the structure-binding relationships where there is a direct correlation between Brønsted acidity and anion binding affinity, though the strength of the correlation also depends on other factors such as hydrogen-bonding geometry and the geometrical distortion that the receptor needs to undergo to bind the anion.
Collapse
Affiliation(s)
- Isolde Sandler
- School of Chemistry, The University of New South Wales, Sydney, New South Wales 2052, Australia
| | - Shaleen Sharma
- School of Chemistry, The University of New South Wales, Sydney, New South Wales 2052, Australia
| | - Bun Chan
- Graduate School of Engineering, Nagasaki University, Bukyo-Machi 1-14, Nagasaki 852-8521, Japan
| | - Junming Ho
- School of Chemistry, The University of New South Wales, Sydney, New South Wales 2052, Australia
| |
Collapse
|
56
|
Kumar A, DeGregorio N, Iyengar SS. Graph-Theory-Based Molecular Fragmentation for Efficient and Accurate Potential Surface Calculations in Multiple Dimensions. J Chem Theory Comput 2021; 17:6671-6690. [PMID: 34623129 DOI: 10.1021/acs.jctc.1c00065] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
We present a multitopology molecular fragmentation approach, based on graph theory, to calculate multidimensional potential energy surfaces in agreement with post-Hartree-Fock levels of theory but at the density functional theory cost. A molecular assembly is coarse-grained into a set of graph-theoretic nodes that are then connected with edges to represent a collection of locally interacting subsystems up to an arbitrary order. Each of the subsystems is treated at two levels of electronic structure theory, the result being used to construct many-body expansions that are embedded within an ONIOM scheme. These expansions converge rapidly with the many-body order (or graphical rank) of subsystems and capture many-body interactions accurately and efficiently. However, multiple graphs, and hence multiple fragmentation topologies, may be defined in molecular configuration space that may arise during conformational sampling or from reactive, bond breaking and bond formation, events. Obtaining the resultant potential surfaces is an exponential scaling proposition, given the number of electronic structure computations needed. We utilize a family of graph-theoretic representations within a variational scheme to obtain multidimensional potential surfaces at a reduced cost. The fast convergence of the graph-theoretic expansion with increasing order of many-body interactions alleviates the exponential scaling cost for computing potential surfaces, with the need to only use molecular fragments that contain a fewer number of quantum nuclear degrees of freedom compared to the full system. This is because the dimensionality of the conformational space sampled by the fragment subsystems is much smaller than the full molecular configurational space. Additionally, we also introduce a multidimensional clustering algorithm, based on physically defined criteria, to reduce the number of energy calculations by orders of magnitude. The molecular systems benchmarked include coupled proton motion in protonated water wires. The potential energy surfaces and multidimensional nuclear eigenstates obtained are shown to be in very good agreement with those from explicit post-Hartree-Fock calculations that become prohibitive as the number of quantum nuclear dimensions grows. The developments here provide a rigorous and efficient alternative to this important chemical physics problem.
Collapse
Affiliation(s)
- Anup Kumar
- Department of Chemistry and Department of Physics, Indiana University, 800 E. Kirkwood Avenue, Bloomington, Indiana 47405, United States
| | - Nicole DeGregorio
- Department of Chemistry and Department of Physics, Indiana University, 800 E. Kirkwood Avenue, Bloomington, Indiana 47405, United States
| | - Srinivasan S Iyengar
- Department of Chemistry and Department of Physics, Indiana University, 800 E. Kirkwood Avenue, Bloomington, Indiana 47405, United States
| |
Collapse
|
57
|
Hellmers J, König C. A unified and flexible formulation of molecular fragmentation schemes. J Chem Phys 2021; 155:164105. [PMID: 34717347 DOI: 10.1063/5.0059598] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022] Open
Abstract
We present a flexible formulation for energy-based molecular fragmentation schemes. This framework does not only incorporate the majority of existing fragmentation expansions but also allows for flexible formulation of novel schemes. We further illustrate its application in multi-level approaches and for electronic interaction energies. For the examples of small water clusters, a small protein, and protein-protein interaction energies, we show how this flexible setup can be exploited to generate a well-suited multi-level fragmentation expansion for the given case. With such a setup, we reproduce the electronic protein-protein interaction energy of ten different structures of a neurotensin and an extracellular loop of its receptor with a mean absolute deviation to the respective super-system calculations below 1 kJ/mol.
Collapse
Affiliation(s)
- Janine Hellmers
- Institute of Physical Chemistry and Electrochemistry, Leibniz University Hannover, Hannover, Germany
| | - Carolin König
- Institute of Physical Chemistry and Electrochemistry, Leibniz University Hannover, Hannover, Germany
| |
Collapse
|
58
|
Westheimer BM, Gordon MS. Scalable ab initio fragmentation methods based on a truncated expansion of the non-orthogonal molecular orbital model. J Chem Phys 2021; 155:154101. [PMID: 34686043 DOI: 10.1063/5.0064864] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
An alternative formulation of the non-orthogonal molecular orbital model of electronic structure theory is developed based on the expansion of the inverse molecular orbital overlap matrix. From this model, a hierarchy of ab initio fragment-based quantum chemistry methods, referred to as the nth-order expanded non-orthogonal molecular orbital methods, are developed using a minimal number of approximations, each of which is frequently employed in intermolecular interaction theory. These novel methods are compared to existing fragment-based quantum chemistry methods, and the implications of those significant differences, where they exist, between the methods developed herein and those already existing methods are examined in detail. Computational complexities and theoretical scaling are also analyzed and discussed. Future extensions for the hierarchy of methods, to account for additional intrafragment and interfragment interactions, are outlined.
Collapse
Affiliation(s)
| | - Mark S Gordon
- Department of Chemistry, Iowa State University, Ames, Iowa 50011, USA
| |
Collapse
|
59
|
Sun Z, Liu Z. BAR‐Based Multi‐Dimensional Nonequilibrium Pulling for Indirect Construction of QM/MM Free Energy Landscapes: Varying the QM Region. ADVANCED THEORY AND SIMULATIONS 2021. [DOI: 10.1002/adts.202100185] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Affiliation(s)
- Zhaoxi Sun
- Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering Peking University Beijing 100871 China
| | - Zhirong Liu
- Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering Peking University Beijing 100871 China
| |
Collapse
|
60
|
Xu S, Wang QD, Sun MM, Yin G, Liang J. Benchmark calculations for bond dissociation energies and enthalpy of formation of chlorinated and brominated polycyclic aromatic hydrocarbons. RSC Adv 2021; 11:29690-29701. [PMID: 35479574 PMCID: PMC9040899 DOI: 10.1039/d1ra05391d] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2021] [Accepted: 08/31/2021] [Indexed: 01/22/2023] Open
Abstract
Thermodynamic properties, i.e., bond dissociation energies and enthalpy of formation, of chlorinated and brominated polycyclic aromatic hydrocarbons play a fundamental role in understanding their formation mechanisms and reactivity. Computational electronic structure calculations routinely used to predict thermodynamic properties of various species are limited for these compounds due to large computational cost to obtain accurate results by employing high-level wave function theory methods. In this work, a number of composite model chemistry methods (CBS-QB3, G3MP2, G3, and G4) are used to compute bond dissociation energies and enthalpies of formation of small to medium-size chlorinated and brominated polycyclic aromatic hydrocarbon compounds. The enthalpy of formation is derived via the atomization method and compared against the recommended values. Statistical analysis indicates that G4 is the best method. For comparison, three commonly used density functional theory (DFT) methods (M06-2X, ωB97X-D and B2PLYP-D3) with various basis sets including 6-311++G(d, p), cc-pVTZ, and cc-pVQZ in the prediction of bond dissociation energies and enthalpies of formation have been tested using the optimized geometries at the same M06-2X/6-311++G(d, p) level of theory. It is found that ωB97X-D/6-311++G(d, p) shows the best performance in computing the bond dissociation energies, while ωB97X-D/cc-pVTZ exhibits the best prediction in enthalpy of formation of the studied reaction systems. The structural effect on the bond dissociation energies and enthalpy of formation of chlorinated and brominated polycyclic aromatic hydrocarbons are then systematically analyzed. Based on comparisons of the various methods, reliable DFT methods are recommended for future theoretical studies on large chlorinated and brominated polycyclic aromatic hydrocarbons considering both accuracy and computational cost. This work, to the authors' knowledge, is the first to systematically benchmark theoretical methods for the accurate prediction of thermodynamic properties for chlorinated and brominated polycyclic aromatic hydrocarbons. Benchmark calculations using state-of-the-art DFT functionals and composite methods for bond dissociation energy and enthalpy of formation of halogenated polycyclic aromatic hydrocarbons are performed.![]()
Collapse
Affiliation(s)
- Shenying Xu
- Faculty of Materials and Chemical Engineering, Yibin University Yibin Sichuan 644000 People's Republic of China
| | - Quan-De Wang
- Faculty of Materials and Chemical Engineering, Yibin University Yibin Sichuan 644000 People's Republic of China .,Low Carbon Energy Institute and School of Chemical Engineering, China University of Mining and Technology Xuzhou 221008 People's Republic of China
| | - Mao-Mao Sun
- Low Carbon Energy Institute and School of Chemical Engineering, China University of Mining and Technology Xuzhou 221008 People's Republic of China
| | - Guoliang Yin
- Faculty of Materials and Chemical Engineering, Yibin University Yibin Sichuan 644000 People's Republic of China
| | - Jinhu Liang
- Faculty of Materials and Chemical Engineering, Yibin University Yibin Sichuan 644000 People's Republic of China .,School of Environment and Safety Engineering, North University of China Taiyuan 030051 People's Republic of China
| |
Collapse
|
61
|
Chan B, Karton A. Polycyclic aromatic hydrocarbons: from small molecules through nano-sized species towards bulk graphene. Phys Chem Chem Phys 2021; 23:17713-17723. [PMID: 34378574 DOI: 10.1039/d1cp01659h] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We have examined the use of systematic bond-separation reactions and purposely constructed chemistry-preserving isodesmic reactions for the thermochemical calculation of aromatic hydrocarbon species. The bond-separation approach yields somewhat disappointing accuracy even when the reaction energies are obtained with generally robust composite and double-hybrid (DH) density functional theory (DFT) methods. In contrast, for the purposely constructed reactions, we find a dramatic improvement in the accuracy for energies calculated with all methods examined. Notably, for medium-sized aromatic hydrocarbons, we find that an effective approach for formulating a well-balanced reaction is to split the target species into two halves with an aromatic overlapping region. Overall, the G4(MP2)-XK, MPW2PLYP, MN15, PBE, and DC-DFTB3 methods are reasonable within their respective classes of methods for the calculation of bond-separation as well as chemistry-preserving isodesmic reactions. We have further computed per-carbon atomization energy (AE) for a series of D6h benzene-type molecules, and thus obtained a formula for extrapolation to the graphene limit [AEn = 711.5 × (1 - 1/n0.640) kJ mol-1, where n = number of carbons]. It suggests that nano-graphene with a length larger than 10 nm would resemble properties of bulk graphene, and conversely, downsizing a nano-graphene beyond this point may lead to considerably altered properties from the bulk.
Collapse
Affiliation(s)
- Bun Chan
- Graduate School of Engineering, Nagasaki University, Bunkyo 1-14, Nagasaki-Shi, Nagasaki 852-8521, Japan.
| | | |
Collapse
|
62
|
Chen J, Kato J, Harper JB, Shao Y, Ho J. On the Accuracy of QM/MM Models: A Systematic Study of Intramolecular Proton Transfer Reactions of Amino Acids in Water. J Phys Chem B 2021; 125:9304-9316. [PMID: 34355564 DOI: 10.1021/acs.jpcb.1c04876] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
This work presents a systematic assessment of QM/QM' and QM/MM models with respect to direct QM calculations for the tautomerization (neutral to zwitterion) reactions of amino acids (glycine, alanine, valine, aspartate, and neutral and protonated histidine) solvated in a 160 water cluster. The effect of varying QM region size and choice of embedding potentials, including fixed-charge and polarizable molecular mechanics force fields (TIP3P and EFP) and various semiempirical QM methods (PM7, GFN2-xTB, DFTBA, DFTB3, HF-3c, and PBEh-3c), on the accuracy of the models was examined. A surprising finding was that molecular mechanics force fields outperformed many of the semiempirical methods. Generally, the errors in the QM/QM' and QM/MM models converge slowly with respect to the QM region size, requiring 50 or more waters to be included in the QM region before the error in the model falls below 1 kcal mol-1 of its pure QM result. Different QM region selection schemes were also compared, and it was found that selection based on Natural Population Analysis (NPA) atomic charges significantly reduced the error in the QM/QM' and QM/MM models particularly if a low-quality embedding potential was used. It is envisaged that these results will be useful for the development of future hybrid QM models.
Collapse
Affiliation(s)
- Junbo Chen
- School of Chemistry, The University of New South Wales, Sydney, NSW 2052, Australia
| | - Jin Kato
- School of Chemistry, The University of New South Wales, Sydney, NSW 2052, Australia
| | - Jason B Harper
- School of Chemistry, The University of New South Wales, Sydney, NSW 2052, Australia
| | - Yihan Shao
- Department of Chemistry and Biochemistry, University of Oklahoma, Norman, Oklahoma 73019, United States
| | - Junming Ho
- School of Chemistry, The University of New South Wales, Sydney, NSW 2052, Australia
| |
Collapse
|
63
|
Collins EM, Raghavachari K. A Fragmentation-Based Graph Embedding Framework for QM/ML. J Phys Chem A 2021; 125:6872-6880. [PMID: 34342449 DOI: 10.1021/acs.jpca.1c06152] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
We introduce a new fragmentation-based molecular representation framework "FragGraph" for QM/ML methods involving embedding fragment-wise fingerprints onto molecular graphs. Our model is specifically designed for delta machine learning (Δ-ML) with the central goal of correcting the deficiencies of approximate methods such as DFT to achieve high accuracy. Our framework is based on a judicious combination of ideas from fragmentation, error cancellation, and a state-of-the-art deep learning architecture. Broadly, we develop a general graph-network framework for molecular machine learning by incorporating the inherent advantages prebuilt into error cancellation methods such as the generalized Connectivity-Based Hierarchy. More specifically, we develop a QM/ML representation through a fragmentation-based attributed graph representation encoded with fragment-wise molecular fingerprints. The utility of our representation is demonstrated through a graph network fingerprint encoder in which a global fingerprint is generated through message passing of local neighborhoods of fragment-wise fingerprints, effectively augmenting standard fingerprints to also include the inbuilt molecular graph structure. On the 130k-GDB9 dataset, our method predicts an out-of-sample mean absolute error significantly lower than 1 kJ/mol compared to target G4(MP2) calculated energies, rivaling current deep learning methods with reduced computational scaling.
Collapse
Affiliation(s)
- Eric M Collins
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| | - Krishnan Raghavachari
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| |
Collapse
|
64
|
Wang Z, Liu W. iOI: An Iterative Orbital Interaction Approach for Solving the Self-Consistent Field Problem. J Chem Theory Comput 2021; 17:4831-4845. [PMID: 34240856 DOI: 10.1021/acs.jctc.1c00445] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
An iterative orbital interaction (iOI) approach is proposed to solve, in a bottom-up fashion, the self-consistent field problem in quantum chemistry. While it belongs grossly to the family of fragment-based quantum chemical methods, iOI is distinctive in that (1) it divides and conquers not only the energy but also the wave function and that (2) the subsystem sizes are automatically determined by successively merging neighboring small subsystems until they are just enough for converging the wave function to a given accuracy. Orthonormal occupied and virtual localized molecular orbitals are obtained in a natural manner, which can be used for all post-SCF purposes.
Collapse
Affiliation(s)
- Zikuan Wang
- Qingdao Institute for Theoretical and Computational Sciences, Institute of Frontier and Interdisciplinary Science, Shandong University, Qingdao, Shandong 266237, P. R. China
| | - Wenjian Liu
- Qingdao Institute for Theoretical and Computational Sciences, Institute of Frontier and Interdisciplinary Science, Shandong University, Qingdao, Shandong 266237, P. R. China
| |
Collapse
|
65
|
Verma P, Huntington L, Coons MP, Kawashima Y, Yamazaki T, Zaribafiyan A. Scaling up electronic structure calculations on quantum computers: The frozen natural orbital based method of increments. J Chem Phys 2021; 155:034110. [DOI: 10.1063/5.0054647] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023] Open
Affiliation(s)
- Prakash Verma
- 1QB Information Technologies (1QBit), 200-1285 W Pender St., Vancouver, British Columbia V6E 4B1, Canada
| | - Lee Huntington
- 1QB Information Technologies (1QBit), 200-1285 W Pender St., Vancouver, British Columbia V6E 4B1, Canada
| | - Marc P. Coons
- Dow, Core R&D, Chemical Science, 1776 Building, Midland, Michigan 48674, USA
| | - Yukio Kawashima
- 1QB Information Technologies (1QBit), 200-1285 W Pender St., Vancouver, British Columbia V6E 4B1, Canada
| | - Takeshi Yamazaki
- 1QB Information Technologies (1QBit), 200-1285 W Pender St., Vancouver, British Columbia V6E 4B1, Canada
| | - Arman Zaribafiyan
- 1QB Information Technologies (1QBit), 200-1285 W Pender St., Vancouver, British Columbia V6E 4B1, Canada
| |
Collapse
|
66
|
Schmitt-Monreal D, Jacob CR. Density-Based Many-Body Expansion as an Efficient and Accurate Quantum-Chemical Fragmentation Method: Application to Water Clusters. J Chem Theory Comput 2021; 17:4144-4156. [PMID: 34196558 DOI: 10.1021/acs.jctc.1c00340] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Fragmentation methods based on the many-body expansion offer an attractive approach for the quantum-chemical treatment of large molecular systems, such as molecular clusters and crystals. Conventionally, the many-body expansion is performed for the total energy, but such an energy-based many-body expansion often suffers from a slow convergence with respect to the expansion order. For systems that show strong polarization effects such as water clusters, this can render the energy-based many-body expansion infeasible. Here, we establish a density-based many-body expansion as a promising alternative approach. By performing the many-body expansion for the electron density instead of the total energy and inserting the resulting total electron density into the total energy functional of density functional theory, one can derive a density-based energy correction, which in principle accounts for all higher-order polarization effects. Here, we systematically assess the accuracy of such a density-based many-body expansion for test sets of water clusters. We show that already a density-based two-body expansion is able to reproduce interaction energies per fragment within chemical accuracy and is able to accurately predict the energetic ordering as well as the relative interaction energies of different isomers of water clusters.
Collapse
Affiliation(s)
- Daniel Schmitt-Monreal
- Institute of Physical and Theoretical Chemistry, Technische Universität Braunschweig, Gaußstr. 17, 38106 Braunschweig, Germany
| | - Christoph R Jacob
- Institute of Physical and Theoretical Chemistry, Technische Universität Braunschweig, Gaußstr. 17, 38106 Braunschweig, Germany
| |
Collapse
|
67
|
Yan Z, Li X, Chung LW. Multiscale Quantum Refinement Approaches for Metalloproteins. J Chem Theory Comput 2021; 17:3783-3796. [PMID: 34032440 DOI: 10.1021/acs.jctc.1c00148] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Abstract
Biomolecules with metal ion(s) (e.g., metalloproteins) play many important biological roles. However, accurate structural determination of metalloproteins, particularly those containing transition metal ion(s), is challenging due to their complicated electronic structure, complex bonding of metal ions, and high number of conformations in biomolecules. Quantum refinement, which was proposed to combine crystallographic data with computational chemistry methods by several groups, can improve the local structures of some proteins. In this study, a quantum refinement method combining several multiscale computational schemes with experimental (X-ray diffraction) information was developed for metalloproteins. Various quantum refinement approaches using different ONIOM (our own N-layered integrated molecular orbital and molecular mechanics) combinations of quantum mechanics (QM), semiempirical (SE), and molecular mechanics (MM) methods were conducted to assess the performance and reliability on the refined local structure in two metalloproteins. The structures for two (Cu- or Zn-containing) metalloproteins were refined by combining two-layer ONIOM2(QM1/QM2) and ONIOM2(QM/MM) and three-layer ONIOM3(QM1/QM2/MM) schemes with experimental data. The accuracy of the quantum-refined metal binding sites was also examined and compared in these multiscale quantum refinement calculations. ONIOM3(QM/SE/MM) schemes were found to give good results with lower computational costs and were proposed to be a good choice for the multiscale computational scheme for quantum refinement calculations of metal binding site(s) in metalloproteins with high efficiency. Additionally, a two-center ONIOM approach was employed to speed up the quantum refinement calculations for the Zn metalloprotein with two remote active sites/ligands. Moreover, a recent quantum-embedding wavefunction-in-density functional theory (WF-in-DFT) method was also adopted as the high-level method in unprecedented ONIOM2(CCSD-in-B3LYP/MM) and ONIOM3(CCSD-in-B3LYP/SE/MM) calculations, which can be regarded as novel pseudo-three- and pseudo-four-layer ONIOM methods, respectively, to refine the key Zn binding site at the coupled-cluster singles and doubles (CCSD) level. These refined results indicate that multiscale quantum refinement schemes can be used to improve the structural accuracy obtained for local metal binding site(s) in metalloproteins with high efficiency.
Collapse
Affiliation(s)
- Zeyin Yan
- Shenzhen Grubbs Institute, Department of Chemistry and Guangdong Provincial Key Laboratory of Catalysis, Southern University of Science and Technology, Shenzhen 518055, China
| | - Xin Li
- Shenzhen Grubbs Institute, Department of Chemistry and Guangdong Provincial Key Laboratory of Catalysis, Southern University of Science and Technology, Shenzhen 518055, China
| | - Lung Wa Chung
- Shenzhen Grubbs Institute, Department of Chemistry and Guangdong Provincial Key Laboratory of Catalysis, Southern University of Science and Technology, Shenzhen 518055, China
| |
Collapse
|
68
|
Zhang JH, Ricard TC, Haycraft C, Iyengar SS. Weighted-Graph-Theoretic Methods for Many-Body Corrections within ONIOM: Smooth AIMD and the Role of High-Order Many-Body Terms. J Chem Theory Comput 2021; 17:2672-2690. [PMID: 33891416 DOI: 10.1021/acs.jctc.0c01287] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
We present a weighted-graph-theoretic approach to adaptively compute contributions from many-body approximations for smooth and accurate post-Hartree-Fock (pHF) ab initio molecular dynamics (AIMD) of highly fluxional chemical systems. This approach is ONIOM-like, where the full system is treated at a computationally feasible quality of treatment (density functional theory (DFT) for the size of systems considered in this publication), which is then improved through a perturbative correction that captures local many-body interactions up to a certain order within a higher level of theory (post-Hartree-Fock in this publication) described through graph-theoretic techniques. Due to the fluxional and dynamical nature of the systems studied here, these graphical representations evolve during dynamics. As a result, energetic "hops" appear as the graphical representation deforms with the evolution of the chemical and physical properties of the system. In this paper, we introduce dynamically weighted, linear combinations of graphs, where the transition between graphical representations is smoothly achieved by considering a range of neighboring graphical representations at a given instant during dynamics. We compare these trajectories with those obtained from a set of trajectories where the range of local many-body interactions considered is increased, sometimes to the maximum available limit, which yields conservative trajectories as the order of interactions is increased. The weighted-graph approach presents improved dynamics trajectories while only using lower-order many-body interaction terms. The methods are compared by computing dynamical properties through time-correlation functions and structural distribution functions. In all cases, the weighted-graph approach provides accurate results at a lower cost.
Collapse
Affiliation(s)
- Juncheng Harry Zhang
- Department of Chemistry and Department of Physics, Indiana University, 800 East Kirkwood Avenue, Bloomington, Indiana 47405, United States
| | - Timothy C Ricard
- Department of Chemistry and Department of Physics, Indiana University, 800 East Kirkwood Avenue, Bloomington, Indiana 47405, United States
| | - Cody Haycraft
- Department of Chemistry and Department of Physics, Indiana University, 800 East Kirkwood Avenue, Bloomington, Indiana 47405, United States
| | - Srinivasan S Iyengar
- Department of Chemistry and Department of Physics, Indiana University, 800 East Kirkwood Avenue, Bloomington, Indiana 47405, United States
| |
Collapse
|
69
|
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] [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
|
70
|
Liu J, Liu Y, Yang J, Zeng XC, He X. Directional Proton Transfer in the Reaction of the Simplest Criegee Intermediate with Water Involving the Formation of Transient H 3O . J Phys Chem Lett 2021; 12:3379-3386. [PMID: 33784110 DOI: 10.1021/acs.jpclett.1c00448] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
The reaction of Criegee intermediates with water vapor has been widely known as a key Criegee reaction in the troposphere. Herein, we investigated the reaction of the smallest Criegee intermediate, CH2OO, with a water cluster through fragment-based ab initio molecular dynamics simulations at the MP2/aug-cc-pVDZ level. Our results show that the CH2OO-water reaction could occur not only at the air/water interface but also inside the water cluster. Moreover, more than one reactive water molecules are required for the CH2OO-water reaction, which is always initiated from the Criegee carbon atom and ends at the terminal Criegee oxygen atom via a directional proton transfer process. The observed reaction pathways include the loop-structure-mediated and stepwise mechanisms, and the latter involves the formation of transient H3O+. The lifetime of transient H3O+ is on the order of a few picoseconds, which may impact the atmospheric budget of the other trace gases in the actual atmosphere.
Collapse
Affiliation(s)
- Jinfeng Liu
- Department of Basic Medicine and Clinical Pharmacy, China Pharmaceutical University, Nanjing, 210009, China
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai, 200062, China
| | - Yanqing Liu
- Department of Basic Medicine and Clinical Pharmacy, China Pharmaceutical University, Nanjing, 210009, China
| | - Jinrong Yang
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai, 200062, China
| | - Xiao Cheng Zeng
- Department of Chemistry, University of Nebraska, Lincoln, Nebraska 68588, United States
| | - Xiao He
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, 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
|
71
|
Szabó PB, Csóka J, Kállay M, Nagy PR. Linear-Scaling Open-Shell MP2 Approach: Algorithm, Benchmarks, and Large-Scale Applications. J Chem Theory Comput 2021; 17:2886-2905. [PMID: 33819030 PMCID: PMC8154337 DOI: 10.1021/acs.jctc.1c00093] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
A linear-scaling
local second-order Møller–Plesset
(MP2) method is presented for high-spin open-shell molecules based
on restricted open-shell (RO) reference functions. The open-shell
local MP2 (LMP2) approach inherits the iteration- and redundancy-free
formulation and the completely integral-direct, OpenMP-parallel, and
memory and disk use economic algorithms of our closed-shell LMP2 implementation.
By utilizing restricted local molecular orbitals for the demanding
integral transformation step and by introducing a novel long-range
spin-polarization approximation, the computational cost of RO-LMP2
approaches that of closed-shell LMP2. Extensive benchmarks were performed
for reactions of radicals, ionization potentials, as well as spin-state
splittings of carbenes and transition-metal complexes. Compared to
the conventional MP2 reference for systems of up to 175 atoms, local
errors of at most 0.1 kcal/mol were found, which are well below the
intrinsic accuracy of MP2. RO-LMP2 computations are presented for
challenging protein models of up to 601 atoms and 11 000 basis
functions, which involve either spin states of a complexed iron ion
or a highly delocalized singly occupied orbital. The corresponding
runtimes of 9–15 h obtained with a single, many-core CPU demonstrate
that MP2, as well as spin-scaled MP2 and double-hybrid density functional
methods, become widely accessible for open-shell systems of unprecedented
size and complexity.
Collapse
Affiliation(s)
- P Bernát Szabó
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - József Csóka
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - Mihály Kállay
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - Péter R Nagy
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| |
Collapse
|
72
|
Uratani H, Yoshikawa T, Nakai H. Trajectory Surface Hopping Approach to Condensed-Phase Nonradiative Relaxation Dynamics Using Divide-and-Conquer Spin-Flip Time-Dependent Density-Functional Tight Binding. J Chem Theory Comput 2021; 17:1290-1300. [PMID: 33577323 DOI: 10.1021/acs.jctc.0c01155] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Nonradiative relaxation of excited molecules is central to many crucial issues in photochemistry. Condensed phases are typical contexts in which such problems are considered, and the nonradiative relaxation dynamics are expected to be significantly affected by interactions with the environment, for example, a solvent. We developed a nonadiabatic molecular dynamics simulation technique that can treat the nonradiative relaxation and explicitly include the environment in the calculations without a heavy computational burden. Specifically, we combined trajectory surface hopping with Tully's fewest-switches algorithm, a tight-binding approximated version of spin-flip time-dependent density-functional theory, and divide-and-conquer (DC) spatial fragmentation scheme. Numerical results showed that this method can treat systems with thousands of atoms within reasonable computational resources, and the error arising from DC fragmentation is negligibly small. Using this method, we obtained molecular insights into the solvent dependence of the photoexcited-state dynamics of trans-azobenzene, which demonstrate the importance of the environment for condensed-phase nonradiative relaxation.
Collapse
Affiliation(s)
- Hiroki Uratani
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan
| | - Takeshi Yoshikawa
- Faculty of Pharmaceutical Sciences, Toho University, 2-2-1 Miyama, Funabashi, Chiba 274-8510, Japan.,Waseda Research Institute for Science and Engineering (WISE), 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan
| | - Hiromi Nakai
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan.,Waseda Research Institute for Science and Engineering (WISE), 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan.,Elements Strategy Initiative for Catalysts and Batteries (ESICB), Kyoto University, Katsura, Kyoto 615-8245, Japan
| |
Collapse
|
73
|
Wolter M, von Looz M, Meyerhenke H, Jacob CR. Systematic Partitioning of Proteins for Quantum-Chemical Fragmentation Methods Using Graph Algorithms. J Chem Theory Comput 2021; 17:1355-1367. [PMID: 33591754 DOI: 10.1021/acs.jctc.0c01054] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Quantum-chemical fragmentation methods offer an efficient approach for the treatment of large proteins, in particular if local target quantities such as protein-ligand interaction energies, enzymatic reaction energies, or spectroscopic properties of embedded chromophores are sought. However, the accuracy that is achievable for such local target quantities intricately depends on how the protein is partitioned into smaller fragments. While the commonly employed naı̈ve approach of using fragments with a fixed size is widely used, it can result in large and unpredictable errors when varying the fragment size. Here, we present a systematic partitioning scheme that aims at minimizing the fragmentation error of a local target quantity for a given maximum fragment size. To this end, we construct a weighted graph representation of the protein, in which the amino acids constitute the nodes. These nodes are connected by edges weighted with an estimate for the fragmentation error that is expected when cutting this edge. This allows us to employ graph partitioning algorithms provided by computer science to determine near-optimal partitions of the protein. We apply this scheme to a test set of six proteins representing various prototypical applications of quantum-chemical fragmentation methods using a simplified molecular fractionation with conjugate caps (MFCC) approach with hydrogen caps. We show that our graph-based scheme consistently improves upon the naı̈ve approach.
Collapse
Affiliation(s)
- Mario Wolter
- Institute of Physical and Theoretical Chemistry, Technische Universität Braunschweig, Gaußstrasse 17, 38106 Braunschweig, Germany
| | - Moritz von Looz
- Department of Computer Science, Humboldt-Universität zu Berlin, Unter den Linden 6, 10099 Berlin, Germany
| | - Henning Meyerhenke
- Department of Computer Science, Humboldt-Universität zu Berlin, Unter den Linden 6, 10099 Berlin, Germany
| | - Christoph R Jacob
- Institute of Physical and Theoretical Chemistry, Technische Universität Braunschweig, Gaußstrasse 17, 38106 Braunschweig, Germany
| |
Collapse
|
74
|
Sandler I, Chen J, Taylor M, Sharma S, Ho J. Accuracy of DLPNO-CCSD(T): Effect of Basis Set and System Size. J Phys Chem A 2021; 125:1553-1563. [PMID: 33560853 DOI: 10.1021/acs.jpca.0c11270] [Citation(s) in RCA: 43] [Impact Index Per Article: 14.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
The DLPNO-CCSD(T) method is designed to study large molecular systems at significantly reduced cost relative to its canonical counterpart. However, the error in this approach is also size-extensive and relies on cancellation of errors for the calculation of relative energies. This work provides a direct comparison of canonical CCSD(T) and TightPNO DLPNO-CCSD(T) calculations of reaction energies and barriers of a broad range of chemical reactions. The dataset includes acidities, anion binding affinities, enolization, Diels-Alder, nucleophilic substitution, and atom transfer reactions and complements existing theoretical datasets in terms of system size as well as new reaction types (e.g., anion binding affinities and chlorine atom transfer reactions). The performance of DLPNO-CCSD(T) was further examined with respect to systematic variation of basis set and system size and amounts of nonbonded interaction present in the system. The errors in the DLPNO-CCSD(T) were found to be relatively insensitive to the choice of basis set for small systems but increase monotonically with system size. Additionally, calculations of barriers appear to be more challenging than reaction energies with errors exceeding 5 kJ mol-1 for many Diels-Alder reactions. Further tests on three realistic organic reactions reveal the impact of the DLPNO approximation in calculating absolute and relative barriers that are important for predictions such as stereoselectivity.
Collapse
Affiliation(s)
- Isolde Sandler
- School of Chemistry, University of New South Wales, Sydney, NSW 2052, Australia
| | - Junbo Chen
- School of Chemistry, University of New South Wales, Sydney, NSW 2052, Australia
| | - Mackenzie Taylor
- School of Chemistry, University of New South Wales, Sydney, NSW 2052, Australia
| | - Shaleen Sharma
- School of Chemistry, University of New South Wales, Sydney, NSW 2052, Australia
| | - Junming Ho
- School of Chemistry, University of New South Wales, Sydney, NSW 2052, Australia
| |
Collapse
|
75
|
Abstract
Computational methods for modeling biochemical processes implemented in GAMESS package are reviewed; in particular, quantum mechanics combined with molecular mechanics (QM/MM), semi-empirical, and fragmentation approaches. A detailed summary of capabilities is provided for the QM/MM implementation in QuanPol program and the fragment molecular orbital (FMO) method. Molecular modeling and visualization packages useful for biochemical simulations with GAMESS are described. GAMESS capabilities with corresponding references are tabulated for reader's convenience.
Collapse
|
76
|
Chandy SK, Thapa B, Raghavachari K. Accurate and cost-effective NMR chemical shift predictions for proteins using a molecules-in-molecules fragmentation-based method. Phys Chem Chem Phys 2020; 22:27781-27799. [PMID: 33244526 DOI: 10.1039/d0cp05064d] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We have developed an efficient protocol using our two-layer Molecules-in-Molecules (MIM2) fragmentation-based quantum chemical method for the prediction of NMR chemical shifts of large biomolecules. To investigate the performance of our fragmentation approach and demonstrate its applicability, MIM-NMR calculations are first calibrated on a test set of six proteins. The MIM2-NMR method yields a mean absolute deviation (MAD) from unfragmented full molecule calculations of 0.01 ppm for 1H and 0.06 ppm for 13C chemical shifts. Thus, the errors from fragmentation are only about 3% of our target accuracy of ∼0.3 ppm for 1H and 2-3 ppm for 13C chemical shifts. To compare with experimental chemical shifts, a standard protocol is first derived using two smaller proteins 2LHY (176 atoms) and 2LI1 (146 atoms) for obtaining an appropriate protein structure for NMR chemical shift calculations. The effect of the solvent environment on the calculated NMR chemical shifts is incorporated through implicit, explicit, or explicit-implicit solvation models. The expensive first solvation shell calculations are replaced by a micro-solvation model in which only the immediate interaction between the protein and the explicit solvation environment is considered. A single explicit water molecule for each amine and amide proton is found to be sufficient to yield accurate results for 1H chemical shifts. The 1H and 13C NMR chemical shifts calculated using our protocol give excellent agreement with experiments for two larger proteins, 2MC5 (the helical part with 265 atoms) and 3UMK (33 residue slice with 547 atoms). Overall, our target accuracy of ∼0.3 ppm for 1H and ∼2-3 ppm for 13C has been achieved for the larger proteins. The proposed MIM-NMR method is accurate and computationally cost-effective and should be applicable to study a wide range of large proteins.
Collapse
Affiliation(s)
- Sruthy K Chandy
- Department of Chemistry, Indiana University, Bloomington, Indiana, USA.
| | | | | |
Collapse
|
77
|
Bao JL, Welch BK, Ulusoy IS, Zhang X, Xu X, Wilson AK, Truhlar DG. Predicting Bond Dissociation Energies and Bond Lengths of Coordinatively Unsaturated Vanadium-Ligand Bonds. J Phys Chem A 2020; 124:9757-9770. [PMID: 33180508 DOI: 10.1021/acs.jpca.0c06519] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Understanding the electronic structure of coordinatively unsaturated transition-metal compounds and predicting their physical properties are of great importance for catalyst design. Bond dissociation energy De and bond length re are two of the fundamental quantities for which good predictions are important for a successful design strategy. In the present work, recent experimentally measured bond energies and bond lengths of VX diatomic molecules (X = C, N, S) are used as a gauge to consider the utility of a number of electronic structure methods. Single-reference methods are one focus because of their efficiency and utility in practical calculations, and multireference configuration interaction (MRCISD) methods and a composite coupled cluster (CCC) method are a second focus because of their potential high accuracy. The comparison is especially challenging because of the large multireference M diagnostics of these molecules, in the range 0.15-0.19. For the single-reference methods, Kohn-Sham density functional theory (KS-DFT) has been tested with a variety of approximate exchange-correlation functionals. Of these, MOHLYP provides the bond dissociation energies in best agreement with experiments, and BLYP provides the bond lengths that are in best agreement with experiments; but by requiring good performance for both the De and re of the vanadium compounds, MOHLYP, MN12-L, MGGA_MS1, MGGA_MS0, O3LYP, and M06-L are the most highly recommended functionals. The CCC calculations include up to connected pentuple excitations for the valence electrons and up to connected quadruple excitations for the core-valence terms; this results in highly accurate dissociation energies and good bond lengths. Averaged over the three molecules, the mean unsigned deviation of CCC bond energies from experimental ones is only 0.4 kcal/mol, demonstrating excellent convergence of theory and experiments.
Collapse
Affiliation(s)
- Junwei Lucas Bao
- Department of Chemistry, Boston College, Chestnut Hill, Massachusetts 02467, United States
| | - Bradley K Welch
- Department of Chemistry, Michigan State University, 578 S Shaw Lane, East Lansing, Michigan 48824-1322, United States
| | - Inga S Ulusoy
- Department of Chemistry, Michigan State University, 578 S Shaw Lane, East Lansing, Michigan 48824-1322, United States.,Theoretical Chemistry, Institute of Physical Chemistry, Heidelberg University, Im Neuenheimer Feld 229, Heidelberg 69120, Germany
| | - Xin Zhang
- State Key Laboratory of Chemical Resource Engineering, Beijing Advanced Innovation Center for Soft Matter Science and Engineering, Beijing University of Chemical Technology, Beijing 100029, People's Republic of China.,Department of Chemistry, Chemical Theory Center, Inorganometallic Catalyst Design Center, and Minnesota Supercomputing Institute, University of Minnesota, 207 Pleasant Street SE, Minneapolis, Minnesota 55455-0431, United States
| | - Xuefei Xu
- Center for Combustion Energy, Department of Energy and Power Engineering, and Key Laboratory for Thermal Science and Power Engineering of Ministry of Education, Tsinghua University, Beijing 100084, China
| | - Angela K Wilson
- Department of Chemistry, Michigan State University, 578 S Shaw Lane, East Lansing, Michigan 48824-1322, United States
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, Inorganometallic Catalyst Design Center, and Minnesota Supercomputing Institute, University of Minnesota, 207 Pleasant Street SE, Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|
78
|
Tribedi S, Chakraborty A, Maitra R. Formulation of a Dressed Coupled-Cluster Method with Implicit Triple Excitations and Benchmark Application to Hydrogen-Bonded Systems. J Chem Theory Comput 2020; 16:6317-6328. [PMID: 32794747 DOI: 10.1021/acs.jctc.0c00736] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
In this paper, we present a coupled-cluster theory based on a double-exponential wave operator ansatz, which is capable of mimicking the effects of connected triple excitations in an iterative manner. The triply excited manifold is spanned via the action of a set of scattering operators on doubly excited determinants, whereas their action annihilates the Hartree-Fock reference determinant. The effect of triple excitations is included at a computational scaling slightly higher than that of conventional coupled-cluster singles and doubles. Furthermore, we demonstrate two approximate schemes, which arise naturally, and argue that both these schemes come equipped with certain renormalization terms capable of handling nonbonding interactions due to robust inclusion of the screened Coulomb interaction. We justify our claims from both a theoretical perspective and a number of numerical applications to prototypical water clusters, in a number of basis functions. Our methods show overall comparable performance to the canonical coupled-cluster theory with singles, doubles, and perturbative triples (CCSD(T)) and allied methods, however, at a lower computational scaling.
Collapse
Affiliation(s)
- Soumi Tribedi
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| | - Anish Chakraborty
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| | - Rahul Maitra
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| |
Collapse
|
79
|
Chan B, Collins E, Raghavachari K. Applications of isodesmic‐type reactions for computational thermochemistry. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2020. [DOI: 10.1002/wcms.1501] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Affiliation(s)
- Bun Chan
- Graduate School of Engineering Nagasaki University Nagasaki Japan
| | - Eric Collins
- Department of Chemistry Indiana University Bloomington Indiana USA
| | | |
Collapse
|
80
|
Sun Z. SAMPL7 TrimerTrip host-guest binding poses and binding affinities from spherical-coordinates-biased simulations. J Comput Aided Mol Des 2020; 35:105-115. [PMID: 32776199 DOI: 10.1007/s10822-020-00335-9] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/04/2020] [Accepted: 08/04/2020] [Indexed: 12/21/2022]
Abstract
Host-guest binding remains a major challenge in modern computational modelling. The newest 7th statistical assessment of the modeling of proteins and ligands (SAMPL) challenge contains a new series of host-guest systems. The TrimerTrip host binds to 16 structurally diverse guests. Previously, we have successfully employed the spherical coordinates as the collective variables coupled with the enhanced sampling technique metadynamics to enhance the sampling of the binding/unbinding event, search for possible binding poses and calculate the binding affinities in all three host-guest binding cases of the 6th SAMPL challenge. In this work, we report a retrospective study on the TrimerTrip host-guest systems by employing the same protocol to investigate the TrimerTrip host in the SAMPL7 challenge. As no binding pose is provided by the SAMPL7 host, our simulations initiate from randomly selected configurations and are proceeded long enough to obtain converged free energy estimates and search for possible binding poses. The calculated binding affinities are in good agreement with the experimental reference, and the obtained binding poses serve as a nice starting point for end-point or alchemical free energy calculations. Note that as the work is performed after the close of the SAMPL7 challenge, we do not participate in the challenge and the results are not formally submitted to the SAMPL7 challenge.
Collapse
Affiliation(s)
- Zhaoxi Sun
- State Key Laboratory of Precision Spectroscopy, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai, 200062, China.
| |
Collapse
|
81
|
Unzueta PA, Beran GJO. Polarizable continuum models provide an effective electrostatic embedding model for fragment-based chemical shift prediction in challenging systems. J Comput Chem 2020; 41:2251-2265. [PMID: 32748418 DOI: 10.1002/jcc.26388] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/12/2020] [Revised: 06/04/2020] [Accepted: 07/04/2020] [Indexed: 12/25/2022]
Abstract
Ab initio nuclear magnetic resonance chemical shift prediction provides an important tool for interpreting and assigning experimental spectra, but it becomes computationally prohibitive in large systems. The computational costs can be reduced considerably by fragmentation of the large system into a series of contributions from many smaller subsystems. However, the presence of charged functional groups and the need to partition the system across covalent bonds create complications in biomolecules that typically require the use of large fragments and careful descriptions of the electrostatic environment. The present work shows how a model that combines chemical shielding contributions from non-overlapping monomer and dimer fragments embedded in a polarizable continuum model provides a simple, easy-to-implement, and computationally inexpensive approach for predicting chemical shifts in complex systems. The model's performance proves rather insensitive to the continuum dielectric constant, making the selection of the optimal embedding dielectric less critical. The PCM-embedded fragment model is demonstrated to perform well across systems ranging from molecular crystals to proteins.
Collapse
Affiliation(s)
- Pablo A Unzueta
- Department of Chemistry, Univeristy of California, Riverside, California, USA
| | - Gregory J O Beran
- Department of Chemistry, Univeristy of California, Riverside, California, USA
| |
Collapse
|
82
|
Chen B, Xu X. XO-PBC: An Accurate and Efficient Method for Molecular Crystals. J Chem Theory Comput 2020; 16:4271-4285. [PMID: 32456429 DOI: 10.1021/acs.jctc.0c00232] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
In this work, we propose the XO-PBC method, which combines the eXtended ONIOM method (XO) with the periodic boundary condition (PBC) for the description of molecular crystals. XO-PBC tries to embed a finite cluster cut out from the solid into the periodic environment, making it feasible to employ advanced molecular quantum chemistry methods, which are usually prohibitively expensive for direct PBC calculations. In particular, XO-PBC utilizes the results from force calculations to design the scheme to fragment the molecule when crystals are made of large molecules and to select cluster model systems automatically consisting of dimer up to tetramer interactions for embedding. By applying an appropriate theory to each model, a satisfactory accuracy for the system under study is ensured, while a high efficiency is achieved with massively parallel computing by distributing model systems onto different processors. A comparison of the XO-PBC calculations with the conventional direct PBC calculations at the B3LYP level demonstrates its accuracy at substantially low cost for the description of molecular crystals. The usefulness of the XO-PBC method is further exemplified, showing that XO-PBC is able to predict the lattice energies of various types of molecular crystals within chemical accuracy (<4 kJ/mol) when the doubly hybrid density functional XYG3 is used as the target high level and the periodic PBE as the basic low level. The XO-PBC method provides a general protocol that brings the great predictive power of advanced electronic structure methods from molecular systems to the extended solids.
Collapse
Affiliation(s)
- Bozhu Chen
- Collaborative Innovation Center of Chemistry for Energy Materials, Shanghai Key Laboratory of Molecular Catalysis and Innovative Materials, Ministry of Education Key Laboratory of Computational Physical Sciences, Department of Chemistry, Fudan University, Shanghai, 200433, China
| | - Xin Xu
- Collaborative Innovation Center of Chemistry for Energy Materials, Shanghai Key Laboratory of Molecular Catalysis and Innovative Materials, Ministry of Education Key Laboratory of Computational Physical Sciences, Department of Chemistry, Fudan University, Shanghai, 200433, China
| |
Collapse
|
83
|
Ricard TC, Iyengar SS. Efficient and Accurate Approach To Estimate Hybrid Functional and Large Basis-Set Contributions to Condensed-Phase Systems and Molecule–Surface Interactions. J Chem Theory Comput 2020; 16:4790-4812. [DOI: 10.1021/acs.jctc.9b01089] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Affiliation(s)
- Timothy C. Ricard
- Department of Chemistry and Department of Physics, Indiana University, 800 E. Kirkwood Avenue, Bloomington, Indiana 47405, United States
| | - Srinivasan S. Iyengar
- Department of Chemistry and Department of Physics, Indiana University, 800 E. Kirkwood Avenue, Bloomington, Indiana 47405, United States
| |
Collapse
|
84
|
Uratani H, Nakai H. Non-adiabatic molecular dynamics with divide-and-conquer type large-scale excited-state calculations. J Chem Phys 2020; 152:224109. [DOI: 10.1063/5.0006831] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Affiliation(s)
- Hiroki Uratani
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1, Okubo, Shinjuku-ku, Tokyo 169-8555, Japan
| | - Hiromi Nakai
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1, Okubo, Shinjuku-ku, Tokyo 169-8555, Japan
- Waseda Research Institute for Science and Engineering (WISE), 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan
- Elements Strategy Initiative for Catalysts and Batteries (ESICB), Kyoto University, Katsura, Kyoto 615-8245, Japan
| |
Collapse
|
85
|
Chen WK, Zhang Y, Jiang B, Fang WH, Cui G. Efficient Construction of Excited-State Hessian Matrices with Machine Learning Accelerated Multilayer Energy-Based Fragment Method. J Phys Chem A 2020; 124:5684-5695. [DOI: 10.1021/acs.jpca.0c04117] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Affiliation(s)
- Wen-Kai Chen
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Yaolong Zhang
- Hefei National Laboratory for Physical Science at the Microscale, Department of Chemical Physics, Key Laboratory of Surface and Interface Chemistry and Energy Catalysis of Anhui Higher Education Institutes, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Bin Jiang
- Hefei National Laboratory for Physical Science at the Microscale, Department of Chemical Physics, Key Laboratory of Surface and Interface Chemistry and Energy Catalysis of Anhui Higher Education Institutes, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Wei-Hai Fang
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Ganglong Cui
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| |
Collapse
|
86
|
Asakereh Z, Zare M, Shakerzadeh E. Exploring enthalpies of formation of imidazolium-, pyridinium-, and pyrrolidinium-based ionic liquids with dicyanamide anion using quantum chemical methods. J Mol Liq 2020. [DOI: 10.1016/j.molliq.2020.113137] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
|
87
|
Liu J, He X. Fragment-based quantum mechanical approach to biomolecules, molecular clusters, molecular crystals and liquids. Phys Chem Chem Phys 2020; 22:12341-12367. [PMID: 32459230 DOI: 10.1039/d0cp01095b] [Citation(s) in RCA: 33] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
To study large molecular systems beyond the system size that the current state-of-the-art ab initio electronic structure methods could handle, fragment-based quantum mechanical (QM) approaches have been developed over the past years, and proved to be efficient in dealing with large molecular systems at various ab initio levels. According to the fragmentation approach, a large molecular system can be divided into subsystems (fragments), and subsequently the property of the whole system can be approximately obtained by taking a proper combination of the corresponding terms of individual fragments. Therefore, the standard QM calculation of a large system could be circumvented by carrying out a series of calculations on small fragments, which significantly promotes computational efficiency. The electrostatically embedded generalized molecular fractionation with conjugate caps (EE-GMFCC) method is one of the fragment-based QM approaches which has been developed by our research group in recent years. This Perspective presents the theoretical framework of this fragmentation method and its applications in biomolecules, molecular clusters, molecular crystals and liquids, including total energy calculation, protein-ligand/protein binding affinity prediction, geometry optimization, vibrational spectrum simulation, ab initio molecular dynamics simulation, and prediction of excited-state properties.
Collapse
Affiliation(s)
- Jinfeng Liu
- Department of Basic Medicine and Clinical Pharmacy, China Pharmaceutical University, Nanjing, 210009, China
| | | |
Collapse
|
88
|
The Impact of Electron Correlation on Describing QM/MM Interactions in the Attendant Molecular Dynamics Simulations of CO in Myoglobin. Sci Rep 2020; 10:8539. [PMID: 32444817 PMCID: PMC7244521 DOI: 10.1038/s41598-020-65475-2] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/22/2019] [Accepted: 05/05/2020] [Indexed: 01/10/2023] Open
Abstract
The impact of the dispersion and electron correlation effects on describing quantum mechanics/molecular mechanics (QM/MM) interactions in QM/MM molecular dynamics (MD) simulations was explored by performing a series of up to 2 ns QM/MM MD simulations on the B states of the myoglobin-carbon monoxide (MbCO) system. The results indicate that both dispersion and electron correlations play significant roles in the simulation of the ratios of two B states (B1/B2), which suggests that the inclusion of the electron correlation effects is essential for accurately modeling the interactions between QM and MM subsystems. We found that the QM/MM interaction energies between the CO and the surroundings statistically present a linear correlation with the electric fields along the CO bond. This indicates that QM/MM interactions can be described by a simple physical model of a dipole with constant moment under the action of the electric fields. The treatment provides us with an accurate and effective approach to account for the electron correlation effects in QM/MM MD simulations.
Collapse
|
89
|
Chen X, Gao J. Fragment Exchange Potential for Realizing Pauli Deformation of Interfragment Interactions. J Phys Chem Lett 2020; 11:4008-4016. [PMID: 32308000 DOI: 10.1021/acs.jpclett.0c00933] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
In fragment-based methods, the lack of explicit short-range exchange interactions between monomers can result in unphysical deformation in charge density. In this study, we describe a fragment exchange potential (XFP) to explicitly account for interfragmental Pauli deformation. In our implementation, a Kohn-Sham exchange potential is adopted along with the Yukawa potential. The method has been validated by comparison of the computed exchange energies using the XFP potential with results obtained from antisymmetrized fragmental orbitals on the S66×8 data set containing 528 bimolecular interactions of equilibrium and arbitrary geometries. It was also found that it is only necessary to deploy numerical grids on atoms within their van der Waals contacts, significantly reducing the small, albeit extra, computational cost. We anticipate that the XFP presented here may be applied to molecular dynamics simulations of macromolecules using a fragment-based quantum mechanical potential with improved SCF convergence and computational accuracy.
Collapse
Affiliation(s)
- Xin Chen
- Laboratory of Theoretical and Computational Chemistry, Institute of Theoretical Chemistry, College of Chemistry, Jilin University, Changchun 130023, P. R. China
| | - Jiali Gao
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen, Guangdong, China
- Department of Chemistry and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, United States
| |
Collapse
|
90
|
Artiukhin DG, Klinting EL, König C, Christiansen O. Adaptive density-guided approach to double incremental potential energy surface construction. J Chem Phys 2020; 152:194105. [PMID: 33687258 DOI: 10.1063/5.0004686] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022] Open
Abstract
We present a combination of the recently developed double incremental expansion of potential energy surfaces with the well-established adaptive density-guided approach to grid construction. This unique methodology is based on the use of an incremental expansion for potential energy surfaces, known as n-mode expansion; an incremental many-body representation of the electronic energy; and an efficient vibrational density-guided approach to automated determination of grid dimensions and granularity. The reliability of the method is validated calculating potential energy surfaces and obtaining fundamental excitation energies for three moderate-size chain-like molecular systems. The use of our methodology leads to considerable computational savings for potential energy surface construction compared to standard approaches while maintaining a high level of accuracy in the resulting potential energy surfaces. Additional investigations indicate that our method can be applied to covalently bound and strongly interacting molecular systems, even though these cases are known to be very unfavorable for fragmentation schemes. We therefore conclude that the presented methodology is a robust and flexible approach to potential energy surface construction, which introduces considerable computational savings without compromising the accuracy of vibrational spectra calculations.
Collapse
Affiliation(s)
- Denis G Artiukhin
- Department of Chemistry, Aarhus Universitet, DK-8000 Aarhus, Denmark
| | | | - Carolin König
- Institute of Physical Chemistry, Christian-Albrechts-University Kiel, Max-Eyth-Straße 1, D-24118 Kiel, Germany
| | - Ove Christiansen
- Department of Chemistry, Aarhus Universitet, DK-8000 Aarhus, Denmark
| |
Collapse
|
91
|
Macetti G, Wieduwilt EK, Assfeld X, Genoni A. Localized Molecular Orbital-Based Embedding Scheme for Correlated Methods. J Chem Theory Comput 2020; 16:3578-3596. [PMID: 32369363 DOI: 10.1021/acs.jctc.0c00084] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Abstract
Embedding strategies currently provide the best compromise between accuracy and computational cost in modeling chemical properties and processes of large and complex systems. In this framework, different methods have been proposed all over the years, from the very popular QM/MM approaches to the more recent and very promising density matrix and density functional embedding techniques. Here, we present a further development of the quantum mechanics/extremely localized molecular orbital technique (QM/ELMO) method, a recently proposed multiscale embedding strategy in which the chemically active region of the investigated system is treated at a fully quantum mechanical level, while the rest is described by frozen extremely localized molecular orbitals previously transferred from proper libraries or tailor-made model molecules. In particular, in this work we discuss and assess in detail the extension of the QM/ELMO approach to density functional theory and post-Hartree-Fock techniques by evaluating its performances when it is used to describe chemical reactions, bond dissociations, and intermolecular interactions. The preliminary test calculations have shown that, in the investigated cases, the new embedding strategy enables the results of the corresponding fully quantum mechanical computations to be reproduced within chemical accuracy in almost all the cases but with a significantly reduced computational cost, especially when correlated post-Hartree-Fock strategies are used to describe the quantum mechanical subsystem. In light of the obtained results, we already envisage the future application of the new correlated QM/ELMO techniques to the investigation of more challenging problems, such as the modeling of enzyme catalysis, the study of excited states of biomolecules, and the refinement of macromolecular X-ray crystal structures.
Collapse
Affiliation(s)
- Giovanni Macetti
- Université de Lorraine & CNRS, Laboratoire de Physique et Chimie Théoriques (LPCT), UMR CNRS 7019, 1 Boulevard Arago, F-57078 Metz, France
| | - Erna K Wieduwilt
- Université de Lorraine & CNRS, Laboratoire de Physique et Chimie Théoriques (LPCT), UMR CNRS 7019, 1 Boulevard Arago, F-57078 Metz, France
| | - Xavier Assfeld
- Université de Lorraine & CNRS, Laboratoire de Physique et Chimie Théoriques (LPCT), UMR CNRS 7019, Boulevard des Aiguilletes, BP 70239, F-54506 Vandoeuvre-lès-Nancy, France
| | - Alessandro Genoni
- Université de Lorraine & CNRS, Laboratoire de Physique et Chimie Théoriques (LPCT), UMR CNRS 7019, 1 Boulevard Arago, F-57078 Metz, France
| |
Collapse
|
92
|
Nishimura Y, Nakai H. Hierarchical parallelization of divide‐and‐conquer density functional tight‐binding molecular dynamics and metadynamics simulations. J Comput Chem 2020; 41:1759-1772. [DOI: 10.1002/jcc.26217] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2020] [Revised: 04/15/2020] [Accepted: 04/20/2020] [Indexed: 11/08/2022]
Affiliation(s)
- Yoshifumi Nishimura
- Waseda Research Institute for Science and Engineering Waseda University Tokyo Japan
| | - Hiromi Nakai
- Waseda Research Institute for Science and Engineering Waseda University Tokyo Japan
- Department of Chemistry and Biochemistry School of Advanced Science and Engineering, Waseda University Tokyo Japan
- Elements Strategy Initiative for Catalysts and Batteries Kyoto University Kyoto Japan
| |
Collapse
|
93
|
Christensen EG, Steele RP. Stepwise Activation of Water by Open-Shell Interactions, Cl(H 2O) n=4–8,17. J Phys Chem A 2020; 124:3417-3437. [DOI: 10.1021/acs.jpca.0c01544] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Elizabeth G. Christensen
- Department of Chemistry and Henry Eyring Center for Theoretical Chemistry, University of Utah, 315 South 1400 East, Salt Lake City, Utah 84112, United States
| | - Ryan P. Steele
- Department of Chemistry and Henry Eyring Center for Theoretical Chemistry, University of Utah, 315 South 1400 East, Salt Lake City, Utah 84112, United States
| |
Collapse
|
94
|
Dawson W, Mohr S, Ratcliff LE, Nakajima T, Genovese L. Complexity Reduction in Density Functional Theory Calculations of Large Systems: System Partitioning and Fragment Embedding. J Chem Theory Comput 2020; 16:2952-2964. [PMID: 32216343 DOI: 10.1021/acs.jctc.9b01152] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
With the development of low order scaling methods for performing Kohn-Sham density functional theory, it is now possible to perform fully quantum mechanical calculations of systems containing tens of thousands of atoms. However, with an increase in the size of the system treated comes an increase in complexity, making it challenging to analyze such large systems and determine the cause of emergent properties. To address this issue, in this paper, we present a systematic complexity reduction methodology which can break down large systems into their constituent fragments and quantify interfragment interactions. The methodology proposed here requires no a priori information or user interaction, allowing a single workflow to be automatically applied to any system of interest. We apply this approach to a variety of different systems and show how it allows for the derivation of new system descriptors, the design of QM/MM partitioning schemes, and the novel application of graph metrics to molecules and materials.
Collapse
Affiliation(s)
- William Dawson
- RIKEN Center for Computational Science, Kobe 650-0047, Japan
| | - Stephan Mohr
- Barcelona Supercomputing Center (BSC), 08034 Barcelona, Spain
| | - Laura E Ratcliff
- Department of Materials, Imperial College London, London SW7 2AZ, United Kingdom
| | | | - Luigi Genovese
- Université Grenoble Alpes, INAC-MEM, L_Sim, Grenoble F-38000, France.,CEA, INAC-MEM, L_Sim, Grenoble F-38000, France
| |
Collapse
|
95
|
Noffke BW, Beckett D, Li LS, Raghavachari K. Aromatic Fragmentation Based on a Ring Overlap Scheme: An Algorithm for Large Polycyclic Aromatic Hydrocarbons Using the Molecules-in-Molecules Fragmentation-Based Method. J Chem Theory Comput 2020; 16:2160-2171. [DOI: 10.1021/acs.jctc.9b00566] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Affiliation(s)
- Benjamin W. Noffke
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| | - Daniel Beckett
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| | - Liang-shi Li
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| | - Krishnan Raghavachari
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| |
Collapse
|
96
|
Förster A, Franchini M, van Lenthe E, Visscher L. A Quadratic Pair Atomic Resolution of the Identity Based SOS-AO-MP2 Algorithm Using Slater Type Orbitals. J Chem Theory Comput 2020; 16:875-891. [PMID: 31930915 PMCID: PMC7027358 DOI: 10.1021/acs.jctc.9b00854] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/27/2019] [Indexed: 01/04/2023]
Abstract
We report a production level implementation of pair atomic resolution of the identity (PARI) based second-order Møller-Plesset perturbation theory (MP2) in the Slater type orbital (STO) based Amsterdam Density Functional (ADF) code. As demonstrated by systematic benchmarks, dimerization and isomerization energies obtained with our code using STO basis sets of triple-ζ-quality show mean absolute deviations from Gaussian type orbital, canonical, basis set limit extrapolated, global density fitting (DF)-MP2 results of less than 1 kcal/mol. Furthermore, we introduce a quadratic scaling atomic orbital based spin-opposite-scaled (SOS)-MP2 approach with a very small prefactor. Due to a worst-case scaling of [Formula: see text], our implementation is very fast already for small systems and shows an exceptionally early crossover to canonical SOS-PARI-MP2. We report computational wall time results for linear as well as for realistic three-dimensional molecules and show that triple-ζ quality calculations on molecules of several hundreds of atoms are only a matter of a few hours on a single compute node, the bottleneck of the computations being the SCF rather than the post-SCF energy correction.
Collapse
Affiliation(s)
- Arno Förster
- Theoretical Chemistry, Vrije
Universiteit, De Boelelaan 1083, NL-1081 HV Amsterdam, The
Netherlands
| | - Mirko Franchini
- Theoretical Chemistry, Vrije
Universiteit, De Boelelaan 1083, NL-1081 HV Amsterdam, The
Netherlands
- Scientific Computing & Modelling
NV, De Boelelaan 1083, NL-1081 HV Amsterdam, The
Netherlands
| | - Erik van Lenthe
- Scientific Computing & Modelling
NV, De Boelelaan 1083, NL-1081 HV Amsterdam, The
Netherlands
| | - Lucas Visscher
- Theoretical Chemistry, Vrije
Universiteit, De Boelelaan 1083, NL-1081 HV Amsterdam, The
Netherlands
| |
Collapse
|
97
|
Jones LO, Mosquera MA, Schatz GC, Ratner MA. Embedding Methods for Quantum Chemistry: Applications from Materials to Life Sciences. J Am Chem Soc 2020; 142:3281-3295. [PMID: 31986877 DOI: 10.1021/jacs.9b10780] [Citation(s) in RCA: 58] [Impact Index Per Article: 14.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
Quantum mechanical embedding methods hold the promise to transform not just the way calculations are performed, but to significantly reduce computational costs and improve scaling for macro-molecular systems containing hundreds if not thousands of atoms. The field of embedding has grown increasingly broad with many approaches of different intersecting flavors. In this perspective, we lay out the methods into two streams: QM:MM and QM:QM, showcasing the advantages and disadvantages of both. We provide a review of the literature, the underpinning theories including our contributions, and we highlight current applications with select examples spanning both materials and life sciences. We conclude with prospects and future outlook on embedding, and our view on the use of universal test case scenarios for cross-comparisons of the many available (and future) embedding theories.
Collapse
Affiliation(s)
- Leighton O Jones
- Department of Chemistry , Northwestern University , Evanston , Illinois 60208 , United States
| | - Martín A Mosquera
- Department of Chemistry , Northwestern University , Evanston , Illinois 60208 , United States
| | - George C Schatz
- Department of Chemistry , Northwestern University , Evanston , Illinois 60208 , United States
| | - Mark A Ratner
- Department of Chemistry , Northwestern University , Evanston , Illinois 60208 , United States
| |
Collapse
|
98
|
Abstract
Basic concepts in the analysis of binding using the fragment molecular orbital method are discussed at length: polarization, desolvation, and interaction. The components in the pair interaction energy decomposition analysis are introduced, and the analysis is illustrated for a water dimer and a protein-ligand complex.
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
|
99
|
Morais PA, Maia FF, Solis-Calero C, Caetano EWS, Freire VN, Carvalho HF. The urokinase plasminogen activator binding to its receptor: a quantum biochemistry description within an in/homogeneous dielectric function framework with application to uPA–uPAR peptide inhibitors. Phys Chem Chem Phys 2020; 22:3570-3583. [DOI: 10.1039/c9cp06530j] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
DFT calculations using the MFCC fragment-based model considering a spatial-dependent dielectric function based on the Poisson–Boltzmann approximation were performed to describe the uPA–uPAR interactions.
Collapse
Affiliation(s)
- Pablo A. Morais
- Instituto Federal de Educação
- Ciência e Tecnologia do Ceará
- Campus Horizonte
- Horizonte
- Brazil
| | - Francisco Franciné Maia
- Departamento de Ciências Naturais
- Matemática e Estatística
- Universidade Federal Rural do Semi-Árido
- Mossoró
- Brazil
| | - Christian Solis-Calero
- Departamento de Biologia Estrutural e Funcional
- Instituto de Biologia
- Universidade Estadual de Campinas
- Campinas
- Brazil
| | | | | | - Hernandes F. Carvalho
- Departamento de Biologia Estrutural e Funcional
- Instituto de Biologia
- Universidade Estadual de Campinas
- Campinas
- Brazil
| |
Collapse
|
100
|
Patel P, Wilson AK. Domain-based local pair natural orbital methods within the correlation consistent composite approach. J Comput Chem 2019; 41:800-813. [PMID: 31891196 DOI: 10.1002/jcc.26129] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/03/2019] [Revised: 10/29/2019] [Accepted: 12/01/2019] [Indexed: 01/15/2023]
Abstract
Ab initio composite approaches have been utilized to model and predict main group thermochemistry within 1 kcal mol-1 , on average, from well-established reliable experiments, primarily for molecules with less than 30 atoms. For molecules of increasing size and complexity, such as biomolecular complexes, composite methodologies have been limited in their application. Therefore, the domain-based local pair natural orbital (DLPNO) methods have been implemented within the correlation consistent composite approach (ccCA) framework, namely DLPNO-ccCA, to reduce the computational cost (disk space, CPU (central processing unit) time, memory) and predict energetic properties such as enthalpies of formation, noncovalent interactions, and conformation energies for organic biomolecular complexes including one of the largest molecules examined via composite strategies, within 1 kcal mol-1 , after calibration with 119 molecules and a set of linear alkanes. © 2019 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Prajay Patel
- Department of Chemistry, Michigan State University, East Lansing, Michigan, 48824
| | - Angela K Wilson
- Department of Chemistry, Michigan State University, East Lansing, Michigan, 48824
| |
Collapse
|