1
|
Mizuide A, Fujii A. Hydrogen bond network structures of protonated dimethylamine clusters H +(DMA) n ( n = 3-7). Phys Chem Chem Phys 2024; 26:19418-19432. [PMID: 38973623 DOI: 10.1039/d4cp01931h] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/09/2024]
Abstract
Infrared spectroscopy of protonated dimethylamine clusters, H+(DMA)n, (n = 3-7), and their Ar-tagged clusters was performed in the NH and CH stretching vibrational region to explore their hydrogen bond network structures. A stable isomer search and vibrational spectral simulations of the clusters were also carried out to support the interpretations of the observed spectra. Weakly hydrogen-bonded NH stretching vibrational bands, which are characteristic of cyclic structures of small-sized protonated clusters, are observed in the spectra of the Ar-tagged clusters of n ≥ 5, while only linear chain type structures are suggested for the Ar-tagged clusters of n = 3-4 and the bare clusters of all the sizes. These results demonstrate that the size and temperature dependence of the hydrogen bond network structures of the protonated dimethylamine clusters is analogous to that of protonated monohydric alcohol clusters.
Collapse
Affiliation(s)
- Atsuya Mizuide
- Department of Chemistry, Graduate School of Science, Tohoku University, Sendai 980-8578, Japan.
| | - Asuka Fujii
- Department of Chemistry, Graduate School of Science, Tohoku University, Sendai 980-8578, Japan.
| |
Collapse
|
2
|
Hino K, Kurashige Y. Encoding a Many-Body Potential Energy Surface into a Grid-Based Matrix Product Operator. J Chem Theory Comput 2024; 20:3839-3849. [PMID: 38647101 DOI: 10.1021/acs.jctc.4c00046] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/25/2024]
Abstract
An efficient algorithm for compressing a given many-body potential energy surface (PES) of molecular systems into a grid-based matrix product operator (MPO) is proposed. The PES is once represented by a full-dimensional or truncated many-body expansion form, which is obtained by ab initio calculations at each grid mesh point, and then all terms in the expansion are compressed and merged into a single MPO while maintaining the bond dimension of the MPO as small as possible. It was shown that the ab initio PES of the H2CO was compressed by more than 2 orders of magnitude in the size of the site operators without loss of accuracy. By the use of grid basis, the tensor rank of the site operators of the MPO is reduced from four to three due to the diagonal nature of the position-dependent operators on grid basis, which significantly reduces the computational cost of the tensor contractions required in the real and imaginary time evolution of the matrix product state (MPS) wave functions with the grid-based MPO (Grid-MPO) Hamiltonian. Similar to other grid-based methods, Grid-MPO is easily applicable to any kinds of potentials of molecular systems, such as analytical empirical model potentials expressed by position operators and ab initio potentials, if the values at the grid points are available. Using the Grid-MPO combined with the MPS, we calculated the time correlation function of the Eigen cation H 3 O + ( H 2 O ) 3 to predict the infrared spectrum and compared with the experimental and the previous theoretical studies. The actual scaling with the size of systems was examined for the multidimensional Henon-Heiles Hamiltonian. It was shown that the method is considerably accelerated by the graphic processing unit (GPU) because the sizes of site operators were kept small and all tensors were able to be stored on the VRAM of a GPU.
Collapse
Affiliation(s)
- Kentaro Hino
- Department of Chemistry, Graduate School of Science, Kyoto University, Kitashirakawa Oiwake-cho, Sakyo-ku, Kyoto 606-8502, Japan
| | - Yuki Kurashige
- Department of Chemistry, Graduate School of Science, Kyoto University, Kitashirakawa Oiwake-cho, Sakyo-ku, Kyoto 606-8502, Japan
- FOREST, JST, Honcho 4-1-8, Kawaguchi, Saitama 332-0012, Japan
- CREST, JST, Honcho 4-1-8, Kawaguchi, Saitama 332-0012, Japan
| |
Collapse
|
3
|
Huang QR, Yano K, Yang Y, Fujii A, Kuo JL. Near-infrared spectroscopy of H 3O +⋯X n (X = Ar, N 2, and CO, n = 1-3). Phys Chem Chem Phys 2024; 26:10757-10768. [PMID: 38516880 DOI: 10.1039/d4cp00458b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/23/2024]
Abstract
Near-infrared (NIR) spectra of H3O+⋯Xn (X = Ar, N2, and CO, n = 1-3) in the first overtone region of OH-stretching vibrations (4800-7000 cm-1) were measured. Not only OH-stretching overtones but also several combination bands are major features in this region, and assignments of these observed bands are not obvious at a glance. High-precision anharmonic vibrational simulations based on the discrete variable representation approach were performed. The simulated spectra show good agreement with the observed ones and provide firm assignments of the observed bands, except in the case of X = CO, in which higher order vibrational mode couplings seem significant. This agreement demonstrates that the present system can be a benchmark for high precision anharmonic vibrational computations of NIR spectra. Band broadening in the observed spectra becomes remarkable with an increase of the interaction with the solvent molecule (X). The origin of the band broadening is explored by rare gas tagging experiments and anharmonic vibrational simulations of hot bands.
Collapse
Affiliation(s)
- Qian-Rui Huang
- Institute of Atomic and Molecular Sciences, Academia Sinica, Taipei, Taiwan.
| | - Kazuyoshi Yano
- Department of Chemistry, Graduate School of Science, Tohoku University, Sendai, Japan.
| | - Yaodi Yang
- Department of Chemistry, Graduate School of Science, Tohoku University, Sendai, Japan.
| | - Asuka Fujii
- Department of Chemistry, Graduate School of Science, Tohoku University, Sendai, Japan.
| | - Jer-Lai Kuo
- Institute of Atomic and Molecular Sciences, Academia Sinica, Taipei, Taiwan.
| |
Collapse
|
4
|
Moscato D, Mandelli G, Bondanza M, Lipparini F, Conte R, Mennucci B, Ceotto M. Unraveling Water Solvation Effects with Quantum Mechanics/Molecular Mechanics Semiclassical Vibrational Spectroscopy: The Case of Thymidine. J Am Chem Soc 2024; 146:8179-8188. [PMID: 38470354 DOI: 10.1021/jacs.3c12700] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/13/2024]
Abstract
We introduce a quantum mechanics/molecular mechanics semiclassical method for studying the solvation process of molecules in water at the nuclear quantum mechanical level with atomistic detail. We employ it in vibrational spectroscopy calculations because this is a tool that is very sensitive to the molecular environment. Specifically, we look at the vibrational spectroscopy of thymidine in liquid water. We find that the C═O frequency red shift and the C═C frequency blue shift, experienced by thymidyne upon solvation, are mainly due to reciprocal polarization effects, that the molecule and the water solvent exert on each other, and nuclear zero-point energy effects. In general, this work provides an accurate and practical tool to study quantum vibrational spectroscopy in solution and condensed phase, incorporating high-level and computationally affordable descriptions of both electronic and nuclear problems.
Collapse
Affiliation(s)
- Davide Moscato
- Dipartimento di Chimica, Università degli Studi di Milano, Via Golgi, 19, 20133 Milano, Italy
| | - Giacomo Mandelli
- Dipartimento di Chimica, Università degli Studi di Milano, Via Golgi, 19, 20133 Milano, Italy
| | - Mattia Bondanza
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via Giuseppe Moruzzi, 13, 56124 Pisa, Italy
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via Giuseppe Moruzzi, 13, 56124 Pisa, Italy
| | - Riccardo Conte
- Dipartimento di Chimica, Università degli Studi di Milano, Via Golgi, 19, 20133 Milano, Italy
| | - Benedetta Mennucci
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via Giuseppe Moruzzi, 13, 56124 Pisa, Italy
| | - Michele Ceotto
- Dipartimento di Chimica, Università degli Studi di Milano, Via Golgi, 19, 20133 Milano, Italy
| |
Collapse
|
5
|
Zhang Y, Wang Y, Xu X, Chen Z, Yang Y. Vibrational Spectra of Highly Anharmonic Water Clusters: Molecular Dynamics and Harmonic Analysis Revisited with Constrained Nuclear-Electronic Orbital Methods. J Chem Theory Comput 2023; 19:9358-9368. [PMID: 38096546 DOI: 10.1021/acs.jctc.3c01037] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2023]
Abstract
Vibrational spectroscopy is widely used to gain insights into structural and dynamic properties of chemical, biological, and materials systems. Thus, an efficient and accurate method to simulate vibrational spectra is desired. In this paper, we justify and employ a microcanonical molecular simulation scheme to calculate the vibrational spectra of three challenging water clusters: the neutral water dimer (H4O2), the protonated water trimer (H7O3+), and the protonated water tetramer (H9O4+). We find that with the accurate description of quantum nuclear delocalization effects through the constrained nuclear-electronic orbital framework, including vibrational mode coupling effects through molecular dynamics simulations can additionally improve the vibrational spectrum calculations. In contrast, without the quantum nuclear delocalization picture, conventional ab initio molecular dynamics may even lead to less accurate results than harmonic analysis.
Collapse
Affiliation(s)
- Yuzhe Zhang
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin-Madison, 1101 University Avenue, Madison, Wisconsin 53706, United States
| | - Yiwen Wang
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin-Madison, 1101 University Avenue, Madison, Wisconsin 53706, United States
| | - Xi Xu
- Center for Advanced Materials Research, Beijing Normal University, Zhuhai 519087, China
| | - Zehua Chen
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin-Madison, 1101 University Avenue, Madison, Wisconsin 53706, United States
| | - Yang Yang
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin-Madison, 1101 University Avenue, Madison, Wisconsin 53706, United States
| |
Collapse
|
6
|
Brünig FN, Hillmann P, Kim WK, Daldrop JO, Netz RR. Proton-transfer spectroscopy beyond the normal-mode scenario. J Chem Phys 2022; 157:174116. [DOI: 10.1063/5.0116686] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022] Open
Abstract
A stochastic theory is developed to predict the spectral signature of proton-transfer processes and is applied to infrared spectra computed from ab initio molecular-dynamics simulations of a single [Formula: see text] cation. By constraining the oxygen atoms to a fixed distance, this system serves as a tunable model for general proton-transfer processes with variable barrier height. Three spectral contributions at distinct frequencies are identified and analytically predicted: the quasi-harmonic motion around the most probable configuration, amenable to normal-mode analysis, the contribution due to transfer paths when the proton moves over the barrier, and a shoulder for low frequencies stemming from the stochastic transfer-waiting-time distribution; the latter two contributions are not captured by normal-mode analysis but exclusively reported on the proton-transfer kinetics. In accordance with reaction rate theory, the transfer-waiting-contribution frequency depends inversely exponentially on the barrier height, whereas the transfer-path-contribution frequency is rather insensitive to the barrier height.
Collapse
Affiliation(s)
- Florian N. Brünig
- Department of Physics, Freie Universität Berlin, 14195 Berlin, Germany
| | - Paul Hillmann
- Department of Physics, Freie Universität Berlin, 14195 Berlin, Germany
| | - Won Kyu Kim
- School of Computational Sciences, Korea Institute for Advanced Study, Seoul 02455, Republic of Korea
| | - Jan O. Daldrop
- Department of Physics, Freie Universität Berlin, 14195 Berlin, Germany
| | - Roland R. Netz
- Department of Physics, Freie Universität Berlin, 14195 Berlin, Germany
| |
Collapse
|
7
|
The coupling of the hydrated proton to its first solvation shell. Nat Commun 2022; 13:6170. [PMID: 36257946 PMCID: PMC9579203 DOI: 10.1038/s41467-022-33650-w] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/22/2022] [Accepted: 09/28/2022] [Indexed: 11/30/2022] Open
Abstract
The Zundel (\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${H}_{5}{O}_{2}^{+}$$\end{document}H5O2+) and Eigen (\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${H}_{9}{O}_{4}^{+}$$\end{document}H9O4+) cations play an important role as intermediate structures for proton transfer processes in liquid water. In the gas phase they exhibit radically different infrared (IR) spectra. The question arises: is there a least common denominator structure that explains the IR spectra of both, the Zundel and Eigen cations, and hence of the solvated proton? Full dimensional quantum simulations of these protonated cations demonstrate that two dynamical water molecules and an excess proton constitute this fundamental subunit. Embedded in the static environment of the parent Eigen cation, this subunit reproduces the positions and broadenings of its main excess-proton bands. In isolation, its spectrum reverts to the well-known Zundel ion. Hence, the dynamics of this subunit polarized by an environment suffice to explain the spectral signatures and anharmonic couplings of the solvated proton in its first solvation shell. The Zundel [H(H2O)2]+ and Eigen [H(H2O)4]+ cations exhibit radicallly different infrared spectra and are the limiting dynamical structures involved in proton mobility in liquid water. Here, the authors find through quantum dynamics simulations that two polarized water molecules and a proton suffice to explain the key spectroscopic features connected to proton mobility for both species.
Collapse
|
8
|
Brünig FN, Rammler M, Adams EM, Havenith M, Netz RR. Spectral signatures of excess-proton waiting and transfer-path dynamics in aqueous hydrochloric acid solutions. Nat Commun 2022; 13:4210. [PMID: 35864099 PMCID: PMC9304333 DOI: 10.1038/s41467-022-31700-x] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/05/2021] [Accepted: 06/22/2022] [Indexed: 11/23/2022] Open
Abstract
The theoretical basis for linking spectral signatures of hydrated excess protons with microscopic proton-transfer mechanisms has so far relied on normal-mode analysis. We introduce trajectory-decomposition techniques to analyze the excess-proton dynamics in ab initio molecular-dynamics simulations of aqueous hydrochloric-acid solutions beyond the normal-mode scenario. We show that the actual proton transfer between two water molecules involves for relatively large water-water separations crossing of a free-energy barrier and thus is not a normal mode, rather it is characterized by two non-vibrational time scales: Firstly, the broadly distributed waiting time for transfer to occur with a mean value of 200–300 fs, which leads to a broad and weak shoulder in the absorption spectrum around 100 cm−1, consistent with our experimental THz spectra. Secondly, the mean duration of a transfer event of about 14 fs, which produces a rather well-defined spectral contribution around 1200 cm−1 and agrees in location and width with previous experimental mid-infrared spectra. The spectroscopic signatures of excess protons in HCl solutions are studied by ab initio simulations and THz experiments. Two contributions beyond the normal-mode scenario are identified that reflect proton-waiting and proton-transfer processes.
Collapse
Affiliation(s)
- Florian N Brünig
- Freie Universität Berlin, Department of Physics, 14195, Berlin, Germany
| | - Manuel Rammler
- Freie Universität Berlin, Department of Physics, 14195, Berlin, Germany
| | - Ellen M Adams
- Ruhr-Universität Bochum, Department of Physical Chemistry II, 44780, Bochum, Germany
| | - Martina Havenith
- Ruhr-Universität Bochum, Department of Physical Chemistry II, 44780, Bochum, Germany
| | - Roland R Netz
- Freie Universität Berlin, Department of Physics, 14195, Berlin, Germany.
| |
Collapse
|
9
|
DiRisio RJ, Finney JM, McCoy AB. Diffusion Monte Carlo approaches for studying nuclear quantum effects in fluxional molecules. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2022. [DOI: 10.1002/wcms.1615] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Affiliation(s)
- Ryan J. DiRisio
- Department of Chemistry University of Washington Seattle Washington USA
| | - Jacob M. Finney
- Department of Chemistry University of Washington Seattle Washington USA
| | - Anne B. McCoy
- Department of Chemistry University of Washington Seattle Washington USA
| |
Collapse
|
10
|
Huchmala RM, McCoy AB. Exploring the Origins of Spectral Signatures of Strong Hydrogen Bonding in Protonated Water Clusters. J Phys Chem A 2022; 126:1360-1368. [PMID: 35171593 DOI: 10.1021/acs.jpca.1c10036] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The effects of anharmonicity on the spectral features of strong ionic hydrogen bonds are explored through reduced dimensional studies of the couplings between the hydrogen bonding OH and the donor-acceptor OO stretching vibrations in protonated water clusters with 2-4 water molecules. Specifically, this study focuses on how the anharmonicities and couplings in these ions are reflected in the vibrational spectra by exploring the intensities of the transitions to states with excitation in both the OH and the OO stretching vibrations and changes in the frequency of the OO stretching vibration when the OH stretching vibration is excited. These questions are addressed through the application of several approximate treatments that are based on an adiabatic separation of the high-frequency OH and low-frequency OO stretching vibrations as well as low-order expansions of the potential and dipole surfaces. While an adiabatic approximation captures most of the trends found in the spectra and from an analysis of the two-dimensional model, a vibrational Franck-Condon approach fails to capture the intensities of these transitions. Of the terms in the expansion of the dipole moment function, those that are proportional to ΔrOH and ΔrOH2 are found to provide the largest contributions to the calculated intensities of the transitions involving excitation of both the OH and the OO stretches. This leads to the conclusion that the intensities of these transitions encode information about the frequency and anharmonicity of the OH stretching vibration and how they are affected by changes in the OO distance. The anharmonicity of the potential also leads to changes in the OO stretching frequency with excitation of the OH stretching vibration. The direction of this change in frequency encodes additional information about the strength of the ionic hydrogen bond.
Collapse
Affiliation(s)
- Rachel M Huchmala
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Anne B McCoy
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| |
Collapse
|
11
|
Boutwell D, Pierre-Jacques D, Cochran O, Dyke J, Salazar D, Tyler C, Kaledin M. Intramolecular Proton Transfer in the Hydrogen Oxalate Anion and the Cooperativity Effects of the Low-Frequency Vibrations: A Driven Molecular Dynamics Study. J Phys Chem A 2022; 126:583-592. [PMID: 35049313 DOI: 10.1021/acs.jpca.1c09686] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We report first-principles molecular dynamics (MD) and dipole-driven molecular dynamics (μ-DMD) simulations of the hydrogen oxalate anion at the MP2/aug-cc-pVDZ level of theory. We examine the role of vibrational coupling between the OH stretching bands, that is, the fundamental and a few combination bands spanning the 2900-3100 cm-1 range, and several of the low-frequency bending and stretching fundamental modes. The low-frequency modes between 300 and 825 cm-1 play a crucial role in the proton-transfer motion. Strong involvement of CO2 and CCO bending and the CC stretching vibrations indicate that these large amplitude motions cause the shortening of the O···O distance and thus promote H+ transfer to the other oxygen by bringing it over the 3.4 kcal/mol barrier. Analysis of resonant μ-DMD trajectories shows that the complex spectral feature near 825 cm-1, closely corresponding to both an overtone of two quanta of 425 cm-1 and a combination band of low-frequency CO2 rocking (300 cm-1) and CCO bending (575 cm-1) modes, is involved in the proton transfer. μ-DMD shows that exciting the system at these mode combinations leads to faster barrier activation than exciting at the OH fundamental mode.
Collapse
Affiliation(s)
- Dalton Boutwell
- Department of Chemistry & Biochemistry, Kennesaw State University, 370 Paulding Ave NW, Box # 1203, Kennesaw, Georgia 30144, United States
| | - Dominick Pierre-Jacques
- Department of Chemistry & Biochemistry, Kennesaw State University, 370 Paulding Ave NW, Box # 1203, Kennesaw, Georgia 30144, United States
| | - Olivia Cochran
- Department of Chemistry & Biochemistry, Kennesaw State University, 370 Paulding Ave NW, Box # 1203, Kennesaw, Georgia 30144, United States
| | - Jason Dyke
- Department of Chemistry & Biochemistry, Kennesaw State University, 370 Paulding Ave NW, Box # 1203, Kennesaw, Georgia 30144, United States
| | - Dayana Salazar
- Department of Chemistry & Biochemistry, Kennesaw State University, 370 Paulding Ave NW, Box # 1203, Kennesaw, Georgia 30144, United States
| | - Ciara Tyler
- Department of Chemistry & Biochemistry, Kennesaw State University, 370 Paulding Ave NW, Box # 1203, Kennesaw, Georgia 30144, United States
| | - Martina Kaledin
- Department of Chemistry & Biochemistry, Kennesaw State University, 370 Paulding Ave NW, Box # 1203, Kennesaw, Georgia 30144, United States
| |
Collapse
|
12
|
Gaigeot MP. Some opinions on MD-based vibrational spectroscopy of gas phase molecules and their assembly: An overview of what has been achieved and where to go. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2021; 260:119864. [PMID: 34052762 DOI: 10.1016/j.saa.2021.119864] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/14/2021] [Revised: 04/13/2021] [Accepted: 04/18/2021] [Indexed: 06/12/2023]
Abstract
We hereby review molecular dynamics simulations for anharmonic gas phase spectroscopy and provide some of our opinions of where the field is heading. With these new directions, the theoretical IR/Raman spectroscopy of large (bio)-molecular systems will be more easily achievable over longer time-scale MD trajectories for an increase in accuracy of the MD-IR and MD-Raman calculated spectra. With the new directions presented here, the high throughput 'decoding' of experimental IR/Raman spectra into 3D-structures should thus be possible, hence advancing e.g. the field of MS-IR for structural characterization by spectroscopy. We also review the assignment of vibrational spectra in terms of anharmonic molecular modes from the MD trajectories, and especially introduce our recent developments based on Graph Theory algorithms. Graph Theory algorithmic is also introduced in this review for the identification of the molecular 3D-structures sampled over MD trajectories.
Collapse
Affiliation(s)
- Marie-Pierre Gaigeot
- Université Paris-Saclay, Univ Evry, CNRS, LAMBE UMR8587, 91025 Evry-Courcouronnes, France.
| |
Collapse
|
13
|
Hanson MD, Readnour JA, Hassanali AA, Corcelli SA. Coupled Local-Mode Approach for the Calculation of Vibrational Spectra: Application to Protonated Water Clusters. J Phys Chem Lett 2021; 12:9226-9232. [PMID: 34529914 DOI: 10.1021/acs.jpclett.1c02254] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
Spectroscopic studies of protonated water clusters (PWCs) have yielded enormous insights into the fundamental nature of the hydrated proton. Here, we introduce a new coupled local-mode (CLM) approach to calculate PWC OH stretch vibrational spectra. The CLM method combines a sampling of representative configurations from density functional theory (DFT)-based ab initio molecular dynamics (AIMD) simulations with DFT calculations of local-mode vibrational frequencies and couplings. Calculations of inhomogeneous OH stretch vibrational spectra for H+(H2O)4 and H+(H2O)21 agree well with experiment and higher-level calculations, and decompositions of the calculated spectra in terms of the coupled modes aids in the interpretation of the spectra. This observation is consistent with the idea that capturing anharmonicity and coupling is as important to accuracy as the underlying level of electronic structure theory. The CLM calculations can easily discern the configuration that dominates the experimental measurement for H+(H2O)5, which can adopt several low-energy conformations.
Collapse
Affiliation(s)
- Matthew D Hanson
- Department of Chemistry and Biochemistry, University of Notre Dame, Notre Dame, Indiana 46556, United States
| | - Janel A Readnour
- Department of Chemistry and Biochemistry, University of Notre Dame, Notre Dame, Indiana 46556, United States
| | - Ali A Hassanali
- Condensed Matter and Statistical Physics, International Centre for Theoretical Physics, Strada Costiera, 11 I - 34151 Trieste, Italy
| | - Steven A Corcelli
- Department of Chemistry and Biochemistry, University of Notre Dame, Notre Dame, Indiana 46556, United States
| |
Collapse
|
14
|
DiRisio RJ, Finney JM, Dzugan LC, Madison LR, McCoy AB. Using Diffusion Monte Carlo Wave Functions to Analyze the Vibrational Spectra of H 7O 3+ and H 9O 4. J Phys Chem A 2021; 125:7185-7197. [PMID: 34433268 DOI: 10.1021/acs.jpca.1c05025] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
An approach for evaluating spectra from ground state probability amplitudes (GSPA) obtained from diffusion Monte Carlo (DMC) simulations is extended to improve the description of excited state energies and allow for coupling among vibrational excited states. This approach is applied to studies of the protonated water trimer and tetramer, and their deuterated analogs. These ions provide models for solvated hydronium, and analysis of these spectra provides insights into spectral signatures of proton transfer in aqueous environments. In this approach, we obtain a separable set of internal coordinates from the DMC ground state probability amplitude. A basis is then developed from products of the DMC ground state wave function and low-order polynomials in these internal coordinates. This approach provides a compact basis in which the Hamiltonian and dipole moment matrix are evaluated and used to obtain the spectrum. The resulting spectra are in good agreement with experiment and in many cases provide comparable agreement to the results obtained using much larger basis sets. In addition, the compact basis allows for interpretation of the spectral features and how they evolve with cluster size and deuteration.
Collapse
Affiliation(s)
- Ryan J DiRisio
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Jacob M Finney
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Laura C Dzugan
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Lindsey R Madison
- Department of Chemistry, Colby College, Waterville, Maine 04901, United States
| | - Anne B McCoy
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| |
Collapse
|
15
|
Opoku E, Pawłowski F, Ortiz JV. Electron binding energies and Dyson orbitals of O nH 2n+1 +,0,- clusters: Double Rydberg anions, Rydberg radicals, and micro-solvated hydronium cations. J Chem Phys 2021; 154:234304. [PMID: 34241254 DOI: 10.1063/5.0053297] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Ab initio electron propagator methods are employed to predict the vertical electron attachment energies (VEAEs) of OH3 +(H2O)n clusters. The VEAEs decrease with increasing n, and the corresponding Dyson orbitals are diffused over exterior, non-hydrogen bonded protons. Clusters formed from OH3 - double Rydberg anions (DRAs) and stabilized by hydrogen bonding or electrostatic interactions between ions and polar molecules are studied through calculations on OH3 -(H2O)n complexes and are compared with more stable H-(H2O)n+1 isomers. Remarkable changes in the geometry of the anionic hydronium-water clusters with respect to their cationic counterparts occur. Rydberg electrons in the uncharged and anionic clusters are held near the exterior protons of the water network. For all values of n, the anion-water complex H-(H2O)n+1 is always the most stable, with large vertical electron detachment energies (VEDEs). OH3 -(H2O)n DRA isomers have well separated VEDEs and may be visible in anion photoelectron spectra. Corresponding Dyson orbitals occupy regions beyond the peripheral O-H bonds and differ significantly from those obtained for the VEAEs of the cations.
Collapse
Affiliation(s)
- Ernest Opoku
- Department of Chemistry and Biochemistry, Auburn University, Auburn, Alabama 36849-5312, USA
| | - Filip Pawłowski
- Department of Chemistry and Biochemistry, Auburn University, Auburn, Alabama 36849-5312, USA
| | - Joseph Vincent Ortiz
- Department of Chemistry and Biochemistry, Auburn University, Auburn, Alabama 36849-5312, USA
| |
Collapse
|
16
|
Pierre-Jacques D, Tyler C, Dyke J, Kaledin AL, Kaledin M. A polarizability driven ab initio molecular dynamics approach to stimulating Raman activity: Application to C20. Mol Phys 2021. [DOI: 10.1080/00268976.2021.1939453] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Affiliation(s)
| | - Ciara Tyler
- Department of Chemistry & Biochemistry, Kennesaw State University, Kennesaw, GA, USA
| | - Jason Dyke
- Department of Chemistry & Biochemistry, Kennesaw State University, Kennesaw, GA, USA
| | - Alexey L. Kaledin
- Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, GA, USA
| | - Martina Kaledin
- Department of Chemistry & Biochemistry, Kennesaw State University, Kennesaw, GA, USA
| |
Collapse
|
17
|
Zeng HJ, Johnson MA. Demystifying the Diffuse Vibrational Spectrum of Aqueous Protons Through Cold Cluster Spectroscopy. Annu Rev Phys Chem 2021; 72:667-691. [PMID: 33646816 DOI: 10.1146/annurev-physchem-061020-053456] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
The ease with which the pH is routinely determined for aqueous solutions masks the fact that the cationic product of Arrhenius acid dissolution, the hydrated proton, or H+(aq), is a remarkably complex species. Here, we review how results obtained over the past 30 years in the study of H+⋅(H2O)n cluster ions isolated in the gas phase shed light on the chemical nature of H+(aq). This effort has also revealed molecular-level aspects of the Grotthuss relay mechanism for positive-charge translocation in water. Recently developed methods involving cryogenic cooling in radiofrequency ion traps and the application of two-color, infrared-infrared (IR-IR) double-resonance spectroscopy have established a clear picture of how local hydrogen-bond topology drives the diverse spectral signatures of the excess proton. This information now enables a new generation of cluster studies designed to unravel the microscopic mechanics underlying the ultrafast relaxation dynamics displayed by H+(aq).
Collapse
Affiliation(s)
- Helen J Zeng
- Sterling Chemistry Laboratory, Yale University, New Haven, Connecticut 06520, USA;
| | - Mark A Johnson
- Sterling Chemistry Laboratory, Yale University, New Haven, Connecticut 06520, USA;
| |
Collapse
|
18
|
Finney JM, DiRisio RJ, McCoy AB. Guided Diffusion Monte Carlo: A Method for Studying Molecules and Ions That Display Large Amplitude Vibrational Motions. J Phys Chem A 2020; 124:9567-9577. [DOI: 10.1021/acs.jpca.0c07181] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Jacob M. Finney
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Ryan J. DiRisio
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Anne B. McCoy
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| |
Collapse
|
19
|
Boutwell D, Okere O, Omodemi O, Toledo A, Barrios A, Olocha M, Kaledin M. Analysis of the Proton Transfer Bands in the Infrared Spectra of Linear N 2H +···OC and N 2D +···OC Complexes Using Electric Field-Driven Classical Trajectories. J Phys Chem A 2020; 124:7549-7558. [PMID: 32808782 DOI: 10.1021/acs.jpca.0c06756] [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/28/2022]
Abstract
In this work, we describe ab initio calculations and assignment of infrared (IR) spectra of hydrogen-bonded ion-molecular complexes that involve a fluxional proton: the linear N2H+···OC and N2D+···OC complexes. Given the challenges of describing fluxional proton dynamics and especially its IR activity, we use electric field-driven classical trajectories, i.e., the driven molecular dynamics (DMD) method that was developed by us in recent years and for similar applications, in conjunction with high-level electronic structure theory. Namely, we present a modified and a numerically efficient implementation of DMD specifically for direct (or "on the fly") calculations, which we carry out at the MP2-F12/AVDZ level of theory for the potential energy surface (PES) and MP2/AVDZ for the dipole moment surfaces (DMSs). Detailed analysis of the PES, DMS, and the time-dependence of the first derivative of the DMS, referred to as the driving force, for the highly fluxional vibrations involving H+/D+ revealed that the strongly non-harmonic PES and non-linear DMS yield remarkably complex vibrational spectra. Interestingly, the classical trajectories reveal a doublet in the proton transfer part of the spectrum with the two peaks at 1800 and 1980 cm-1. We find that their shared intensity is due to a Fermi-like resonance interaction, within the classical limit, of the H+ parallel stretch fundamental and an H+ perpendicular bending overtone. This doublet is also observed in the deuterated species at 1360 and 1460 cm-1.
Collapse
Affiliation(s)
- Dalton Boutwell
- Department of Chemistry & Biochemistry, Kennesaw State University, 370 Paulding Ave NW, Box # 1203, Kennesaw, Georgia 30144, United States
| | - Onyinye Okere
- Department of Chemistry & Biochemistry, Kennesaw State University, 370 Paulding Ave NW, Box # 1203, Kennesaw, Georgia 30144, United States
| | - Oluwaseun Omodemi
- Department of Chemistry & Biochemistry, Kennesaw State University, 370 Paulding Ave NW, Box # 1203, Kennesaw, Georgia 30144, United States
| | - Alexander Toledo
- Department of Chemistry & Biochemistry, Kennesaw State University, 370 Paulding Ave NW, Box # 1203, Kennesaw, Georgia 30144, United States
| | - Antonio Barrios
- Department of Chemistry & Biochemistry, Kennesaw State University, 370 Paulding Ave NW, Box # 1203, Kennesaw, Georgia 30144, United States
| | - Monique Olocha
- Department of Chemistry & Biochemistry, Kennesaw State University, 370 Paulding Ave NW, Box # 1203, Kennesaw, Georgia 30144, United States
| | - Martina Kaledin
- Department of Chemistry & Biochemistry, Kennesaw State University, 370 Paulding Ave NW, Box # 1203, Kennesaw, Georgia 30144, United States
| |
Collapse
|
20
|
Bakels S, Gaigeot MP, Rijs AM. Gas-Phase Infrared Spectroscopy of Neutral Peptides: Insights from the Far-IR and THz Domain. Chem Rev 2020; 120:3233-3260. [PMID: 32073261 PMCID: PMC7146864 DOI: 10.1021/acs.chemrev.9b00547] [Citation(s) in RCA: 34] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
![]()
Gas-phase, double
resonance IR spectroscopy has proven to be an
excellent approach to obtain structural information on peptides ranging
from single amino acids to large peptides and peptide clusters. In
this review, we discuss the state-of-the-art of infrared action spectroscopy
of peptides in the far-IR and THz regime. An introduction to the field
of far-IR spectroscopy is given, thereby highlighting the opportunities
that are provided for gas-phase research on neutral peptides. Current
experimental methods, including spectroscopic schemes, have been reviewed.
Structural information from the experimental far-IR spectra can be
obtained with the help of suitable theoretical approaches such as
dynamical DFT techniques and the recently developed Graph Theory.
The aim of this review is to underline how the synergy between far-IR
spectroscopy and theory can provide an unprecedented picture of the
structure of neutral biomolecules in the gas phase. The far-IR signatures
of the discussed studies are summarized in a far-IR map, in order
to gain insight into the origin of the far-IR localized and delocalized
motions present in peptides and where they can be found in the electromagnetic
spectrum.
Collapse
Affiliation(s)
- Sjors Bakels
- Radboud University, Institute for Molecules and Materials, FELIX Laboratory, Toernooiveld 7-c, 6525 ED Nijmegen, The Netherlands
| | - Marie-Pierre Gaigeot
- LAMBE CNRS UMR8587, Université d'Evry val d'Essonne, Blvd F. Mitterrand, Bât Maupertuis, 91025 Evry, France
| | - Anouk M Rijs
- Radboud University, Institute for Molecules and Materials, FELIX Laboratory, Toernooiveld 7-c, 6525 ED Nijmegen, The Netherlands
| |
Collapse
|
21
|
Chatterjee K, Dopfer O. Protonation of Naphthalene–(Water)n Nanoclusters: Intracluster Proton Transfer to Hydration Shell Revealed by Infrared Photodissociation Spectroscopy. J Phys Chem A 2020; 124:1134-1151. [DOI: 10.1021/acs.jpca.9b11779] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Kuntal Chatterjee
- Institut für Optik und Atomare Physik, Technische Universität Berlin, Hardenbergstrasse 36, 10623 Berlin, Germany
| | - Otto Dopfer
- Institut für Optik und Atomare Physik, Technische Universität Berlin, Hardenbergstrasse 36, 10623 Berlin, Germany
| |
Collapse
|
22
|
Ikeda K, Mahyuddin MH, Shiota Y, Staykov A, Matsumoto T, Ogo S, Yoshizawa K. Computational Study on the Light-Induced Oxidation of Iridium-Aqua Complex to Iridium-Oxo Complex over WO 3(001) Surface. Inorg Chem 2019; 59:415-422. [PMID: 31829576 DOI: 10.1021/acs.inorgchem.9b02704] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
An iridium aqua complex [IrIII(η5-C5Me5){bpy(COOH)2}(H2O)]2+ under visible light irradiation has been experimentally reported to form an iridium-oxo (Ir-oxo) complex [IrV(η5-C5Me5){bpy(COOH)2}(O)]2+, which oxidizes H2O to O2. However, the mechanism for the formation of this Ir-oxo complex remains unclear, due to the difficulties in observing the unstable Ir-oxo complex and computing light-induced systems having different numbers of electrons. In this study, we perform density functional theory (DFT) and time-dependent DFT (TD-DFT) calculations to investigate more in detail our previously proposed deprotonation and light-induced oxidation reactions composing the formation of the Ir-oxo complex. In particular, we discuss effects of light irradiation and WO3 support on the formation of the Ir-oxo complex. We suggest two distinct mechanisms, that is, direct and indirect for the light-induced oxidation. In the direct mechanism electrons are directly transferred from the occupied π* orbitals of IrIII-OH or IrIV=O• to the conduction band of the WO3 surface, whereas in the indirect mechanism electrons are first excited from the valence band to the conduction band of the WO3 surface due to the UV light, and then the resultant electron hole oxidizes the Ir complex. In the direct mechanism, in particular, we found that the lowest energy of the anode's conduction band determines the adsorption wavelength of the light irradiation, enabling us to predict alternative semiconductor anodes for more efficient formation of the Ir-oxo complex.
Collapse
Affiliation(s)
- Kei Ikeda
- Institute for Materials Chemistry and Engineering , Kyushu University , Fukuoka 819-0395 , Japan
| | - Muhammad Haris Mahyuddin
- Institute for Materials Chemistry and Engineering , Kyushu University , Fukuoka 819-0395 , Japan
| | - Yoshihito Shiota
- Institute for Materials Chemistry and Engineering , Kyushu University , Fukuoka 819-0395 , Japan
| | - Aleksandar Staykov
- International Institute for Carbon-Neutral Energy Research (WPI-I2CNER) , Kyushu University , Fukuoka 819-0395 , Japan
| | - Takahiro Matsumoto
- International Institute for Carbon-Neutral Energy Research (WPI-I2CNER) , Kyushu University , Fukuoka 819-0395 , Japan.,Department of Chemistry and Biochemistry, Graduate School of Engineering , Kyushu University , Fukuoka 819-0395 , Japan.,Center for Small Molecule Energy , Kyushu University , Fukuoka 819-0395 , Japan
| | - Seiji Ogo
- International Institute for Carbon-Neutral Energy Research (WPI-I2CNER) , Kyushu University , Fukuoka 819-0395 , Japan.,Department of Chemistry and Biochemistry, Graduate School of Engineering , Kyushu University , Fukuoka 819-0395 , Japan.,Center for Small Molecule Energy , Kyushu University , Fukuoka 819-0395 , Japan
| | - Kazunari Yoshizawa
- Institute for Materials Chemistry and Engineering , Kyushu University , Fukuoka 819-0395 , Japan.,International Institute for Carbon-Neutral Energy Research (WPI-I2CNER) , Kyushu University , Fukuoka 819-0395 , Japan.,Department of Chemistry and Biochemistry, Graduate School of Engineering , Kyushu University , Fukuoka 819-0395 , Japan.,Center for Small Molecule Energy , Kyushu University , Fukuoka 819-0395 , Japan
| |
Collapse
|
23
|
Yu Q, Carpenter WB, Lewis NHC, Tokmakoff A, Bowman JM. High-Level VSCF/VCI Calculations Decode the Vibrational Spectrum of the Aqueous Proton. J Phys Chem B 2019; 123:7214-7224. [PMID: 31361141 DOI: 10.1021/acs.jpcb.9b05723] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
The hydrated excess proton is a common species in aqueous chemistry, which complexes with water in a variety of structures. The infrared spectrum of the aqueous proton is particularly sensitive to this array of structures, which manifests as continuous IR absorption from 1000 to 3000 cm-1 known as the "proton continuum". Because of the extreme breadth of the continuum and strong anharmonicity of the involved vibrational modes, this spectrum has eluded straightforward interpretation and simulation. Using protonated water hexamer clusters from reactive molecular dynamics trajectories, and focusing on their central H+(H2O)2 structures' spectral contribution, we reproduce the linear IR spectrum of the aqueous proton with a high-level local monomer quantum method and highly accurate many-body potential energy surface. The accuracy of this approach is first verified in the vibrational spectra of the two isomers of the protonated water hexamer in the gas phase. We then apply this approach to 800 H+(H2O)6 clusters, also written as [H+(H2O)2](H2O)4, drawn from multistate empirical valence bond simulations of the bulk liquid to calculate the infrared spectrum of the aqueous proton complex. Incorporation of anharmonic effects to the vibrational potential and quantum mechanical treatment of the proton produces a better agreement to the infrared spectrum compared to that of the double-harmonic approximation. We assess the correlation of the proton stretching mode with different atomistic coordinates, finding the best correlation with ⟨ROH⟩, the expectation value of the proton-oxygen distance ROH. We also decompose the IR spectrum based on normal mode vibrations and ⟨ROH⟩ to provide insight on how different frequency regions in the continuum report on different configurations, vibrational modes, and mode couplings.
Collapse
Affiliation(s)
- Qi Yu
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation , Emory University , Atlanta , Georgia 30322 , United States
| | - William B Carpenter
- Department of Chemistry, James Frank Institute, and Institute for Biophysical Dynamics , The University of Chicago , Chicago , Illinois 60637 , United States
| | - Nicholas H C Lewis
- Department of Chemistry, James Frank Institute, and Institute for Biophysical Dynamics , The University of Chicago , Chicago , Illinois 60637 , United States
| | - Andrei Tokmakoff
- Department of Chemistry, James Frank Institute, and Institute for Biophysical Dynamics , The University of Chicago , Chicago , Illinois 60637 , United States
| | - Joel M Bowman
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation , Emory University , Atlanta , Georgia 30322 , United States
| |
Collapse
|
24
|
Hooper R, Boutwell D, Kaledin M. Assignment of Infrared-Active Combination Bands in the Vibrational Spectra of Protonated Molecular Clusters Using Driven Classical Trajectories: Application to N 4H + and N 4D . J Phys Chem A 2019; 123:5613-5620. [PMID: 31244119 DOI: 10.1021/acs.jpca.9b04480] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We investigate the utility of the driven molecular dynamics (DMD) approach to complex molecular vibrations by applying it to linear clusters with several degenerate vibrational modes and infrared (IR) intense combination bands. Here, the prominent features in N4H+ and N4D+ IR spectra, reported and described by others previously, have been characterized for the first time by DMD using recently published high-level potential and dipole moment surfaces. Namely, the calculations closely correlate the parallel proton stretch vibration in N4H+, at 750 cm-1, with the one observed experimentally at 743 cm-1. Second, the intense IR-active combination bands found in experimental spectra within 900-1100 cm-1 have been properly recovered by DMD at 950 cm-1 as strongly IR-active and confirmed as consisting of H+ asymmetric stretch and N2···N2 intermolecular symmetric stretch modes. Furthermore, we show that certain combination bands involving overtone transitions may be recovered by DMD using a hard-driving regime, such as the 1409 cm-1 band measured in N4H+, revealed by DMD at 1375 cm-1, and assigned to a progressive combination of the parallel H+ stretch and two quanta of N2···N2 stretch, in agreement with quantum mechanical studies reported previously by others.
Collapse
Affiliation(s)
- Reagan Hooper
- Department of Chemistry & Biochemistry , Kennesaw State University , 370 Paulding Avenue NW , Box # 1203, Kennesaw , Georgia 30144 , Unites States
| | - Dalton Boutwell
- Department of Chemistry & Biochemistry , Kennesaw State University , 370 Paulding Avenue NW , Box # 1203, Kennesaw , Georgia 30144 , Unites States
| | - Martina Kaledin
- Department of Chemistry & Biochemistry , Kennesaw State University , 370 Paulding Avenue NW , Box # 1203, Kennesaw , Georgia 30144 , Unites States
| |
Collapse
|
25
|
Quantum structural fluctuations of protonated water clusters (H2O) H+ (n = 1 − 4) studied by variational molecular dynamics method. J Mol Liq 2019. [DOI: 10.1016/j.molliq.2019.03.170] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
26
|
Yu Q, Bowman JM. Classical, Thermostated Ring Polymer, and Quantum VSCF/VCI Calculations of IR Spectra of H7O3+ and H9O4+ (Eigen) and Comparison with Experiment. J Phys Chem A 2019; 123:1399-1409. [DOI: 10.1021/acs.jpca.8b11603] [Citation(s) in RCA: 25] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Affiliation(s)
- Qi Yu
- Cherry L. Emerson Center for Scientific Computation and Department of Chemistry, Emory University, Atlanta, Georgia 30322, United States
| | - Joel M. Bowman
- Cherry L. Emerson Center for Scientific Computation and Department of Chemistry, Emory University, Atlanta, Georgia 30322, United States
| |
Collapse
|
27
|
Qu C, Bowman JM. Quantum approaches to vibrational dynamics and spectroscopy: is ease of interpretation sacrificed as rigor increases? Phys Chem Chem Phys 2019; 21:3397-3413. [DOI: 10.1039/c8cp04990d] [Citation(s) in RCA: 25] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
The subject of this Perspective is quantum approaches, beyond the harmonic approximation, to vibrational dynamics and IR spectroscopy.
Collapse
Affiliation(s)
- Chen Qu
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University
- Atlanta
- USA
| | - Joel M. Bowman
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University
- Atlanta
- USA
| |
Collapse
|
28
|
Chatterjee K, Dopfer O. Intracluster proton transfer in protonated benzonitrile–(H2O)n≤6 nanoclusters: hydrated hydronium core for n ≥ 2. Phys Chem Chem Phys 2019; 21:25226-25246. [DOI: 10.1039/c9cp05042f] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Infrared spectroscopy and density functional theory calculations of protonated benzonitrile–(H2O)n clusters reveal proton transfer to solvent for n ≥ 2 and the drastic effects of the aromatic dopant molecule on the network of H+(H2O)n+1.
Collapse
Affiliation(s)
- Kuntal Chatterjee
- Institut für Optik und Atomare Physik
- Technische Universität Berlin
- 10623 Berlin
- Germany
| | - Otto Dopfer
- Institut für Optik und Atomare Physik
- Technische Universität Berlin
- 10623 Berlin
- Germany
| |
Collapse
|
29
|
Galimberti DR, Bougueroua S, Mahé J, Tommasini M, Rijs AM, Gaigeot MP. Conformational assignment of gas phase peptides and their H-bonded complexes using far-IR/THz: IR-UV ion dip experiment, DFT-MD spectroscopy, and graph theory for mode assignment. Faraday Discuss 2019; 217:67-97. [DOI: 10.1039/c8fd00211h] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Graph theory based vibrational modes as new entities for vibrational THz spectroscopy.
Collapse
Affiliation(s)
| | | | - Jérôme Mahé
- LAMBE UMR8587
- Univ Evry
- Université Paris-Saclay
- CNRS
- 91025 Evry
| | - Matteo Tommasini
- Department of Chemistry, Materials, Chemical Engineering “G. Natta” Politecnico di Milano
- 20133 Milano
- Italy
| | - Anouk M. Rijs
- Radboud University
- Institute for Molecules and Materials
- FELIX Laboratory
- 6525 ED Nijmegen
- The Netherlands
| | | |
Collapse
|
30
|
Lewis NHC, Fournier JA, Carpenter WB, Tokmakoff A. Direct Observation of Ion Pairing in Aqueous Nitric Acid Using 2D Infrared Spectroscopy. J Phys Chem B 2018; 123:225-238. [DOI: 10.1021/acs.jpcb.8b10019] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Affiliation(s)
- Nicholas H. C. Lewis
- Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, United States
| | - Joseph A. Fournier
- Department of Chemistry, Washington University in St. Louis, St. Louis, Missouri 63130, United States
| | - William B. Carpenter
- Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, United States
| | - Andrei Tokmakoff
- Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, United States
| |
Collapse
|
31
|
Duong CH, Yang N, Kelleher PJ, Johnson MA, DiRisio RJ, McCoy AB, Yu Q, Bowman JM, Henderson BV, Jordan KD. Tag-Free and Isotopomer-Selective Vibrational Spectroscopy of the Cryogenically Cooled H9O4+ Cation with Two-Color, IR–IR Double-Resonance Photoexcitation: Isolating the Spectral Signature of a Single OH Group in the Hydronium Ion Core. J Phys Chem A 2018; 122:9275-9284. [DOI: 10.1021/acs.jpca.8b08507] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Affiliation(s)
- Chinh H. Duong
- Sterling Chemistry Laboratory, Yale University, New Haven, Connecticut 06520, United States
| | - Nan Yang
- Sterling Chemistry Laboratory, Yale University, New Haven, Connecticut 06520, United States
| | - Patrick J. Kelleher
- Sterling Chemistry Laboratory, Yale University, New Haven, Connecticut 06520, United States
| | - Mark A. Johnson
- Sterling Chemistry Laboratory, Yale University, New Haven, Connecticut 06520, United States
| | - Ryan J. DiRisio
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Anne B. McCoy
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Qi Yu
- Department of Chemistry and Cherry L. Emerson Center for Computational Science, Emory University, Atlanta, Georgia 30322, United States
| | - Joel M. Bowman
- Department of Chemistry and Cherry L. Emerson Center for Computational Science, Emory University, Atlanta, Georgia 30322, United States
| | - Bryan V. Henderson
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, United States
| | - Kenneth D. Jordan
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, United States
| |
Collapse
|
32
|
McDonald DC, Wagner JP, McCoy AB, Duncan MA. Near-Infrared Spectroscopy and Anharmonic Theory of Protonated Water Clusters: Higher Elevations in the Hydrogen Bonding Landscape. J Phys Chem Lett 2018; 9:5664-5671. [PMID: 30205006 DOI: 10.1021/acs.jpclett.8b02499] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
Near-infrared spectroscopy measurements are presented for protonated water clusters, H+(H2O) n, in the size range of n = 1-8. Clusters are produced in a pulsed-discharge supersonic expansion, mass selected, and studied with infrared laser photodissociation spectroscopy in the regions of 3600-4550 and 4850-7350 cm-1. Although there is some variation with cluster size, the main features of these spectra are a broad absorption near 5300 cm-1, a sharp doublet near 7200 cm-1, as well as a structured absorption near 4100 cm-1 for n ≥ 2. The vibrational patterns measured for the hydronium, Zundel, and Eigen ions are compared to those predicted by different forms of anharmonic theory. Second-order vibrational perturbation theory (VPT2) and a local mode treatment of the OH stretches both capture key aspects of the spectra but suffer understandable deficiencies in the quantitative description of band positions and intensities.
Collapse
Affiliation(s)
- D C McDonald
- Department of Chemistry , University of Georgia , Athens , Georgia 30602 , United States
| | - J P Wagner
- Department of Chemistry , University of Georgia , Athens , Georgia 30602 , United States
| | - A B McCoy
- Department of Chemistry , University of Washington , Seattle , Washington 98195 , United States
| | - M A Duncan
- Department of Chemistry , University of Georgia , Athens , Georgia 30602 , United States
| |
Collapse
|
33
|
|
34
|
Heindel JP, Yu Q, Bowman JM, Xantheas SS. Benchmark Electronic Structure Calculations for H3O+(H2O)n, n = 0–5, Clusters and Tests of an Existing 1,2,3-Body Potential Energy Surface with a New 4-Body Correction. J Chem Theory Comput 2018; 14:4553-4566. [DOI: 10.1021/acs.jctc.8b00598] [Citation(s) in RCA: 27] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Affiliation(s)
- Joseph P. Heindel
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Qi Yu
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| | - Joel M. Bowman
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| | - Sotiris S. Xantheas
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
- Advanced Computing, Mathematics and Data Division, Pacific Northwest National Laboratory, 902 Battelle Boulevard, P.O. Box
999, MS K1-83, Richland, Washington 99352, United States
| |
Collapse
|
35
|
Qu C, Bowman JM. IR Spectra of (HCOOH) 2 and (DCOOH) 2: Experiment, VSCF/VCI, and Ab Initio Molecular Dynamics Calculations Using Full-Dimensional Potential and Dipole Moment Surfaces. J Phys Chem Lett 2018; 9:2604-2610. [PMID: 29709189 DOI: 10.1021/acs.jpclett.8b00447] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
We report quantum VSCF/VCI and ab initio molecular dynamics (AIMD) calculations of the IR spectra of (HCOOH)2 and (DCOOH)2, using full-dimensional, ab initio potential energy and dipole moment surfaces (PES and DMS). These surfaces are fits, using permutationally invariant polynomials, to 13 475 ab initio CCSD(T)-F12a electronic energies and MP2 dipole moments. Here "AIMD" means using these ab initio potential and dipole moment surfaces in the MD calculations. The VSCF/VCI calculations use all (24) normal modes for coupling, with a four-mode representation of the potential. The quantum spectra align well with jet-cooled and room-temperature experimental spectra over the spectral range 600-3600 cm-1. Analyses of the complex O-H and C-H stretch bands are made based on the mixing of the VSCF/VCI basis functions. The comparisons of the AIMD IR spectra with both experimental and VSCF/VCI ones provide tests of the accuracy of the AIMD approach. These indicate good accuracy for simple bands but not for the complex O-H stretch band, which is upshifted from experimental and VSCF/VCI bands by roughly 300 cm-1. In addition to testing the AIMD approach, the PES, DMS, and VSCF/VCI calculations for formic acid dimer provide opportunities for testing other methods to represent high-dimensional data and other methods that perform postharmonic vibrational calculations.
Collapse
Affiliation(s)
- Chen Qu
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation , Emory University , Atlanta , Georgia 30322 , United States
| | - Joel M Bowman
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation , Emory University , Atlanta , Georgia 30322 , United States
| |
Collapse
|
36
|
Qu C, Bowman J. Quantum and classical IR spectra of (HCOOH)2, (DCOOH)2 and (DCOOD)2 using ab initio potential energy and dipole moment surfaces. Faraday Discuss 2018; 212:33-49. [DOI: 10.1039/c8fd00077h] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Abstract
Full-dimensional (24 modes) quantum calculation of the IR spectrum of (DCOOD)2, and comparison with classical MD one.
Collapse
Affiliation(s)
- Chen Qu
- Department of Chemistry
- Cherry L. Emerson Center for Scientific Computations
- Emory University
- Atlanta
- USA
| | - Joel M. Bowman
- Department of Chemistry
- Cherry L. Emerson Center for Scientific Computations
- Emory University
- Atlanta
- USA
| |
Collapse
|