1
|
Cao Y, Balduf T, Beachy MD, Bennett MC, Bochevarov AD, Chien A, Dub PA, Dyall KG, Furness JW, Halls MD, Hughes TF, Jacobson LD, Kwak HS, Levine DS, Mainz DT, Moore KB, Svensson M, Videla PE, Watson MA, Friesner RA. Quantum chemical package Jaguar: A survey of recent developments and unique features. J Chem Phys 2024; 161:052502. [PMID: 39092934 DOI: 10.1063/5.0213317] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/09/2024] [Accepted: 07/12/2024] [Indexed: 08/04/2024] Open
Abstract
This paper is dedicated to the quantum chemical package Jaguar, which is commercial software developed and distributed by Schrödinger, Inc. We discuss Jaguar's scientific features that are relevant to chemical research as well as describe those aspects of the program that are pertinent to the user interface, the organization of the computer code, and its maintenance and testing. Among the scientific topics that feature prominently in this paper are the quantum chemical methods grounded in the pseudospectral approach. A number of multistep workflows dependent on Jaguar are covered: prediction of protonation equilibria in aqueous solutions (particularly calculations of tautomeric stability and pKa), reactivity predictions based on automated transition state search, assembly of Boltzmann-averaged spectra such as vibrational and electronic circular dichroism, as well as nuclear magnetic resonance. Discussed also are quantum chemical calculations that are oriented toward materials science applications, in particular, prediction of properties of optoelectronic materials and organic semiconductors, and molecular catalyst design. The topic of treatment of conformations inevitably comes up in real world research projects and is considered as part of all the workflows mentioned above. In addition, we examine the role of machine learning methods in quantum chemical calculations performed by Jaguar, from auxiliary functions that return the approximate calculation runtime in a user interface, to prediction of actual molecular properties. The current work is second in a series of reviews of Jaguar, the first having been published more than ten years ago. Thus, this paper serves as a rare milestone on the path that is being traversed by Jaguar's development in more than thirty years of its existence.
Collapse
Affiliation(s)
- Yixiang Cao
- Schrödinger, Inc., 1540 Broadway, Floor 24, New York, New York 10036, USA
| | - Ty Balduf
- Schrödinger, Inc., 1540 Broadway, Floor 24, New York, New York 10036, USA
| | - Michael D Beachy
- Schrödinger, Inc., 1540 Broadway, Floor 24, New York, New York 10036, USA
| | - M Chandler Bennett
- Schrödinger, Inc., 1540 Broadway, Floor 24, New York, New York 10036, USA
| | - Art D Bochevarov
- Schrödinger, Inc., 1540 Broadway, Floor 24, New York, New York 10036, USA
| | - Alan Chien
- Schrödinger, Inc., 1540 Broadway, Floor 24, New York, New York 10036, USA
| | - Pavel A Dub
- Schrödinger, Inc., 9868 Scranton Road, Suite 3200, San Diego, California 92121, USA
| | - Kenneth G Dyall
- Schrödinger, Inc., 101 SW Main St., Suite 1300, Portland, Oregon 97204, USA
| | - James W Furness
- Schrödinger, Inc., 1540 Broadway, Floor 24, New York, New York 10036, USA
| | - Mathew D Halls
- Schrödinger, Inc., 9868 Scranton Road, Suite 3200, San Diego, California 92121, USA
| | - Thomas F Hughes
- Schrödinger, Inc., 1540 Broadway, Floor 24, New York, New York 10036, USA
| | - Leif D Jacobson
- Schrödinger, Inc., 101 SW Main St., Suite 1300, Portland, Oregon 97204, USA
| | - H Shaun Kwak
- Schrödinger, Inc., 101 SW Main St., Suite 1300, Portland, Oregon 97204, USA
| | - Daniel S Levine
- Schrödinger, Inc., 1540 Broadway, Floor 24, New York, New York 10036, USA
| | - Daniel T Mainz
- Schrödinger, Inc., 1540 Broadway, Floor 24, New York, New York 10036, USA
| | - Kevin B Moore
- Schrödinger, Inc., 1540 Broadway, Floor 24, New York, New York 10036, USA
| | - Mats Svensson
- Schrödinger, Inc., 1540 Broadway, Floor 24, New York, New York 10036, USA
| | - Pablo E Videla
- Schrödinger, Inc., 1540 Broadway, Floor 24, New York, New York 10036, USA
| | - Mark A Watson
- Schrödinger, Inc., 1540 Broadway, Floor 24, New York, New York 10036, USA
| | - Richard A Friesner
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, USA
| |
Collapse
|
2
|
Cui Q, Peng J, Xu C, Lan Z. Automatic Approach to Explore the Multireaction Mechanism for Medium-Sized Bimolecular Reactions via Collision Dynamics Simulations and Transition State Searches. J Chem Theory Comput 2022; 18:910-924. [DOI: 10.1021/acs.jctc.1c00795] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Qinghai Cui
- Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety and MOE Key Laboratory of Environmental Theoretical Chemistry, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, P. R. China
| | - Jiawei Peng
- Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety and MOE Key Laboratory of Environmental Theoretical Chemistry, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, P. R. China
| | - Chao Xu
- Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety and MOE Key Laboratory of Environmental Theoretical Chemistry, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, P. R. China
- Key Laboratory of Theoretical Chemistry of Environment, Ministry of Education; School of Chemistry, South China Normal University, Guangzhou 510006, P. R. China
| | - Zhenggang Lan
- Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety and MOE Key Laboratory of Environmental Theoretical Chemistry, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, P. R. China
| |
Collapse
|
3
|
Coley CW, Eyke NS, Jensen KF. Autonomous Discovery in the Chemical Sciences Part I: Progress. Angew Chem Int Ed Engl 2020; 59:22858-22893. [DOI: 10.1002/anie.201909987] [Citation(s) in RCA: 100] [Impact Index Per Article: 25.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/06/2019] [Indexed: 01/05/2023]
Affiliation(s)
- Connor W. Coley
- Department of Chemical Engineering Massachusetts Institute of Technology Cambridge MA 02139 USA
| | - Natalie S. Eyke
- Department of Chemical Engineering Massachusetts Institute of Technology Cambridge MA 02139 USA
| | - Klavs F. Jensen
- Department of Chemical Engineering Massachusetts Institute of Technology Cambridge MA 02139 USA
| |
Collapse
|
4
|
Coley CW, Eyke NS, Jensen KF. Autonome Entdeckung in den chemischen Wissenschaften, Teil I: Fortschritt. Angew Chem Int Ed Engl 2020. [DOI: 10.1002/ange.201909987] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Affiliation(s)
- Connor W. Coley
- Department of Chemical Engineering Massachusetts Institute of Technology Cambridge MA 02139 USA
| | - Natalie S. Eyke
- Department of Chemical Engineering Massachusetts Institute of Technology Cambridge MA 02139 USA
| | - Klavs F. Jensen
- Department of Chemical Engineering Massachusetts Institute of Technology Cambridge MA 02139 USA
| |
Collapse
|
5
|
Lindgren P, Kastlunger G, Peterson AA. Scaled and Dynamic Optimizations of Nudged Elastic Bands. J Chem Theory Comput 2019; 15:5787-5793. [PMID: 31600078 DOI: 10.1021/acs.jctc.9b00633] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We present a modified nudged elastic band routine that can reduce the number of force calls by more than 50% for bands with nonuniform convergence. The method, which we call "dyNEB", dynamically and selectively optimizes images on the basis of the perpendicular PES-derived forces and parallel spring forces acting on that region of the band. The convergence criteria are scaled to focus on the region of interest, i.e., the saddle point, while maintaining continuity of the band and avoiding truncation. We show that this method works well for solid state reaction barriers-nonelectrochemical in general and electrochemical in particular-and that the number of force calls can be significantly reduced without loss of resolution at the saddle point.
Collapse
Affiliation(s)
- Per Lindgren
- School of Engineering , Brown University , Providence , Rhode Island 02912 , United States
| | - Georg Kastlunger
- School of Engineering , Brown University , Providence , Rhode Island 02912 , United States
| | - Andrew A Peterson
- School of Engineering , Brown University , Providence , Rhode Island 02912 , United States
| |
Collapse
|
6
|
Meyer R, Schmuck KS, Hauser AW. Machine Learning in Computational Chemistry: An Evaluation of Method Performance for Nudged Elastic Band Calculations. J Chem Theory Comput 2019; 15:6513-6523. [PMID: 31553610 DOI: 10.1021/acs.jctc.9b00708] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The localization of transition states and the calculation of reaction pathways are routine tasks of computational chemists but often very CPU-intense problems, in particular for large systems. The standard algorithm for this purpose is the nudged elastic band method, but it has become obvious that an "intelligent" selection of points to be evaluated on the potential energy surface can improve its convergence significantly. This article summarizes, compares, and extends known strategies that have been heavily inspired by the machine learning developments of recent years. It presents advantages and disadvantages and provides an unbiased comparison of neural network based approaches, Gaussian process regression in Cartesian coordinates, and Gaussian approximation potentials. We test their performance on two example reactions, the ethane rotation and the activation of carbon dioxide on a metal catalyst, and provide a clear ranking in terms of usability for future implementations.
Collapse
Affiliation(s)
- Ralf Meyer
- Graz University of Technology , Institute of Experimental Physics , Petersgasse 16 , 8010 Graz , Austria
| | - Klemens S Schmuck
- Graz University of Technology , Institute of Experimental Physics , Petersgasse 16 , 8010 Graz , Austria
| | - Andreas W Hauser
- Graz University of Technology , Institute of Experimental Physics , Petersgasse 16 , 8010 Graz , Austria
| |
Collapse
|
7
|
Plessow PN. Efficient Transition State Optimization of Periodic Structures through Automated Relaxed Potential Energy Surface Scans. J Chem Theory Comput 2018; 14:981-990. [DOI: 10.1021/acs.jctc.7b01070] [Citation(s) in RCA: 35] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Philipp N. Plessow
- Institute of Catalysis Research and Technology (IKFT), Hermann-von-Helmholtz-Platz 1, D-76344 Eggenstein-Leopoldshafen, Germany
| |
Collapse
|
8
|
Celik FE, Peters B, Coppens MO, McCormick A, Hicks RF, Ekerdt J. A Career in Catalysis: Alexis T. Bell. ACS Catal 2017. [DOI: 10.1021/acscatal.7b03218] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Fuat E. Celik
- Chemical and Biochemical Engineering, Rutgers University, Piscataway, New Jersey 08854, United States
| | - Baron Peters
- Chemical Engineering, University of California, Santa Barbara, Santa Barbara, California, 93106, United States
| | - Marc-Olivier Coppens
- Chemical Engineering, University College London, London WC1E 7JE, United Kingdom
| | - Alon McCormick
- Chemical Engineering and Materials Science, University of Minnesota, Minneapolis, Minneapolis 55455, United States
| | - Robert F. Hicks
- Chemical and Biomolecular Engineering, University of California, Los Angeles, Los Angeles, California 90095, United States
| | - John Ekerdt
- McKetta Department of Chemical Engineering, University of Texas, Austin, Texas 78712, United States
| |
Collapse
|
9
|
Fast exploration of an optimal path on the multidimensional free energy surface. PLoS One 2017; 12:e0177740. [PMID: 28542475 PMCID: PMC5436793 DOI: 10.1371/journal.pone.0177740] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/24/2016] [Accepted: 05/02/2017] [Indexed: 11/29/2022] Open
Abstract
In a reaction, determination of an optimal path with a high reaction rate (or a low free energy barrier) is important for the study of the reaction mechanism. This is a complicated problem that involves lots of degrees of freedom. For simple models, one can build an initial path in the collective variable space by the interpolation method first and then update the whole path constantly in the optimization. However, such interpolation method could be risky in the high dimensional space for large molecules. On the path, steric clashes between neighboring atoms could cause extremely high energy barriers and thus fail the optimization. Moreover, performing simulations for all the snapshots on the path is also time-consuming. In this paper, we build and optimize the path by a growing method on the free energy surface. The method grows a path from the reactant and extends its length in the collective variable space step by step. The growing direction is determined by both the free energy gradient at the end of the path and the direction vector pointing at the product. With fewer snapshots on the path, this strategy can let the path avoid the high energy states in the growing process and save the precious simulation time at each iteration step. Applications show that the presented method is efficient enough to produce optimal paths on either the two-dimensional or the twelve-dimensional free energy surfaces of different small molecules.
Collapse
|
10
|
Quantum chemical approaches to [NiFe] hydrogenase. Essays Biochem 2017; 61:293-303. [PMID: 28487405 DOI: 10.1042/ebc20160079] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/15/2016] [Revised: 02/22/2017] [Accepted: 03/01/2017] [Indexed: 11/17/2022]
Abstract
The mechanism by which [NiFe] hydrogenase catalyses the oxidation of molecular hydrogen is a significant yet challenging topic in bioinorganic chemistry. With far-reaching applications in renewable energy and carbon mitigation, significant effort has been invested in the study of these complexes. In particular, computational approaches offer a unique perspective on how this enzyme functions at an electronic and atomistic level. In this article, we discuss state-of-the art quantum chemical methods and how they have helped deepen our comprehension of [NiFe] hydrogenase. We outline the key strategies that can be used to compute the (i) geometry, (ii) electronic structure, (iii) thermodynamics and (iv) kinetic properties associated with the enzymatic activity of [NiFe] hydrogenase and other bioinorganic complexes.
Collapse
|
11
|
Jafari M, Zimmerman PM. Reliable and efficient reaction path and transition state finding for surface reactions with the growing string method. J Comput Chem 2017; 38:645-658. [PMID: 28130776 DOI: 10.1002/jcc.24720] [Citation(s) in RCA: 31] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/19/2016] [Revised: 10/20/2016] [Accepted: 12/10/2016] [Indexed: 11/08/2022]
Abstract
The computational challenge of fast and reliable transition state and reaction path optimization requires new methodological strategies to maintain low cost, high accuracy, and systematic searching capabilities. The growing string method using internal coordinates has proven to be highly effective for the study of molecular, gas phase reactions, but difficulties in choosing a suitable coordinate system for periodic systems has prevented its use for surface chemistry. New developments are therefore needed, and presented herein, to handle surface reactions which include atoms with large coordination numbers that cannot be treated using standard internal coordinates. The double-ended and single-ended growing string methods are implemented using a hybrid coordinate system, then benchmarked for a test set of 43 elementary reactions occurring on surfaces. These results show that the growing string method is at least 45% faster than the widely used climbing image-nudged elastic band method, which also fails to converge in several of the test cases. Additionally, the surface growing string method has a unique single-ended search method which can move outward from an initial structure to find the intermediates, transition states, and reaction paths simultaneously. This powerful explorative feature of single ended-growing string method is demonstrated to uncover, for the first time, the mechanism for atomic layer deposition of TiN on Cu(111) surface. This reaction is found to proceed through multiple hydrogen-transfer and ligand-exchange events, while formation of H-bonds stabilizes intermediates of the reaction. Purging gaseous products out of the reaction environment is the driving force for these reactions. © 2017 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Mina Jafari
- Department of Chemistry, University of Michigan, 930 N. University Ave, Ann Arbor, Michigan, 48109
| | - Paul M Zimmerman
- Department of Chemistry, University of Michigan, 930 N. University Ave, Ann Arbor, Michigan, 48109
| |
Collapse
|
12
|
Lee J, Small DW, Epifanovsky E, Head-Gordon M. Coupled-Cluster Valence-Bond Singles and Doubles for Strongly Correlated Systems: Block-Tensor Based Implementation and Application to Oligoacenes. J Chem Theory Comput 2017; 13:602-615. [DOI: 10.1021/acs.jctc.6b01092] [Citation(s) in RCA: 50] [Impact Index Per Article: 7.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Affiliation(s)
- Joonho Lee
- Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - David W. Small
- Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Evgeny Epifanovsky
- Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Martin Head-Gordon
- Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
13
|
Jha SK, Brown K, Todde G, Subramanian G. A mechanochemical study of the effects of compression on a Diels-Alder reaction. J Chem Phys 2016; 145:074307. [DOI: 10.1063/1.4960955] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Sanjiv K. Jha
- School of Polymers and High Performance Materials, University of Southern Mississippi, Hattiesburg, Mississippi 39402, USA
| | - Katie Brown
- Department of Polymer and Fiber Engineering, Auburn University, Auburn, Alabama 36849, USA
| | - Guido Todde
- School of Polymers and High Performance Materials, University of Southern Mississippi, Hattiesburg, Mississippi 39402, USA
| | - Gopinath Subramanian
- School of Polymers and High Performance Materials, University of Southern Mississippi, Hattiesburg, Mississippi 39402, USA
| |
Collapse
|
14
|
Kearns FL, Hudson PS, Boresch S, Woodcock HL. Methods for Efficiently and Accurately Computing Quantum Mechanical Free Energies for Enzyme Catalysis. Methods Enzymol 2016; 577:75-104. [PMID: 27498635 DOI: 10.1016/bs.mie.2016.05.053] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Abstract
Enzyme activity is inherently linked to free energies of transition states, ligand binding, protonation/deprotonation, etc.; these free energies, and thus enzyme function, can be affected by residue mutations, allosterically induced conformational changes, and much more. Therefore, being able to predict free energies associated with enzymatic processes is critical to understanding and predicting their function. Free energy simulation (FES) has historically been a computational challenge as it requires both the accurate description of inter- and intramolecular interactions and adequate sampling of all relevant conformational degrees of freedom. The hybrid quantum mechanical molecular mechanical (QM/MM) framework is the current tool of choice when accurate computations of macromolecular systems are essential. Unfortunately, robust and efficient approaches that employ the high levels of computational theory needed to accurately describe many reactive processes (ie, ab initio, DFT), while also including explicit solvation effects and accounting for extensive conformational sampling are essentially nonexistent. In this chapter, we will give a brief overview of two recently developed methods that mitigate several major challenges associated with QM/MM FES: the QM non-Boltzmann Bennett's acceptance ratio method and the QM nonequilibrium work method. We will also describe usage of these methods to calculate free energies associated with (1) relative properties and (2) along reaction paths, using simple test cases with relevance to enzymes examples.
Collapse
Affiliation(s)
- F L Kearns
- University of South Florida, Tampa, FL, United States
| | - P S Hudson
- University of South Florida, Tampa, FL, United States
| | - S Boresch
- Faculty of Chemistry, University of Vienna, Vienna, Austria.
| | - H L Woodcock
- University of South Florida, Tampa, FL, United States.
| |
Collapse
|
15
|
Maximova T, Moffatt R, Ma B, Nussinov R, Shehu A. Principles and Overview of Sampling Methods for Modeling Macromolecular Structure and Dynamics. PLoS Comput Biol 2016; 12:e1004619. [PMID: 27124275 PMCID: PMC4849799 DOI: 10.1371/journal.pcbi.1004619] [Citation(s) in RCA: 132] [Impact Index Per Article: 16.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023] Open
Abstract
Investigation of macromolecular structure and dynamics is fundamental to understanding how macromolecules carry out their functions in the cell. Significant advances have been made toward this end in silico, with a growing number of computational methods proposed yearly to study and simulate various aspects of macromolecular structure and dynamics. This review aims to provide an overview of recent advances, focusing primarily on methods proposed for exploring the structure space of macromolecules in isolation and in assemblies for the purpose of characterizing equilibrium structure and dynamics. In addition to surveying recent applications that showcase current capabilities of computational methods, this review highlights state-of-the-art algorithmic techniques proposed to overcome challenges posed in silico by the disparate spatial and time scales accessed by dynamic macromolecules. This review is not meant to be exhaustive, as such an endeavor is impossible, but rather aims to balance breadth and depth of strategies for modeling macromolecular structure and dynamics for a broad audience of novices and experts.
Collapse
Affiliation(s)
- Tatiana Maximova
- Department of Computer Science, George Mason University, Fairfax, Virginia, United States of America
| | - Ryan Moffatt
- Department of Computer Science, George Mason University, Fairfax, Virginia, United States of America
| | - Buyong Ma
- Basic Science Program, Leidos Biomedical Research, Inc. Cancer and Inflammation Program, National Cancer Institute, Frederick, Maryland, United States of America
| | - Ruth Nussinov
- Basic Science Program, Leidos Biomedical Research, Inc. Cancer and Inflammation Program, National Cancer Institute, Frederick, Maryland, United States of America
- Sackler Institute of Molecular Medicine, Department of Human Genetics and Molecular Medicine, Sackler School of Medicine, Tel Aviv University, Tel Aviv, Israel
| | - Amarda Shehu
- Department of Computer Science, George Mason University, Fairfax, Virginia, United States of America
- Department of Biongineering, George Mason University, Fairfax, Virginia, United States of America
- School of Systems Biology, George Mason University, Manassas, Virginia, United States of America
| |
Collapse
|
16
|
Abstract
SUMMARYEvidence is emerging that the role of protein structure in disease needs to be rethought. Sequence mutations in proteins are often found to affect the rate at which a protein switches between structures. Modeling structural transitions in wildtype and variant proteins is central to understanding the molecular basis of disease. This paper investigates an efficient algorithmic realization of the stochastic roadmap simulation framework to model structural transitions in wildtype and variants of proteins implicated in human disorders. Our results indicate that the algorithm is able to extract useful information on the impact of mutations on protein structure and function.
Collapse
|
17
|
|
18
|
Chen C, Huang Y, Jiang X, Xiao Y. A fast tomographic method for searching the minimum free energy path. J Chem Phys 2014; 141:154109. [DOI: 10.1063/1.4897983] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022] Open
Affiliation(s)
- Changjun Chen
- Biomolecular Physics and Modeling Group, School of Physics, Huazhong University of Science and Technology, Wuhan 430074, Hubei, China
| | - Yanzhao Huang
- Biomolecular Physics and Modeling Group, School of Physics, Huazhong University of Science and Technology, Wuhan 430074, Hubei, China
| | - Xuewei Jiang
- School of Fashion, Wuhan Textile University, Wuhan 430073, Hubei, China
| | - Yi Xiao
- Biomolecular Physics and Modeling Group, School of Physics, Huazhong University of Science and Technology, Wuhan 430074, Hubei, China
| |
Collapse
|
19
|
Sun YG, Guo Y, Sun D, Zhu MC, Ding F, Liu YN, Gao EJ, Wang SJ, Xiong G, Dragutan I, Dragutan V. Palladium(II) and Platinum(II) Complexes Containing Six-Membered N-Heterocyclic Ligands: Synthesis, Characterization, Interaction with DNA, DFT Calculation, and Cytotoxicity. Eur J Inorg Chem 2014. [DOI: 10.1002/ejic.201402532] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
|
20
|
Hudson PS, White JK, Kearns FL, Hodoscek M, Boresch S, Lee Woodcock H. Efficiently computing pathway free energies: New approaches based on chain-of-replica and Non-Boltzmann Bennett reweighting schemes. Biochim Biophys Acta Gen Subj 2014; 1850:944-953. [PMID: 25239198 DOI: 10.1016/j.bbagen.2014.09.016] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2014] [Revised: 09/09/2014] [Accepted: 09/10/2014] [Indexed: 11/25/2022]
Abstract
BACKGROUND Accurately modeling condensed phase processes is one of computation's most difficult challenges. Include the possibility that conformational dynamics may be coupled to chemical reactions, where multiscale (i.e., QM/MM) methods are needed, and this task becomes even more daunting. METHODS Free energy simulations (i.e., molecular dynamics), multiscale modeling, and reweighting schemes. RESULTS Herein, we present two new approaches for mitigating the aforementioned challenges. The first is a new chain-of-replica method (off-path simulations, OPS) for computing potentials of mean force (PMFs) along an easily defined reaction coordinate. This development is coupled with a new distributed, highly-parallel replica framework (REPDstr) within the CHARMM package. Validation of these new schemes is carried out on two processes that undergo conformational changes. First is the simple torsional rotation of butane, while a much more challenging glycosidic rotation (in vacuo and solvated) is the second. Additionally, a new approach that greatly improves (i.e., possibly an order of magnitude) the efficiency of computing QM/MM PMFs is introduced and compared to standard schemes. Our efforts are grounded in the recently developed method for efficiently computing QM-based free energies (i.e., QM-Non-Boltzmann Bennett, QM-NBB). Again, we validate this new technique by computing the QM/MM PMF of butane's torsional rotation. CONCLUSIONS The OPS-REPDstr method is a promising new approach that overcomes many limitations of standard pathway simulations in CHARMM. The combination of QM-NBB with pathway techniques is very promising as it offers significant advantages over current procedures. GENERAL SIGNIFICANCE Efficiently computing potentials of mean force is a major, unresolved, area of interest. This article is part of a Special Issue entitled Recent developments of molecular dynamics.
Collapse
Affiliation(s)
- Phillip S Hudson
- Department of Chemistry, University of South Florida, 4202 E. Fowler Ave., CHE205, Tampa, FL 33620-5250, USA
| | - Justin K White
- Department of Chemistry, University of South Florida, 4202 E. Fowler Ave., CHE205, Tampa, FL 33620-5250, USA
| | - Fiona L Kearns
- Department of Chemistry, University of South Florida, 4202 E. Fowler Ave., CHE205, Tampa, FL 33620-5250, USA
| | - Milan Hodoscek
- Center for Molecular Modeling, National Institute of Chemistry, Hajdrihova 19, SI-1000 Ljubljana, Slovenia
| | - Stefan Boresch
- Department of Computational Biological Chemistry, Faculty of Chemistry, University of Vienna, Währingerstraße 17, A-1090 Vienna, Austria
| | - H Lee Woodcock
- Department of Chemistry, University of South Florida, 4202 E. Fowler Ave., CHE205, Tampa, FL 33620-5250, USA.
| |
Collapse
|
21
|
Schmidt TC, Paasche A, Grebner C, Ansorg K, Becker J, Lee W, Engels B. QM/MM investigations of organic chemistry oriented questions. Top Curr Chem (Cham) 2014; 351:25-101. [PMID: 22392477 DOI: 10.1007/128_2011_309] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/31/2023]
Abstract
About 35 years after its first suggestion, QM/MM became the standard theoretical approach to investigate enzymatic structures and processes. The success is due to the ability of QM/MM to provide an accurate atomistic picture of enzymes and related processes. This picture can even be turned into a movie if nuclei-dynamics is taken into account to describe enzymatic processes. In the field of organic chemistry, QM/MM methods are used to a much lesser extent although almost all relevant processes happen in condensed matter or are influenced by complicated interactions between substrate and catalyst. There is less importance for theoretical organic chemistry since the influence of nonpolar solvents is rather weak and the effect of polar solvents can often be accurately described by continuum approaches. Catalytic processes (homogeneous and heterogeneous) can often be reduced to truncated model systems, which are so small that pure quantum-mechanical approaches can be employed. However, since QM/MM becomes more and more efficient due to the success in software and hardware developments, it is more and more used in theoretical organic chemistry to study effects which result from the molecular nature of the environment. It is shown by many examples discussed in this review that the influence can be tremendous, even for nonpolar reactions. The importance of environmental effects in theoretical spectroscopy was already known. Due to its benefits, QM/MM can be expected to experience ongoing growth for the next decade.In the present chapter we give an overview of QM/MM developments and their importance in theoretical organic chemistry, and review applications which give impressions of the possibilities and the importance of the relevant effects. Since there is already a bunch of excellent reviews dealing with QM/MM, we will discuss fundamental ingredients and developments of QM/MM very briefly with a focus on very recent progress. For the applications we follow a similar strategy.
Collapse
Affiliation(s)
- Thomas C Schmidt
- Institut für Phys. und Theor. Chemie, Emil-Fischer-Strasse 42, Campus Hubland Nord, 97074, Würzburg, Germany
| | | | | | | | | | | | | |
Collapse
|
22
|
Jung J, Re S, Sugita Y, Ten-no S. Improved constrained optimization method for reaction-path determination in the generalized hybrid orbital quantum mechanical/molecular mechanical calculations. J Chem Phys 2013; 138:044106. [DOI: 10.1063/1.4775812] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022] Open
|
23
|
Nikodem A, Matveev AV, Zheng BX, Rösch N. Efficient Two-Step Procedures for Locating Transition States of Surface Reactions. J Chem Theory Comput 2012; 9:588-99. [DOI: 10.1021/ct300728a] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Affiliation(s)
- Astrid Nikodem
- Department Chemie and Catalysis
Research Center, Technische Universität München, 85747
Garching, Germany
| | - Alexei V. Matveev
- Department Chemie and Catalysis
Research Center, Technische Universität München, 85747
Garching, Germany
| | - Bo-Xiao Zheng
- Department Chemie and Catalysis
Research Center, Technische Universität München, 85747
Garching, Germany
- Department
of Chemistry and
Key Laboratory of Organic Optoelectronics and Molecular Engineering
of Ministry of Education, Tsinghua University, Beijing 100084, China
| | - Notker Rösch
- Department Chemie and Catalysis
Research Center, Technische Universität München, 85747
Garching, Germany
- Institute of High Performance
Computing, 1 Fusionopolis Way, #16-16 Connexis, Singapore 138632,
Singapore
| |
Collapse
|
24
|
Mallikarjun Sharada S, Zimmerman PM, Bell AT, Head-Gordon M. Automated Transition State Searches without Evaluating the Hessian. J Chem Theory Comput 2012; 8:5166-74. [PMID: 26593206 DOI: 10.1021/ct300659d] [Citation(s) in RCA: 51] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Accurate and speedy determination of transition structures (TSs) is essential for computational studies on reaction pathways, particularly when the process involves expensive electronic structure calculations. Many search algorithms require a good initial guess of the TS geometry, as well as a Hessian input that possesses a structure consistent with the desired saddle point. Among the double-ended interpolation methods for generation of the guess for the TS, the freezing string method (FSM) is proven to be far less expensive compared to its predecessor, the growing string method (GSM). In this paper, it is demonstrated that the efficiency of this technique can be improved further by replacing the conjugate gradient optimization step (FSM-CG) with a quasi-Newton line search coupled with a BFGS Hessian update (FSM-BFGS). A second crucial factor that affects the speed with which convergence to the TS is achieved is the quality and cost of the Hessian of the energy for the guessed TS. For electronic structure calculations, the cost of calculating an exact Hessian increases more rapidly with system size than the energy and gradient. Therefore, to sidestep calculation of the exact Hessian, an approximate Hessian is constructed, using the tangent direction and local curvature at the TS guess. It is demonstrated that the partitioned-rational function optimization algorithm for locating TSs with this approximate Hessian input performs at least as well as with an exact Hessian input in most test cases. The two techniques, FSM and approximate Hessian construction, therefore can significantly reduce costs associated with finding TSs.
Collapse
Affiliation(s)
- Shaama Mallikarjun Sharada
- Department of Chemical and Biomolecular Engineering and ‡Department of Chemistry, University of California, Berkeley , Berkeley, California 94720, United States
| | - Paul M Zimmerman
- Department of Chemical and Biomolecular Engineering and ‡Department of Chemistry, University of California, Berkeley , Berkeley, California 94720, United States
| | - Alexis T Bell
- Department of Chemical and Biomolecular Engineering and ‡Department of Chemistry, University of California, Berkeley , Berkeley, California 94720, United States
| | - Martin Head-Gordon
- Department of Chemical and Biomolecular Engineering and ‡Department of Chemistry, University of California, Berkeley , Berkeley, California 94720, United States
| |
Collapse
|
25
|
Tao P, Hodošček M, Larkin JD, Shao Y, Brooks BR. Comparison of Three Chain-of-States Methods: Nudged Elastic Band and Replica Path with Restraints or Constraints. J Chem Theory Comput 2012; 8:5035-5051. [PMID: 23526888 PMCID: PMC3604905 DOI: 10.1021/ct3006248] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Chain-of-state methods are becoming important tools in studying the chemical reaction mechanisms, especially for biomacromolecules. In this article, three chain-of-state methods, nudged elastic band (NEB) method and the replica path method with restraints or constraints, were tested and compared using three model systems with various sizes and at different levels of theory: alanine dipeptide isomerization, β-alanine intramolecular condensation, and the matrix metalloproteinase 2 inhibition mechanism. The levels of theory used to describe the three model systems include molecular mechanics (MM), quantum mechanics (QM), and combined quantum mechanics and molecular mechanics (QM/MM). All three methods could correctly determine a reaction path with reasonable estimation of reaction barriers in most cases. The RMSD measurement with additional weighting schemes provides practically infinite choices of reaction coordinates to describe the reaction progress. These findings demonstrate that the chain-of-state methods are powerful tools when being used carefully to generate a plausible reaction mechanism with full pathway for complex systems at an affordable computational cost.
Collapse
Affiliation(s)
- Peng Tao
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, United States
| | - Milan Hodošček
- Center for Molecular Modeling, National Institute of Chemistry, Hajdrihova 19, SI-1000 Ljubljana, Slovenia
| | - Joseph D. Larkin
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, United States
| | - Yihan Shao
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, United States
- Q-Chem Inc., 5001 Baum Boulevard, Suite 690, Pittsburgh, Pennsylvania 15213, United States
| | - Bernard R. Brooks
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, United States
| |
Collapse
|
26
|
Shiga M, Fujisaki H. A quantum generalization of intrinsic reaction coordinate using path integral centroid coordinates. J Chem Phys 2012; 136:184103. [DOI: 10.1063/1.4709723] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022] Open
|
27
|
Abstract
The evolution of many systems is dominated by rare activated events that occur on timescale
ranging from nanoseconds to the hour or more. For such systems, simulations must leave aside the
full thermal description to focus specifically on mechanisms that generate a configurational change.
We present here the activation relaxation technique (ART), an open-ended saddle point search
algorithm, and a series of recent improvements to ART nouveau and kinetic ART, an ART-based
on-the-fly off-lattice self-learning kinetic Monte Carlo method.
Collapse
|
28
|
Chaffey-Millar H, Nikodem A, Matveev AV, Krüger S, Rösch N. Improving Upon String Methods for Transition State Discovery. J Chem Theory Comput 2012; 8:777-86. [DOI: 10.1021/ct200639w] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
Affiliation(s)
- Hugh Chaffey-Millar
- Molekulare Katalyse, Department
Chemie, Technische Universität München, 85747 Garching,
Germany
| | - Astrid Nikodem
- Department Chemie & Catalysis Research Center, Technische Universität München, 85747 Garching, Germany
| | - Alexei V. Matveev
- Department Chemie & Catalysis Research Center, Technische Universität München, 85747 Garching, Germany
| | - Sven Krüger
- Department Chemie & Catalysis Research Center, Technische Universität München, 85747 Garching, Germany
| | - Notker Rösch
- Department Chemie & Catalysis Research Center, Technische Universität München, 85747 Garching, Germany
| |
Collapse
|
29
|
Behn A, Zimmerman PM, Bell AT, Head-Gordon M. Efficient exploration of reaction paths via a freezing string method. J Chem Phys 2011; 135:224108. [DOI: 10.1063/1.3664901] [Citation(s) in RCA: 128] [Impact Index Per Article: 9.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
30
|
Behn A, Zimmerman PM, Bell AT, Head-Gordon M. Incorporating Linear Synchronous Transit Interpolation into the Growing String Method: Algorithm and Applications. J Chem Theory Comput 2011; 7:4019-25. [DOI: 10.1021/ct200654u] [Citation(s) in RCA: 33] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Affiliation(s)
- Andrew Behn
- Department of Chemical and Biomolecular Engineering and ‡Department of Chemistry, University of California, Berkeley, California 94720-1462, United States
| | - Paul M. Zimmerman
- Department of Chemical and Biomolecular Engineering and ‡Department of Chemistry, University of California, Berkeley, California 94720-1462, United States
| | - Alexis T. Bell
- Department of Chemical and Biomolecular Engineering and ‡Department of Chemistry, University of California, Berkeley, California 94720-1462, United States
| | - Martin Head-Gordon
- Department of Chemical and Biomolecular Engineering and ‡Department of Chemistry, University of California, Berkeley, California 94720-1462, United States
| |
Collapse
|
31
|
Abstract
We present modifications for the method recently developed by Granot and Baer [J. Chem. Phys. 128, 184111 (2008)]. These modifications significantly enhance the efficiency and reliability of the method. In addition, we discuss some specific features of this method. These features provide important flexibilities which are crucial for a double-ended saddle point search method in order to be applicable to complex reaction mechanisms. Furthermore, it is discussed under what circumstances this methods might fail to find the transition state and remedies to avoid such situations are provided. We demonstrate the performance of the enhanced splined saddle method on several examples with increasing complexity, isomerization of ammonia, ethane and cyclopropane molecules, tautomerization of cytosine, the ring opening of cyclobutene, the Stone-Wales transformation of the C(60) fullerene, and finally rolling a small NaCl cube on NaCl(001) surface. All of these calculations are based on density functional theory. The efficiency of the method is remarkable in regard to the reduction of the total computational time.
Collapse
Affiliation(s)
- S Alireza Ghasemi
- Department of Physics, Universität Basel, Klingelbergstr. 82, 4056 Basel, Switzerland.
| | | |
Collapse
|
32
|
Chen JL, Hu WP. Theoretical Prediction on the Thermal Stability of Cyclic Ozone and Strong Oxygen Tunneling. J Am Chem Soc 2011; 133:16045-53. [DOI: 10.1021/ja203428x] [Citation(s) in RCA: 30] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Jien-Lian Chen
- Department of Chemistry and Biochemistry, National Chung Cheng University, Chia-Yi, 621, Taiwan
| | - Wei-Ping Hu
- Department of Chemistry and Biochemistry, National Chung Cheng University, Chia-Yi, 621, Taiwan
| |
Collapse
|
33
|
Fournier R, Bulusu S, Chen S, Tung J. Using swarm intelligence for finding transition states and reaction paths. J Chem Phys 2011; 135:104117. [DOI: 10.1063/1.3633515] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
34
|
Machado-Charry E, Béland LK, Caliste D, Genovese L, Deutsch T, Mousseau N, Pochet P. Optimized energy landscape exploration using the ab initio based activation-relaxation technique. J Chem Phys 2011; 135:034102. [DOI: 10.1063/1.3609924] [Citation(s) in RCA: 74] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
35
|
Affiliation(s)
- Alexis T. Bell
- Department of Chemical and Biomolecular Engineering, University of California, Berkeley, California 94720
| | - Martin Head-Gordon
- Department of Chemistry, University of California, Berkeley, California 94720; ,
| |
Collapse
|
36
|
Fuchigami S, Fujisaki H, Matsunaga Y, Kidera A. Protein Functional Motions: Basic Concepts and Computational Methodologies. ADVANCING THEORY FOR KINETICS AND DYNAMICS OF COMPLEX, MANY-DIMENSIONAL SYSTEMS: CLUSTERS AND PROTEINS 2011. [DOI: 10.1002/9781118087817.ch2] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
|
37
|
Siebert MR, Zhang J, Addepalli SV, Tantillo DJ, Hase WL. The need for enzymatic steering in abietic acid biosynthesis: gas-phase chemical dynamics simulations of carbocation rearrangements on a bifurcating potential energy surface. J Am Chem Soc 2011; 133:8335-43. [PMID: 21548620 DOI: 10.1021/ja201730y] [Citation(s) in RCA: 66] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Abietic acid, a constituent of pine resin, is naturally derived from abietadiene --a process that requires four enzymes: one (abietadiene synthase) for conversion of the acyclic, achiral geranylgeranyl diphosphate to the polycyclic, chiral abietadiene (a complex process involving the copalyl diphosphate intermediate) and then three to oxidize a single methyl group of abietadiene to the corresponding carboxylic acid. In previous work (Nature Chem.2009, 1, 384), electronic structure calculations on carbocation rearrangements leading to abietadienyl cation revealed an interesting potential energy surface with a bifurcating reaction pathway (two transition-state structures connected directly with no intervening minimum), which links two products--one natural and one not yet isolated from Nature. Herein we describe direct dynamics simulations of the key step in the formation of abietadiene (in the gas phase and in the absence of the enzyme). The simulations reveal that abietadiene synthase must intervene in order to produce abietadiene selectively, in essence steering this reaction to avoid the generation of byproducts with different molecular architectures.
Collapse
Affiliation(s)
- Matthew R Siebert
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409, USA
| | | | | | | | | |
Collapse
|
38
|
Schlegel HB. Geometry optimization. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2011. [DOI: 10.1002/wcms.34] [Citation(s) in RCA: 186] [Impact Index Per Article: 14.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
|
39
|
Liu YN, Liang WZ, Sang XG, Huo YQ, Sze-to L, Yung KF, Liu XX. Syntheses, characterizations and theoretical calculations of rhodium(III) 1,2-naphthoquinone-1-oxime complexes. Inorganica Chim Acta 2010. [DOI: 10.1016/j.ica.2009.12.033] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
40
|
Aguilar-Mogas A, Giménez X, Bofill JM. Implementation of an algorithm based on the Runge-Kutta-Fehlberg technique and the potential energy as a reaction coordinate to locate intrinsic reaction paths. J Comput Chem 2010; 31:2510-25. [DOI: 10.1002/jcc.21539] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
|