1
|
Yu Q, Bowman JM. Fully Quantum Simulation of Polaritonic Vibrational Spectra of Large Cavity-Molecule System. J Chem Theory Comput 2024; 20:4278-4287. [PMID: 38717309 DOI: 10.1021/acs.jctc.4c00129] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/29/2024]
Abstract
The formation of molecular vibrational polaritons, arising from the interplay between molecular vibrations and infrared cavity modes, is a quantum phenomenon necessitating accurate quantum dynamical simulations. Here, we introduce the cavity vibrational self-consistent field/virtual state configuration interaction method, enabling quantum simulation of the vibrational spectra of many-molecule systems within the optical cavity. Focusing on the representative (H2O)21 system, we showcase this parameter-free quantum approach's ability to capture both linear and nonlinear vibrational spectral features. Our findings highlight the growing prominence of molecular couplings among OH stretches and bending excited bands with increased light-matter interaction, revealing distinctive nonlinear spectral features induced by vibrational strong coupling.
Collapse
Affiliation(s)
- Qi Yu
- Department of Chemistry, Fudan University, Shanghai 200438, P. R. China
| | - Joel M Bowman
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| |
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
|
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
|
5
|
Dean JLS, Winkler VS, Boyer MA, Sibert EL, Fournier JA. Investigating Intramolecular H Atom Transfer Dynamics in β-Diketones with Ultrafast Infrared Spectroscopies and Theoretical Modeling. J Phys Chem A 2023; 127:9258-9272. [PMID: 37882618 DOI: 10.1021/acs.jpca.3c05417] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2023]
Abstract
The vibrational signatures and ultrafast dynamics of the intramolecular H-bond in a series of β-diketones are investigated with 2D IR spectroscopy and computational modeling. The chosen β-diketones exhibit a range of H atom donor-acceptor distances and asymmetry along the H atom transfer coordinate that tunes the intramolecular H-bond strength. The species with the strongest H-bonds are calculated to have very soft H atom potentials, resulting in highly red-shifted OH stretch fundamental frequencies and dislocation of the H atom upon vibrational excitation. These soft potentials lead to significant coupling to the other normal mode coordinates and give rise to the very broad vibrational signatures observed experimentally. The 2D IR spectra in both the OH and OD stretch regions of the light and deuterated isotopologues reveal broadened and long-lived ground-state bleach signatures of the vibrationally hot molecules. Polarization-sensitive transient absorption measurements in the OH and OD stretch regions reveal notable isotopic differences in orientational dynamics. Orientational relaxation was measured to occur on ∼600 fs and ∼2 ps time scales for the light and deuterated isotopologues, respectively. The orientational dynamics are interpreted in terms of activated H/D atom transfer events driven by collective intramolecular structural rearrangements.
Collapse
Affiliation(s)
- Jessika L S Dean
- Department of Chemistry, Washington University in St. Louis, St. Louis, Missouri 63130, United States
| | - Valerie S Winkler
- Department of Chemistry, Washington University in St. Louis, St. Louis, Missouri 63130, United States
| | - Mark A Boyer
- Department of Chemistry and Theoretical Chemistry Institute, University of Wisconsin Madison, Madison, Wisconsin 53706, United States
| | - Edwin L Sibert
- Department of Chemistry and Theoretical Chemistry Institute, University of Wisconsin Madison, Madison, Wisconsin 53706, United States
| | - Joseph A Fournier
- Department of Chemistry, Washington University in St. Louis, St. Louis, Missouri 63130, United States
| |
Collapse
|
6
|
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
|
7
|
Muñoz-Santiburcio D. Accurate diffusion coefficients of the excess proton and hydroxide in water via extensive ab initio simulations with different schemes. J Chem Phys 2022; 157:024504. [PMID: 35840376 DOI: 10.1063/5.0093958] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Despite its simple molecular formula, obtaining an accurate in silico description of water is far from straightforward. Many of its very peculiar properties are quite elusive, and in particular, obtaining good estimations of the diffusion coefficients of the solvated proton and hydroxide at a reasonable computational cost has been an unsolved challenge until now. Here, I present extensive results of several unusually long ab initio molecular dynamics (MD) simulations employing different combinations of the Born-Oppenheimer and second-generation Car-Parrinello MD propagation methods with different ensembles (NVE and NVT) and thermostats, which show that these methods together with the RPBE-D3 functional provide a very accurate estimation of the diffusion coefficients of the solvated H3O+ and OH- ions, together with an extremely accurate description of several properties of neutral water (such as the structure of the liquid and its diffusion and shear viscosity coefficients). In addition, I show that the estimations of DH3O+ and DOH- depend dramatically on the simulation length, being necessary to reach timescales in the order of hundreds of picoseconds to obtain reliable results.
Collapse
Affiliation(s)
- Daniel Muñoz-Santiburcio
- CIC nanoGUNE BRTA, Tolosa Hiribidea 76, 20018 San Sebastián, Spain and Instituto de Fusión Nuclear "Guillermo Velarde," Universidad Politécnica de Madrid, C/ José Gutiérrez Abascal 2, 28006 Madrid, Spain
| |
Collapse
|
8
|
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
|
9
|
Meuwly M. Atomistic Simulations for Reactions and Vibrational Spectroscopy in the Era of Machine Learning─ Quo Vadis?. J Phys Chem B 2022; 126:2155-2167. [PMID: 35286087 DOI: 10.1021/acs.jpcb.2c00212] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Atomistic simulations using accurate energy functions can provide molecular-level insight into functional motions of molecules in the gas and in the condensed phase. This Perspective delineates the present status of the field from the efforts of others and some of our own work and discusses open questions and future prospects. The combination of physics-based long-range representations using multipolar charge distributions and kernel representations for the bonded interactions is shown to provide realistic models for the exploration of the infrared spectroscopy of molecules in solution. For reactions, empirical models connecting dedicated energy functions for the reactant and product states allow statistically meaningful sampling of conformational space whereas machine-learned energy functions are superior in accuracy. The future combination of physics-based models with machine-learning techniques and integration into all-purpose molecular simulation software provides a unique opportunity to bring such dynamics simulations closer to reality.
Collapse
Affiliation(s)
- Markus Meuwly
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, 4056 Basel, Switzerland
| |
Collapse
|
10
|
Theoretical Description of Water from Single-Molecule to Condensed Phase: a Review of Recent Progress on Potential Energy Surfaces and Molecular Dynamics. CHINESE J CHEM PHYS 2022. [DOI: 10.1063/1674-0068/cjcp2201005] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022]
|
11
|
Towards complete assignment of the infrared spectrum of the protonated water cluster H +(H 2O) 21. Nat Commun 2021; 12:6141. [PMID: 34686665 PMCID: PMC8536673 DOI: 10.1038/s41467-021-26284-x] [Citation(s) in RCA: 27] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/24/2021] [Accepted: 09/22/2021] [Indexed: 11/08/2022] Open
Abstract
The spectroscopic features of protonated water species in dilute acid solutions have been long sought after for understanding the microscopic behavior of the proton in water with gas-phase water clusters H+(H2O)n extensively studied as bottom-up model systems. We present a new protocol for the calculation of the infrared (IR) spectra of complex systems, which combines the fragment-based Coupled Cluster method and anharmonic vibrational quasi-degenerate perturbation theory, and demonstrate its accuracy towards the complete and accurate assignment of the IR spectrum of the H+(H2O)21 cluster. The site-specific IR spectral signatures reveal two distinct structures for the internal and surface four-coordinated water molecules, which are ice-like and liquid-like, respectively. The effect of inter-molecular interaction between water molecules is addressed, and the vibrational resonance is found between the O-H stretching fundamental and the bending overtone of the nearest neighboring water molecule. The revelation of the spectral signature of the excess proton offers deeper insight into the nature of charge accommodation in the extended hydrogen-bonding network underpinning this aqueous cluster. Protonated water species have been the subject of numerous experimental and computational studies. Here the authors provide a nearly complete assignment of the experimental IR spectrum of the H+(H2O)21 water cluster based on high-level wavefunction theory and anharmonic vibrational quasi-degenerate perturbation theory.
Collapse
|
12
|
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
|
13
|
Inakollu VSS, Yu H. Comparative studies of IR spectra of deprotonated serine with classical and thermostated ring polymer molecular dynamics simulations. STRUCTURAL DYNAMICS (MELVILLE, N.Y.) 2021; 8:054101. [PMID: 34549074 PMCID: PMC8443303 DOI: 10.1063/4.0000124] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 06/17/2021] [Accepted: 08/30/2021] [Indexed: 06/13/2023]
Abstract
Here we report the vibrational spectra of deprotonated serine calculated from the classical molecular dynamics (MD) simulations and thermostated ring-polymer molecular dynamics (TRPMD) simulation with third-order density-functional tight-binding. In our earlier study [Inakollu and Yu, "A systematic benchmarking of computational vibrational spectroscopy with DFTB3: Normal mode analysis and fast Fourier transform dipole autocorrelation function," J. Comput. Chem. 39, 2067 (2018)] of deprotonated serine, we observed a significant difference in the vibrational spectra with the classical MD simulations compared to the infrared multiple photon dissociation spectra. It was postulated that this is due to neglecting the nuclear quantum effects (NQEs). In this work, NQEs are considered in spectral calculation using the TRPMD simulations. With the help of potential of mean force calculations, the conformational space of deprotonated serine is analyzed and used to understand the difference in the spectra of classical MD and TRPMD simulations at 298.15 and 100 K. The high-frequency vibrational bands in the spectra are characterized using Fourier transform localized vibrational mode (FT-νN AC) and interatomic distance histograms. At room temperature, the quantum effects are less significant, and the free energy profiles in the classical MD and the TRPMD simulations are very similar. However, the hydrogen bond between the hydroxyl-carboxyl bond is slightly stronger in TRPMD simulations. At 100 K, the quantum effects are more prominent, especially in the 2600-3600 cm-1, and the free energy profile slightly differs between the classical MD and TRPMD simulations. Using the FT-νN AC and the interatomic distance histograms, the high-frequency vibrational bands are discussed in detail.
Collapse
Affiliation(s)
| | - Haibo Yu
- Author to whom correspondence should be addressed:
| |
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
|
Pavošević F, Tao Z, Hammes-Schiffer S. Multicomponent Coupled Cluster Singles and Doubles with Density Fitting: Protonated Water Tetramers with Quantized Protons. J Phys Chem Lett 2021; 12:1631-1637. [PMID: 33555187 DOI: 10.1021/acs.jpclett.0c03771] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
Nuclear quantum effects such as zero-point energy are important for describing a wide range of chemical properties. The nuclear-electronic orbital (NEO) approach incorporates such effects into quantum chemistry calculations by treating specified nuclei, typically protons, quantum mechanically on the same level as electrons. Herein, both the traditional and t1-transformed NEO coupled cluster with singles and doubles (NEO-CCSD) methods are implemented with a density fitting (DF) scheme for approximating the four-center two-particle integrals. The enhanced computational efficiency enables calculations on larger molecules with multiple quantum protons. The NEO-DF-CCSD method predicts proton affinities within chemical accuracy. Its application to protonated water tetramers with all nine protons treated quantum mechanically produces the qualitatively correct ordering of the isomer energies, which are strongly influenced by the zero-point energy contributions inherently included in NEO energy calculations. This work showcases the capabilities of the NEO-DF-CCSD method and provides the foundation for future developments and applications.
Collapse
Affiliation(s)
- Fabijan Pavošević
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, United States
| | - Zhen Tao
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, United States
| | - Sharon Hammes-Schiffer
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, United States
| |
Collapse
|
17
|
Kwan V, Consta S. Molecular Characterization of the Surface Excess Charge Layer in Droplets. JOURNAL OF THE AMERICAN SOCIETY FOR MASS SPECTROMETRY 2021; 32:33-45. [PMID: 32597645 DOI: 10.1021/jasms.0c00053] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
The surface excess charge layer (SECL) in droplets has often been associated with distinct chemistry. We examine the effect of the nature of ions in the composition and structure of SECL by using molecular dynamics. We find that in the presence of simple ions the thickness of SECL is invariant not only with respect to droplet size but also with respect to the nature of the ions. In the presence of simple ions, this layer has a thickness of ∼1.5-1.7 nm but in the presence of macroions it may extend to ∼2.0 nm. The proportion of ions contained in SECL depends on the nature of the ions and the droplet size. For the same droplet size, I- and model H3O+ ions show considerably higher concentration than Na+ and Cl- ions. We identify the maximum ion concentration region, which, in nanodrops, may partially overlap with SECL. As the relative shape fluctuations decrease when microdrop size is approached, the overlap between SECL and maximum ion concentration region increases. We suggest the extension of the bilayer droplet structure assumed in the equilibrium partitioning model of Enke to include the maximum ion concentration region that may not coincide with SECL in nanodrops. We compute the ion concentrations in SECL, which are those that should enter the kinetic equation in the ion-evaporation mechanism, instead of the overall drop ion concentration that has been used.
Collapse
Affiliation(s)
- Victor Kwan
- Department of Chemistry, The University of Western Ontario, London, Ontario, Canada N6A 5B7
| | - Styliani Consta
- Department of Chemistry, The University of Western Ontario, London, Ontario, Canada N6A 5B7
| |
Collapse
|
18
|
Huang QR, Li YC, Nishigori T, Katada M, Fujii A, Kuo JL. Vibrational Coupling in Solvated H 3O +: Interplay between Fermi Resonance and Combination Band. J Phys Chem Lett 2020; 11:10067-10072. [PMID: 33179938 DOI: 10.1021/acs.jpclett.0c03059] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
Complex vibrational features of solvated hydronium ion, H3O+, in 3 μm enable us to look into the vibrational coupling among O-H stretching modes and other degrees of freedom. Two anharmonic coupling schemes have often been engaged to explain observed spectra: coupling with the OH bending overtone, known as Fermi resonance (FR), has been proposed to account for the splitting of the OH stretch band at ∼3300 cm-1 in H3O+···Ar3, but an additional peak in H3O+···(N2)3 at the similar frequency region has been assigned to a combination band (CB) with the low-frequency intermolecular stretches. While even stronger vibrational coupling is expected in H3O+···(H2O)3, such pronounced peaks are absent. In the present study, vibrational spectra of H3O+···Kr3 and H3O+···(CO)3 are measured to complement the existing spectra. Using ab initio anharmonic algorithms, we are able to assign the observed complex spectral features, to resolve seemingly contradictory notions in the interpretations, and to reveal simple pictures of the interplay between FR and CB.
Collapse
Affiliation(s)
- Qian-Rui Huang
- Institute of Atomic and Molecular Sciences, Academia Sinica, Taipei, Taiwan
| | - Ying-Cheng Li
- Institute of Atomic and Molecular Sciences, Academia Sinica, Taipei, Taiwan
| | - Tomoki Nishigori
- Department of Chemistry, Graduate School of Science, Tohoku University, Sendai, Japan
| | - Marusu Katada
- 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
|
19
|
Carpenter WB, Yu Q, Hack JH, Dereka B, Bowman JM, Tokmakoff A. Decoding the 2D IR spectrum of the aqueous proton with high-level VSCF/VCI calculations. J Chem Phys 2020; 153:124506. [PMID: 33003749 DOI: 10.1063/5.0020279] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023] Open
Abstract
The aqueous proton is a common and long-studied species in chemistry, yet there is currently intense interest devoted to understanding its hydration structure and transport dynamics. Typically described in terms of two limiting structures observed in gas-phase clusters, the Zundel H5O2 + and Eigen H9O4 + ions, the aqueous structure is less clear due to the heterogeneity of hydrogen bonding environments and room-temperature structural fluctuations in water. The linear infrared (IR) spectrum, which reports on structural configurations, is challenging to interpret because it appears as a continuum of absorption, and the underlying vibrational modes are strongly anharmonically coupled to each other. Recent two-dimensional IR (2D IR) experiments presented strong evidence for asymmetric Zundel-like motifs in solution, but true structure-spectrum correlations are missing and complicated by the anharmonicity of the system. In this study, we employ high-level vibrational self-consistent field/virtual state configuration interaction calculations to demonstrate that the 2D IR spectrum reports on a broad distribution of geometric configurations of the aqueous proton. We find that the diagonal 2D IR spectrum around 1200 cm-1 is dominated by the proton stretch vibrations of Zundel-like and intermediate geometries, broadened by the heterogeneity of aqueous configurations. There is a wide distribution of multidimensional potential shapes for the proton stretching vibration with varying degrees of potential asymmetry and confinement. Finally, we find specific cross peak patterns due to aqueous Zundel-like species. These studies provide clarity on highly debated spectral assignments and stringent spectroscopic benchmarks for future simulations.
Collapse
Affiliation(s)
- William B Carpenter
- Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, USA
| | - Qi Yu
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| | - John H Hack
- Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, USA
| | - Bogdan Dereka
- Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, USA
| | - Joel M Bowman
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| | - Andrei Tokmakoff
- Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, USA
| |
Collapse
|
20
|
Tyrode E, Sengupta S, Sthoer A. Identifying Eigen-like hydrated protons at negatively charged interfaces. Nat Commun 2020; 11:493. [PMID: 31980619 PMCID: PMC6981112 DOI: 10.1038/s41467-020-14370-5] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/23/2019] [Accepted: 01/06/2020] [Indexed: 11/24/2022] Open
Abstract
Despite the importance of the hydrogen ion in a wide range of biological, chemical, and physical processes, its molecular structure in solution remains lively debated. Progress has been primarily hampered by the extreme diffuse nature of the vibrational signatures of hydrated protons in bulk solution. Using the inherently surface-specific vibrational sum frequency spectroscopy technique, we show that at selected negatively charged interfaces, a resolved spectral feature directly linked to the H3O+ core in an Eigen-like species can be readily identified in a biologically compatible pH range. Centered at ~2540 cm−1, the band is seen to shift to ~1875 cm−1 when forming D3O+ upon isotopic substitution. The results offer the possibility of tracking and understanding from a molecular perspective the behavior of hydrated protons at charged interfaces. Hydrated protons are always present in aqueous solution, but their molecular structure remains under debate. Here the authors use vibrational sum frequency spectroscopy to show that at negatively charged liquid–vapor interfaces, protons adopt a specific configuration characteristic of Eigen-like species.
Collapse
Affiliation(s)
- Eric Tyrode
- Department of Chemistry, KTH Royal Institute of Technology, SE-10044, Stockholm, Sweden.
| | - Sanghamitra Sengupta
- Department of Chemistry, KTH Royal Institute of Technology, SE-10044, Stockholm, Sweden
| | - Adrien Sthoer
- Department of Chemistry, KTH Royal Institute of Technology, SE-10044, Stockholm, Sweden
| |
Collapse
|
21
|
Yu Q, Bowman JM. Tracking Hydronium/Water Stretches in Magic H3O+(H2O)20 Clusters through High-level Quantum VSCF/VCI Calculations. J Phys Chem A 2020; 124:1167-1175. [DOI: 10.1021/acs.jpca.9b11983] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/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
|
22
|
Lebedev AV. The H3O+(H2O)n Reagent Ion: Calculations of the Structure, Thermodynamic Parameters of Hydration, Equilibrium Composition, and Mobility. JOURNAL OF ANALYTICAL CHEMISTRY 2019. [DOI: 10.1134/s1061934819130082] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
23
|
Bertaina G, Di Liberto G, Ceotto M. Reduced rovibrational coupling Cartesian dynamics for semiclassical calculations: Application to the spectrum of the Zundel cation. J Chem Phys 2019; 151:114307. [PMID: 31542046 DOI: 10.1063/1.5114616] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We study the vibrational spectrum of the protonated water dimer, by means of a divide-and-conquer semiclassical initial value representation of the quantum propagator, as a first step in the study of larger protonated water clusters. We use the potential energy surface from the work of Huang et al. [J. Chem. Phys. 122, 044308 (2005)]. To tackle such an anharmonic and floppy molecule, we employ fully Cartesian dynamics and carefully reduce the coupling to global rotations in the definition of normal modes. We apply the time-averaging filter and obtain clean power spectra relative to suitable reference states that highlight the spectral peaks corresponding to the fundamental excitations of the system. Our trajectory-based approach allows for the physical interpretation of the very challenging proton transfer modes. We find that it is important, for such a floppy molecule, to selectively avoid initially exciting lower energy modes, in order to obtain cleaner spectra. The estimated vibrational energies display a mean absolute error (MAE) of ∼29 cm-1 with respect to available multiconfiguration time-dependent Hartree calculations and MAE ∼ 14 cm-1 when compared to the optically active experimental excitations of the Ne-tagged Zundel cation. The reasonable scaling in the number of trajectories for Monte Carlo convergence is promising for applications to higher dimensional protonated cluster systems.
Collapse
Affiliation(s)
- G Bertaina
- Dipartimento di Chimica, Università degli Studi di Milano, via C. Golgi 19, 20133 Milano, Italy
| | - G Di Liberto
- Dipartimento di Chimica, Università degli Studi di Milano, via C. Golgi 19, 20133 Milano, Italy
| | - M Ceotto
- Dipartimento di Chimica, Università degli Studi di Milano, via C. Golgi 19, 20133 Milano, Italy
| |
Collapse
|
24
|
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
|
25
|
|
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
|
Dahms F, Kundu A, Pines E, Fingerhut BP, Nibbering ETJ, Elsaesser T. Ultrafast dynamics of hydrated excess protons in CH 3CN:H 2O mixtures. EPJ WEB OF CONFERENCES 2019. [DOI: 10.1051/epjconf/201920509003] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
In a combined experimental and theoretical 2D-IR and pump-probe study we determine how ultrafast solvent motions govern the vibrational dynamics of the hydrated proton and the key role played by the underlying proton potential.
Collapse
|
29
|
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
|
30
|
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
|
31
|
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
|
32
|
Shi R, Li K, Su Y, Tang L, Huang X, Sai L, Zhao J. Revisit the landscape of protonated water clusters H +(H 2O) n with n = 10-17: An ab initio global search. J Chem Phys 2018; 148:174305. [PMID: 29739201 DOI: 10.1063/1.5026383] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023] Open
Abstract
Using a genetic algorithm incorporated with density functional theory, we explore the ground state structures of protonated water clusters H+(H2O)n with n = 10-17. Then we re-optimize the isomers at B97-D/aug-cc-pVDZ level of theory. The extra proton connects with a H2O molecule to form a H3O+ ion in all H+(H2O)10-17 clusters. The lowest-energy structures adopt a monocage form at n = 10-16 and core-shell structure at n = 17 based on the MP2/aug-cc-pVTZ//B97-D/aug-cc-pVDZ+ZPE single-point-energy calculation. Using second-order vibrational perturbation theory, we further calculate the infrared spectra with anharmonic correction for the ground state structures of H+(H2O)10-17 clusters at the PBE0/aug-cc-pVDZ level. The anharmonic correction to the spectra is crucial since it reproduces the experimental results quite well. The extra proton weakens the O-H bond strength in the H3O+ ion since the Wiberg bond order of the O-H bond in the H3O+ ion is smaller than that in H2O molecules, which causes a red shift of the O-H stretching mode in the H3O+ ion.
Collapse
Affiliation(s)
- Ruili Shi
- Key Laboratory of Materials Modification by Laser, Ion and Electron Beams (Dalian University of Technology), Ministry of Education, Dalian 116024, China
| | - Keyao Li
- Key Laboratory of Materials Modification by Laser, Ion and Electron Beams (Dalian University of Technology), Ministry of Education, Dalian 116024, China
| | - Yan Su
- Key Laboratory of Materials Modification by Laser, Ion and Electron Beams (Dalian University of Technology), Ministry of Education, Dalian 116024, China
| | - Lingli Tang
- College of Science, Dalian Minzu University, Dalian 116600, China
| | - Xiaoming Huang
- School of Ocean Science and Technology, Dalian University of Technology, Panjin Campus, Panjin 124221, China
| | - Linwei Sai
- Department of Mathematics and Physics, Hohai University, Changzhou 213022, China
| | - Jijun Zhao
- Key Laboratory of Materials Modification by Laser, Ion and Electron Beams (Dalian University of Technology), Ministry of Education, Dalian 116024, China
| |
Collapse
|
33
|
Affiliation(s)
- Lorenz S. Cederbaum
- Theoretische Chemie, Physikalisch–Chemisches Institut, Heidelberg University, Heidelberg, Germany
| |
Collapse
|
34
|
Affiliation(s)
- Chen Qu
- Department of Chemistry, Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| | - Qi Yu
- Department of Chemistry, Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| | - Joel M. Bowman
- Department of Chemistry, Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| |
Collapse
|
35
|
Di Liberto G, Conte R, Ceotto M. “Divide-and-conquer” semiclassical molecular dynamics: An application to water clusters. J Chem Phys 2018; 148:104302. [DOI: 10.1063/1.5023155] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Giovanni Di Liberto
- Dipartimento di Chimica, Università degli Studi di Milano, Via C. Golgi 19, 20133 Milano, Italy
| | - Riccardo Conte
- Dipartimento di Chimica, Università degli Studi di Milano, Via C. Golgi 19, 20133 Milano, Italy
| | - Michele Ceotto
- Dipartimento di Chimica, Università degli Studi di Milano, Via C. Golgi 19, 20133 Milano, Italy
| |
Collapse
|
36
|
Esser TK, Knorke H, Asmis KR, Schöllkopf W, Yu Q, Qu C, Bowman JM, Kaledin M. Deconstructing Prominent Bands in the Terahertz Spectra of H 7O 3+ and H 9O 4+: Intermolecular Modes in Eigen Clusters. J Phys Chem Lett 2018; 9:798-803. [PMID: 29360366 DOI: 10.1021/acs.jpclett.7b03395] [Citation(s) in RCA: 32] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
We report experimental vibrational action spectra (210-2200 cm-1) and calculated IR spectra, using recent ab initio potential energy and dipole moment surfaces, of H7O3+ and H9O4+. We focus on prominent far-IR bands, which postharmonic analyses show, arise from two types of intermolecular motions, i.e., hydrogen bond stretching and terminal water wagging modes, that are similar in both clusters. The good agreement between experiment and theory further validates the accuracy of the potential and dipole moment surfaces, which was used in a recent theoretical study that concluded that infrared photodissociation spectra of the cold clusters correspond to the Eigen isomer. The comparison between theory and experiment adds further confirmation of the need of postharmonic analysis for these clusters.
Collapse
Affiliation(s)
- Tim K Esser
- Wilhelm-Ostwald-Institut für Physikalische und Theoretische Chemie, Universität Leipzig , Linnéstraße 2, 04103 Leipzig, Germany
| | - Harald Knorke
- Wilhelm-Ostwald-Institut für Physikalische und Theoretische Chemie, Universität Leipzig , Linnéstraße 2, 04103 Leipzig, Germany
| | - Knut R Asmis
- Wilhelm-Ostwald-Institut für Physikalische und Theoretische Chemie, Universität Leipzig , Linnéstraße 2, 04103 Leipzig, Germany
| | - Wieland Schöllkopf
- Fritz-Haber-Institut der Max-Planck-Gesellschaft , Faradayweg 4-6, D-14195 Berlin, Germany
| | - Qi Yu
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University , Atlanta, Georgia 30322, United States
| | - 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
| | - Martina Kaledin
- Department of Chemistry and Biochemistry, Kennesaw State University , Kennesaw, Georgia 30144, United States
| |
Collapse
|
37
|
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
|
38
|
DeGregorio N, Iyengar SS. Efficient and Adaptive Methods for Computing Accurate Potential Surfaces for Quantum Nuclear Effects: Applications to Hydrogen-Transfer Reactions. J Chem Theory Comput 2017; 14:30-47. [DOI: 10.1021/acs.jctc.7b00927] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Affiliation(s)
- Nicole DeGregorio
- Department of Chemistry and
Department of Physics, Indiana University, 800 E. Kirkwood Avenue, Bloomington, Indiana 47405, United States
| | - Srinivasan S. Iyengar
- Department of Chemistry and
Department of Physics, Indiana University, 800 E. Kirkwood Avenue, Bloomington, Indiana 47405, United States
| |
Collapse
|