1
|
Tao Y, Zou W, Nanayakkara S, Freindorf M, Kraka E. A revised formulation of the generalized subsystem vibrational analysis (GSVA). Theor Chem Acc 2021; 140:31. [PMID: 33716564 PMCID: PMC7942689 DOI: 10.1007/s00214-021-02727-y] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/09/2020] [Accepted: 02/08/2021] [Indexed: 11/30/2022]
Abstract
In this work, a simplified formulation of our recently developed generalized subsystem vibrational analysis (GSVA) for obtaining intrinsic fragmental vibrations (J Chem Theory Comput 14:2558, 2018) is presented. In contrast to the earlier implementation, which requires the explicit definition of a non-redundant set of internal coordinate parameters to be constructed for the subsystem, the new implementation circumvents this process by employing massless Eckart conditions to the subsystem fragment paired with a Gram-Schmidt orthogonalization to span the same internal vibration space indirectly. This revised version of GSVA (rev-GSVA) can be applied to equilibrium structure as well as transition state structure, and it has been incorporated into the open-source package UniMoVib (https://github.com/zorkzou/UniMoVib). Supplementary Information The online version contains supplementary material available at 10.1007/s00214-021-02727-y.
Collapse
Affiliation(s)
- Yunwen Tao
- Department of Chemistry, Southern Methodist University, 3215 Daniel Ave, Dallas, TX 75275-0314 USA
| | - Wenli Zou
- Institute of Modern Physics, Northwest University, and Shaanxi Key Laboratory for Theoretical Physics Frontiers, Xi'an, 710127 Shaanxi People's Republic of China
| | - Sadisha Nanayakkara
- Department of Chemistry, Southern Methodist University, 3215 Daniel Ave, Dallas, TX 75275-0314 USA
| | - Marek Freindorf
- Department of Chemistry, Southern Methodist University, 3215 Daniel Ave, Dallas, TX 75275-0314 USA
| | - Elfi Kraka
- Department of Chemistry, Southern Methodist University, 3215 Daniel Ave, Dallas, TX 75275-0314 USA
| |
Collapse
|
2
|
Kraka E, Zou W, Tao Y. Decoding chemical information from vibrational spectroscopy data: Local vibrational mode theory. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2020. [DOI: 10.1002/wcms.1480] [Citation(s) in RCA: 43] [Impact Index Per Article: 10.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Affiliation(s)
- Elfi Kraka
- Department of Chemistry Southern Methodist University Dallas Texas USA
| | - Wenli Zou
- Institute of Modern Physics Northwest University and Shaanxi Key Laboratory for Theoretical Physics Frontiers, Xi'an Shaanxi PR China
| | - Yunwen Tao
- Department of Chemistry Southern Methodist University Dallas Texas USA
| |
Collapse
|
3
|
Tao Y, Tian C, Verma N, Zou W, Wang C, Cremer D, Kraka E. Recovering Intrinsic Fragmental Vibrations Using the Generalized Subsystem Vibrational Analysis. J Chem Theory Comput 2018; 14:2558-2569. [DOI: 10.1021/acs.jctc.7b01171] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Yunwen Tao
- Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, Texas 75275-0314, United States
| | - Chuan Tian
- Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, United States
| | - Niraj Verma
- Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, Texas 75275-0314, United States
| | - Wenli Zou
- Institute of Modern Physics, Northwest University, Xi’an, Shaanxi 710127, P. R. China
| | - Chao Wang
- Institute of Nuclear Energy Safety Technology, Chinese Academy of Sciences, Hefei, Anhui 230031, P. R. China
| | - Dieter Cremer
- Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, Texas 75275-0314, United States
| | - Elfi Kraka
- Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, Texas 75275-0314, United States
| |
Collapse
|
4
|
Hahn S. Effective representation of amide III, II, I, and A modes on local vibrational modes: Analysis of ab initio quantum calculation results. J Chem Phys 2016; 145:164113. [PMID: 27802648 DOI: 10.1063/1.4965958] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The Hamiltonian matrix for the first excited vibrational states of a protein can be effectively represented by local vibrational modes constituting amide III, II, I, and A modes to simulate various vibrational spectra. Methods for obtaining the Hamiltonian matrix from ab initio quantum calculation results are discussed, where the methods consist of three steps: selection of local vibrational mode coordinates, calculation of a reduced Hessian matrix, and extraction of the Hamiltonian matrix from the Hessian matrix. We introduce several methods for each step. The methods were assessed based on the density functional theory calculation results of 24 oligopeptides with four different peptide lengths and six different secondary structures. The completeness of a Hamiltonian matrix represented in the reduced local mode space is improved by adopting a specific atom group for each amide mode and reducing the effect of ignored local modes. The calculation results are also compared to previous models using C=O stretching vibration and transition dipole couplings. We found that local electric transition dipole moments of the amide modes are mainly bound on the local peptide planes. Their direction and magnitude are well conserved except amide A modes, which show large variation. Contrary to amide I modes, the vibrational coupling constants of amide III, II, and A modes obtained by analysis of a dipeptide are not transferable to oligopeptides with the same secondary conformation because coupling constants are affected by the surrounding atomic environment.
Collapse
Affiliation(s)
- Seungsoo Hahn
- Department of Chemistry, Chung-Ang University, 84 Heukseok-ro, Dongjak-gu 156-756, Seoul, South Korea
| |
Collapse
|
5
|
Jose KVJ, Beckett D, Raghavachari K. Vibrational Circular Dichroism Spectra for Large Molecules through Molecules-in-Molecules Fragment-Based Approach. J Chem Theory Comput 2015; 11:4238-47. [DOI: 10.1021/acs.jctc.5b00647] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- K. V. Jovan Jose
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| | - Daniel Beckett
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| | - Krishnan Raghavachari
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| |
Collapse
|
6
|
Nguyen CM, Reyniers MF, Marin GB. Adsorption thermodynamics of C1–C4 alcohols in H-FAU, H-MOR, H-ZSM-5, and H-ZSM-22. J Catal 2015. [DOI: 10.1016/j.jcat.2014.11.013] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
|
7
|
Correcting for the free energy costs of bond or angle constraints in molecular dynamics simulations. Biochim Biophys Acta Gen Subj 2014; 1850:932-943. [PMID: 25218695 DOI: 10.1016/j.bbagen.2014.09.001] [Citation(s) in RCA: 37] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2014] [Revised: 08/28/2014] [Accepted: 09/01/2014] [Indexed: 11/22/2022]
Abstract
BACKGROUND Free energy simulations are an important tool in the arsenal of computational biophysics, allowing the calculation of thermodynamic properties of binding or enzymatic reactions. This paper introduces methods to increase the accuracy and precision of free energy calculations by calculating the free energy costs of constraints during post-processing. The primary purpose of employing constraints for these free energy methods is to increase the phase space overlap between ensembles, which is required for accuracy and convergence. METHODS The free energy costs of applying or removing constraints are calculated as additional explicit steps in the free energy cycle. The new techniques focus on hard degrees of freedom and use both gradients and Hessian estimation. Enthalpy, vibrational entropy, and Jacobian free energy terms are considered. RESULTS We demonstrate the utility of this method with simple classical systems involving harmonic and anharmonic oscillators, four-atomic benchmark systems, an alchemical mutation of ethane to methanol, and free energy simulations between alanine and serine. The errors for the analytical test cases are all below 0.0007kcal/mol, and the accuracy of the free energy results of ethane to methanol is improved from 0.15 to 0.04kcal/mol. For the alanine to serine case, the phase space overlaps of the unconstrained simulations range between 0.15 and 0.9%. The introduction of constraints increases the overlap up to 2.05%. On average, the overlap increases by 94% relative to the unconstrained value and precision is doubled. CONCLUSIONS The approach reduces errors arising from constraints by about an order of magnitude. Free energy simulations benefit from the use of constraints through enhanced convergence and higher precision. GENERAL SIGNIFICANCE The primary utility of this approach is to calculate free energies for systems with disparate energy surfaces and bonded terms, especially in multi-scale molecular mechanics/quantum mechanics simulations. This article is part of a Special Issue entitled Recent developments of molecular dynamics.
Collapse
|
8
|
Sweet JC, Nowling RJ, Cickovski T, Sweet CR, Pande VS, Izaguirre JA. Long Timestep Molecular Dynamics on the Graphical Processing Unit. J Chem Theory Comput 2013; 9:3267-3281. [PMID: 24436689 DOI: 10.1021/ct400331r] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
Molecular dynamics (MD) simulations now play a key role in many areas of theoretical chemistry, biology, physics, and materials science. In many cases, such calculations are significantly limited by the massive amount of computer time needed to perform calculations of interest. Herein, we present Long Timestep Molecular Dynamics (LTMD), a method to significantly speed MD simulations. In particular, we discuss new methods to calculate the needed terms in LTMD as well as issues germane to a GPU implementation. The resulting code, implemented in the OpenMM MD library, can achieve a significant 6-fold speed increase, leading to MD simulations on the order of 5 μs/day using implicit solvent models.
Collapse
Affiliation(s)
- James C Sweet
- Department of Computer Science and Engineering, University of Notre Dame, Notre Dame, IN 46556, USA, Department of Computer Science, Eckerd College, Saint Petersburg, FL 33712, USA, and Department of Chemistry, Stanford University, Stanford, CA 94305, USA
| | - Ronald J Nowling
- Department of Computer Science and Engineering, University of Notre Dame, Notre Dame, IN 46556, USA, Department of Computer Science, Eckerd College, Saint Petersburg, FL 33712, USA, and Department of Chemistry, Stanford University, Stanford, CA 94305, USA
| | - Trevor Cickovski
- Department of Computer Science and Engineering, University of Notre Dame, Notre Dame, IN 46556, USA, Department of Computer Science, Eckerd College, Saint Petersburg, FL 33712, USA, and Department of Chemistry, Stanford University, Stanford, CA 94305, USA
| | - Christopher R Sweet
- Department of Computer Science and Engineering, University of Notre Dame, Notre Dame, IN 46556, USA, Department of Computer Science, Eckerd College, Saint Petersburg, FL 33712, USA, and Department of Chemistry, Stanford University, Stanford, CA 94305, USA
| | - Vijay S Pande
- Department of Computer Science and Engineering, University of Notre Dame, Notre Dame, IN 46556, USA, Department of Computer Science, Eckerd College, Saint Petersburg, FL 33712, USA, and Department of Chemistry, Stanford University, Stanford, CA 94305, USA
| | - Jesús A Izaguirre
- Department of Computer Science and Engineering, University of Notre Dame, Notre Dame, IN 46556, USA, Department of Computer Science, Eckerd College, Saint Petersburg, FL 33712, USA, and Department of Chemistry, Stanford University, Stanford, CA 94305, USA
| |
Collapse
|
9
|
Nakata H, Nagata T, Fedorov DG, Yokojima S, Kitaura K, Nakamura S. Analytic second derivatives of the energy in the fragment molecular orbital method. J Chem Phys 2013; 138:164103. [DOI: 10.1063/1.4800990] [Citation(s) in RCA: 37] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
10
|
Ghysels A, Miller BT, Pickard FC, Brooks BR. Comparing normal modes across different models and scales: Hessian reductionversuscoarse-graining. J Comput Chem 2012; 33:2250-75. [DOI: 10.1002/jcc.23076] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/23/2011] [Revised: 05/09/2012] [Accepted: 06/24/2012] [Indexed: 12/24/2022]
|
11
|
Hemelsoet K, Ghysels A, Mores D, De Wispelaere K, Van Speybroeck V, Weckhuysen BM, Waroquier M. Experimental and theoretical IR study of methanol and ethanol conversion over H-SAPO-34. Catal Today 2011. [DOI: 10.1016/j.cattod.2011.05.040] [Citation(s) in RCA: 37] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
12
|
Bieler NS, Haag MP, Jacob CR, Reiher M. Analysis of the Cartesian Tensor Transfer Method for Calculating Vibrational Spectra of Polypeptides. J Chem Theory Comput 2011; 7:1867-81. [DOI: 10.1021/ct2001478] [Citation(s) in RCA: 42] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Noah S. Bieler
- ETH Zurich, Laboratorium für Physikalische Chemie, Wolfgang-Pauli-Strasse 10, 8093 Zurich, Switzerland
| | - Moritz P. Haag
- ETH Zurich, Laboratorium für Physikalische Chemie, Wolfgang-Pauli-Strasse 10, 8093 Zurich, Switzerland
| | - Christoph R. Jacob
- Karlsruhe Institute of Technology (KIT), Center for Functional Nanostructures, Wolfgang-Gaede-Str. 1a, 76131 Karlsruhe, Germany
| | - Markus Reiher
- ETH Zurich, Laboratorium für Physikalische Chemie, Wolfgang-Pauli-Strasse 10, 8093 Zurich, Switzerland
| |
Collapse
|
13
|
De Moor BA, Ghysels A, Reyniers MF, Van Speybroeck V, Waroquier M, Marin GB. Normal Mode Analysis in Zeolites: Toward an Efficient Calculation of Adsorption Entropies. J Chem Theory Comput 2011; 7:1090-101. [DOI: 10.1021/ct1005505] [Citation(s) in RCA: 85] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Affiliation(s)
- Bart A. De Moor
- Laboratory for Chemical Technology, Ghent University, Krijgslaan 281 S5, 9000 Ghent, Belgium
| | - An Ghysels
- Center for Molecular Modeling, Ghent University, Technologiepark 903, 9052 Zwijnaarde, Belgium
| | | | | | - Michel Waroquier
- Center for Molecular Modeling, Ghent University, Technologiepark 903, 9052 Zwijnaarde, Belgium
| | - Guy B. Marin
- Laboratory for Chemical Technology, Ghent University, Krijgslaan 281 S5, 9000 Ghent, Belgium
| |
Collapse
|
14
|
Ghysels A, Woodcock HL, Larkin JD, Miller BT, Shao Y, Kong J, Neck DV, Speybroeck VV, Waroquier M, Brooks BR. Efficient Calculation of QM/MM Frequencies with the Mobile Block Hessian. J Chem Theory Comput 2011; 7:496-514. [DOI: 10.1021/ct100473f] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Affiliation(s)
- An Ghysels
- Center for Molecular Modeling, Ghent University, Technologiepark 903, 9052 Zwijnaarde, Belgium, Department of Chemistry, University of South Florida, 4202 E. Fowler Avenue, CHE 205, Tampa, Florida 33620-5240, United States, Laboratory of Computational Biology, National Heart Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, United States, and Q-Chem Inc., 5001 Baum Blvd, Suite 690, Pittsburgh, Pennsylvania 15213, United States
| | - H. Lee Woodcock
- Center for Molecular Modeling, Ghent University, Technologiepark 903, 9052 Zwijnaarde, Belgium, Department of Chemistry, University of South Florida, 4202 E. Fowler Avenue, CHE 205, Tampa, Florida 33620-5240, United States, Laboratory of Computational Biology, National Heart Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, United States, and Q-Chem Inc., 5001 Baum Blvd, Suite 690, Pittsburgh, Pennsylvania 15213, United States
| | - Joseph D. Larkin
- Center for Molecular Modeling, Ghent University, Technologiepark 903, 9052 Zwijnaarde, Belgium, Department of Chemistry, University of South Florida, 4202 E. Fowler Avenue, CHE 205, Tampa, Florida 33620-5240, United States, Laboratory of Computational Biology, National Heart Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, United States, and Q-Chem Inc., 5001 Baum Blvd, Suite 690, Pittsburgh, Pennsylvania 15213, United States
| | - Benjamin T. Miller
- Center for Molecular Modeling, Ghent University, Technologiepark 903, 9052 Zwijnaarde, Belgium, Department of Chemistry, University of South Florida, 4202 E. Fowler Avenue, CHE 205, Tampa, Florida 33620-5240, United States, Laboratory of Computational Biology, National Heart Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, United States, and Q-Chem Inc., 5001 Baum Blvd, Suite 690, Pittsburgh, Pennsylvania 15213, United States
| | - Yihan Shao
- Center for Molecular Modeling, Ghent University, Technologiepark 903, 9052 Zwijnaarde, Belgium, Department of Chemistry, University of South Florida, 4202 E. Fowler Avenue, CHE 205, Tampa, Florida 33620-5240, United States, Laboratory of Computational Biology, National Heart Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, United States, and Q-Chem Inc., 5001 Baum Blvd, Suite 690, Pittsburgh, Pennsylvania 15213, United States
| | - Jing Kong
- Center for Molecular Modeling, Ghent University, Technologiepark 903, 9052 Zwijnaarde, Belgium, Department of Chemistry, University of South Florida, 4202 E. Fowler Avenue, CHE 205, Tampa, Florida 33620-5240, United States, Laboratory of Computational Biology, National Heart Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, United States, and Q-Chem Inc., 5001 Baum Blvd, Suite 690, Pittsburgh, Pennsylvania 15213, United States
| | - Dimitri Van Neck
- Center for Molecular Modeling, Ghent University, Technologiepark 903, 9052 Zwijnaarde, Belgium, Department of Chemistry, University of South Florida, 4202 E. Fowler Avenue, CHE 205, Tampa, Florida 33620-5240, United States, Laboratory of Computational Biology, National Heart Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, United States, and Q-Chem Inc., 5001 Baum Blvd, Suite 690, Pittsburgh, Pennsylvania 15213, United States
| | - Veronique Van Speybroeck
- Center for Molecular Modeling, Ghent University, Technologiepark 903, 9052 Zwijnaarde, Belgium, Department of Chemistry, University of South Florida, 4202 E. Fowler Avenue, CHE 205, Tampa, Florida 33620-5240, United States, Laboratory of Computational Biology, National Heart Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, United States, and Q-Chem Inc., 5001 Baum Blvd, Suite 690, Pittsburgh, Pennsylvania 15213, United States
| | - Michel Waroquier
- Center for Molecular Modeling, Ghent University, Technologiepark 903, 9052 Zwijnaarde, Belgium, Department of Chemistry, University of South Florida, 4202 E. Fowler Avenue, CHE 205, Tampa, Florida 33620-5240, United States, Laboratory of Computational Biology, National Heart Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, United States, and Q-Chem Inc., 5001 Baum Blvd, Suite 690, Pittsburgh, Pennsylvania 15213, United States
| | - Bernard R. Brooks
- Center for Molecular Modeling, Ghent University, Technologiepark 903, 9052 Zwijnaarde, Belgium, Department of Chemistry, University of South Florida, 4202 E. Fowler Avenue, CHE 205, Tampa, Florida 33620-5240, United States, Laboratory of Computational Biology, National Heart Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, United States, and Q-Chem Inc., 5001 Baum Blvd, Suite 690, Pittsburgh, Pennsylvania 15213, United States
| |
Collapse
|
15
|
Ghysels A, Verstraelen T, Hemelsoet K, Waroquier M, Van Speybroeck V. TAMkin: A Versatile Package for Vibrational Analysis and Chemical Kinetics. J Chem Inf Model 2010; 50:1736-50. [DOI: 10.1021/ci100099g] [Citation(s) in RCA: 131] [Impact Index Per Article: 9.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- An Ghysels
- Center for Molecular Modeling, QCMM Alliance Ghent-Brussels, Ghent University, Technologiepark 903, 9052 Zwijnaarde, Belgium
| | - Toon Verstraelen
- Center for Molecular Modeling, QCMM Alliance Ghent-Brussels, Ghent University, Technologiepark 903, 9052 Zwijnaarde, Belgium
| | - Karen Hemelsoet
- Center for Molecular Modeling, QCMM Alliance Ghent-Brussels, Ghent University, Technologiepark 903, 9052 Zwijnaarde, Belgium
| | - Michel Waroquier
- Center for Molecular Modeling, QCMM Alliance Ghent-Brussels, Ghent University, Technologiepark 903, 9052 Zwijnaarde, Belgium
| | - Veronique Van Speybroeck
- Center for Molecular Modeling, QCMM Alliance Ghent-Brussels, Ghent University, Technologiepark 903, 9052 Zwijnaarde, Belgium
| |
Collapse
|
16
|
Burger SK, Ayers PW. Quasi-Newton parallel geometry optimization methods. J Chem Phys 2010; 133:034116. [DOI: 10.1063/1.3455719] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022] Open
|
17
|
Ghysels A, Van Speybroeck V, Pauwels E, Catak S, Brooks BR, Van Neck D, Waroquier M. Comparative study of various normal mode analysis techniques based on partial Hessians. J Comput Chem 2010; 31:994-1007. [PMID: 19813181 DOI: 10.1002/jcc.21386] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Abstract
Standard normal mode analysis becomes problematic for complex molecular systems, as a result of both the high computational cost and the excessive amount of information when the full Hessian matrix is used. Several partial Hessian methods have been proposed in the literature, yielding approximate normal modes. These methods aim at reducing the computational load and/or calculating only the relevant normal modes of interest in a specific application. Each method has its own (dis)advantages and application field but guidelines for the most suitable choice are lacking. We have investigated several partial Hessian methods, including the Partial Hessian Vibrational Analysis (PHVA), the Mobile Block Hessian (MBH), and the Vibrational Subsystem Analysis (VSA). In this article, we focus on the benefits and drawbacks of these methods, in terms of the reproduction of localized modes, collective modes, and the performance in partially optimized structures. We find that the PHVA is suitable for describing localized modes, that the MBH not only reproduces localized and global modes but also serves as an analysis tool of the spectrum, and that the VSA is mostly useful for the reproduction of the low frequency spectrum. These guidelines are illustrated with the reproduction of the localized amine-stretch, the spectrum of quinine and a bis-cinchona derivative, and the low frequency modes of the LAO binding protein.
Collapse
Affiliation(s)
- An Ghysels
- Center for Molecular Modeling, Ghent University, Proeftuinstraat 86, 9000 Gent, Belgium.
| | | | | | | | | | | | | |
Collapse
|