1
|
Tripathy V, Raghavachari K. Fragment-based models for dissociation of strong acids in water: Electrostatic embedding minimizes the dependence on the fragmentation schemes. J Chem Phys 2023; 159:124106. [PMID: 38127382 DOI: 10.1063/5.0164089] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/20/2023] [Accepted: 08/28/2023] [Indexed: 12/23/2023] Open
Abstract
Fragmentation methods such as MIM (Molecules-in-Molecules) provide a route to accurately model large systems and have been successful in predicting their structures, energies, and spectroscopic properties. However, their use is often limited to systems at equilibrium due to the inherent complications in the choice of fragments in systems away from equilibrium. Furthermore, the presence of charges resulting from any heterolytic bond breaking may increase the fragmentation error. We have previously suggested EE-MIM (Electrostatically Embedded Molecules-In-Molecules) as a method to mitigate the errors resulting from the missing long-range interactions in molecular clusters in equilibrium. Here, we show that the same method can be applied to improve the performance of MIM to solve the longstanding problem of dependency of the fragmentation energy error on the choice of the fragmentation scheme. We chose four widely used acid dissociation reactions (HCl, HClO4, HNO3, and H2SO4) as test cases due to their importance in chemical processes and complex reaction potential energy surfaces. Electrostatic embedding improves the performance at both one and two-layer MIM as shown by lower EE-MIM1 and EE-MIM2 errors. The EE-MIM errors are also demonstrated to be less dependent on the choice of the fragmentation scheme by analyzing the variation in fragmentation energy at the points with more than one possible fragmentation scheme (points where the fragmentation scheme changes). EE-MIM2 with M06-2X as the low-level resulted in a variation of less than 1 kcal/mol for all the cases and 1 kJ/mol for all but three cases, rendering our method fragmentation scheme-independent for acid dissociation processes.
Collapse
Affiliation(s)
- Vikrant Tripathy
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, USA
| | | |
Collapse
|
2
|
Galvez Vallejo JL, Snowdon C, Stocks R, Kazemian F, Yan Yu FC, Seidl C, Seeger Z, Alkan M, Poole D, Westheimer BM, Basha M, De La Pierre M, Rendell A, Izgorodina EI, Gordon MS, Barca GMJ. Toward an extreme-scale electronic structure system. J Chem Phys 2023; 159:044112. [PMID: 37497819 DOI: 10.1063/5.0156399] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/29/2023] [Accepted: 07/03/2023] [Indexed: 07/28/2023] Open
Abstract
Electronic structure calculations have the potential to predict key matter transformations for applications of strategic technological importance, from drug discovery to material science and catalysis. However, a predictive physicochemical characterization of these processes often requires accurate quantum chemical modeling of complex molecular systems with hundreds to thousands of atoms. Due to the computationally demanding nature of electronic structure calculations and the complexity of modern high-performance computing hardware, quantum chemistry software has historically failed to operate at such large molecular scales with accuracy and speed that are useful in practice. In this paper, novel algorithms and software are presented that enable extreme-scale quantum chemistry capabilities with particular emphasis on exascale calculations. This includes the development and application of the multi-Graphics Processing Unit (GPU) library LibCChem 2.0 as part of the General Atomic and Molecular Electronic Structure System package and of the standalone Extreme-scale Electronic Structure System (EXESS), designed from the ground up for scaling on thousands of GPUs to perform high-performance accurate quantum chemistry calculations at unprecedented speed and molecular scales. Among various results, we report that the EXESS implementation enables Hartree-Fock/cc-pVDZ plus RI-MP2/cc-pVDZ/cc-pVDZ-RIFIT calculations on an ionic liquid system with 623 016 electrons and 146 592 atoms in less than 45 min using 27 600 GPUs on the Summit supercomputer with a 94.6% parallel efficiency.
Collapse
Affiliation(s)
| | - Calum Snowdon
- School of Computing, Australian National University, Canberra 2601, ACT, Australia
| | - Ryan Stocks
- School of Computing, Australian National University, Canberra 2601, ACT, Australia
| | - Fazeleh Kazemian
- School of Computing, Australian National University, Canberra 2601, ACT, Australia
| | - Fiona Chuo Yan Yu
- School of Computing, Australian National University, Canberra 2601, ACT, Australia
| | - Christopher Seidl
- School of Computing, Australian National University, Canberra 2601, ACT, Australia
| | - Zoe Seeger
- School of Chemistry, Monash University, Clayton 3800, VIC, Australia
| | - Melisa Alkan
- Department of Chemistry, Iowa State University, Ames, Iowa 50011-3111, USA
| | - David Poole
- School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332, USA
| | - Bryce M Westheimer
- Department of Chemistry, Iowa State University, Ames, Iowa 50011-3111, USA
| | - Mehaboob Basha
- Pawsey Supercomputing Research Centre, Kensington, WA 6151, Australia
| | | | - Alistair Rendell
- College of Science and Engineering, Flinders University, Adelaide, SA 5042, Australia
| | | | | | - Giuseppe M J Barca
- School of Computing, Australian National University, Canberra 2601, ACT, Australia
| |
Collapse
|
3
|
Vornweg JR, Wolter M, Jacob CR. A simple and consistent quantum-chemical fragmentation scheme for proteins that includes two-body contributions. J Comput Chem 2023; 44:1634-1644. [PMID: 37171574 DOI: 10.1002/jcc.27114] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/27/2023] [Revised: 03/28/2023] [Accepted: 03/30/2023] [Indexed: 05/13/2023]
Abstract
The Molecular Fractionation with Conjugate Caps (MFCC) method is a popular fragmentation method for the quantum-chemical treatment of proteins. However, it does not account for interactions between the amino acid fragments, such as intramolecular hydrogen bonding. Here, we present a combination of the MFCC fragmentation scheme with a second-order many-body expansion (MBE) that consistently accounts for all fragment-fragment, fragment-cap, and cap-cap interactions, while retaining the overall simplicity of the MFCC scheme with its chemically meaningful fragments. We show that with the resulting MFCC-MBE(2) scheme, the errors in the total energies of selected polypeptides and proteins can be reduced by up to one order of magnitude and relative energies of different protein conformers can be predicted accurately.
Collapse
Affiliation(s)
- Johannes R Vornweg
- Institute of Physical and Theoretical Chemistry, Technische Universität Braunschweig, Braunschweig, Germany
| | - Mario Wolter
- Institute of Physical and Theoretical Chemistry, Technische Universität Braunschweig, Braunschweig, Germany
| | - Christoph R Jacob
- Institute of Physical and Theoretical Chemistry, Technische Universität Braunschweig, Braunschweig, Germany
| |
Collapse
|
4
|
Chen WK, Fang WH, Cui G. Extending multi-layer energy-based fragment method for excited-state calculations of large covalently bonded fragment systems. J Chem Phys 2023; 158:044110. [PMID: 36725521 DOI: 10.1063/5.0129458] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023] Open
Abstract
Recently, we developed a low-scaling Multi-Layer Energy-Based Fragment (MLEBF) method for accurate excited-state calculations and nonadiabatic dynamics simulations of nonbonded fragment systems. In this work, we extend the MLEBF method to treat covalently bonded fragment ones. The main idea is cutting a target system into many fragments according to chemical properties. Fragments with dangling bonds are first saturated by chemical groups; then, saturated fragments, together with the original fragments without dangling bonds, are grouped into different layers. The accurate total energy expression is formulated with the many-body energy expansion theory, in combination with the inclusion-exclusion principle that is used to delete the contribution of chemical groups introduced to saturate dangling bonds. Specifically, in a two-layer MLEBF model, the photochemically active and inert layers are calculated with high-level and efficient electronic structure methods, respectively. Intralayer and interlayer energies can be truncated at the two- or three-body interaction level. Subsequently, through several systems, including neutral and charged covalently bonded fragment systems, we demonstrate that MLEBF can provide accurate ground- and excited-state energies and gradients. Finally, we realize the structure, conical intersection, and path optimizations by combining our MLEBF program with commercial and free packages, e.g., ASE and SciPy. These developments make MLEBF a practical and reliable tool for studying complex photochemical and photophysical processes of large nonbonded and bonded fragment systems.
Collapse
Affiliation(s)
- Wen-Kai Chen
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, 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
|
5
|
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
|
6
|
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
|
7
|
Liao K, Wang S, Li W, Li S. Generalized energy-based fragmentation approach for calculations of solvation energies of large systems. Phys Chem Chem Phys 2021; 23:19394-19401. [PMID: 34490874 DOI: 10.1039/d1cp02814f] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Abstract
A generalized energy-based fragmentation (GEBF) approach has been combined with a universal solvation model based on solute electron density (SMD) to compute the solvation energies of general large systems (such as protein molecules) in solutions. In the GEBF-SMD method, the solvation energy of a target system could be combined by the corresponding solvation energies of various subsystems, each of which is embedded in the background point charges and surface charges on the surface of solute cavity at the positions of its atoms and neighbouring atoms outside of the subsystem. Our results show that the GEBF-SMD model could reproduce the conventional SMD solvation energies quite well for various proteins in solutions, and could significantly reduce the computational costs for the SMD calculations of large proteins. In addition, the GEBF-SMD approach is almost independent of the basis sets and the types of solvents (including protic, polar, and nonpolar ones). Also, the GEBF-SMD approach could reproduce the relative energies of various conformers of large systems in solutions. Therefore, the GEBF-SMD method is expected to be applicable for computing the solvation energies of a broad range of large systems.
Collapse
Affiliation(s)
- Kang Liao
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic, Chemistry of Ministry of Education, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing, 210023, P. R. China.
| | - Shirong Wang
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic, Chemistry of Ministry of Education, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing, 210023, P. R. China.
| | - Wei Li
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic, Chemistry of Ministry of Education, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing, 210023, P. R. China.
| | - Shuhua Li
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic, Chemistry of Ministry of Education, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing, 210023, P. R. China.
| |
Collapse
|
8
|
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
|
9
|
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
|
10
|
Paz ASP, Glover WJ. Diabatic Many-Body Expansion: Development and Application to Charge-Transfer Reactions. J Chem Theory Comput 2021; 17:1497-1511. [DOI: 10.1021/acs.jctc.0c01231] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Amiel S. P. Paz
- NYU Shanghai, 1555 Century Avenue, Shanghai 200122, China
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, 3663 Zhongshang Road North, Shanghai 200062, China
- Department of Chemistry, New York University, New York, New York 10003, United States
| | - William J. Glover
- NYU Shanghai, 1555 Century Avenue, Shanghai 200122, China
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, 3663 Zhongshang Road North, Shanghai 200062, China
- Department of Chemistry, New York University, New York, New York 10003, United States
| |
Collapse
|
11
|
Fu F, Liao K, Liu Z, Hong D, Yang H, Tian Y, Wei W, Liu C, Li S, Ma J, Li W. Controlled Fluorescence Enhancement of DNA-Binding Dye Through Chain Length Match between Oligoguanine and TOTO. J Phys Chem B 2021; 125:518-527. [PMID: 33426891 DOI: 10.1021/acs.jpcb.0c09611] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
Abstract
Fluorescent DNA-binding dyes are extensively employed as probe and biosensing in biological detection and imaging. Experiments and theoretical calculations of thiazole orange homodimeric (TOTO) dye binding to a single-strand DNA (ssDNA), poly(dG)n (n = 2, 4, 6, 8), reveal that the n = 6 complex shows about 300-fold stronger fluorescence than n = 2, 4 and a slightly stronger one than n = 8 complexes, which is benefited from the length match between TOTO and poly(dG)6. The machine learning, based on molecular dynamics trajectories, indicates that TOTO is featured by the dihedral angle along its backbone and its end-to-end distance, in which the latter one defines the stretch and hairpin structures of TOTO, respectively. The time-dependent density functional theory calculations on the low-lying excited states show that the stretched TOTO with π-π end-stacking binding mode can bring about strong fluorescence with localized π-π* transitions. For the n = 2, 4, and 8 complexes, the linear scaling quantum mechanics calculations indicate that the dominant hairpin TOTO with intercalative binding modes have relatively larger binding energies, leading to fluorescence quenching by intramolecular charge transfer. Our results may provide an insight for modulating the DNA-dye binding modes to tune the degree of charge transfer and designing fluorescent probes for the recognition of specific DNA sequences.
Collapse
Affiliation(s)
- Fangjia Fu
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing 210023, P. R. China
| | - Kang Liao
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing 210023, P. R. China
| | - Ziteng Liu
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing 210023, P. R. China
| | - Daocheng Hong
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing 210023, P. R. China
| | - Haitang Yang
- Jiangsu Engineering Laboratory of Smart Carbon-Rich Materials and Device, Jiangsu Province Hi-Tech Key Laboratory for Bio-medical Research, School of Chemistry and Chemical Engineering, Southeast University, Nanjing 211189, China
| | - Yuxi Tian
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing 210023, P. R. China
| | - Wei Wei
- Jiangsu Engineering Laboratory of Smart Carbon-Rich Materials and Device, Jiangsu Province Hi-Tech Key Laboratory for Bio-medical Research, School of Chemistry and Chemical Engineering, Southeast University, Nanjing 211189, China
| | - Chungen Liu
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing 210023, P. R. China
| | - Shuhua Li
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing 210023, P. R. China
| | - Jing Ma
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing 210023, P. R. China
| | - Wei Li
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing 210023, P. R. China
| |
Collapse
|
12
|
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
|
13
|
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
|
14
|
Wang Z, Han Y, Li J, He X. Combining the Fragmentation Approach and Neural Network Potential Energy Surfaces of Fragments for Accurate Calculation of Protein Energy. J Phys Chem B 2020; 124:3027-3035. [PMID: 32208716 DOI: 10.1021/acs.jpcb.0c01370] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Accurate and efficient all-atom quantum mechanical (QM) calculations for biomolecules still present a challenge to computational physicists and chemists. In this study, an extensible generalized molecular fractionation with a conjugate caps method combined with neural networks (NN-GMFCC) is developed for efficient QM calculation of protein energy. In the NN-GMFCC scheme, the total energy of a given protein is calculated by taking a proper combination of the high-precision neural network potential energies of all capped residues and overlapping conjugate caps. In addition, the two-body interaction energies of residue pairs are calculated by molecular mechanics (MM). With reference to the GMFCC/MM calculation at the ωB97XD/6-31G* level, the overall mean unsigned errors of the energy deviations and atomic force root-mean-squared errors calculated by NN-GMFCC are only 2.01 kcal/mol and 0.68 kcal/mol/Å, respectively, for 14 proteins (containing up to 13,728 atoms). Meanwhile, the NN-GMFCC approach is about 4 orders of magnitude faster than the GMFCC/MM method. The NN-GMFCC method could be systematically improved by inclusion of two-body QM interaction and multibody electronic polarization effect. Moreover, the NN-GMFCC approach can also be applied to other macromolecular systems such as DNA/RNA, and it is capable of providing a powerful and efficient approach for exploration of structures and functions of proteins with QM accuracy.
Collapse
Affiliation(s)
- Zhilong Wang
- Key Laboratory of Thin Film and Micro Fabrication, Ministry of Education, Department of Micro/Nano-electronics, Shanghai Jiao Tong University, Shanghai 200240, China
| | - Yanqiang Han
- Key Laboratory of Thin Film and Micro Fabrication, Ministry of Education, Department of Micro/Nano-electronics, Shanghai Jiao Tong University, Shanghai 200240, China
| | - Jinjin Li
- Key Laboratory of Thin Film and Micro Fabrication, Ministry of Education, Department of Micro/Nano-electronics, Shanghai Jiao Tong University, Shanghai 200240, China
| | - 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
|
15
|
Abstract
Since the introduction of the fragment molecular orbital method 20 years ago, fragment-based approaches have occupied a small but growing niche in quantum chemistry. These methods decompose a large molecular system into subsystems small enough to be amenable to electronic structure calculations, following which the subsystem information is reassembled in order to approximate an otherwise intractable supersystem calculation. Fragmentation sidesteps the steep rise (with respect to system size) in the cost of ab initio calculations, replacing it with a distributed cost across numerous computer processors. Such methods are attractive, in part, because they are easily parallelizable and therefore readily amenable to exascale computing. As such, there has been hope that distributed computing might offer the proverbial "free lunch" in quantum chemistry, with the entrée being high-level calculations on very large systems. While fragment-based quantum chemistry can count many success stories, there also exists a seedy underbelly of rarely acknowledged problems. As these methods begin to mature, it is time to have a serious conversation about what they can and cannot be expected to accomplish in the near future. Both successes and challenges are highlighted in this Perspective.
Collapse
Affiliation(s)
- John M Herbert
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| |
Collapse
|
16
|
Liu J, Rana B, Liu KY, Herbert JM. Variational Formulation of the Generalized Many-Body Expansion with Self-Consistent Charge Embedding: Simple and Correct Analytic Energy Gradient for Fragment-Based ab Initio Molecular Dynamics. J Phys Chem Lett 2019; 10:3877-3886. [PMID: 31251619 DOI: 10.1021/acs.jpclett.9b01214] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
The many-body expansion (MBE) and its extension to overlapping fragments, the generalized (G)MBE, constitute the theoretical basis for most fragment-based approaches for large-scale quantum chemistry. We reformulate the GMBE for use with embedding charges determined self-consistently from the fragment wave functions, in a manner that preserves the variational nature of the underlying self-consistent field method. As a result, the analytic gradient retains the simple "sum of fragment gradients" form that is often assumed in practice, sometimes incorrectly. This obviates (without approximation) the need to solve coupled-perturbed equations, and we demonstrate stable, fragment-based ab initio molecular dynamics simulations using this technique. Energy conservation fails when charge-response contributions to the Fock matrix are neglected, even while geometry optimizations and vibrational frequency calculations may yet be accurate. Stable simulations can be recovered by means of straightforward modifications introduced here, providing a general paradigm for fragment-based ab initio molecular dynamics.
Collapse
Affiliation(s)
- Jie Liu
- Department of Chemistry and Biochemistry , The Ohio State University , Columbus , Ohio 43210 , United States
| | - Bhaskar Rana
- Department of Chemistry and Biochemistry , The Ohio State University , Columbus , Ohio 43210 , United States
| | - Kuan-Yu Liu
- Department of Chemistry and Biochemistry , The Ohio State University , Columbus , Ohio 43210 , United States
| | - John M Herbert
- Department of Chemistry and Biochemistry , The Ohio State University , Columbus , Ohio 43210 , United States
| |
Collapse
|
17
|
Chi Y, You X. Kinetics of Hydrogen Abstraction Reactions of Methyl Palmitate and Octadecane by Hydrogen Atoms. J Phys Chem A 2019; 123:3058-3067. [PMID: 30893997 DOI: 10.1021/acs.jpca.8b08802] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Hydrogen abstractions play a crucial role in the consumption of fuel molecules during fuel pyrolysis and combustion processes. In this study, a generalized energy-based fragmentation approach was used to obtain CCSD(T)-F12a/cc-pVTZ energy barriers of hydrogen abstraction reactions by hydrogen atoms from methyl palmitate (C15H31COOCH3), a key component of biodiesel. The accuracy of M06-2X/6-311++G(d,p) for obtaining the energy barriers was evaluated against the CCSD(T) results. Based on the quantum chemical results, the high-pressure-limit rate constants for C15H31COOCH3 + H were calculated and compared with those of octadecane ( n-C18H38) reacting with H. The treatment of hindered internal rotations for such long-chain molecules was discussed and the rate rules for different abstraction sites were summarized. The results show that in the C15H31COOCH3 + H system, the α hydrogen abstraction no longer plays a dominant role as in small methyl esters, and the hydrogen atoms of CH2 groups far away from the ester group are more easily abstracted than those near the ester group.
Collapse
Affiliation(s)
- Yawei Chi
- Center for Combustion Energy , Tsinghua University , Beijing 100084 , China.,Key Laboratory for Thermal Science and Power Engineering of Ministry of Education , Tsinghua University , Beijing 100084 , China
| | - Xiaoqing You
- Center for Combustion Energy , Tsinghua University , Beijing 100084 , China.,Key Laboratory for Thermal Science and Power Engineering of Ministry of Education , Tsinghua University , Beijing 100084 , China
| |
Collapse
|
18
|
Fu F, Liao K, Ma J, Cheng Z, Zheng D, Gao L, Liu C, Li S, Li W. How intermolecular interactions influence electronic absorption spectra: insights from the molecular packing of uracil in condensed phases. Phys Chem Chem Phys 2019; 21:4072-4081. [DOI: 10.1039/c8cp06152a] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
Intermolecular interactions in terms of molecular packing are crucial for the investigation of the absorption spectra of uracil in different environments.
Collapse
Affiliation(s)
- Fangjia Fu
- School of Chemistry and Chemical Engineering
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education
- Institute of Theoretical and Computational Chemistry
- Nanjing University
- Nanjing
| | - Kang Liao
- School of Chemistry and Chemical Engineering
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education
- Institute of Theoretical and Computational Chemistry
- Nanjing University
- Nanjing
| | - Jing Ma
- School of Chemistry and Chemical Engineering
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education
- Institute of Theoretical and Computational Chemistry
- Nanjing University
- Nanjing
| | - Zheng Cheng
- School of Chemistry and Chemical Engineering
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education
- Institute of Theoretical and Computational Chemistry
- Nanjing University
- Nanjing
| | - Dong Zheng
- School of Chemistry and Chemical Engineering
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education
- Institute of Theoretical and Computational Chemistry
- Nanjing University
- Nanjing
| | - Liuzhou Gao
- School of Chemistry and Chemical Engineering
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education
- Institute of Theoretical and Computational Chemistry
- Nanjing University
- Nanjing
| | - Chungen Liu
- School of Chemistry and Chemical Engineering
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education
- Institute of Theoretical and Computational Chemistry
- Nanjing University
- Nanjing
| | - Shuhua Li
- School of Chemistry and Chemical Engineering
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education
- Institute of Theoretical and Computational Chemistry
- Nanjing University
- Nanjing
| | - Wei Li
- School of Chemistry and Chemical Engineering
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education
- Institute of Theoretical and Computational Chemistry
- Nanjing University
- Nanjing
| |
Collapse
|
19
|
Meitei OR, Heßelmann A. Geometry optimizations with the incremental molecular fragmentation method. JOURNAL OF THEORETICAL & COMPUTATIONAL CHEMISTRY 2018. [DOI: 10.1142/s0219633618500372] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
Nuclear energy gradients for the incremental molecular fragmentation (IMF) method presented in our previous work [Meitei OR, Heßelmann A, Molecular energies from an incremental fragmentation method, J Chem Phys 144(8):084109, 2016] have been derived. Using the second-order Møller–Plesset perturbation theory method to describe the bonded and nonbonded energy and gradient contributions and the uncorrelated Hartree–Fock method to describe the correction increment, it is shown that the IMF gradient can be easily computed by a sum of the underlying individual derivatives of the energy contributions. The performance of the method has been compared against the supermolecular method by optimizing the structures of a range of polyglycine molecules with up to 36 glycine residues in the chain. It is shown that with a sensible set of parameters used in the fragmentation the supermolecular structures can be fairly well reproduced. In a few cases the optimization with the IMF method leads to structures that differ from the supermolecular ones. It was found, however, that these are more stable geometries also on the supermolecular potential energy surface.
Collapse
Affiliation(s)
- Oinam Romesh Meitei
- Department Chemie und Pharmazie, Lehrstuhl für Theoretische Chemie, Friedrich-Alexander Universität Erlangen-Nürnberg, Egerlandstr. 3, D-91058 Erlangen, Germany
| | - Andreas Heßelmann
- Department Chemie und Pharmazie, Lehrstuhl für Theoretische Chemie, Friedrich-Alexander Universität Erlangen-Nürnberg, Egerlandstr. 3, D-91058 Erlangen, Germany
| |
Collapse
|
20
|
Thapa B, Beckett D, Jovan Jose KV, Raghavachari K. Assessment of Fragmentation Strategies for Large Proteins Using the Multilayer Molecules-in-Molecules Approach. J Chem Theory Comput 2018; 14:1383-1394. [DOI: 10.1021/acs.jctc.7b01198] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Affiliation(s)
- Bishnu Thapa
- Department of Chemistry, Indiana University, Bloomington 47405, Indiana, United States
| | - Daniel Beckett
- Department of Chemistry, Indiana University, Bloomington 47405, Indiana, United States
| | - K. V. Jovan Jose
- Department of Chemistry, Indiana University, Bloomington 47405, Indiana, United States
| | - Krishnan Raghavachari
- Department of Chemistry, Indiana University, Bloomington 47405, Indiana, United States
| |
Collapse
|
21
|
Sun CL, Liu LP, Tian F, Ding F, Wang LW. Charge-patching method for the calculation of electronic structure of polypeptides. Phys Chem Chem Phys 2018; 20:23301-23310. [DOI: 10.1039/c8cp01803k] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Based on the CPM method, the charge densities of polypeptides can be generated and their electronic structure can be further calculated.
Collapse
Affiliation(s)
- Chang-Liang Sun
- Center of Physical Chemistry Test
- Shenyang University of Chemical Technology
- Shenyang 110142
- People's Republic of China
- Materials Science Division
| | - Li-Ping Liu
- Materials Science Division
- Lawrence Berkeley National Laboratory
- Berkeley
- USA
- School of Physics
| | - Fubo Tian
- Materials Science Division
- Lawrence Berkeley National Laboratory
- Berkeley
- USA
- College of Physics
| | - Fu Ding
- Center of Physical Chemistry Test
- Shenyang University of Chemical Technology
- Shenyang 110142
- People's Republic of China
| | - Lin-Wang Wang
- Materials Science Division
- Lawrence Berkeley National Laboratory
- Berkeley
- USA
| |
Collapse
|
22
|
Li Y, Yuan D, Wang Q, Li W, Li S. Accurate prediction of the structure and vibrational spectra of ionic liquid clusters with the generalized energy-based fragmentation approach: critical role of ion-pair-based fragmentation. Phys Chem Chem Phys 2018; 20:13547-13557. [DOI: 10.1039/c8cp00513c] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The GEBF method with the ion-pair-based fragmentation has been developed to facilitate ab initio calculations of general ionic liquid clusters.
Collapse
Affiliation(s)
- Yunzhi Li
- School of Chemistry and Chemical Engineering
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education
- Institute of Theoretical and Computational Chemistry
- Nanjing University
- Nanjing
| | - Dandan Yuan
- School of Chemistry and Chemical Engineering
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education
- Institute of Theoretical and Computational Chemistry
- Nanjing University
- Nanjing
| | - Qingchun Wang
- School of Chemistry and Chemical Engineering
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education
- Institute of Theoretical and Computational Chemistry
- Nanjing University
- Nanjing
| | - Wei Li
- School of Chemistry and Chemical Engineering
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education
- Institute of Theoretical and Computational Chemistry
- Nanjing University
- Nanjing
| | - Shuhua Li
- School of Chemistry and Chemical Engineering
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education
- Institute of Theoretical and Computational Chemistry
- Nanjing University
- Nanjing
| |
Collapse
|
23
|
Zhang L, Li W, Fang T, Li S. Accurate Relative Energies and Binding Energies of Large Ice–Liquid Water Clusters and Periodic Structures. J Phys Chem A 2017; 121:4030-4038. [DOI: 10.1021/acs.jpca.7b03376] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Lei Zhang
- Institute of Theoretical
and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry
of MOE, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Wei Li
- Institute of Theoretical
and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry
of MOE, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Tao Fang
- Institute of Theoretical
and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry
of MOE, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Shuhua Li
- Institute of Theoretical
and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry
of MOE, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| |
Collapse
|
24
|
Yuan D, Li Y, Ni Z, Pulay P, Li W, Li S. Benchmark Relative Energies for Large Water Clusters with the Generalized Energy-Based Fragmentation Method. J Chem Theory Comput 2017; 13:2696-2704. [DOI: 10.1021/acs.jctc.7b00284] [Citation(s) in RCA: 29] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Affiliation(s)
- Dandan Yuan
- School
of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic
Chemistry of Ministry of Education, Institute of Theoretical and Computational
Chemistry, Nanjing University, Nanjing 210023, P. R. China
| | - Yunzhi Li
- School
of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic
Chemistry of Ministry of Education, Institute of Theoretical and Computational
Chemistry, Nanjing University, Nanjing 210023, P. R. China
| | - Zhigang Ni
- School
of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic
Chemistry of Ministry of Education, Institute of Theoretical and Computational
Chemistry, Nanjing University, Nanjing 210023, P. R. China
| | - Peter Pulay
- Department
of Chemistry and Biochemistry, University of Arkansas, Fayetteville, Arkansas 72701, United States
| | - Wei Li
- School
of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic
Chemistry of Ministry of Education, Institute of Theoretical and Computational
Chemistry, Nanjing University, Nanjing 210023, P. R. China
| | - Shuhua Li
- School
of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic
Chemistry of Ministry of Education, Institute of Theoretical and Computational
Chemistry, Nanjing University, Nanjing 210023, P. R. China
| |
Collapse
|
25
|
Chen GD, Weng J, Song G, Li ZH. Generalized Switch Functions in the Multilevel Many-Body Expansion Method and Its Application to Water Clusters. J Chem Theory Comput 2017; 13:2010-2020. [DOI: 10.1021/acs.jctc.7b00144] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Affiliation(s)
- Guo Dong Chen
- Collaborative Innovation Center of Chemistry for Energy Material, Shanghai Key Laboratory of Molecular Catalysis & Innovative Materials, Department of Chemistry, Fudan University, Shanghai 200433, China
| | - Jingwei Weng
- Collaborative Innovation Center of Chemistry for Energy Material, Shanghai Key Laboratory of Molecular Catalysis & Innovative Materials, Department of Chemistry, Fudan University, Shanghai 200433, China
| | - Guoliang Song
- Collaborative Innovation Center of Chemistry for Energy Material, Shanghai Key Laboratory of Molecular Catalysis & Innovative Materials, Department of Chemistry, Fudan University, Shanghai 200433, China
| | - Zhen Hua Li
- Collaborative Innovation Center of Chemistry for Energy Material, Shanghai Key Laboratory of Molecular Catalysis & Innovative Materials, Department of Chemistry, Fudan University, Shanghai 200433, China
| |
Collapse
|
26
|
Liu J, Qi LW, Zhang JZH, He X. Fragment Quantum Mechanical Method for Large-Sized Ion–Water Clusters. J Chem Theory Comput 2017; 13:2021-2034. [DOI: 10.1021/acs.jctc.7b00149] [Citation(s) in RCA: 43] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Affiliation(s)
- Jinfeng Liu
- State
Key Laboratory of Natural Medicines, Department of Basic Medicine
and Clinical Pharmacy, China Pharmaceutical University, Nanjing, 210009, China
| | - Lian-Wen Qi
- State
Key Laboratory of Natural Medicines, Department of Basic Medicine
and Clinical Pharmacy, China Pharmaceutical University, Nanjing, 210009, China
| | - John Z. H. Zhang
- School
of Chemistry and Molecular Engineering, East China Normal University, Shanghai, 200062, China
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai, 200062, China
- Department
of Chemistry, New York University, New York, New York 10003, United States
| | - Xiao He
- School
of Chemistry and Molecular Engineering, East China Normal University, Shanghai, 200062, China
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai, 200062, China
| |
Collapse
|
27
|
Jin X, Zhang JZH, He X. Full QM Calculation of RNA Energy Using Electrostatically Embedded Generalized Molecular Fractionation with Conjugate Caps Method. J Phys Chem A 2017; 121:2503-2514. [PMID: 28264557 DOI: 10.1021/acs.jpca.7b00859] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
In this study, the electrostatically embedded generalized molecular fractionation with conjugate caps (concaps) method (EE-GMFCC) was employed for efficient linear-scaling quantum mechanical (QM) calculation of total energies of RNAs. In the EE-GMFCC approach, the total energy of RNA is calculated by taking a proper combination of the QM energy of each nucleotide-centric fragment with large caps or small caps (termed EE-GMFCC-LC and EE-GMFCC-SC, respectively) deducted by the energies of concaps. The two-body QM interaction energy between non-neighboring ribonucleotides which are spatially in close contact are also taken into account for the energy calculation. Numerical studies were carried out to calculate the total energies of a number of RNAs using the EE-GMFCC-LC and EE-GMFCC-SC methods at levels of the Hartree-Fock (HF) method, density functional theory (DFT), and second-order many-body perturbation theory (MP2), respectively. The results show that the efficiency of the EE-GMFCC-SC method is about 3 times faster than the EE-GMFCC-LC method with minimal accuracy sacrifice. The EE-GMFCC-SC method is also applied for relative energy calculations of 20 different conformers of two RNA systems using HF and DFT, respectively. Both single-point and relative energy calculations demonstrate that the EE-GMFCC method has deviations from the full system results of only a few kcal/mol.
Collapse
Affiliation(s)
- Xinsheng Jin
- School of Chemistry and Molecular Engineering, East China Normal University , Shanghai 200062, China
| | - John Z H Zhang
- School of Chemistry and Molecular Engineering, East China Normal University , Shanghai 200062, China.,NYU-ECNU Center for Computational Chemistry at NYU Shanghai , Shanghai 200062, China.,Department of Chemistry, New York University , New York, New York 10003, United States
| | - Xiao He
- School of Chemistry and Molecular Engineering, East China Normal University , Shanghai 200062, China.,NYU-ECNU Center for Computational Chemistry at NYU Shanghai , Shanghai 200062, China
| |
Collapse
|
28
|
Fang T, Li Y, Li S. Generalized energy‐based fragmentation approach for modeling condensed phase systems. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2017. [DOI: 10.1002/wcms.1297] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Affiliation(s)
- Tao Fang
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute of Theoretical and Computational ChemistryNanjing University Nanjing P. R. China
| | - Yunzhi Li
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute of Theoretical and Computational ChemistryNanjing University Nanjing P. R. China
| | - Shuhua Li
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute of Theoretical and Computational ChemistryNanjing University Nanjing P. R. China
| |
Collapse
|
29
|
Lao KU, Liu KY, Richard RM, Herbert JM. Understanding the many-body expansion for large systems. II. Accuracy considerations. J Chem Phys 2017; 144:164105. [PMID: 27131529 DOI: 10.1063/1.4947087] [Citation(s) in RCA: 57] [Impact Index Per Article: 8.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022] Open
Abstract
To complement our study of the role of finite precision in electronic structure calculations based on a truncated many-body expansion (MBE, or "n-body expansion"), we examine the accuracy of such methods in the present work. Accuracy may be defined either with respect to a supersystem calculation computed at the same level of theory as the n-body calculations, or alternatively with respect to high-quality benchmarks. Both metrics are considered here. In applications to a sequence of water clusters, (H2O)N=6-55 described at the B3LYP/cc-pVDZ level, we obtain mean absolute errors (MAEs) per H2O monomer of ∼1.0 kcal/mol for two-body expansions, where the benchmark is a B3LYP/cc-pVDZ calculation on the entire cluster. Three- and four-body expansions exhibit MAEs of 0.5 and 0.1 kcal/mol/monomer, respectively, without resort to charge embedding. A generalized many-body expansion truncated at two-body terms [GMBE(2)], using 3-4 H2O molecules per fragment, outperforms all of these methods and affords a MAE of ∼0.02 kcal/mol/monomer, also without charge embedding. GMBE(2) requires significantly fewer (although somewhat larger) subsystem calculations as compared to MBE(4), reducing problems associated with floating-point roundoff errors. When compared to high-quality benchmarks, we find that error cancellation often plays a critical role in the success of MBE(n) calculations, even at the four-body level, as basis-set superposition error can compensate for higher-order polarization interactions. A many-body counterpoise correction is introduced for the GMBE, and its two-body truncation [GMBCP(2)] is found to afford good results without error cancellation. Together with a method such as ωB97X-V/aug-cc-pVTZ that can describe both covalent and non-covalent interactions, the GMBE(2)+GMBCP(2) approach provides an accurate, stable, and tractable approach for large systems.
Collapse
Affiliation(s)
- Ka Un Lao
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Kuan-Yu Liu
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Ryan M Richard
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - John M Herbert
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| |
Collapse
|
30
|
Li W, Li Y, Lin R, Li S. Generalized Energy-Based Fragmentation Approach for Localized Excited States of Large Systems. J Phys Chem A 2016; 120:9667-9677. [DOI: 10.1021/acs.jpca.6b11193] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Wei Li
- Institute of Theoretical
and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry
of MOE, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, People’s Repubic of China
| | - Yunzhi Li
- Institute of Theoretical
and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry
of MOE, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, People’s Repubic of China
| | - Ruochen Lin
- Institute of Theoretical
and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry
of MOE, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, People’s Repubic of China
| | - Shuhua Li
- Institute of Theoretical
and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry
of MOE, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, People’s Repubic of China
| |
Collapse
|
31
|
Collins MA. Can Systematic Molecular Fragmentation Be Applied to Direct Ab Initio Molecular Dynamics? J Phys Chem A 2016; 120:9281-9291. [DOI: 10.1021/acs.jpca.6b08739] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Michael A. Collins
- Research School of Chemistry, Australian National University, Canberra, ACT 0200, Australia
| |
Collapse
|
32
|
Saha A, Raghavachari K. Analysis of Different Fragmentation Strategies on a Variety of Large Peptides: Implementation of a Low Level of Theory in Fragment-Based Methods Can Be a Crucial Factor. J Chem Theory Comput 2016; 11:2012-23. [PMID: 26574406 DOI: 10.1021/ct501045s] [Citation(s) in RCA: 32] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023]
Abstract
We have investigated the performance of two classes of fragmentation methods developed in our group (Molecules-in-Molecules (MIM) and Many-Overlapping-Body (MOB) expansion), to reproduce the unfragmented MP2 energies on a test set composed of 10 small to large biomolecules. They have also been assessed to recover the relative energies of different motifs of the acetyl(ala)18NH2 system. Performance of different bond-cutting environments and the use of Hartree-Fock and different density functionals (as a low level of theory) in conjunction with the fragmentation strategies have been analyzed. Our investigation shows that while a low level of theory (for recovering long-range interactions) may not be necessary for small peptides, it provides a very effective strategy to accurately reproduce the total and relative energies of larger peptides such as the different motifs of the acetyl(ala)18NH2 system. Employing M06-2X as the low level of theory, the calculated mean total energy deviation (maximum deviation) in the total MP2 energies for the 10 molecules in the test set at MIM(d=3.5Å), MIM(η=9), and MOB(d=5Å) are 1.16 (2.31), 0.72 (1.87), and 0.43 (2.02) kcal/mol, respectively. The excellent performance suggests that such fragment-based methods should be of general use for the computation of accurate energies of large biomolecular systems.
Collapse
Affiliation(s)
- Arjun Saha
- Department of Chemistry, Indiana University , Bloomington, Indiana 47405, United States
| | - Krishnan Raghavachari
- Department of Chemistry, Indiana University , Bloomington, Indiana 47405, United States
| |
Collapse
|
33
|
Fang T, Jia J, Li S. Vibrational Spectra of Molecular Crystals with the Generalized Energy-Based Fragmentation Approach. J Phys Chem A 2016; 120:2700-11. [DOI: 10.1021/acs.jpca.5b10927] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Tao Fang
- School of Chemistry and Chemical
Engineering, Key Laboratory of Mesoscopic Chemistry of Ministry of
Education, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing 210093, People’s Republic of China
| | - Junteng Jia
- School of Chemistry and Chemical
Engineering, Key Laboratory of Mesoscopic Chemistry of Ministry of
Education, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing 210093, People’s Republic of China
| | - Shuhua Li
- School of Chemistry and Chemical
Engineering, Key Laboratory of Mesoscopic Chemistry of Ministry of
Education, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing 210093, People’s Republic of China
| |
Collapse
|
34
|
Zhang L, Li W, Fang T, Li S. Ab initio molecular dynamics with intramolecular noncovalent interactions for unsolvated polypeptides. Theor Chem Acc 2016. [DOI: 10.1007/s00214-015-1799-z] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023]
|
35
|
Liu J, Herbert JM. Pair-Pair Approximation to the Generalized Many-Body Expansion: An Alternative to the Four-Body Expansion for ab Initio Prediction of Protein Energetics via Molecular Fragmentation. J Chem Theory Comput 2016; 12:572-84. [PMID: 26730608 DOI: 10.1021/acs.jctc.5b00955] [Citation(s) in RCA: 43] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We introduce a "pair-pair" approximation to the generalized many-body expansion (pp-GMBE) as an approximation to a traditional four-body expansion, the latter of which is accurate but quickly becomes numerically unstable and ultimately intractable as the number of "bodies" (fragments) increases. The pp-GMBE method achieves a good balance between accuracy and efficiency by defining significant fragment pairs and then fragment quartets. An efficient fragmentation scheme is introduced for proteins such that the largest subsystems contain about 60 atoms. Application of the pp-GMBE method to proteins with as many as 70 amino acids (1142 atoms) reveals that pp-GMBE energies are quite faithful to those obtained when the same level of density functional theory is applied to the entire macromolecule. When combined with embedding charges obtained from natural population analysis, the pp-GMBE approach affords absolute energies that differ by 1-3 kcal/mol from full supersystem results, but it yields conformational energy profiles that are practically indistinguishable from the supersystem calculation at the same level of theory.
Collapse
Affiliation(s)
- Jie Liu
- Department of Chemistry and Biochemistry, The Ohio State University , Columbus, Ohio 43210, United States
| | - John M Herbert
- Department of Chemistry and Biochemistry, The Ohio State University , Columbus, Ohio 43210, United States
| |
Collapse
|
36
|
Yuan D, Shen X, Li W, Li S. Are fragment-based quantum chemistry methods applicable to medium-sized water clusters? Phys Chem Chem Phys 2016; 18:16491-500. [DOI: 10.1039/c6cp01931e] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023]
Abstract
The GEBF method is demonstrated to be more accurate than the EE-MB method for medium-sized water clusters.
Collapse
Affiliation(s)
- Dandan Yuan
- School of Chemistry and Chemical Engineering
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education
- Institute of Theoretical and Computational Chemistry
- Nanjing University
- Nanjing
| | - Xiaoling Shen
- School of Chemistry and Chemical Engineering
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education
- Institute of Theoretical and Computational Chemistry
- Nanjing University
- Nanjing
| | - Wei Li
- School of Chemistry and Chemical Engineering
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education
- Institute of Theoretical and Computational Chemistry
- Nanjing University
- Nanjing
| | - Shuhua Li
- School of Chemistry and Chemical Engineering
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education
- Institute of Theoretical and Computational Chemistry
- Nanjing University
- Nanjing
| |
Collapse
|
37
|
Liu J, Zhang JZH, He X. Fragment quantum chemical approach to geometry optimization and vibrational spectrum calculation of proteins. Phys Chem Chem Phys 2016; 18:1864-75. [DOI: 10.1039/c5cp05693d] [Citation(s) in RCA: 52] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Geometry optimization and vibrational spectra (infrared and Raman spectra) calculations of proteins are carried out by a quantum chemical approach using the EE-GMFCC (electrostatically embedded generalized molecular fractionation with conjugate caps) method (J. Phys. Chem. A, 2013, 117, 7149).
Collapse
Affiliation(s)
- Jinfeng Liu
- State Key Laboratory of Precision Spectroscopy
- Institute of Theoretical and Computational Science
- College of Chemistry and Molecular Engineering
- East China Normal University
- Shanghai
| | - John Z. H. Zhang
- State Key Laboratory of Precision Spectroscopy
- Institute of Theoretical and Computational Science
- College of Chemistry and Molecular Engineering
- East China Normal University
- Shanghai
| | - Xiao He
- State Key Laboratory of Precision Spectroscopy
- Institute of Theoretical and Computational Science
- College of Chemistry and Molecular Engineering
- East China Normal University
- Shanghai
| |
Collapse
|
38
|
Wen J, Li W, Chen S, Ma J. Simulations of molecular self-assembled monolayers on surfaces: packing structures, formation processes and functions tuned by intermolecular and interfacial interactions. Phys Chem Chem Phys 2016; 18:22757-71. [DOI: 10.1039/c6cp01049k] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
Simulations using QM and MM methods guide the rational design of functionalized SAMs on surfaces.
Collapse
Affiliation(s)
- Jin Wen
- Institute of Theoretical and Computational Chemistry
- Key Laboratory of Mesoscopic Chemistry of MOE
- School of Chemistry and Chemical Engineering
- Nanjing University
- Nanjing 210023
| | - Wei Li
- Institute of Theoretical and Computational Chemistry
- Key Laboratory of Mesoscopic Chemistry of MOE
- School of Chemistry and Chemical Engineering
- Nanjing University
- Nanjing 210023
| | - Shuang Chen
- Kuang Yaming Honors School
- Nanjing University
- Nanjing
- P. R. China
| | - Jing Ma
- Institute of Theoretical and Computational Chemistry
- Key Laboratory of Mesoscopic Chemistry of MOE
- School of Chemistry and Chemical Engineering
- Nanjing University
- Nanjing 210023
| |
Collapse
|
39
|
Collins MA. Molecular forces, geometries, and frequencies by systematic molecular fragmentation including embedded charges. J Chem Phys 2015; 141:094108. [PMID: 25194365 DOI: 10.1063/1.4894185] [Citation(s) in RCA: 37] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The accuracy of energies, energy gradients, and hessians evaluated by systematic molecular fragmentation is examined for a wide range of neutral molecules, zwitterions, and ions. A protocol is established that may employ embedded charges in conjunction with fragmentation to provide accurate evaluation of minimum energy geometries and vibrational frequencies in an automated procedure.
Collapse
Affiliation(s)
- Michael A Collins
- Research School of Chemistry, Australian National University, Canberra, ACT, Australia
| |
Collapse
|
40
|
Akimov AV, Prezhdo OV. Large-Scale Computations in Chemistry: A Bird’s Eye View of a Vibrant Field. Chem Rev 2015; 115:5797-890. [DOI: 10.1021/cr500524c] [Citation(s) in RCA: 159] [Impact Index Per Article: 17.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Affiliation(s)
- Alexey V. Akimov
- Department
of Chemistry, University of South California, Los Angeles, California 90089, United States
| | - Oleg V. Prezhdo
- Department
of Chemistry, University of South California, Los Angeles, California 90089, United States
| |
Collapse
|
41
|
Raghavachari K, Saha A. Accurate Composite and Fragment-Based Quantum Chemical Models for Large Molecules. Chem Rev 2015; 115:5643-77. [PMID: 25849163 DOI: 10.1021/cr500606e] [Citation(s) in RCA: 183] [Impact Index Per Article: 20.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Affiliation(s)
- Krishnan Raghavachari
- Department of Chemistry, Indiana University, 800 E. Kirkwood Avenue, Bloomington, Indiana 47405, United States
| | - Arjun Saha
- Department of Chemistry, Indiana University, 800 E. Kirkwood Avenue, Bloomington, Indiana 47405, United States
| |
Collapse
|
42
|
Affiliation(s)
- Michael A Collins
- †Research School of Chemistry, Australian National University, Canberra, ACT 0200, Australia
| | - Ryan P A Bettens
- ‡Department of Chemistry, National University of Singapore, 3 Science Drive 3, Singapore 117543, Singapore
| |
Collapse
|
43
|
Lao KU, Herbert JM. Accurate and Efficient Quantum Chemistry Calculations for Noncovalent Interactions in Many-Body Systems: The XSAPT Family of Methods. J Phys Chem A 2014; 119:235-52. [DOI: 10.1021/jp5098603] [Citation(s) in RCA: 72] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Affiliation(s)
- Ka Un Lao
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| | - John M. Herbert
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| |
Collapse
|
44
|
Fang T, Li W, Gu F, Li S. Accurate Prediction of Lattice Energies and Structures of Molecular Crystals with Molecular Quantum Chemistry Methods. J Chem Theory Comput 2014; 11:91-8. [DOI: 10.1021/ct500833k] [Citation(s) in RCA: 36] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023]
Affiliation(s)
- Tao Fang
- School of Chemistry and Chemical
Engineering, Key Laboratory of Mesoscopic Chemistry of Ministry of
Education, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing, 210093, People’s Republic of China
| | - Wei Li
- School of Chemistry and Chemical
Engineering, Key Laboratory of Mesoscopic Chemistry of Ministry of
Education, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing, 210093, People’s Republic of China
| | - Fangwei Gu
- School of Chemistry and Chemical
Engineering, Key Laboratory of Mesoscopic Chemistry of Ministry of
Education, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing, 210093, People’s Republic of China
| | - Shuhua Li
- School of Chemistry and Chemical
Engineering, Key Laboratory of Mesoscopic Chemistry of Ministry of
Education, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing, 210093, People’s Republic of China
| |
Collapse
|
45
|
Timm MJ, Matta CF, Massa L, Huang L. The Localization–Delocalization Matrix and the Electron-Density-Weighted Connectivity Matrix of a Finite Graphene Nanoribbon Reconstructed from Kernel Fragments. J Phys Chem A 2014; 118:11304-16. [DOI: 10.1021/jp508490p] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Matthew J. Timm
- Department
of Chemistry and Physics, Mount Saint Vincent University, Halifax, Nova Scotia B3M2J6, Canada
| | - Chérif F. Matta
- Department
of Chemistry and Physics, Mount Saint Vincent University, Halifax, Nova Scotia B3M2J6, Canada
- Department
of Chemistry, Dalhousie University, Halifax, Nova Scotia B3H4J3, Canada
- Department
of Chemistry, Saint Mary’s University, Halifax, Nova Scotia B3H3C3, Canada
| | - Lou Massa
- Hunter
College and the Graduate School, City University of New York, New York, New York 10065, United States
| | - Lulu Huang
- Center
for Computational Materials Science, Naval Research Laboratory, Washington, D.C. 20375-5341, United States
| |
Collapse
|
46
|
Brorsen KR, Zahariev F, Nakata H, Fedorov DG, Gordon MS. Analytic Gradient for Density Functional Theory Based on the Fragment Molecular Orbital Method. J Chem Theory Comput 2014; 10:5297-307. [DOI: 10.1021/ct500808p] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Kurt R. Brorsen
- Ames
Laboratory, U.S. Department of Energy (US-DOE), Department of Chemistry, Iowa State University, Ames, Iowa 50011, United States
| | - Federico Zahariev
- Ames
Laboratory, U.S. Department of Energy (US-DOE), Department of Chemistry, Iowa State University, Ames, Iowa 50011, United States
| | - Hiroya Nakata
- Department
of Biomolecular Engineering, Tokyo Institute of Technology, 4259 Nagatsuta-cho, Midori-ku, Kanagawa, Yokohama 226-8501, Japan
- Nakamura
Laboratory, RIKEN Research Cluster for Innovation, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
- Japan Society
for the Promotion of Science, Kojimachi
Business Center Building, 5-3-1 Kojimachi, Chiyoda-ku, Tokyo 102-0083, Japan
| | - Dmitri G. Fedorov
- NRI, National
Institute of Advanced Industrial Science and Technology (AIST), 1-1-1 Umezono, Tsukuba, Ibaraki 305-8568, Japan
| | - Mark S. Gordon
- Ames
Laboratory, U.S. Department of Energy (US-DOE), Department of Chemistry, Iowa State University, Ames, Iowa 50011, United States
| |
Collapse
|
47
|
Shrestha K, González-Delgado JM, Blew JH, Jakubikova E. Electronic Structure of Covalently Linked Zinc Bacteriochlorin Molecular Arrays: Insights into Molecular Design for NIR Light Harvesting. J Phys Chem A 2014; 118:9901-13. [DOI: 10.1021/jp507749c] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Affiliation(s)
- Kushal Shrestha
- Department
of Chemistry, North Carolina State University, Raleigh, North Carolina 27695, United States
| | - Jessica M. González-Delgado
- Department
of Chemistry, University of Puerto Rico, Rio Piedras Campus, San Juan, Puerto Rico 00931, United States
| | - James H. Blew
- Department
of Chemistry, North Carolina State University, Raleigh, North Carolina 27695, United States
| | - Elena Jakubikova
- Department
of Chemistry, North Carolina State University, Raleigh, North Carolina 27695, United States
| |
Collapse
|
48
|
Li S, Li W, Ma J. Generalized energy-based fragmentation approach and its applications to macromolecules and molecular aggregates. Acc Chem Res 2014; 47:2712-20. [PMID: 24873495 DOI: 10.1021/ar500038z] [Citation(s) in RCA: 125] [Impact Index Per Article: 12.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Conspectus The generalized energy-based fragmentation (GEBF) approach provides a very simple way of approximately evaluating the ground-state energy or properties of a large system in terms of ground-state energies of various small "electrostatically embedded" subsystems, which can be calculated with any traditional ab initio quantum chemistry (X) method (X = Hartree-Fock, density functional theory, and so on). Due to its excellent parallel efficiency, the GEBF approach at the X theory level (GEBF-X) allows full quantum mechanical (QM) calculations to be accessible for systems with hundreds and even thousands of atoms on ordinary workstations. The implementation of the GEBF approach at various theoretical levels can be easily done with existing quantum chemistry programs. This Account reviews the methodology, implementation, and applications of the GEBF-X approach. This method has been successfully applied to optimize the structures of various large systems including molecular clusters, polypeptides, proteins, and foldamers. Such investigations could allow us to elucidate the origin and nature of the cooperative interaction in secondary structures of long peptides or the driving force of the self-assembly processes of aromatic oligoamides. These GEBF-based QM calculations reveal that the structures and stability of various complex systems result from a subtle balance of many types of noncovalent interactions such as hydrogen bonding and van der Waals interactions. The GEBF-based ab initio molecular dynamics (AIMD) method also allows the investigation of dynamic behaviors of large systems on the order of tens of picoseconds. It was demonstrated that the conformational dynamics of two model peptides predicted by GEBF-based AIMD are noticeably different from those predicted by the classical force field MD method. With the target of extending QM calculations to molecular aggregates in the condensed phase, we have implemented the GEBF-based multilayer hybrid models, which could provide satisfactory descriptions of the binding energies between a solute molecule and its surrounding waters and the chain-length dependence of the conformational changes of oligomers in aqueous solutions. A coarse-grained polarizable molecular mechanics model, furnished with GEBF-X dipole moments of subsystems, exhibits some advantages of treating the electrostatic polarization with reduced computational costs. We anticipate that the GEBF approach will continue to develop with the ultimate goal of studying complicated phenomena at mesoscopic scales and serve as a practical tool to elucidate the structure and dynamics of chemical and biological systems.
Collapse
Affiliation(s)
- Shuhua Li
- School of Chemistry and Chemical
Engineering, Key Laboratory of Mesoscopic Chemistry of MOE, Institute
of Theoretical and Computational Chemistry, Nanjing University, Nanjing 210093, People’s Republic of China
| | - Wei Li
- School of Chemistry and Chemical
Engineering, Key Laboratory of Mesoscopic Chemistry of MOE, Institute
of Theoretical and Computational Chemistry, Nanjing University, Nanjing 210093, People’s Republic of China
| | - Jing Ma
- School of Chemistry and Chemical
Engineering, Key Laboratory of Mesoscopic Chemistry of MOE, Institute
of Theoretical and Computational Chemistry, Nanjing University, Nanjing 210093, People’s Republic of China
| |
Collapse
|
49
|
Richard RM, Lao KU, Herbert JM. Understanding the many-body expansion for large systems. I. Precision considerations. J Chem Phys 2014; 141:014108. [DOI: 10.1063/1.4885846] [Citation(s) in RCA: 58] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Ryan M. Richard
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Ka Un Lao
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - John M. Herbert
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| |
Collapse
|
50
|
Guo Y, Li W, Li S. Improved cluster-in-molecule local correlation approach for electron correlation calculation of large systems. J Phys Chem A 2014; 118:8996-9004. [PMID: 24963784 DOI: 10.1021/jp501976x] [Citation(s) in RCA: 36] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
An improved cluster-in-molecule (CIM) local correlation approach is developed to allow electron correlation calculations of large systems more accurate and faster. We have proposed a refined strategy of constructing virtual LMOs of various clusters, which is suitable for basis sets of various types. To recover medium-range electron correlation, which is important for quantitative descriptions of large systems, we find that a larger distance threshold (ξ) is necessary for highly accurate results. Our illustrative calculations show that the present CIM-MP2 (second-order Møller-Plesser perturbation theory, MP2) or CIM-CCSD (coupled cluster singles and doubles, CCSD) scheme with a suitable ξ value is capable of recovering more than 99.8% correlation energies for a wide range of systems at different basis sets. Furthermore, the present CIM-MP2 scheme can provide reliable relative energy differences as the conventional MP2 method for secondary structures of polypeptides.
Collapse
Affiliation(s)
- Yang Guo
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of MOE, Institute of Theoretical and Computational Chemistry, Nanjing University , Nanjing 210093, P. R. China
| | | | | |
Collapse
|