1
|
Pitsevich GA, Malevich AE, Kamnev AA. A convenient set of vibrational coordinates for 2D calculation of the tunneling splittings of the ground state and some excited vibrational states for the inversion motion in H 3O +, H 3O -, and H 3O. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2023; 296:122660. [PMID: 37043918 DOI: 10.1016/j.saa.2023.122660] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/05/2022] [Revised: 03/10/2023] [Accepted: 03/19/2023] [Indexed: 06/19/2023]
Abstract
Splitting of the ground state and some excited symmetric bending vibrational states due to inversion tunneling of the oxygen atom in the H3O+, H3O- ions and in the H3O radical are analyzed by numerically solving the vibrational Schrödinger equation of restricted (2D) dimensionality. As two vibrational coordinates, we used 1) the distance of the oxygen atom from the plane of a regular triangle formed by three hydrogen atoms and 2) a symmetry coordinate composed of three distances between chemically non-bonded hydrogen atoms. The kinetic energy operator in this case takes the simplest form. The 2D potential energy surface (PES) in the given coordinates was calculated for H3O+ at the CCSD(T)/aug-cc-pVTZ and CCSD(T)-F12/cc-pVTZ-F12 levels of theory. The same 2D PES for the H3O- anion and H3O radical were calculated at the CCSD(T)/aug-cc-pVQZ, CCSD(T)/d-aug-cc-pVQZ and UCCSD(T)/aug-cc-pVQZ, UCCSD(T)/d-aug-cc-pVQZ levels of theory, respectively. The tunneling splittings were calculated for the cations H316O+, D316O+, T316O+, H318O+, D318O+, T318O+. The tunneling splittings for the H3O-, D3O-, T3O- anions and H3O, D3O, T3O radicals were calculated for the first time. The results of calculations demonstrate good agreement with experimental values of the tunneling splittings in the ground state and in some excited vibrational states of the H316O+ and D316O+ cations.
Collapse
Affiliation(s)
- George A Pitsevich
- Department of Physical Optics and Applied Informatics, Faculty of Physics, Belarusian State University, Nezavisimosti Ave., 4, 220030 Minsk, Belarus.
| | - Alex E Malevich
- Department of Differential Equations and System Analysis, Faculty of Mechanics and Mathematics, Belarusian State University, Nezavisimosti Ave., 4, 220030 Minsk, Belarus
| | - Alexander A Kamnev
- Laboratory of Biochemistry, Institute of Biochemistry and Physiology of Plants and Microorganisms - Subdivision of the Federal State Budgetary Research Institution Saratov Federal Scientific Center of the Russian Academy of Sciences, Prosp. Entuziastov, 13, 410049 Saratov, Russia.
| |
Collapse
|
2
|
Sahu N, Richardson JO, Berger R. Instanton calculations of tunneling splittings in chiral molecules. J Comput Chem 2021; 42:210-221. [PMID: 33259074 DOI: 10.1002/jcc.26447] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/29/2020] [Revised: 08/20/2020] [Accepted: 08/27/2020] [Indexed: 11/11/2022]
Abstract
We report the ground state tunneling splittings (ΔE± ) of a number of axially chiral molecules using the ring-polymer instanton (RPI) method (J. Chem. Phys., 2011, 134, 054109). The list includes isotopomers of hydrogen dichalcogenides H2 X2 (X = O, S, Se, Te, and Po), hydrogen thioperoxide HSOH and dichlorodisulfane S2 Cl2 . Ab initio electronic-structure calculations have been performed on the level of second-order Møller-Plesset perturbation (MP2) theory either with split-valance basis sets or augmented correlation-consistent basis sets on H, O, S, and Cl atoms. Energy-consistent pseudopotential and corresponding triple zeta basis sets of the Stuttgart group are used on Se, Te, and Po atoms. The results are further improved using single point calculations performed at the coupled cluster level with iterative singles and doubles and perturbative triples amplitudes. When available for comparison, our computed values of ΔE± are found to lie within the same order of magnitude as values reported in the literature, although RPI also provides predictions for H2 Po2 and S2 Cl2 , which have not previously been directly calculated. Since RPI is a single-shot method which does not require detailed prior knowledge of the optimal tunneling path, it offers an effective way for estimating the tunneling dynamics of more complex chiral molecules, and especially those with small tunneling splittings.
Collapse
Affiliation(s)
- Nityananda Sahu
- Fachbereich Chemie, Theoretische Chemie, Philipps Universität Marburg, Marburg, Germany
| | | | - Robert Berger
- Fachbereich Chemie, Theoretische Chemie, Philipps Universität Marburg, Marburg, Germany
| |
Collapse
|
3
|
Eraković M, Cvitaš MT. Tunneling splittings of vibrationally excited states using general instanton paths. J Chem Phys 2020; 153:134106. [PMID: 33032414 DOI: 10.1063/5.0024210] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023] Open
Abstract
A multidimensional semiclassical method for calculating tunneling splittings in vibrationally excited states of molecules using Cartesian coordinates is developed. It is an extension of the theory by Mil'nikov and Nakamura [J. Chem. Phys. 122, 124311 (2005)] to asymmetric paths that are necessary for calculating tunneling splitting patterns in multi-well systems, such as water clusters. Additionally, new terms are introduced in the description of the semiclassical wavefunction that drastically improves the splitting estimates for certain systems. The method is based on the instanton theory and builds the semiclassical wavefunction of the vibrationally excited states from the ground-state instanton wavefunction along the minimum action path and its harmonic neighborhood. The splittings of excited states are thus obtained at a negligible added numerical effort. The cost is concentrated, as for the ground-state splittings, in the instanton path optimization and the hessian evaluation along the path. The method can thus be applied without modification to many mid-sized molecules in full dimensionality and in combination with on-the-fly evaluation of electronic potentials. The tests were performed on several model potentials and on the water dimer.
Collapse
Affiliation(s)
- Mihael Eraković
- Department of Physical Chemistry, Ruder Bošković Institute, Bijenička Cesta 54, 10000 Zagreb, Croatia
| | - Marko T Cvitaš
- Department of Physical Chemistry, Ruder Bošković Institute, Bijenička Cesta 54, 10000 Zagreb, Croatia
| |
Collapse
|
4
|
Nejad A, Crittenden DL. On the separability of large-amplitude motions in anharmonic frequency calculations. Phys Chem Chem Phys 2020; 22:20588-20601. [PMID: 32966420 DOI: 10.1039/d0cp03515g] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Nuclear vibrational theories based upon the Watson Hamiltonian are ubiquitous in quantum chemistry, but are generally unable to model systems in which the wavefunction can delocalise over multiple energy minima, i.e. molecules that have low-energy torsion and inversion barriers. In a 2019 Chemical Reviews article, Puzzarini et al. note that a common workaround is to simply decouple these problematic modes from all other vibrations in the system during anharmonic frequency calculations. They also point out that this approximation can be "ill-suited", but do not quantify the errors introduced. In this work, we present the first systematic investigation into how separating out or constraining torsion and inversion vibrations within potential energy surface (PES) expansions affects the accuracy of computed fundamental wavenumbers for the remaining vibrational modes, using a test set of 19 tetratomic molecules for which high quality analytic potential energy surfaces and fully-coupled anharmonic reference fundamental frequencies are available. We find that the most effective and efficient strategy is to remove the mode in question from the PES expansion entirely. This introduces errors of up to +10 cm-1 in stretching fundamentals that would otherwise couple to the dropped mode, and ±5 cm-1 in all other fundamentals. These errors are approximately commensurate with, but not necessarily additional to, errors due to the choice of electronic structure model used in constructing spectroscopically accurate PES.
Collapse
Affiliation(s)
- Arman Nejad
- Institute of Physical Chemistry, University of Göttingen, Tammannstr. 6, D-37077 Göttingen, Germany.
| | - Deborah L Crittenden
- School of Physical and Chemical Sciences, University of Canterbury, Christchurch, New Zealand
| |
Collapse
|
5
|
Werner HJ, Knowles PJ, Manby FR, Black JA, Doll K, Heßelmann A, Kats D, Köhn A, Korona T, Kreplin DA, Ma Q, Miller TF, Mitrushchenkov A, Peterson KA, Polyak I, Rauhut G, Sibaev M. The Molpro quantum chemistry package. J Chem Phys 2020; 152:144107. [PMID: 32295355 DOI: 10.1063/5.0005081] [Citation(s) in RCA: 441] [Impact Index Per Article: 110.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Abstract
Molpro is a general purpose quantum chemistry software package with a long development history. It was originally focused on accurate wavefunction calculations for small molecules but now has many additional distinctive capabilities that include, inter alia, local correlation approximations combined with explicit correlation, highly efficient implementations of single-reference correlation methods, robust and efficient multireference methods for large molecules, projection embedding, and anharmonic vibrational spectra. In addition to conventional input-file specification of calculations, Molpro calculations can now be specified and analyzed via a new graphical user interface and through a Python framework.
Collapse
Affiliation(s)
- Hans-Joachim Werner
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Peter J Knowles
- School of Chemistry, Cardiff University, Main Building, Park Place, Cardiff CF10 3AT, United Kingdom
| | - Frederick R Manby
- School of Chemistry, University of Bristol, Cantock's Close, Bristol BS8 1TS, United Kingdom
| | - Joshua A Black
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Klaus Doll
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Andreas Heßelmann
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Daniel Kats
- Max-Planck Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | - Andreas Köhn
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Tatiana Korona
- Faculty of Chemistry, University of Warsaw, L. Pasteura 1 St., 02-093 Warsaw, Poland
| | - David A Kreplin
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Qianli Ma
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Thomas F Miller
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | | | - Kirk A Peterson
- Department of Chemistry, Washington State University, Pullman, Washington 99164-4630, USA
| | - Iakov Polyak
- School of Chemistry, Cardiff University, Main Building, Park Place, Cardiff CF10 3AT, United Kingdom
| | - Guntram Rauhut
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Marat Sibaev
- School of Chemistry, Cardiff University, Main Building, Park Place, Cardiff CF10 3AT, United Kingdom
| |
Collapse
|
6
|
Petrenko T, Rauhut G. Account of non-Condon effects in time-independent Raman wavefunction theory: Calculation of the S 1 ← S 0 vibronic absorption spectrum of formaldehyde. J Chem Phys 2020; 152:114109. [DOI: 10.1063/5.0003272] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Taras Petrenko
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Guntram Rauhut
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| |
Collapse
|
7
|
Meisner J, Hallmen PP, Kästner J, Rauhut G. Vibrational analysis of methyl cation—Rare gas atom complexes: CH3+—Rg (Rg = He, Ne, Ar, Kr). J Chem Phys 2019; 150:084306. [DOI: 10.1063/1.5084100] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Jan Meisner
- Institute for Theoretical Chemistry, University of Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, USA
- SLAC National Accelerator Laboratory, Menlo Park, California 94025, USA
| | - Philipp P. Hallmen
- Institute for Theoretical Chemistry, University of Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
| | - Johannes Kästner
- Institute for Theoretical Chemistry, University of Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
| | - Guntram Rauhut
- Institute for Theoretical Chemistry, University of Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
| |
Collapse
|
8
|
Cvitaš MT. Quadratic String Method for Locating Instantons in Tunneling Splitting Calculations. J Chem Theory Comput 2018; 14:1487-1500. [DOI: 10.1021/acs.jctc.7b00881] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Marko T. Cvitaš
- Department of Physical Chemistry, Ruđer Bošković Institute, Bijenička Cesta 54, 10000 Zagreb, Croatia
| |
Collapse
|
9
|
Sibaev M, Crittenden DL. An efficient and numerically stable procedure for generating sextic force fields in normal mode coordinates. J Chem Phys 2016; 144:214107. [DOI: 10.1063/1.4953080] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- M. Sibaev
- Department of Chemistry, University of Canterbury, Christchurch, New Zealand
| | - D. L. Crittenden
- Department of Chemistry, University of Canterbury, Christchurch, New Zealand
| |
Collapse
|
10
|
Cvitaš MT, Althorpe SC. Locating Instantons in Calculations of Tunneling Splittings: The Test Case of Malonaldehyde. J Chem Theory Comput 2016; 12:787-803. [PMID: 26756608 DOI: 10.1021/acs.jctc.5b01073] [Citation(s) in RCA: 32] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The recently developed ring-polymer instanton (RPI) method [J. Chem. Phys. 2011, 134, 054109] is an efficient technique for calculating approximate tunneling splittings in high-dimensional molecular systems. The key step is locating the instanton tunneling-path at zero temperature. Here, we show that techniques previously designed for locating instantons in finite-temperature rate calculations can be adapted to the RPI method, where they become extremely efficient, reducing the number of potential energy calls by 2 orders of magnitude. We investigate one technique that employs variable time steps to minimize the action integral, and two that employ equally spaced position steps to minimize the abbreviated (i.e., Jacobi) action integral, using respectively the nudged elastic band (NEB) and string methods. We recommend use of the latter because it is parameter-free, but all three methods give comparable efficiency savings. Having located the instanton pathway, we then interpolate the instanton path onto a fine grid of imaginary time points, allowing us to compute the fluctuation prefactor. The crucial modification needed to the original finite-temperature algorithms is to allow the end points of the zero-temperature instanton path to describe overall rotations, which is done using a standard quaternion algorithm. These approaches will allow the RPI method to be combined effectively with expensive potential energy surfaces or on-the-fly electronic structure methods.
Collapse
Affiliation(s)
- Marko T Cvitaš
- Department of Physical Chemistry, Ruđer Bošković Institute , Bijenička Cesta 54, 10000 Zagreb, Croatia
| | - Stuart C Althorpe
- Department of Chemistry, University of Cambridge , Lensfield Road, Cambridge CB2 1EW, United Kingdom
| |
Collapse
|
11
|
Fletcher TL, Davie SJ, Popelier PLA. Prediction of Intramolecular Polarization of Aromatic Amino Acids Using Kriging Machine Learning. J Chem Theory Comput 2015; 10:3708-19. [PMID: 26588516 DOI: 10.1021/ct500416k] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Present computing power enables novel ways of modeling polarization. Here we show that the machine learning method kriging accurately captures the way the electron density of a topological atom responds to a change in the positions of the surrounding atoms. The success of this method is demonstrated on the four aromatic amino acids histidine, phenylalanine, tryptophan, and tyrosine. A new technique of varying training set sizes to vastly reduce training times while maintaining accuracy is described and applied to each amino acid. Each amino acid has its geometry distorted via normal modes of vibration over all local energy minima in the Ramachandran map. These geometries are then used to train the kriging models. Total electrostatic energies predicted by the kriging models for previously unseen geometries are compared to the true energies, yielding mean absolute errors of 2.9, 5.1, 4.2, and 2.8 kJ mol(-1) for histidine, phenylalanine, tryptophan, and tyrosine, respectively.
Collapse
Affiliation(s)
- Timothy L Fletcher
- Manchester Institute of Biotechnology (MIB) , 131 Princess Street, Manchester, M1 7DN, Great Britain.,School of Chemistry, University of Manchester , Oxford Road, Manchester, M13 9PL, Great Britain
| | - Stuart J Davie
- Manchester Institute of Biotechnology (MIB) , 131 Princess Street, Manchester, M1 7DN, Great Britain.,School of Chemistry, University of Manchester , Oxford Road, Manchester, M13 9PL, Great Britain
| | - Paul L A Popelier
- Manchester Institute of Biotechnology (MIB) , 131 Princess Street, Manchester, M1 7DN, Great Britain.,School of Chemistry, University of Manchester , Oxford Road, Manchester, M13 9PL, Great Britain
| |
Collapse
|
12
|
Isotope effects of ammonia umbrella flip using semiclassical instanton calculations based on discretized path integrals. Chem Phys Lett 2015. [DOI: 10.1016/j.cplett.2015.05.072] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
|
13
|
Kawatsu T, Miura S. Efficient algorithms for semiclassical instanton calculations based on discretized path integrals. J Chem Phys 2014; 141:024101. [DOI: 10.1063/1.4885437] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022] Open
Affiliation(s)
- Tsutomu Kawatsu
- Institute for Molecular Science, National Institute of Natural Science, 38 Nishigonaka, Myodaiji, Okazaki 222-8585, Japan
- School of Mathematics and Physics, Kanazawa University, Kanazawa 920-1192, Japan
| | - Shinichi Miura
- School of Mathematics and Physics, Kanazawa University, Kanazawa 920-1192, Japan
| |
Collapse
|