51
|
Schmitz G, Klinting EL, Christiansen O. A Gaussian process regression adaptive density guided approach for potential energy surface construction. J Chem Phys 2020; 153:064105. [DOI: 10.1063/5.0015344] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Gunnar Schmitz
- Department of Chemistry, Aarhus Universitet, DK-8000 Aarhus, Denmark
| | | | - Ove Christiansen
- Department of Chemistry, Aarhus Universitet, DK-8000 Aarhus, Denmark
| |
Collapse
|
52
|
Ollitrault PJ, Baiardi A, Reiher M, Tavernelli I. Hardware efficient quantum algorithms for vibrational structure calculations. Chem Sci 2020; 11:6842-6855. [PMID: 32874524 PMCID: PMC7448527 DOI: 10.1039/d0sc01908a] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/02/2020] [Accepted: 06/02/2020] [Indexed: 11/21/2022] Open
Abstract
We introduce a framework for the calculation of ground and excited state energies of bosonic systems suitable for near-term quantum devices and apply it to molecular vibrational anharmonic Hamiltonians. Our method supports generic reference modal bases and Hamiltonian representations, including the ones that are routinely used in classical vibrational structure calculations. We test different parametrizations of the vibrational wavefunction, which can be encoded in quantum hardware, based either on heuristic circuits or on the bosonic Unitary Coupled Cluster Ansatz. In particular, we define a novel compact heuristic circuit and demonstrate that it provides a good compromise in terms of circuit depth, optimization costs, and accuracy. We evaluate the requirements, number of qubits and circuit depth, for the calculation of vibrational energies on quantum hardware and compare them with state-of-the-art classical vibrational structure algorithms for molecules with up to seven atoms.
Collapse
Affiliation(s)
- Pauline J Ollitrault
- IBM Quantum , IBM Research - Zurich , Säumerstrasse 4 , 8803 Rüschlikon , Switzerland .
- Laboratory of Physical Chemistry , ETH Zürich , Vladimir-Prelog-Weg 2 , 8093 Zürich , Switzerland .
| | - Alberto Baiardi
- Laboratory of Physical Chemistry , ETH Zürich , Vladimir-Prelog-Weg 2 , 8093 Zürich , Switzerland .
| | - Markus Reiher
- Laboratory of Physical Chemistry , ETH Zürich , Vladimir-Prelog-Weg 2 , 8093 Zürich , Switzerland .
| | - Ivano Tavernelli
- IBM Quantum , IBM Research - Zurich , Säumerstrasse 4 , 8803 Rüschlikon , Switzerland .
| |
Collapse
|
53
|
Erfort S, Tschöpe M, Rauhut G. Toward a fully automated calculation of rovibrational infrared intensities for semi-rigid polyatomic molecules. J Chem Phys 2020; 152:244104. [DOI: 10.1063/5.0011832] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Sebastian Erfort
- Institute for Theoretical Chemistry, University of Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Martin Tschöpe
- Institute for Theoretical Chemistry, University of Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Guntram Rauhut
- Institute for Theoretical Chemistry, University of Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| |
Collapse
|
54
|
Artiukhin DG, Klinting EL, König C, Christiansen O. Adaptive density-guided approach to double incremental potential energy surface construction. J Chem Phys 2020; 152:194105. [PMID: 33687258 DOI: 10.1063/5.0004686] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022] Open
Abstract
We present a combination of the recently developed double incremental expansion of potential energy surfaces with the well-established adaptive density-guided approach to grid construction. This unique methodology is based on the use of an incremental expansion for potential energy surfaces, known as n-mode expansion; an incremental many-body representation of the electronic energy; and an efficient vibrational density-guided approach to automated determination of grid dimensions and granularity. The reliability of the method is validated calculating potential energy surfaces and obtaining fundamental excitation energies for three moderate-size chain-like molecular systems. The use of our methodology leads to considerable computational savings for potential energy surface construction compared to standard approaches while maintaining a high level of accuracy in the resulting potential energy surfaces. Additional investigations indicate that our method can be applied to covalently bound and strongly interacting molecular systems, even though these cases are known to be very unfavorable for fragmentation schemes. We therefore conclude that the presented methodology is a robust and flexible approach to potential energy surface construction, which introduces considerable computational savings without compromising the accuracy of vibrational spectra calculations.
Collapse
Affiliation(s)
- Denis G Artiukhin
- Department of Chemistry, Aarhus Universitet, DK-8000 Aarhus, Denmark
| | | | - Carolin König
- Institute of Physical Chemistry, Christian-Albrechts-University Kiel, Max-Eyth-Straße 1, D-24118 Kiel, Germany
| | - Ove Christiansen
- Department of Chemistry, Aarhus Universitet, DK-8000 Aarhus, Denmark
| |
Collapse
|
55
|
Klinting EL, Lauvergnat D, Christiansen O. Vibrational Coupled Cluster Computations in Polyspherical Coordinates with the Exact Analytical Kinetic Energy Operator. J Chem Theory Comput 2020; 16:4505-4520. [DOI: 10.1021/acs.jctc.0c00261] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Affiliation(s)
| | - David Lauvergnat
- Université Paris-Saclay, CNRS, Institut de Chimie Physique, UMR8000, 91405 Orsay, France
| | - Ove Christiansen
- Department of Chemistry, Aarhus University, DK-8000 Aarhus C, Denmark
| |
Collapse
|
56
|
Erfort S, Tschöpe M, Rauhut G, Zeng X, Tew DP. Ab initio calculation of rovibrational states for non-degenerate double-well potentials: cis-trans isomerization of HOPO. J Chem Phys 2020; 152:174306. [PMID: 32384829 DOI: 10.1063/5.0005497] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The rovibrational spectra of metaphosphorous acid, HOPO, and its deuterated isotopologue have been studied by vibrational configuration interaction calculations, relying on the internal coordinate path Hamiltonian and the Watson Hamiltonian. Tunneling effects for the overtones of the torsional mode, which gives rise to the cis-trans isomerization, and its rovibrational transitions have been investigated in detail. Due to strong matrix effects, comparison with experimental data is hindered, and thus, the calculations provide accurate estimates for the fundamental modes of these species.
Collapse
Affiliation(s)
- Sebastian Erfort
- Institute for Theoretical Chemistry, University of Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Martin Tschöpe
- Institute for Theoretical Chemistry, University of Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Guntram Rauhut
- Institute for Theoretical Chemistry, University of Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Xiaoqing Zeng
- Department of Chemistry, Fudan University, Shanghai 200433, China
| | - David P Tew
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany
| |
Collapse
|
57
|
|
58
|
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: 514] [Impact Index Per Article: 128.5] [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
|
59
|
Le Bris V, Odunlami M, Bégué D, Baraille I, Coulaud O. Using computed infrared intensities for the reduction of vibrational configuration interaction bases. Phys Chem Chem Phys 2020; 22:7021-7030. [PMID: 32191243 DOI: 10.1039/d0cp00593b] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The Adaptive Vibrational Configuration Interaction (A-VCI) algorithm is an iterative process that computes the anharmonic spectrum of a molecule using nested bases to discretize the Hamiltonian operator. For large molecular systems, the size of the discretization space and the computation time quickly become prohibitive. It is therefore necessary to develop new methods to further limit the number of basis functions. Most of the time, the interpretation of an experimental infrared spectrum does not require the calculation of all eigenvalues but only those corresponding to vibrational states with significant intensity. In this paper, a technique that uses infrared intensities is introduced to select a subset of eigenvalues to be precisely calculated. Thus, we build smaller nested bases and reduce both the memory footprint and the computational time. We validate the advantages of this new approach on a well-studied 7-atom molecular system (C2H4O), and we apply it on a larger 10-atom molecule (C4H4N2).
Collapse
Affiliation(s)
- Vincent Le Bris
- Universite de Pau et des Pays de l'Adour, E2S UPPA, CNRS, IPREM, Institut des Sciences Analytiques et de Physicochimie pour l'Environnement et les Matériaux, UMR5254, Pau, France.
| | - Marc Odunlami
- Universite de Pau et des Pays de l'Adour, E2S UPPA, CNRS, IPREM, Institut des Sciences Analytiques et de Physicochimie pour l'Environnement et les Matériaux, UMR5254, Pau, France.
| | - Didier Bégué
- Universite de Pau et des Pays de l'Adour, E2S UPPA, CNRS, IPREM, Institut des Sciences Analytiques et de Physicochimie pour l'Environnement et les Matériaux, UMR5254, Pau, France.
| | - Isabelle Baraille
- Universite de Pau et des Pays de l'Adour, E2S UPPA, CNRS, IPREM, Institut des Sciences Analytiques et de Physicochimie pour l'Environnement et les Matériaux, UMR5254, Pau, France.
| | - Olivier Coulaud
- HiePACS project-team, Inria Bordeaux Sud-Ouest, 200, avenue de la Vieille Tour, F-33405 Talence Cedex, France
| |
Collapse
|
60
|
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
|
61
|
Klinting EL, Christiansen O, König C. Toward Accurate Theoretical Vibrational Spectra: A Case Study for Maleimide. J Phys Chem A 2020; 124:2616-2627. [DOI: 10.1021/acs.jpca.9b11915] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Affiliation(s)
| | - Ove Christiansen
- Department of Chemistry, Aarhus University, DK-8000 Aarhus C, Denmark
| | - Carolin König
- Institute of Physical Chemistry, Kiel University, Max-Eyth-Straße 1, D-24118 Kiel, Germany
| |
Collapse
|
62
|
Menezes da Silva VH, Ornellas FR. Characterizing structures, energetics, and spectra of species on the 1,3[H, C, As] potential energy surfaces: A high-level theoretical contribution. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2020; 226:117578. [PMID: 31670030 DOI: 10.1016/j.saa.2019.117578] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/08/2019] [Revised: 09/23/2019] [Accepted: 09/28/2019] [Indexed: 06/10/2023]
Abstract
The ground and the low lying electronic states of structures on the 1,3[H, C, As] potential energy surfaces were investigated with the highly correlated theoretical approaches CCSD(T), CCSD(T)-F12b, and CASSCF/MRCI along with the series of correlation consistent (aug-cc-pVnZ, n = D, T, Q, 5) basis sets. Energetic and spectroscopic parameters were obtained at the complete basis set limit, and the effect of core-valence correlation on these properties evaluated. Fundamental frequencies were also computed with the variational configuration interaction (VCI) approach. Heats of formation at 0 and 298.15 K were estimated for HCAs and CH, AsH, CAs, and HCAs, as well as the calculation of ionization potentials for HCAs. Comparisons of the present results with literature ones for the systems HCN/HNC, HCP/HPC highlight similarities and differences among these systems. Altogether, this investigation provides a very reliable characterization of the species on the surfaces and should guide future experimental studies on these systems.
Collapse
Affiliation(s)
- Vitor H Menezes da Silva
- Universidade de São Paulo, Instituto de Química, Departamento de Química Fundamental, Av. Lineu Prestes, 748, São Paulo, São Paulo, 05508-000, Brazil
| | - Fernando R Ornellas
- Universidade de São Paulo, Instituto de Química, Departamento de Química Fundamental, Av. Lineu Prestes, 748, São Paulo, São Paulo, 05508-000, Brazil.
| |
Collapse
|
63
|
Dinu DF, Podewitz M, Grothe H, Loerting T, Liedl KR. Decomposing anharmonicity and mode-coupling from matrix effects in the IR spectra of matrix-isolated carbon dioxide and methane. Phys Chem Chem Phys 2020; 22:17932-17947. [DOI: 10.1039/d0cp02121k] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
A combined experimental and computational approach revealed similarities and differences in the vibrational signature of matrix-isolated carbon dioxide and methane.
Collapse
Affiliation(s)
- Dennis F. Dinu
- Institute of General
- Inorganic and Theoretical Chemistry
- University of Innsbruck
- A-6020 Innsbruck
- Austria
| | - Maren Podewitz
- Institute of General
- Inorganic and Theoretical Chemistry
- University of Innsbruck
- A-6020 Innsbruck
- Austria
| | - Hinrich Grothe
- Institute of Materials Chemistry
- TU Wien
- A-1060 Vienna
- Austria
| | - Thomas Loerting
- Institute of Physical Chemistry
- University of Innsbruck
- A-6020 Innsbruck
- Austria
| | - Klaus R. Liedl
- Institute of General
- Inorganic and Theoretical Chemistry
- University of Innsbruck
- A-6020 Innsbruck
- Austria
| |
Collapse
|
64
|
Bergmann TG, Welzel MO, Jacob CR. Towards theoretical spectroscopy with error bars: systematic quantification of the structural sensitivity of calculated spectra. Chem Sci 2019; 11:1862-1877. [PMID: 34123280 PMCID: PMC8148348 DOI: 10.1039/c9sc05103a] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Abstract
Molecular spectra calculated with quantum-chemical methods are subject to a number of uncertainties (e.g., errors introduced by the computational methodology) that hamper the direct comparison of experiment and computation. Judging these uncertainties is crucial for drawing reliable conclusions from the interplay of experimental and theoretical spectroscopy, but largely relies on subjective judgment. Here, we explore the application of methods from uncertainty quantification to theoretical spectroscopy, with the ultimate goal of providing systematic error bars for calculated spectra. As a first target, we consider distortions of the underlying molecular structure as one important source of uncertainty. We show that by performing a principal component analysis, the most influential collective distortions can be identified, which allows for the construction of surrogate models that are amenable to a statistical analysis of the propagation of uncertainties in the molecular structure to uncertainties in the calculated spectrum. This is applied to the calculation of X-ray emission spectra of iron carbonyl complexes, of the electronic excitation spectrum of a coumarin dye, and of the infrared spectrum of alanine. We show that with our approach it becomes possible to obtain error bars for calculated spectra that account for uncertainties in the molecular structure. This is an important first step towards systematically quantifying other relevant sources of uncertainty in theoretical spectroscopy. Uncertainty quantification is applied in theoretical spectroscopy to obtain error bars accounting for the structural sensitivity of calculated spectra.![]()
Collapse
Affiliation(s)
- Tobias G Bergmann
- Technische Universität Braunschweig, Institute of Physical and Theoretical Chemistry Gaußstraße 17 38106 Braunschweig Germany
| | - Michael O Welzel
- Technische Universität Braunschweig, Institute of Physical and Theoretical Chemistry Gaußstraße 17 38106 Braunschweig Germany
| | - Christoph R Jacob
- Technische Universität Braunschweig, Institute of Physical and Theoretical Chemistry Gaußstraße 17 38106 Braunschweig Germany
| |
Collapse
|
65
|
Ku J, Kamath A, Carrington T, Manzhos S. Machine Learning Optimization of the Collocation Point Set for Solving the Kohn–Sham Equation. J Phys Chem A 2019; 123:10631-10642. [DOI: 10.1021/acs.jpca.9b09732] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Jonas Ku
- Department of Mechanical Engineering, National University of Singapore, Block EA #07-08, 9 Engineering Drive 1, Singapore 117576, Singapore
| | - Aditya Kamath
- Department of Mechanical Engineering, National University of Singapore, Block EA #07-08, 9 Engineering Drive 1, Singapore 117576, Singapore
| | - Tucker Carrington
- Chemistry Department, Queen’s University, Kingston, Ontario K7L 3N6, Canada
| | - Sergei Manzhos
- Centre Énergie Matériaux Télécommunications, Institut National de la Recherche Scientifique, 1650, boulevard Lionel-Boulet, Varennes QC J3X 1S2, Canada
| |
Collapse
|
66
|
Dinu DF, Podewitz M, Grothe H, Liedl KR, Loerting T. Toward Elimination of Discrepancies between Theory and Experiment: Anharmonic Rotational-Vibrational Spectrum of Water in Solid Noble Gas Matrices. J Phys Chem A 2019; 123:8234-8242. [PMID: 31433184 PMCID: PMC6767348 DOI: 10.1021/acs.jpca.9b07221] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/03/2022]
Abstract
![]()
Rotational–vibrational
spectroscopy of water in solid noble
gas matrices has been studied for many decades. Despite that, discrepancies
persist in the literature about the assignment of specific bands.
We tackle the involved rotational–vibrational spectrum of the
water isotopologues H216O, HD16O,
and D216O with an unprecedented combination
of experimental high-resolution matrix isolation infrared (MI-IR)
spectroscopy and computational anharmonic vibrational spectroscopy
by vibrational configuration interaction (VCI) on high-level ab initio
potential energy surfaces. With VCI, the average deviation to gas-phase
experiments is reduced from >100 to ≈1 cm–1 when compared to harmonic vibrational spectra. Discrepancies between
MI-IR and VCI spectra are identified as matrix effects rather than
missing anharmonicity in the theoretical approach. Matrix effects
are small in Ne (≈1.5 cm–1) and a bit larger
in Ar (≈10 cm–1). Controversial assignments
in Ne MI-IR spectra are resolved, for example, concerning the ν3 triad in HDO. We identify new transitions, for example, the
ν2 101 ← 110 transition
in D2O and H2O or the ν3 000 ← 101 transition in D2O, and
reassign bands, for example, the band at 3718.9 cm–1 that is newly assigned as the 110 ← 111 transition. The identification and solution of discrepancies for
a well-studied benchmark system such as water prove the importance
of an iterative and one-hand combination of theory and experiment
in the field of high-resolution infrared spectroscopy of single molecules.
As the computational costs involved in the VCI approach are reasonably
low, such combined experimental/theoretical studies can be extended
to molecules larger than triatomics.
Collapse
Affiliation(s)
- Dennis F Dinu
- Institute of Materials Chemistry , TU Wien , A-1060 Vienna , Austria
| | | | - Hinrich Grothe
- Institute of Materials Chemistry , TU Wien , A-1060 Vienna , Austria
| | | | | |
Collapse
|
67
|
Panadés-Barrueta RL, Martínez-Núñez E, Peláez D. Specific Reaction Parameter Multigrid POTFIT (SRP-MGPF): Automatic Generation of Sum-of-Products Form Potential Energy Surfaces for Quantum Dynamical Calculations. Front Chem 2019; 7:576. [PMID: 31475138 PMCID: PMC6702682 DOI: 10.3389/fchem.2019.00576] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/04/2019] [Accepted: 07/30/2019] [Indexed: 11/13/2022] Open
Abstract
We present Specific Reaction Parameter Multigrid POTFIT (SRP-MGPF), an automated methodology for the generation of global potential energy surfaces (PES), molecular properties surfaces, e.g., dipole, polarizabilities, etc. using a single random geometry as input. The SRP-MGPF workflow integrates: (i) a fully automated procedure for the global topographical characterization of a (intermolecular) PES based on the Transition State Search Using Chemical Dynamical Simulations (TSSCDS) family of methods;i (ii) the global optimization of the parameters of a semiempirical Hamiltonian in order to reproduce a given level of electronic structure theory; and (iii) a tensor decomposition algorithm which turns the resulting SRP-PES into sum of products (Tucker) form with the Multigrid POTFIT algorithm. The latter is necessary for quantum dynamical studies within the Multiconfiguration Time-Dependent Hartree (MCTDH) quantum dynamics method. To demonstrate our approach, we have applied our methodology to the cis-trans isomerization reaction in HONO in full dimensionality (6D). The resulting SRP-PES has been validated through the computation of classical on-the-fly dynamical calculations as well as calculations of the lowest vibrational eigenstates of HONO as well as high-energy wavepacket propagations.
Collapse
Affiliation(s)
- Ramón L. Panadés-Barrueta
- Laboratoire de Physique des Lasers, Atomes et Molécules (PhLAM), Université de Lille, Villeneuve-d'Ascq, France
| | - Emilio Martínez-Núñez
- Departamento de Química Física, Facultade de Química, Universidade de Santiago de Compostela, Santiago de Compostela, Spain
| | - Daniel Peláez
- Laboratoire de Physique des Lasers, Atomes et Molécules (PhLAM), Université de Lille, Villeneuve-d'Ascq, France
| |
Collapse
|
68
|
Schmitz G, Godtliebsen IH, Christiansen O. Machine learning for potential energy surfaces: An extensive database and assessment of methods. J Chem Phys 2019; 150:244113. [DOI: 10.1063/1.5100141] [Citation(s) in RCA: 32] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022] Open
Affiliation(s)
- Gunnar Schmitz
- Department of Chemistry, Aarhus Universitet, DK-8000 Aarhus, Denmark
| | | | - Ove Christiansen
- Department of Chemistry, Aarhus Universitet, DK-8000 Aarhus, Denmark
| |
Collapse
|
69
|
Ziegler B, Rauhut G. Localized Normal Coordinates in Accurate Vibrational Structure Calculations: Benchmarks for Small Molecules. J Chem Theory Comput 2019; 15:4187-4196. [DOI: 10.1021/acs.jctc.9b00381] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Benjamin Ziegler
- Institute for Theoretical Chemistry, University of Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Guntram Rauhut
- Institute for Theoretical Chemistry, University of Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| |
Collapse
|
70
|
Welch BK, Dawes R, Bross DH, Ruscic B. An Automated Thermochemistry Protocol Based on Explicitly Correlated Coupled-Cluster Theory: The Methyl and Ethyl Peroxy Families. J Phys Chem A 2019; 123:5673-5682. [DOI: 10.1021/acs.jpca.9b04381] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Bradley K. Welch
- Department of Chemistry, Missouri University of Science and Technology, Rolla, Missouri 65409, United States
| | - Richard Dawes
- Department of Chemistry, Missouri University of Science and Technology, Rolla, Missouri 65409, United States
| | - David H. Bross
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Argonne, Illinois 60439, United States
| | - Branko Ruscic
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Argonne, Illinois 60439, United States
- Consortium for Advanced Science and Engineering, University of Chicago, Chicago, Illinois 60637, United States
| |
Collapse
|
71
|
Erba A, Maul J, Ferrabone M, Dovesi R, Rérat M, Carbonnière P. Anharmonic Vibrational States of Solids from DFT Calculations. Part II: Implementation of the VSCF and VCI Methods. J Chem Theory Comput 2019; 15:3766-3777. [DOI: 10.1021/acs.jctc.9b00294] [Citation(s) in RCA: 26] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Alessandro Erba
- Dipartimento di Chimica, Universitá di Torino, via Giuria 5, 10125 Torino, Italy
| | - Jefferson Maul
- Dipartimento di Chimica, Universitá di Torino, via Giuria 5, 10125 Torino, Italy
| | - Matteo Ferrabone
- Dipartimento di Chimica, Universitá di Torino, via Giuria 5, 10125 Torino, Italy
| | - Roberto Dovesi
- Dipartimento di Chimica, Universitá di Torino, via Giuria 5, 10125 Torino, Italy
| | - Michel Rérat
- IPREM, Université de Pau et des Pays de l’Adour, IPREM-CAPT UMR CNRS 5254, Hélioparc Pau Pyrénées, 2 avenue du Président Angot, 64053 PAU CEDEX 9, Pau, France
| | - Philippe Carbonnière
- IPREM, Université de Pau et des Pays de l’Adour, IPREM-CAPT UMR CNRS 5254, Hélioparc Pau Pyrénées, 2 avenue du Président Angot, 64053 PAU CEDEX 9, Pau, France
| |
Collapse
|
72
|
Kuenzer U, Hofer TS. A four-dimensional Numerov approach and its application to the vibrational eigenstates of linear triatomic molecules – The interplay between anharmonicity and inter-mode coupling. Chem Phys 2019. [DOI: 10.1016/j.chemphys.2019.01.007] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
|
73
|
Ziegler B, Rauhut G. Accurate Vibrational Configuration Interaction Calculations on Diborane and Its Isotopologues. J Phys Chem A 2019; 123:3367-3373. [DOI: 10.1021/acs.jpca.9b01604] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Benjamin Ziegler
- University of Stuttgart, Institute for Theoretical Chemistry, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Guntram Rauhut
- University of Stuttgart, Institute for Theoretical Chemistry, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| |
Collapse
|
74
|
Panek PT, Hoeske AA, Jacob CR. On the choice of coordinates in anharmonic theoretical vibrational spectroscopy: Harmonic vs. anharmonic coupling in vibrational configuration interaction. J Chem Phys 2019; 150:054107. [PMID: 30736699 DOI: 10.1063/1.5083186] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023] Open
Abstract
By a suitable choice of coordinates, the computational effort required for calculations of anharmonic vibrational spectra can be reduced significantly. By using suitable localized-mode coordinates obtained from an orthogonal transformation of the conventionally used normal-mode coordinates, anharmonic couplings between modes can be significantly reduced. However, such a transformation introduces harmonic couplings between the localized modes. To elucidate the role of these harmonic couplings, we consider the vibrational self-consistent field (VSCF)/vibrational configuration interaction (VCI) calculations for both few-mode model systems and for ethene as a molecular test case. We show that large harmonic couplings can result in significant errors in localized-mode L-VSCF/L-VCI calculations and study the convergence with respect to the size of the VCI excitation space. To further elucidate the errors introduced by harmonic couplings, we discuss the connection between L-VSCF/L-VCI and vibrational exciton models. With the help of our results, we propose an algorithm for the localization of normal modes in suitable subsets that are chosen to strictly limit the errors introduced by the harmonic couplings while still leading to maximally localized modes.
Collapse
Affiliation(s)
- Paweł T Panek
- Technische Universität Braunschweig, Institute of Physical and Theoretical Chemistry, Gaußstraße 17, 38106 Braunschweig, Germany
| | - Adrian A Hoeske
- Technische Universität Braunschweig, Institute of Physical and Theoretical Chemistry, Gaußstraße 17, 38106 Braunschweig, Germany
| | - Christoph R Jacob
- Technische Universität Braunschweig, Institute of Physical and Theoretical Chemistry, Gaußstraße 17, 38106 Braunschweig, Germany
| |
Collapse
|
75
|
Yagi K, Yamada K, Kobayashi C, Sugita Y. Anharmonic Vibrational Analysis of Biomolecules and Solvated Molecules Using Hybrid QM/MM Computations. J Chem Theory Comput 2019; 15:1924-1938. [PMID: 30730746 PMCID: PMC8864611 DOI: 10.1021/acs.jctc.8b01193] [Citation(s) in RCA: 22] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
![]()
Quantum
mechanics/molecular mechanics (QM/MM) calculations are
applied for anharmonic vibrational analyses of biomolecules and solvated
molecules. The QM/MM method is implemented into a molecular dynamics
(MD) program, GENESIS, by interfacing with external electronic structure
programs. Following the geometry optimization and the harmonic normal-mode
analysis based on a partial Hessian, the anharmonic potential energy
surface (PES) is generated from QM/MM energies and gradients calculated
at grid points. The PES is used for vibrational self-consistent field
(VSCF) and post-VSCF calculations to compute the vibrational spectrum.
The method is first applied to a phosphate ion in solution. With both
the ion and neighboring water molecules taken as a QM region, IR spectra
of representative hydration structures are calculated by the second-order
vibrational quasi-degenerate perturbation theory (VQDPT2) at the level
of B3LYP/cc-pVTZ and TIP3P force field. A weight-average of IR spectra
over the structures reproduces the experimental spectrum with a mean
absolute deviation of 16 cm–1. Then, the method
is applied to an enzyme, P450 nitric oxide reductase (P450nor), with
the NO molecule bound to a ferric (FeIII) heme. Starting
from snapshot structures obtained from MD simulations of P450nor in
solution, QM/MM calculations have been carried out at the level of
B3LYP-D3/def2-SVP(D). The spin state of FeIII(NO) is likely
a closed-shell singlet state based on a ratio of N–O and Fe–NO
stretching frequencies (νN–O and νFe–NO) calculated for closed- and open-shell singlet
states. The calculated νN–O and νFe–NO overestimate the experimental ones by 120 and
75 cm–1, respectively. The electronic structure
and solvation of FeIII(NO) affect the structure around
the heme of P450nor leading to an increase in νN–O and νFe–NO.
Collapse
Affiliation(s)
- Kiyoshi Yagi
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
| | - Kenta Yamada
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
| | - Chigusa Kobayashi
- Computational Biophysics Research Team, RIKEN Center for Computational Science, 7-1-26 Minatojima-Minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
| | - Yuji Sugita
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
- Computational Biophysics Research Team, RIKEN Center for Computational Science, 7-1-26 Minatojima-Minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
- Laboratory for Biomolecular Function Simulation, RIKEN Center for Biosystems Dynamics Research, 1-6-5 Minatojima-Minamimachi,
Chuo-ku, Kobe, Hyogo 650-0047, Japan
| |
Collapse
|
76
|
Tan JA, Kuo JL. Multilevel Approach for Direct VSCF/VCI MULTIMODE Calculations with Applications to Large “Zundel” Cations. J Chem Theory Comput 2018; 14:6405-6416. [DOI: 10.1021/acs.jctc.8b00679] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Jake A. Tan
- Institute of Atomic and Molecular Sciences, Academia Sinica, No. 1 Roosevelt Road, Section 4, Taipei 10617, Taiwan (ROC)
| | - Jer-Lai Kuo
- Institute of Atomic and Molecular Sciences, Academia Sinica, No. 1 Roosevelt Road, Section 4, Taipei 10617, Taiwan (ROC)
| |
Collapse
|
77
|
Spieler S, Dinu DF, Jusko P, Bastian B, Simpson M, Podewitz M, Liedl KR, Schlemmer S, Brünken S, Wester R. Low frequency vibrational anharmonicity and nuclear spin effects of Cl -(H 2) and Cl -(D 2). J Chem Phys 2018; 149:174310. [PMID: 30409015 DOI: 10.1063/1.5049680] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Low frequency combination bands of 35Cl-(H2) and 35Cl-(D2) have been measured in the region between 600 and 1100 cm-1 by infrared predissociation spectroscopy in a cryogenic 22-pole ion trap using a free electron laser at the FELIX Laboratory as a tunable light source. The 35Cl-(H2) (35Cl-(D2)) spectrum contains three bands at 773 cm-1 (620 cm-1), 889 cm-1 (692 cm-1), and 978 cm-1 (750 cm-1) with decreasing intensity toward higher photon energies. Comparison of the experimentally determined transition frequencies with anharmonic vibrational self-consistent field and vibrational configuration interaction calculations suggests the assignment of the combination bands v1 + v2, 2v1 + v2, and 3v1 + v2 for 35Cl-(H2) and 2v1 + v2, 3v1 + v2, and 4v1 + v2 for 35Cl-(D2), where v1 is the 35Cl-⋯H2 stretching fundamental and v2 is the Cl-(H2) bend. The observed asymmetric temperature dependent line shape of the v1 + v2 transition can be modeled by a series of ∑+-∏ ro-vibrational transitions, when substantially decreasing the rotational constant in the vibrationally excited state by 35%. The spectrum of 35Cl-(D2) shows a splitting of 7 cm-1 for the strongest band which can be attributed to the tunneling of the ortho/para states of D2.
Collapse
Affiliation(s)
- Steffen Spieler
- Institut für Ionenphysik und Angewandte Physik, Universität Innsbruck, Technikerstraße 25, 6020 Innsbruck, Austria
| | - Dennis F Dinu
- Institut für Allgemeine, Anorganische und Theoretische Chemie, und Center for Molecular Biosciences Innsbruck (CMBI), Universität Innsbruck, Innrain 80-82, A-6020 Innsbruck, Austria
| | - Pavol Jusko
- Institut de Recherche en Astrophysique et Planétologie (IRAP), Université de Toulouse (UPS), CNRS, CNES, 9 Av. du Colonel Roche, 31028 Toulouse Cedex 4, France
| | - Björn Bastian
- Institut für Ionenphysik und Angewandte Physik, Universität Innsbruck, Technikerstraße 25, 6020 Innsbruck, Austria
| | - Malcolm Simpson
- Institut für Ionenphysik und Angewandte Physik, Universität Innsbruck, Technikerstraße 25, 6020 Innsbruck, Austria
| | - Maren Podewitz
- Institut für Allgemeine, Anorganische und Theoretische Chemie, und Center for Molecular Biosciences Innsbruck (CMBI), Universität Innsbruck, Innrain 80-82, A-6020 Innsbruck, Austria
| | - Klaus R Liedl
- Institut für Allgemeine, Anorganische und Theoretische Chemie, und Center for Molecular Biosciences Innsbruck (CMBI), Universität Innsbruck, Innrain 80-82, A-6020 Innsbruck, Austria
| | - Stephan Schlemmer
- I. Physikalisches Institut, Universität zu Köln, Zülpicher Str. 77, 50937 Köln, Germany
| | - Sandra Brünken
- Institute for Molecules and Materials, FELIX Laboratory, Radboud University, Toernooiveld 7c, 6525 ED Nijmegen, The Netherlands
| | - Roland Wester
- Institut für Ionenphysik und Angewandte Physik, Universität Innsbruck, Technikerstraße 25, 6020 Innsbruck, Austria
| |
Collapse
|
78
|
Zhao J, Wang Q, Yu W, Huang T, Wang X. M-S Multiple Bond in HMSH, H 2MS, and HMS Molecules (M = B, Al, Ga): Matrix Infrared Spectra and Theoretical Calculations. J Phys Chem A 2018; 122:8626-8635. [PMID: 30335387 DOI: 10.1021/acs.jpca.8b08266] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Reaction products of B, Al, and Ga atoms with H2S have been identified in solid argon using matrix isolation infrared absorption spectroscopy. The results show that the ground state B atom reaction with H2S gives the H2BS molecule while for the Al atom the HAlSH molecule forms first, which then further isomerizes to H2AlS upon >500 nm irradiation. The reaction of the Ga atom with H2S only takes place upon photolysis to produce HGaSH in the matrix. The assignments of the major modes for these products were confirmed by appropriate 10B, 11B, D2S, and H234S isotopic shifts and theoretical frequency calculations. Topological analysis of the electron density suggests that both HBSH and H2BS molecules possess covalent B-S bond with significant double bond character, while the M-S bond in the heavier group 13 homologues (Al, Ga) was characterized as a polar covalent strong interaction.
Collapse
Affiliation(s)
- Jie Zhao
- School of Chemistry and Chemical Engineering , Guizhou University , Guiyang 550025 , Guizhou China.,School of Chemical Science and Engineering , Tongji University , Shanghai 200092 , China
| | - Qiang Wang
- State Key Laboratory of Coal Conversion, Institute of Coal Chemistry , Chinese Academy of Sciences , Taiyuan 030001 , Shanxi China
| | - Wenjie Yu
- School of Chemical Science and Engineering , Tongji University , Shanghai 200092 , China
| | - Tengfei Huang
- School of Chemical Science and Engineering , Tongji University , Shanghai 200092 , China
| | - Xuefeng Wang
- School of Chemical Science and Engineering , Tongji University , Shanghai 200092 , China
| |
Collapse
|
79
|
Ziegler B, Rauhut G. Rigorous use of symmetry within the construction of multidimensional potential energy surfaces. J Chem Phys 2018; 149:164110. [DOI: 10.1063/1.5047912] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023] Open
Affiliation(s)
- Benjamin Ziegler
- 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
|
80
|
Roy TK, Kopysov V, Pereverzev A, Šebek J, Gerber RB, Boyarkin OV. Intrinsic structure of pentapeptide Leu-enkephalin: geometry optimization and validation by comparison of VSCF-PT2 calculations with cold ion spectroscopy. Phys Chem Chem Phys 2018; 20:24894-24901. [PMID: 30234204 DOI: 10.1039/c8cp03989e] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
The intrinsic structure of an opioid peptide [Ala2, Leu5]-leucine enkephalin (ALE) has been investigated using first-principles based vibrational self-consistent field (VSCF) theory and cold ion spectroscopy. IR-UV double resonance spectroscopy revealed the presence of only one highly abundant conformer of the singly protonated ALE, isolated and cryogenically cooled in the gas phase. High-level quantum mechanical calculations of electronic structures in conjunction with a systematic conformational search allowed for finding a few low-energy candidate structures. In order to identify the observed structure, we computed vibrational spectra of the candidate structures and employed the theory at the semi-empirically scaled harmonic level and at the first-principles based anharmonic VSCF levels. The best match between the calculated "anharmonic" and the measured spectra appeared, indeed, for the most stable candidate. An average of two spectra calculated with different quantum mechanical potentials is proposed for the best match with experiment. The match thus validates the calculated intrinsic structure of ALE and demonstrates the predictive power of first-principles theory for solving structures of such large molecules.
Collapse
Affiliation(s)
- Tapta Kanchan Roy
- Department of Chemistry and Chemical Sciences, Central University of Jammu, Rahya-Suchani (Bagla), Jammu 181143, India
| | | | | | | | | | | |
Collapse
|
81
|
Kuenzer U, Klotz M, Hofer TS. Probing vibrational coupling via a grid-based quantum approach-an efficient strategy for accurate calculations of localized normal modes in solid-state systems. J Comput Chem 2018; 39:2196-2209. [PMID: 30341952 PMCID: PMC6767160 DOI: 10.1002/jcc.25533] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/16/2018] [Revised: 06/08/2018] [Accepted: 06/29/2018] [Indexed: 11/19/2022]
Abstract
In this work an approach to investigate the properties of strongly localized vibrational modes of functional groups in bulk material and on solid-state surfaces is presented. The associated normal mode vectors are approximated solely on the basis of structural information and obtained via diagonalization of a reduced Hessian. The grid-based Numerov procedure in one and two dimensions is then applied to an adequate scan of the respective potential surface yielding the associated vibrational wave functions and energy eigenvalues. This not only provides a detailed description of anharmonic effects but also an accurate inclusion of the coupling between the investigated vibrational states on a quantum mechanical level. All results obtained for the constructed normal modes are benchmarked against their analytical counterparts obtained from the diagonalization of the total Hessian of the entire system. Three increasingly complex systems treated at quantum chemical level of theory have been considered, namely the symmetric and asymmetric stretch vibrations of an isolated water molecule, hydroxyl groups bound to the surface of GeO2 (001), α-quartz(001) and Rutil (001) as well as crystalline Li2 NH serving as an example for a bulk material. While the data obtained for the individual systems verify the applicability of the proposed methodology, comparison to experimental data demonstrates the accuracy of this methodology despite the restriction to limit this methodology to a few selected vibrational modes. The possibility to investigate vibrational phenomena of localized normal modes without the requirement of executing costly harmonic frequency calculations of the entire system enables the application of this method to cases in which the determination of normal modes is prohibitively expensive or not available for a particular level of theory. © 2018 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Ulrich Kuenzer
- University of Innsbruck, Institute of General, Inorganic and Theoretical Chemistry, Innrain 80‐826020 InnsbruckAustria
| | - Martin Klotz
- University of Innsbruck, Institute of General, Inorganic and Theoretical Chemistry, Innrain 80‐826020 InnsbruckAustria
| | - Thomas S. Hofer
- University of Innsbruck, Institute of General, Inorganic and Theoretical Chemistry, Innrain 80‐826020 InnsbruckAustria
| |
Collapse
|
82
|
Roy S, Maiti KS. Structural sensitivity of CH vibrational band in methyl benzoate. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2018; 196:289-294. [PMID: 29459159 DOI: 10.1016/j.saa.2018.02.031] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/05/2017] [Accepted: 02/08/2018] [Indexed: 06/08/2023]
Abstract
The CH vibrational bands of methyl benzoate are studied to understand its coupling pattern with other vibrational bands of the biological molecule. This will facilitate to understand the biological structure and dynamics in spectroscopic as well as in microscopic study. Due to the congested spectroscopic pattern, near degeneracy, and strong anharmonicity of the CH stretch vibrations, assignment of the CH vibrational frequencies are often misleading. Anharmonic vibrational frequency calculation with multidimensional potential energy surface interprets the CH vibrational spectra more accurately. In this article we have presented the importance of multidimensional potential energy surface in anharmonic vibrational frequency calculation and discuss the unexpected red shift of asymmetric CH stretch vibration of methyl group. The CD stretch vibrational band which is splitted to double peaks due to the Fermi resonance is also discussed here.
Collapse
Affiliation(s)
- Susmita Roy
- Die Klinik und Poliklinik für Orthopädie und Sportorthopädie, Technische Universität München, Ismaninger Straße 22, München D81675, Germany
| | - Kiran Sankar Maiti
- Lehrstuhl für Experimental Physik, Ludwig-Maximilians-Universität München, Am Couombwall 1, Garching 85748, Germany; Max-Planck-Institut für Quantenoptik, Hans-Kopfermann-Straße 1, Garching 85748, Germany.
| |
Collapse
|
83
|
Abstract
Transition metal complexes bearing terminal oxido ligands are quite common, yet group 11 terminal oxo complexes remain elusive. Here we show that excited coinage metal atoms M (M = Au, Ag, Cu) react with OF2 to form hypofluorites FOMF and group 11 oxygen metal fluorides OMF2, OAuF and OAgF. These compounds have been characterized by IR matrix-isolation spectroscopy in conjunction with state-of-the-art quantum-chemical calculations. The oxygen fluorides are formed by photolysis of the initially prepared hypofluorites. The linear molecules OAgF and OAuF have a 3Σ − ground state with a biradical character. Two unpaired electrons are located mainly at the oxygen ligand in antibonding O−M π* orbitals. For the 2B2 ground state of the OMIIIF2 compounds only an O−M single bond arises and a significant spin-density contribution was found at the oxygen atom as well. While transition metal complexes bearing terminal oxido ligands are common, those of group 11 elements have yet to be experimentally observed. Here, Riedel and colleagues synthesise molecular oxygen fluorides of copper, silver and gold, and show that the oxo ligands possess radical character.
Collapse
|
84
|
Morgan WJ, Matthews DA, Ringholm M, Agarwal J, Gong JZ, Ruud K, Allen WD, Stanton JF, Schaefer HF. Geometric Energy Derivatives at the Complete Basis Set Limit: Application to the Equilibrium Structure and Molecular Force Field of Formaldehyde. J Chem Theory Comput 2018; 14:1333-1350. [DOI: 10.1021/acs.jctc.7b01138] [Citation(s) in RCA: 27] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- W. James Morgan
- Center for Computational Quantum Chemistry (CCQC), University of Georgia, Athens, Georgia 30602, United States
| | - Devin A. Matthews
- Institute for Computational Engineering and Sciences (ICES), University of Texas at Austin, Austin, Texas 78712, United States
| | - Magnus Ringholm
- Hylleraas Centre for Quantum Molecular Science, Department of Chemistry, University of Tromsø − The Arctic University of Norway, N-9037 Tromsø, Norway
| | - Jay Agarwal
- Center for Computational Quantum Chemistry (CCQC), University of Georgia, Athens, Georgia 30602, United States
| | - Justin Z. Gong
- Department of Chemistry, University of Texas at Austin, Austin, Texas 78712, United States
| | - Kenneth Ruud
- Hylleraas Centre for Quantum Molecular Science, Department of Chemistry, University of Tromsø − The Arctic University of Norway, N-9037 Tromsø, Norway
| | - Wesley D. Allen
- Center for Computational Quantum Chemistry (CCQC), University of Georgia, Athens, Georgia 30602, United States
| | - John F. Stanton
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611, United States
| | - Henry F. Schaefer
- Center for Computational Quantum Chemistry (CCQC), University of Georgia, Athens, Georgia 30602, United States
| |
Collapse
|
85
|
Klinting EL, Thomsen B, Godtliebsen IH, Christiansen O. Employing general fit-bases for construction of potential energy surfaces with an adaptive density-guided approach. J Chem Phys 2018; 148:064113. [DOI: 10.1063/1.5016259] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
| | - Bo Thomsen
- Department of Chemistry, Aarhus University, DK-8000 Aarhus C, Denmark
| | | | - Ove Christiansen
- Department of Chemistry, Aarhus University, DK-8000 Aarhus C, Denmark
| |
Collapse
|
86
|
Ravichandran L, Banik S. Anomalous description of the anharmonicity of bending motions of carbon–carbon double bonded molecules with the MP2 method: ethylene as a case study. Phys Chem Chem Phys 2018; 20:27329-27341. [DOI: 10.1039/c8cp04682d] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
In this work, we report anomalous descriptions of bending modes of the smallest carbon–carbon double bonded molecule ethylene with the MP2 method.
Collapse
Affiliation(s)
| | - Subrata Banik
- Department of Chemistry and Chemical Sciences
- Central University of Jammu
- Samba-181143
- India
| |
Collapse
|
87
|
Madsen D, Christiansen O, König C. Anharmonic vibrational spectra from double incremental potential energy and dipole surfaces. Phys Chem Chem Phys 2018; 20:3445-3456. [DOI: 10.1039/c7cp07190f] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Using incremental approaches, size limitations for property surface generations are pushed significantly, enabling accurate large molecule anharmonic vibrational spectra calculations.
Collapse
Affiliation(s)
- Diana Madsen
- Department of Chemistry
- Aarhus University
- DK-8000 Aarhus C
- Denmark
| | | | - Carolin König
- Division of Theoretical Chemistry & Biology
- Royal Institute of Technology
- SE-106 91 Stockholm
- Sweden
| |
Collapse
|
88
|
Petrenko T, Rauhut G. A General Approach for Calculating Strongly Anharmonic Vibronic Spectra with a High Density of States: The X̃2B1 ← X̃1A1 Photoelectron Spectrum of Difluoromethane. J Chem Theory Comput 2017; 13:5515-5527. [DOI: 10.1021/acs.jctc.7b00468] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Taras Petrenko
- Institute for Theoretical Chemistry, Pfaffenwaldring
55, 70569 Stuttgart, Germany
| | - Guntram Rauhut
- Institute for Theoretical Chemistry, Pfaffenwaldring
55, 70569 Stuttgart, Germany
| |
Collapse
|
89
|
Belinassi AR, Ornellas FR. A theoretical exploration of new species on the 1[H, Se, I] potential energy surface: Energetics, structures, IR spectra, and heats of formation. COMPUT THEOR CHEM 2017. [DOI: 10.1016/j.comptc.2017.06.008] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
90
|
Gronowski M, Eluszkiewicz P, Custer T. Structure and Spectroscopy of C2HNO Isomers. J Phys Chem A 2017; 121:3263-3273. [DOI: 10.1021/acs.jpca.6b12609] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Marcin Gronowski
- Institute of Physical Chemistry, Polish Academy of Sciences, Kasprzaka 44/52, 01-224 Warsaw, Poland
| | - Piotr Eluszkiewicz
- Institute of Physical Chemistry, Polish Academy of Sciences, Kasprzaka 44/52, 01-224 Warsaw, Poland
| | - Thomas Custer
- Institute of Physical Chemistry, Polish Academy of Sciences, Kasprzaka 44/52, 01-224 Warsaw, Poland
| |
Collapse
|
91
|
Petrenko T, Rauhut G. A new efficient method for the calculation of interior eigenpairs and its application to vibrational structure problems. J Chem Phys 2017; 146:124101. [DOI: 10.1063/1.4978581] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.4] [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
|
92
|
Carrington T. Perspective: Computing (ro-)vibrational spectra of molecules with more than four atoms. J Chem Phys 2017; 146:120902. [DOI: 10.1063/1.4979117] [Citation(s) in RCA: 67] [Impact Index Per Article: 9.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Tucker Carrington
- Chemistry Department, Queen’s University, Kingston, Ontario
K7L 3N6, Canada
| |
Collapse
|
93
|
Panek PT, Jacob CR. On the benefits of localized modes in anharmonic vibrational calculations for small molecules. J Chem Phys 2017; 144:164111. [PMID: 27131535 DOI: 10.1063/1.4947213] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Anharmonic vibrational calculations can already be computationally demanding for relatively small molecules. The main bottlenecks lie in the construction of the potential energy surface and in the size of the excitation space in the vibrational configuration interaction (VCI) calculations. To address these challenges, we use localized-mode coordinates to construct potential energy surfaces and perform vibrational self-consistent field and L-VCI calculations [P. T. Panek and C. R. Jacob, ChemPhysChem 15, 3365 (2014)] for all vibrational modes of two prototypical test cases, the ethene and furan molecules. We find that the mutual coupling between modes is reduced when switching from normal-mode coordinates to localized-mode coordinates. When using such localized-mode coordinates, we observe a faster convergence of the n-mode expansion of the potential energy surface. This makes it possible to neglect higher-order contributions in the n-mode expansion of the potential energy surface or to approximate higher-order contributions in hybrid potential energy surfaces, which reduced the computational effort for the construction of the anharmonic potential energy surface significantly. Moreover, we find that when using localized-mode coordinates, the convergence with respect to the VCI excitation space proceeds more smoothly and that the error at low orders is reduced significantly. This makes it possible to devise low-cost models for obtaining a first approximation of anharmonic corrections. This demonstrates that the use of localized-mode coordinates can be beneficial already in anharmonic vibrational calculations of small molecules and provides a possible avenue for enabling such accurate calculations also for larger molecules.
Collapse
Affiliation(s)
- Paweł T Panek
- TU Braunschweig, Institute of Physical and Theoretical Chemistry, Hans-Sommer-Str. 10, 38106 Braunschweig, Germany
| | - Christoph R Jacob
- TU Braunschweig, Institute of Physical and Theoretical Chemistry, Hans-Sommer-Str. 10, 38106 Braunschweig, Germany
| |
Collapse
|
94
|
Chakraborty S, Banik S, Das PK. Anharmonicity in the Vibrational Spectra of Naphthalene and Naphthalene-d8: Experiment and Theory. J Phys Chem A 2016; 120:9707-9718. [DOI: 10.1021/acs.jpca.6b09034] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Shubhadip Chakraborty
- Department
of Inorganic and Physical Chemistry, Indian Institute of Science, Bangalore, India
| | - Subrata Banik
- Advanced
Centre for Research in High Energy Materials and School of Chemistry, University of Hyderabad, Hyderabad, India
| | - Puspendu K. Das
- Department
of Inorganic and Physical Chemistry, Indian Institute of Science, Bangalore, India
| |
Collapse
|
95
|
Partanen L, Hänninen V, Halonen L. Effects of Global and Local Anharmonicities on the Thermodynamic Properties of Sulfuric Acid Monohydrate. J Chem Theory Comput 2016; 12:5511-5524. [PMID: 27662456 DOI: 10.1021/acs.jctc.6b00683] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
We use state-of-the-art electronic structure calculation methods and large basis sets to obtain reliable values for the thermodynamic properties of sulfuric acid monohydrate and study the effects of vibrational anharmonicity on these properties. We distinguish between two forms of vibrational anharmonicity: local anharmonicity, which refers to the anharmonicity of the vibrational modes of a given cluster conformer, and global anharmonicity, which originates from accounting for the presence of different conformers in the first place. In our most accurate approach, we solve the nuclear Schrödinger equation variationally for the intermolecular large-amplitude motions, thus quantum-mechanically accounting for the presence of higher-energy conformers for both reactants and products, while using the standard vibrational perturbational approach for the other vibrational modes. This results in a value of -11.0 kJ/mol for the reaction Gibbs free energy at 298.15 K. When standard vibrational perturbational approaches are employed, the effects of local anharmonicity depend heavily on the choice of the electronic structure calculation basis set. In fact, better results can often be achieved by combining a simple harmonic treatment for the vibrational partition function with a statistical mechanical accounting of global anharmonicity. Thus, we recommend that future studies that intend to include anharmonicity start by accounting for the presence of higher-energy conformers and only then consider whether local anharmonicity calculations are feasible and necessary.
Collapse
Affiliation(s)
- Lauri Partanen
- Laboratory of Physical Chemistry, Department of Chemistry, University of Helsinki , P.O. Box 55, A.I. Virtasen aukio 1, FIN-00014 Helsinki, Finland
| | - Vesa Hänninen
- Laboratory of Physical Chemistry, Department of Chemistry, University of Helsinki , P.O. Box 55, A.I. Virtasen aukio 1, FIN-00014 Helsinki, Finland
| | - Lauri Halonen
- Laboratory of Physical Chemistry, Department of Chemistry, University of Helsinki , P.O. Box 55, A.I. Virtasen aukio 1, FIN-00014 Helsinki, Finland
| |
Collapse
|
96
|
Fingerhut BP, Costard R, Elsaesser T. Predominance of short range Coulomb forces in phosphate-water interactions—a theoretical analysis. J Chem Phys 2016. [DOI: 10.1063/1.4962755] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023] Open
Affiliation(s)
- Benjamin P. Fingerhut
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie, D-12489 Berlin, Germany
| | - Rene Costard
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie, D-12489 Berlin, Germany
| | - Thomas Elsaesser
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie, D-12489 Berlin, Germany
| |
Collapse
|
97
|
Ostrowski L, Ziegler B, Rauhut G. Tensor decomposition in potential energy surface representations. J Chem Phys 2016; 145:104103. [DOI: 10.1063/1.4962368] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022] Open
Affiliation(s)
- Lukas Ostrowski
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Benjamin Ziegler
- 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
|
98
|
Kowalewski M, Larsson E, Heryudono A. An adaptive interpolation scheme for molecular potential energy surfaces. J Chem Phys 2016; 145:084104. [DOI: 10.1063/1.4961148] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Markus Kowalewski
- Department of Information Technology, Uppsala University, Box 337, SE-751 05 Uppsala, Sweden
| | - Elisabeth Larsson
- Department of Information Technology, Uppsala University, Box 337, SE-751 05 Uppsala, Sweden
| | - Alfa Heryudono
- Department of Mathematics, University of Massachusetts Dartmouth, Dartmouth, Massachusetts 02747, USA
| |
Collapse
|
99
|
König C, Christiansen O. Linear-scaling generation of potential energy surfaces using a double incremental expansion. J Chem Phys 2016. [DOI: 10.1063/1.4960189] [Citation(s) in RCA: 35] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Carolin König
- Department of Chemistry, Aarhus University, DK-8000 Aarhus C, Denmark
| | - Ove Christiansen
- Department of Chemistry, Aarhus University, DK-8000 Aarhus C, Denmark
| |
Collapse
|
100
|
On the choice electronic structure method to calculate the quartic potential energy surface for the vibrational calculation of polyatomic molecules. Theor Chem Acc 2016. [DOI: 10.1007/s00214-016-1962-1] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|