1
|
Spencer RJ, Zhanserkeev AA, Yang EL, Steele RP. The Near-Sightedness of Many-Body Interactions in Anharmonic Vibrational Couplings. J Am Chem Soc 2024; 146:15376-15392. [PMID: 38771156 DOI: 10.1021/jacs.4c03198] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/22/2024]
Abstract
Couplings between vibrational motions are driven by electronic interactions, and these couplings carry special significance in vibrational energy transfer, multidimensional spectroscopy experiments, and simulations of vibrational spectra. In this investigation, the many-body contributions to these couplings are analyzed computationally in the context of clathrate-like alkali metal cation hydrates, including Cs+(H2O)20, Rb+(H2O)20, and K+(H2O)20, using both analytic and quantum-chemistry potential energy surfaces. Although the harmonic spectra and one-dimensional anharmonic spectra depend strongly on these many-body interactions, the mode-pair couplings were, perhaps surprisingly, found to be dominated by one-body effects, even in cases of couplings to low-frequency modes that involved the motion of multiple water molecules. The origin of this effect was traced mainly to geometric distortion within water monomers and cancellation of many-body effects in differential couplings, and the effect was also shown to be agnostic to the identity of the ion. These outcomes provide new understanding of vibrational couplings and suggest the possibility of improved computational methods for the simulation of infrared and Raman spectra.
Collapse
Affiliation(s)
- Ryan J Spencer
- Department of Chemistry and Henry Eyring Center for Theoretical Chemistry, University of Utah, Salt Lake City, Utah 84112, United States
| | - Asylbek A Zhanserkeev
- Department of Chemistry and Henry Eyring Center for Theoretical Chemistry, University of Utah, Salt Lake City, Utah 84112, United States
| | - Emily L Yang
- Department of Chemistry and Henry Eyring Center for Theoretical Chemistry, University of Utah, Salt Lake City, Utah 84112, United States
| | - Ryan P Steele
- Department of Chemistry and Henry Eyring Center for Theoretical Chemistry, University of Utah, Salt Lake City, Utah 84112, United States
| |
Collapse
|
2
|
Hunter KM, Paesani F. Monitoring water harvesting in metal-organic frameworks, one water molecule at a time. Chem Sci 2024; 15:5303-5310. [PMID: 38577368 PMCID: PMC10988614 DOI: 10.1039/d3sc06162k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/17/2023] [Accepted: 03/05/2024] [Indexed: 04/06/2024] Open
Abstract
Metal-organic frameworks (MOFs) have gained prominence as potential materials for atmospheric water harvesting, a vital solution for arid regions and areas experiencing severe water shortages. However, the molecular factors influencing the performance of MOFs in capturing water from the air remain elusive. Among all MOFs, Ni2X2BTDD (X = F, Cl, Br) stands out as a promising water harvester due to its ability to adsorb substantial amounts of water at low relative humidity (RH). Here, we use advanced molecular dynamics simulations carried out with the state-of-the-art MB-pol data-driven many-body potential to monitor water adsorption in the three Ni2X2BTDD variants as a function of RH. Our simulations reveal that the type of halide atom in the three Ni2X2BTDD frameworks significantly influences the corresponding molecular mechanisms of water adsorption: while water molecules form strong hydrogen bonds with the fluoride atoms in Ni2F2BTDD, they tend to form hydrogen bonds with the nitrogen atoms of the triazolate linkers in Ni2Cl2BTDD and Ni2Br2BTDD. Importantly, the large size of the bromide atoms reduces the void volume in the Ni2Br2BTDD pores, which enable water molecules to initiate an extended hydrogen-bond network at lower RH. These findings not only underscore the prospect for precisely tuning structural and chemical modifications of the frameworks to optimize their interaction with water, but also highlight the predictive power of simulations with the MB-pol data-driven many-body potential. By providing a realistic description of water under different thermodynamic conditions and environments, these simulations yield unique, molecular-level insights that can guide the design and optimization of energy-efficient water harvesting materials.
Collapse
Affiliation(s)
- Kelly M Hunter
- Department of Chemistry and Biochemistry, University of California La Jolla San Diego California 92093 USA
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, University of California La Jolla San Diego California 92093 USA
- Materials Science and Engineering, University of California La Jolla San Diego California 92093 USA
- Halicioğlu Data Science Institute, University of California La Jolla San Diego California 92093 USA
- San Diego Supercomputer Center, University of California La Jolla San Diego California 92093 USA
| |
Collapse
|
3
|
Streu K, Hunsberger S, Patel J, Wan X, Daly CA. Development of a universal method for vibrational analysis of the terminal alkyne C≡C stretch. J Chem Phys 2024; 160:074106. [PMID: 38364010 DOI: 10.1063/5.0185580] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2023] [Accepted: 01/07/2024] [Indexed: 02/18/2024] Open
Abstract
The terminal alkyne C≡C stretch has a large Raman scattering cross section in the "silent" region for biomolecules. This has led to many Raman tag and probe studies using this moiety to study biomolecular systems. A computational investigation of these systems is vital to aid in the interpretation of these results. In this work, we develop a method for computing terminal alkyne vibrational frequencies and isotropic transition polarizabilities that can easily and accurately be applied to any terminal alkyne molecule. We apply the discrete variable representation method to a localized version of the C≡C stretch normal mode. The errors of (1) vibrational localization to the terminal alkyne moiety, (2) anharmonic normal mode isolation, and (3) discretization of the Born-Oppenheimer potential energy surface are quantified and found to be generally small and cancel each other. This results in a method with low error compared to other anharmonic vibrational methods like second-order vibrational perturbation theory and to experiments. Several density functionals are tested using the method, and TPSS-D3, an inexpensive nonempirical density functional with dispersion corrections, is found to perform surprisingly well. Diffuse basis functions are found to be important for the accuracy of computed frequencies. Finally, the computation of vibrational properties like isotropic transition polarizabilities and the universality of the localized normal mode for terminal alkynes are demonstrated.
Collapse
Affiliation(s)
- Kristina Streu
- Department of Chemistry, Haverford College, 370 Lancaster Ave., Haverford, Pennsylvania 19041, USA
| | - Sara Hunsberger
- Department of Chemistry, Haverford College, 370 Lancaster Ave., Haverford, Pennsylvania 19041, USA
| | - Jeanette Patel
- Department of Chemistry, Haverford College, 370 Lancaster Ave., Haverford, Pennsylvania 19041, USA
| | - Xiang Wan
- Department of Mathematics and Statistics, Loyola University Chicago, 1032 W. Sheridan Rd., Chicago, Illinois 60660, USA
| | - Clyde A Daly
- Department of Chemistry, Haverford College, 370 Lancaster Ave., Haverford, Pennsylvania 19041, USA
| |
Collapse
|
4
|
Oenen K, Dinu DF, Liedl KR. Determining internal coordinate sets for optimal representation of molecular vibration. J Chem Phys 2024; 160:014104. [PMID: 38180253 DOI: 10.1063/5.0180657] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/11/2023] [Accepted: 12/10/2023] [Indexed: 01/06/2024] Open
Abstract
Arising from the harmonic approximation in solving the vibrational Schrödinger equation, normal modes dissect molecular vibrations into distinct degrees of freedom. Normal modes are widely used as they give rise to descriptive vibrational notations and are convenient for expanding anharmonic potential energy surfaces as an alternative to higher-order Taylor series representations. Usually, normal modes are expressed in Cartesian coordinates, which bears drawbacks that can be overcome by switching to internal coordinates. Considering vibrational notations, normal modes with delocalized characters are difficult to denote, but internal coordinates offer a route to clearer notations. Based on the Hessian, normal mode decomposition schemes for a given set of internal coordinates can describe a normal mode by its contributions from internal coordinates. However, choosing a set of internal coordinates is not straightforward. While the Hessian provides unique sets of normal modes, various internal coordinate sets are possible for a given system. In the present work, we employ a normal mode decomposition scheme to choose an optimal set. Therefore, we screen reasonable sets based on topology and symmetry considerations and rely on a metric that minimizes coupling between internal coordinates. Ultimately, the Nomodeco toolkit presented here generates internal coordinate sets to find an optimal set for representing molecular vibrations. The resulting contribution tables can be used to clarify vibrational notations. We test our scheme on small to mid-sized molecules, showing how the space of definable internal coordinate sets can significantly be reduced.
Collapse
Affiliation(s)
- Kemal Oenen
- Department of General, Inorganic and Theoretical Chemistry, University of Innsbruck, Innrain 80, 6020 Innsbruck, Austria
| | - Dennis F Dinu
- Department of General, Inorganic and Theoretical Chemistry, University of Innsbruck, Innrain 80, 6020 Innsbruck, Austria
| | - Klaus R Liedl
- Department of General, Inorganic and Theoretical Chemistry, University of Innsbruck, Innrain 80, 6020 Innsbruck, Austria
| |
Collapse
|
5
|
Glaser N, Baiardi A, Reiher M. Flexible DMRG-Based Framework for Anharmonic Vibrational Calculations. J Chem Theory Comput 2023; 19:9329-9343. [PMID: 38060309 PMCID: PMC10753801 DOI: 10.1021/acs.jctc.3c00902] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/16/2023] [Revised: 10/16/2023] [Accepted: 10/23/2023] [Indexed: 12/08/2023]
Abstract
We present a novel formulation of the vibrational density matrix renormalization group (vDMRG) algorithm tailored to strongly anharmonic molecules described by general, high-dimensional model representations of potential energy surfaces. For this purpose, we extend the vDMRG framework to support vibrational Hamiltonians expressed in the so-called n-mode second-quantization formalism. The resulting n-mode vDMRG method offers full flexibility with respect to both the functional form of the PES and the choice of the single-particle basis set. We leverage this framework to apply, for the first time, vDMRG based on an anharmonic modal basis set optimized with the vibrational self-consistent field algorithm on an on-the-fly constructed PES. We also extend the n-mode vDMRG framework to include excited-state-targeting algorithms in order to efficiently calculate anharmonic transition frequencies. We demonstrate the capabilities of our novel n-mode vDMRG framework for methyloxirane, a challenging molecule with 24 coupled vibrational modes.
Collapse
Affiliation(s)
- Nina Glaser
- Department of Chemistry
and Applied Biosciences, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Alberto Baiardi
- Department of Chemistry
and Applied Biosciences, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Markus Reiher
- Department of Chemistry
and Applied Biosciences, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| |
Collapse
|
6
|
Boyer MA, Sibert EL. A general expression for vibrational Hamiltonians expressed in oblique coordinates. J Chem Phys 2023; 159:234108. [PMID: 38108478 DOI: 10.1063/5.0181135] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/14/2023] [Accepted: 11/29/2023] [Indexed: 12/19/2023] Open
Abstract
We examine the properties of oblique coordinates. The coordinates, introduced by Zúñiga et al. [J. Phys. B: At., Mol. Opt. Phys. 52, 055101, (2019)], reduce vibrational mode-mixing and enhance the quality of vibrational assignments in quantum mechanical investigations of two-dimensional model Hamiltonians. Oblique coordinates are obtained by making non-orthogonal rotations of the original coordinates that convert the matrix representation of the quadratic Hamiltonian operator into a block-diagonal matrix where the blocks are distinguished by the total quanta of vibrational excitation. Using techniques for the polar decomposition of matrices, we present a scheme for finding these coordinates for systems of arbitrary dimensions. Several molecular examples are presented that highlight the advantages of these coordinates.
Collapse
Affiliation(s)
- Mark A Boyer
- Department of Chemistry and Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA
| | - Edwin L Sibert
- Department of Chemistry and Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA
| |
Collapse
|
7
|
Yang EL, Talbot JJ, Spencer RJ, Steele RP. Pitfalls in the n-mode representation of vibrational potentials. J Chem Phys 2023; 159:204104. [PMID: 38010326 DOI: 10.1063/5.0176612] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/15/2023] [Accepted: 10/19/2023] [Indexed: 11/29/2023] Open
Abstract
Simulations of anharmonic vibrational motion rely on computationally expedient representations of the governing potential energy surface. The n-mode representation (n-MR)-effectively a many-body expansion in the space of molecular vibrations-is a general and efficient approach that is often used for this purpose in vibrational self-consistent field (VSCF) calculations and correlated analogues thereof. In the present analysis, a lack of convergence in many VSCF calculations is shown to originate from negative and unbound potentials at truncated orders of the n-MR expansion. For cases of strong anharmonic coupling between modes, the n-MR can both dip below the true global minimum of the potential surface and lead to effective single-mode potentials in VSCF that do not correspond to bound vibrational problems, even for bound total potentials. The present analysis serves mainly as a pathology report of this issue. Furthermore, this insight into the origin of VSCF non-convergence provides a simple, albeit ad hoc, route to correct the problem by "painting in" the full representation of groups of modes that exhibit these negative potentials at little additional computational cost. Somewhat surprisingly, this approach also reasonably approximates the results of the next-higher n-MR order and identifies groups of modes with particularly strong coupling. The method is shown to identify and correct problematic triples of modes-and restore SCF convergence-in two-mode representations of challenging test systems, including the water dimer and trimer, as well as protonated tropine.
Collapse
Affiliation(s)
- Emily L Yang
- Department of Chemistry, The University of Utah, 315 S 1400 E, Salt Lake City, Utah 84112, USA
- Henry Eyring Center for Theoretical Chemistry, The University of Utah, Salt Lake City, Utah 84112, USA
| | - Justin J Talbot
- Department of Chemistry, University of California-Berkeley, 420 Latimer Hall, Berkeley, California 94720, USA
| | - Ryan J Spencer
- Department of Chemistry, The University of Utah, 315 S 1400 E, Salt Lake City, Utah 84112, USA
- Henry Eyring Center for Theoretical Chemistry, The University of Utah, Salt Lake City, Utah 84112, USA
| | - Ryan P Steele
- Department of Chemistry, The University of Utah, 315 S 1400 E, Salt Lake City, Utah 84112, USA
- Henry Eyring Center for Theoretical Chemistry, The University of Utah, Salt Lake City, Utah 84112, USA
| |
Collapse
|
8
|
Zhanserkeev AA, Yang EL, Steele RP. Accelerating Anharmonic Spectroscopy Simulations via Local-Mode, Multilevel Methods. J Chem Theory Comput 2023; 19:5572-5585. [PMID: 37555634 DOI: 10.1021/acs.jctc.3c00589] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/10/2023]
Abstract
Ab initio computer simulations of anharmonic vibrational spectra provide nuanced insight into the vibrational behavior of molecules and complexes. The computational bottleneck in such simulations, particularly for ab initio potentials, is often the generation of mode-coupling potentials. Focusing specifically on two-mode couplings in this analysis, the combination of a local-mode representation and multilevel methods is demonstrated to be particularly symbiotic. In this approach, a low-level quantum chemistry method is employed to predict the pairwise couplings that should be included at the target level of theory in vibrational self-consistent field (and similar) calculations. Pairs that are excluded by this approach are "recycled" at the low level of theory. Furthermore, because this low-level pre-screening will eventually become the computational bottleneck for sufficiently large chemical systems, distance-based truncation is applied to these low-level predictions without substantive loss of accuracy. This combination is demonstrated to yield sub-wavenumber fidelity with reference vibrational transitions when including only a small fraction of target-level couplings; the overhead of predicting these couplings, particularly when employing distance-based, local-mode cutoffs, is a trivial added cost. This combined approach is assessed on a series of test cases, including ethylene, hexatriene, and the alanine dipeptide. Vibrational self-consistent field (VSCF) spectra were obtained with an RI-MP2/cc-pVTZ potential for the dipeptide, at approximately a 5-fold reduction in computational cost. Considerable optimism for increased accelerations for larger systems and higher-order couplings is also justified, based on this investigation.
Collapse
Affiliation(s)
- Asylbek A Zhanserkeev
- Department of Chemistry and Henry Eyring Center for Theoretical Chemistry, University of Utah, 315 South 1400 East, Salt Lake City, Utah 84112, United States
| | - Emily L Yang
- Department of Chemistry and Henry Eyring Center for Theoretical Chemistry, University of Utah, 315 South 1400 East, Salt Lake City, Utah 84112, United States
| | - Ryan P Steele
- Department of Chemistry and Henry Eyring Center for Theoretical Chemistry, University of Utah, 315 South 1400 East, Salt Lake City, Utah 84112, United States
| |
Collapse
|
9
|
Foley CD, Allen CD, Au K, Lee C, Rempe SB, Ren P, Sibert EL, Zwier TS. Molecular Cage Reports on Its Contents: Spectroscopic Signatures of Cryo-Cooled K +- and Ba 2+-Benzocryptand Complexes. J Phys Chem A 2023. [PMID: 37478410 DOI: 10.1021/acs.jpca.3c03457] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/23/2023]
Abstract
UV photofragment spectroscopy and IR-UV double resonance methods are used to determine the structure and spectroscopic responses of a three-dimensional [2.2.2]-benzocryptand cage to the incorporation of a single K+ or Ba2+ imbedded inside it (labeled as K+-BzCrypt, Ba2+-BzCrypt). We studied the isolated ion-cryptand complex under cryo-cooled conditions, brought into the gas phase by nano-electrospray ionization. Incorporation of a phenyl ring in place of the central ethyl group in one of the three N-CH2-CH2-O-CH2-CH2-O-CH2-CH2-N chains provides a UV chromophore whose S0-S1 transition we probe. K+-BzCrypt and Ba2+-BzCrypt have their S0-S1 origin transitions at 35,925 and 36,446 cm-1, respectively, blue-shifted by 174 and 695 cm-1 from that of 1,2-dimethoxybenzene. These origins are used to excite a single conformation of each complex selectively and record their IR spectra using IR-UV dip spectroscopy. The alkyl CH stretch region (2800-3000 cm-1) is surprisingly sensitive to the presence and nature of the encapsulated ion. We carried out an exhaustive conformational search of cage conformations for K+-BzCrypt and Ba2+-BzCrypt, identifying two conformations (A and B) that lie below all others in energy. We extend our local mode anharmonic model of the CH stretch region to these strongly bound ion-cage complexes to predict conformation-specific alkyl CH stretch spectra, obtaining quantitative agreement with experiment for conformer A, the gas-phase global minimum. The large electrostatic effect of the charge on the O- and N-lone pairs affects the local mode frequencies of the CH2 groups adjacent to these atoms. The localized CH2 scissors modes are pushed up in frequency by the adjacent O/N-atoms so that their overtones have little effect on the alkyl CH stretch region. However, the localized CH2 wags are nearly degenerate and strongly coupled to one another, producing an array of delocalized wag normal modes, whose highest frequency members reach up above 1400 cm-1. As such, their overtones mix significantly with the CH stretch modes, most notably involving the CH2 symmetric stretch fundamentals of the central ethyl groups in the all-alkyl chains and the CH stretches adjacent to the N-atoms and antiperiplanar to the nitrogen lone pair.
Collapse
Affiliation(s)
- Casey D Foley
- Gas Phase Chemical Physics, Sandia National Laboratories, Livermore, California 94550, United States
| | - Cole D Allen
- Department of Biomedical Engineering, University of Texas-Austin, Austin, Texas 78712, United States
| | - Kendrew Au
- Gas Phase Chemical Physics, Sandia National Laboratories, Livermore, California 94550, United States
| | - Chin Lee
- Gas Phase Chemical Physics, Sandia National Laboratories, Livermore, California 94550, United States
| | - Susan B Rempe
- Center for Integrated Nanotechnologies, Sandia National Laboratories, Albuquerque, New Mexico 87185, United States
| | - Pengyu Ren
- Department of Biomedical Engineering, University of Texas-Austin, Austin, Texas 78712, United States
| | - Edwin L Sibert
- Department of Chemistry, University of Wisconsin-Madison, Madison, Wisconsin 53706, United States
| | - Timothy S Zwier
- Gas Phase Chemical Physics, Sandia National Laboratories, Livermore, California 94550, United States
| |
Collapse
|
10
|
Artiukhin DG, Godtliebsen IH, Schmitz G, Christiansen O. Gaussian process regression adaptive density-guided approach: Toward calculations of potential energy surfaces for larger molecules. J Chem Phys 2023; 159:024102. [PMID: 37428042 DOI: 10.1063/5.0152367] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/29/2023] [Accepted: 06/15/2023] [Indexed: 07/11/2023] Open
Abstract
We present a new program implementation of the Gaussian process regression adaptive density-guided approach [Schmitz et al., J. Chem. Phys. 153, 064105 (2020)] for automatic and cost-efficient potential energy surface construction in the MidasCpp program. A number of technical and methodological improvements made allowed us to extend this approach toward calculations of larger molecular systems than those previously accessible and maintain the very high accuracy of constructed potential energy surfaces. On the methodological side, improvements were made by using a Δ-learning approach, predicting the difference against a fully harmonic potential, and employing a computationally more efficient hyperparameter optimization procedure. We demonstrate the performance of this method on a test set of molecules of growing size and show that up to 80% of single point calculations could be avoided, introducing a root mean square deviation in fundamental excitations of about 3 cm-1. A much higher accuracy with errors below 1 cm-1 could be achieved with tighter convergence thresholds still reducing the number of single point computations by up to 68%. We further support our findings with a detailed analysis of wall times measured while employing different electronic structure methods. Our results demonstrate that GPR-ADGA is an effective tool, which could be applied for cost-efficient calculations of potential energy surfaces suitable for highly accurate vibrational spectra simulations.
Collapse
Affiliation(s)
- Denis G Artiukhin
- Institut für Chemie und Biochemie, Freie Universität Berlin, Arnimallee 22, 14195 Berlin, Germany
| | - Ian H Godtliebsen
- Department of Chemistry, Aarhus Universitet, DK-8000 Aarhus, Denmark
| | - Gunnar Schmitz
- Lehrstuhl für Theoretische Chemie II, Ruhr-Universität Bochum, Universitätstraße 150, 44801 Bochum, Germany
| | - Ove Christiansen
- Department of Chemistry, Aarhus Universitet, DK-8000 Aarhus, Denmark
| |
Collapse
|
11
|
Brüggemann J, Wolter M, Jacob CR. Quantum-chemical calculation of two-dimensional infrared spectra using localized-mode VSCF/VCI. J Chem Phys 2022; 157:244107. [PMID: 36586972 DOI: 10.1063/5.0135273] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022] Open
Abstract
Computational protocols for the simulation of two-dimensional infrared (2D IR) spectroscopy usually rely on vibrational exciton models which require an empirical parameterization. Here, we present an efficient quantum-chemical protocol for predicting static 2D IR spectra that does not require any empirical parameters. For the calculation of anharmonic vibrational energy levels and transition dipole moments, we employ the localized-mode vibrational self-consistent field (L-VSCF)/vibrational configuration interaction (L-VCI) approach previously established for (linear) anharmonic theoretical vibrational spectroscopy [P. T. Panek and C. R. Jacob, ChemPhysChem 15, 3365-3377 (2014)]. We demonstrate that with an efficient expansion of the potential energy surface using anharmonic one-mode potentials and harmonic two-mode potentials, 2D IR spectra of metal carbonyl complexes and dipeptides can be predicted reliably. We further show how the close connection between L-VCI and vibrational exciton models can be exploited to extract the parameters of such models from those calculations. This provides a novel route to the fully quantum-chemical parameterization of vibrational exciton models for predicting 2D IR spectra.
Collapse
Affiliation(s)
- Julia Brüggemann
- Institute of Physical and Theoretical Chemistry, Technische Universität Braunschweig, Gaußstraße 17, 38106 Braunschweig, Germany
| | - Mario Wolter
- Institute of Physical and Theoretical Chemistry, Technische Universität Braunschweig, Gaußstraße 17, 38106 Braunschweig, Germany
| | - Christoph R Jacob
- Institute of Physical and Theoretical Chemistry, Technische Universität Braunschweig, Gaußstraße 17, 38106 Braunschweig, Germany
| |
Collapse
|
12
|
Li SC, Lin YC, Li YP. Comparative Analysis of Uncoupled Mode Approximations for Molecular Thermochemistry and Kinetics. J Chem Theory Comput 2022; 18:6866-6877. [PMID: 36269729 DOI: 10.1021/acs.jctc.2c00664] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The accurate prediction of thermochemistry and kinetic parameters is an important task for reaction modeling. Unfortunately, the commonly used harmonic oscillator model is often not accurate enough due to the absence of anharmonic effects. In this work, we improve the representation of an anharmonic potential energy surface (PES) using uncoupled mode (UM) approximations, which model the full-dimensional PES as a sum of one-dimensional potentials of each mode. We systematically analyze different PES sampling schemes and coordinate systems for constructing the one-dimensional potentials, and benchmark the performance of UM methods on data sets of molecular thermochemistry and kinetic properties. The results show that the accuracy of the UM approach strongly depends on how the one-dimensional potentials are defined. If one-dimensional potentials are constructed by sampling along normal mode directions (UM-N) or along the directions that minimize intermode coupling (E- and E'-optimized), the accuracies of the predicted properties are not significantly improved compared to the harmonic oscillator model. However, significant improvements can be achieved by sampling the torsional modes separately from the vibrational modes (UM-T and UM-VT). In this work, three types of coordinate systems are examined, including redundant internal coordinates (RIC), hybrid internal coordinates (HIC), and translation-rotation-internal coordinates (TRIC). The HIC and TRIC coordinate systems can outperform RIC since transition state species may contain large-amplitude interfragmentary motions that regular internal coordinates can not describe adequately. Among all the methods we examined, the activation energies and pre-exponential factors calculated using UM-VT with either TRIC or HIC best agree with the reference values. Since UM-VT requires only a number of additional single point energy calculations for each independent mode, the scaling of computational costs of UM-VT is the same as that of the standard harmonic oscillator model, making UM-VT an appealing way of calculating the thermochemistry and kinetic properties for large-size systems.
Collapse
Affiliation(s)
- Shih-Cheng Li
- Department of Chemical Engineering, National Taiwan University, No. 1, Sec. 4, Roosevelt Road, Taipei 10617, Taiwan
| | - Yen-Chun Lin
- Department of Chemical Engineering, National Taiwan University, No. 1, Sec. 4, Roosevelt Road, Taipei 10617, Taiwan
| | - Yi-Pei Li
- Department of Chemical Engineering, National Taiwan University, No. 1, Sec. 4, Roosevelt Road, Taipei 10617, Taiwan.,Taiwan International Graduate Program on Sustainable Chemical Science and Technology (TIGP-SCST), Academia Sinica, No. 128, Sec. 2, Academia Road, Taipei 11529, Taiwan
| |
Collapse
|
13
|
Hanson-Heine MWD. Static Electron Correlation in Anharmonic Molecular Vibrations: A Hybrid TAO-DFT Study. J Phys Chem A 2022; 126:7273-7282. [PMID: 36164938 PMCID: PMC9574917 DOI: 10.1021/acs.jpca.2c05881] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
Hybrid thermally-assisted-occupation density functional theory is used to examine the effects of static electron correlation on the prediction of a benchmark set of experimentally observed molecular vibrational frequencies. The B3LYP and B97-1 thermally-assisted-occupation measure of static electron correlation is important for describing the vibrations of many of the molecules that make up several popular test sets of experimental data. Shifts are seen for known multireference systems and for many molecules containing atoms from the second row of the periodic table of elements. Several molecules only show significant shifts in select vibrational modes, and significant improvements are seen for the prediction of hydrogen stretching frequencies throughout the test set.
Collapse
|
14
|
Sibert EL, Bernath PF. A local mode study of ring puckering effects in the infrared spectra of cyclopentane. J Chem Phys 2022; 156:214305. [DOI: 10.1063/5.0095010] [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/14/2022] Open
Abstract
We report and interpret recently recorded high-resolution infrared spectra for the fundamentals of the CH2 scissors and CH stretches of gas phase cyclopentane at −26.1 and −50 °C, respectively. We extend previous theoretical studies of this molecule, which is known to undergo barrierless pseudorotation due to ring puckering, by constructing local mode Hamiltonians of the stretching and scissor vibrations for which the frequencies, couplings, and linear dipoles are calculated as functions of the pseudorotation angle using B3LYP/6-311++(d,p) and MP2/cc-pVTZ levels of theory. Symmetrization ( D5 h) of the vibrational basis sets leads to simple vibration/pseudorotation Hamiltonians whose solutions lead to good agreement with the experiment at medium resolution, but which miss interesting line fractionation when compared to the high-resolution spectra. In contrast to the scissor motion, pseudorotation leads to significant state mixing of the CH stretches, which themselves are Fermi coupled to the scissor overtones.
Collapse
Affiliation(s)
- Edwin L. Sibert
- Department of Chemistry and Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA
| | - Peter F. Bernath
- Department of Chemistry and Biochemistry, Old Dominion University, Norfolk, Virginia 23529, USA
- Department of Physics, Old Dominion University, Norfolk, Virginia 23529, USA
| |
Collapse
|
15
|
Hino K, Kurashige Y. Matrix Product State Formulation of the MCTDH Theory in Local Mode Representations for Anharmonic Potentials. J Chem Theory Comput 2022; 18:3347-3356. [PMID: 35606892 DOI: 10.1021/acs.jctc.2c00243] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The matrix product state formulation of the multiconfiguration time-dependent Hartree theory, MPS-MCTDH, reported previously [Kurashige, J. Chem. Phys. 2018, 19, 194114] is extended to realistic anharmonic potentials with n-mode representations beyond the linear vibronic coupling model. For realistic vibrational potentials, the local mode representation should give a more compact representation of the potentials, i.e., lowering the dimensionality of the entanglements, than the normal coordinates, and the MPS-MCTDH formulation should work more efficiently and maintain the accuracy with a small bond dimension of the MPS ansatz. In fact, it was confirmed that the use of the local coordinates made the interaction matrices diagonal dominant and the number of terms in the n-body expansion of the potentials was significantly reduced. The method was applied to the IR spectrum of the CH2O molecule, the zero-point energies, and the vibrational energy redistribution dynamics of polyenes C2nH2n+2. The results showed that the efficiency of the MPS-MCTDH method is significantly accelerated by the use of local coordinates even if the long-range interactions are included in the potential.
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
| |
Collapse
|
16
|
McNary CP, Demireva M, Martens J, Berden G, Oomens J, Hamlow LA, Rodgers MT, Armentrout PB. Infrared multiple photon dissociation action spectroscopy of protonated unsymmetrical dimethylhydrazine and proton-bound dimers of hydrazine and unsymmetrical dimethylhydrazine. Phys Chem Chem Phys 2021; 23:25877-25885. [PMID: 34766618 DOI: 10.1039/d1cp03781a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
Abstract
The gas-phase structures of protonated unsymmetrical 1,1-dimethylhydrazine (UDMH) and the proton-bound dimers of UDMH and hydrazine are examined by infrared multiple photon dissociation (IRMPD) action spectroscopy utilizing light generated by a free electron laser and an optical parametric oscillator laser system. To identify the structures present in the experimental studies, the measured IRMPD spectra are compared to spectra calculated at the B3LYP-GD3BJ/6-311+G(d,p) level of theory. These comparisons show that protonated UDMH binds the proton at the methylated nitrogen atom (α) with two low-lying α conformers probably being populated. For (UDMH)2H+, the proton is shared between the methylated nitrogen atoms with several low-lying α conformers likely to be populated. Higher-lying conformers of (UDMH)2H+ in which the proton is shared between α and β (unmethylated) nitrogen atoms cannot be ruled out on the basis of the IRPMD spectrum. For (N2H4)2H+, there are four low-lying conformers that all reproduce the IRMPD spectrum reasonably well. As hydrazine and UDMH see usage as fuels for rocket engines, such spectra are potentially useful as a means of remotely monitoring rocket launches, especially in cases of unsuccessful launches where environmental hazards need to be assessed.
Collapse
Affiliation(s)
| | - Maria Demireva
- Department of Chemistry, University of Utah, Salt Lake City, UT 84112, USA.
| | - Jonathan Martens
- Radboud University, Institute for Molecules and Materials, FELIX Laboratory, Toernooiveld 7, 6525ED Nijmegen, The Netherlands
| | - Giel Berden
- Radboud University, Institute for Molecules and Materials, FELIX Laboratory, Toernooiveld 7, 6525ED Nijmegen, The Netherlands
| | - Jos Oomens
- Radboud University, Institute for Molecules and Materials, FELIX Laboratory, Toernooiveld 7, 6525ED Nijmegen, The Netherlands.,Van't Hoff Institute for Molecular Sciences, University of Amsterdam, Science Park 904, 1098XH Amsterdam, The Netherlands
| | - L A Hamlow
- Department of Chemistry, Wayne State University, Detroit, Michigan 48202, USA
| | - M T Rodgers
- Department of Chemistry, Wayne State University, Detroit, Michigan 48202, USA
| | - P B Armentrout
- Department of Chemistry, University of Utah, Salt Lake City, UT 84112, USA.
| |
Collapse
|
17
|
Yagi K, Sugita Y. Anharmonic Vibrational Calculations Based on Group-Localized Coordinates: Applications to Internal Water Molecules in Bacteriorhodopsin. J Chem Theory Comput 2021; 17:5007-5020. [PMID: 34296615 PMCID: PMC10986902 DOI: 10.1021/acs.jctc.1c00060] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
An efficient anharmonic vibrational method is developed exploiting the locality of molecular vibration. Vibrational coordinates localized to a group of atoms are employed to divide the potential energy surface (PES) of a system into intra- and inter-group contributions. Then, the vibrational Schrödinger equation is solved based on a PES, in which the inter-group coupling is truncated at the harmonic level while accounting for the intra-group anharmonicity. The method is applied to a pentagonal hydrogen bond network (HBN) composed of internal water molecules and charged residues in a membrane protein, bacteriorhodopsin. The PES is calculated by the quantum mechanics/molecular mechanics (QM/MM) calculation at the level of B3LYP-D3/aug-cc-pVDZ. The infrared (IR) spectrum is computed using a set of coordinates localized to each water molecule and amino acid residue by second-order vibrational quasi-degenerate perturbation theory (VQDPT2). Benchmark calculations show that the proposed method yields the N-D/O-D stretching frequencies with an error of 7 cm-1 at the cost reduced by more than five times. In contrast, the harmonic approximation results in a severe error of 150 cm-1. Furthermore, the size of QM regions is carefully assessed to find that the QM regions should include not only the pentagonal HBN itself but also its HB partners. VQDPT2 calculations starting from transient structures obtained by molecular dynamics simulations have shown that the structural sampling has a significant impact on the calculated IR spectrum. The incorporation of anharmonicity, sufficiently large QM regions, and structural samplings are of essential importance to reproduce the experimental IR spectrum. The computational spectrum paves the way for decoding the IR signal of strong HBNs and helps elucidate their functional roles in biomolecules.
Collapse
Affiliation(s)
- Kiyoshi Yagi
- Theoretical
Molecular Science Laboratory, RIKEN Cluster
for Pioneering Research, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
| | - Yuji Sugita
- Theoretical
Molecular Science Laboratory, RIKEN Cluster
for Pioneering Research, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
- Computational
Biophysics Research Team, RIKEN Center for
Computational Science, 7-1-26 Minatojima-Minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
- Laboratory
for Biomolecular Function Simulation, RIKEN
Center for Biosystems Dynamics Research, 1-6-5 Minatojima-Minamimachi,
Chuo-ku, Kobe, Hyogo 650-0047, Japan
| |
Collapse
|
18
|
Hanson-Heine MWD. Reduced Two-Electron Interactions in Anharmonic Molecular Vibrational Calculations Involving Localized Normal Coordinates. J Chem Theory Comput 2021; 17:4383-4391. [PMID: 34087068 DOI: 10.1021/acs.jctc.1c00314] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Spatially localized vibrational normal mode coordinates are shown to reduce the importance of calculating the full set of two-electron terms in the molecular electronic Schrödinger equation. Electron correlation and dispersion interactions become less significant in (E,E)-1,3,5,7-octatetraene vibrational self-consistent field calculations when displacing remote atoms along multiple coordinates. Electron correlation interactions between spatially remote modes are also found to be less important compared to their corresponding uncorrelated interaction terms. Attenuation of the Coulomb operator indicates that the two-electron terms between remote electrons become less important for accurately describing the strongly contributing mode-coupling terms between sets of localized vibrational modes.
Collapse
|
19
|
Bhatty AU, Brorsen KR. An alternative formulation of vibrational heat-bath configuration interaction. Mol Phys 2021. [DOI: 10.1080/00268976.2021.1936250] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
Affiliation(s)
- Abuzar U. Bhatty
- Department of Chemistry, University of Missouri, Columbia, MO, USA
| | - Kurt R. Brorsen
- Department of Chemistry, University of Missouri, Columbia, MO, USA
| |
Collapse
|
20
|
Development of an Electron-Atom Compton Scattering Apparatus Using a Picosecond Pulsed Electron Gun. ATOMS 2021. [DOI: 10.3390/atoms9010019] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022] Open
Abstract
An apparatus has been developed for electron-atom Compton scattering experiments that can employ a pulsed laser and a picosecond pulsed electron beam in a pump-and-probe scheme. The design and technical details of the apparatus are described. Furthermore, experimental results on the Xe atom in its ground state are presented to illustrate the performance of the pulsed electron gun and the detection and spectrometric capabilities for scattered electrons. The scope of future application is also discussed, involving real-time measurement of intramolecular force acting on each constituent atom with different mass numbers, in a transient, evolving system during a molecular reaction.
Collapse
|
21
|
Mathea T, Petrenko T, Rauhut G. VCI Calculations Based on Canonical and Localized Normal Coordinates for Non-Abelian Molecules: Accurate Assignment of the Vibrational Overtones of Allene. J Phys Chem A 2021; 125:990-998. [DOI: 10.1021/acs.jpca.0c10429] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Tina Mathea
- Institute for Theoretical Chemistry, University of Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Taras Petrenko
- Institute for Theoretical Chemistry, University of Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Guntram Rauhut
- Institute for Theoretical Chemistry, University of Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| |
Collapse
|
22
|
Gloaguen E, Mons M, Schwing K, Gerhards M. Neutral Peptides in the Gas Phase: Conformation and Aggregation Issues. Chem Rev 2020; 120:12490-12562. [PMID: 33152238 DOI: 10.1021/acs.chemrev.0c00168] [Citation(s) in RCA: 38] [Impact Index Per Article: 9.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
Combined IR and UV laser spectroscopic techniques in molecular beams merged with theoretical approaches have proven to be an ideal tool to elucidate intrinsic structural properties on a molecular level. It offers the possibility to analyze structural changes, in a controlled molecular environment, when successively adding aggregation partners. By this, it further makes these techniques a valuable starting point for a bottom-up approach in understanding the forces shaping larger molecular systems. This bottom-up approach was successfully applied to neutral amino acids starting around the 1990s. Ever since, experimental and theoretical methods developed further, and investigations could be extended to larger peptide systems. Against this background, the review gives an introduction to secondary structures and experimental methods as well as a summary on theoretical approaches. Vibrational frequencies being characteristic probes of molecular structure and interactions are especially addressed. Archetypal biologically relevant secondary structures investigated by molecular beam spectroscopy are described, and the influences of specific peptide residues on conformational preferences as well as the competition between secondary structures are discussed. Important influences like microsolvation or aggregation behavior are presented. Beyond the linear α-peptides, the main results of structural analysis on cyclic systems as well as on β- and γ-peptides are summarized. Overall, this contribution addresses current aspects of molecular beam spectroscopy on peptides and related species and provides molecular level insights into manifold issues of chemical and biochemical relevance.
Collapse
Affiliation(s)
- Eric Gloaguen
- CEA, CNRS, Université Paris-Saclay, CEA Paris-Saclay, Bât 522, 91191 Gif-sur-Yvette, France
| | - Michel Mons
- CEA, CNRS, Université Paris-Saclay, CEA Paris-Saclay, Bât 522, 91191 Gif-sur-Yvette, France
| | - Kirsten Schwing
- TU Kaiserslautern & Research Center Optimas, Erwin-Schrödinger-Straße 52, D-67663 Kaiserslautern, Germany
| | - Markus Gerhards
- TU Kaiserslautern & Research Center Optimas, Erwin-Schrödinger-Straße 52, D-67663 Kaiserslautern, Germany
| |
Collapse
|
23
|
Aieta C, Micciarelli M, Bertaina G, Ceotto M. Anharmonic quantum nuclear densities from full dimensional vibrational eigenfunctions with application to protonated glycine. Nat Commun 2020; 11:4348. [PMID: 32859910 PMCID: PMC7455743 DOI: 10.1038/s41467-020-18211-3] [Citation(s) in RCA: 23] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/23/2020] [Accepted: 07/29/2020] [Indexed: 11/24/2022] Open
Abstract
The interpretation of molecular vibrational spectroscopic signals in terms of atomic motion is essential to understand molecular mechanisms and for chemical characterization. The signals are usually assigned after harmonic normal mode analysis, even if molecular vibrations are known to be anharmonic. Here we obtain the quantum anharmonic vibrational eigenfunctions of the 11-atom protonated glycine molecule and we calculate the density distribution of its nuclei and its geometry parameters, for both the ground and the O-H stretch excited states, using our semiclassical method based on ab initio molecular dynamics trajectories. Our quantum mechanical results describe a molecule elongated and more flexible with respect to what previously thought. More importantly, our method is able to assign each spectral peak in vibrational spectroscopy by showing quantitatively how normal modes involving different functional groups cooperate to originate that spectroscopic signal. The method will possibly allow for a better rationalization of experimental spectroscopy.
Collapse
Affiliation(s)
- Chiara Aieta
- Dipartimento di Chimica, Università degli Studi di Milano, via C. Golgi 19, 20133, Milano, Italy
| | - Marco Micciarelli
- Dipartimento di Chimica, Università degli Studi di Milano, via C. Golgi 19, 20133, Milano, Italy
| | - Gianluca Bertaina
- Dipartimento di Chimica, Università degli Studi di Milano, via C. Golgi 19, 20133, Milano, Italy
- Istituto Nazionale di Ricerca Metrologica, Strada delle Cacce 91, 10135, Torino, Italy
| | - Michele Ceotto
- Dipartimento di Chimica, Università degli Studi di Milano, via C. Golgi 19, 20133, Milano, Italy.
| |
Collapse
|
24
|
Baiz CR, Błasiak B, Bredenbeck J, Cho M, Choi JH, Corcelli SA, Dijkstra AG, Feng CJ, Garrett-Roe S, Ge NH, Hanson-Heine MWD, Hirst JD, Jansen TLC, Kwac K, Kubarych KJ, Londergan CH, Maekawa H, Reppert M, Saito S, Roy S, Skinner JL, Stock G, Straub JE, Thielges MC, Tominaga K, Tokmakoff A, Torii H, Wang L, Webb LJ, Zanni MT. Vibrational Spectroscopic Map, Vibrational Spectroscopy, and Intermolecular Interaction. Chem Rev 2020; 120:7152-7218. [PMID: 32598850 PMCID: PMC7710120 DOI: 10.1021/acs.chemrev.9b00813] [Citation(s) in RCA: 165] [Impact Index Per Article: 41.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
Vibrational spectroscopy is an essential tool in chemical analyses, biological assays, and studies of functional materials. Over the past decade, various coherent nonlinear vibrational spectroscopic techniques have been developed and enabled researchers to study time-correlations of the fluctuating frequencies that are directly related to solute-solvent dynamics, dynamical changes in molecular conformations and local electrostatic environments, chemical and biochemical reactions, protein structural dynamics and functions, characteristic processes of functional materials, and so on. In order to gain incisive and quantitative information on the local electrostatic environment, molecular conformation, protein structure and interprotein contacts, ligand binding kinetics, and electric and optical properties of functional materials, a variety of vibrational probes have been developed and site-specifically incorporated into molecular, biological, and material systems for time-resolved vibrational spectroscopic investigation. However, still, an all-encompassing theory that describes the vibrational solvatochromism, electrochromism, and dynamic fluctuation of vibrational frequencies has not been completely established mainly due to the intrinsic complexity of intermolecular interactions in condensed phases. In particular, the amount of data obtained from the linear and nonlinear vibrational spectroscopic experiments has been rapidly increasing, but the lack of a quantitative method to interpret these measurements has been one major obstacle in broadening the applications of these methods. Among various theoretical models, one of the most successful approaches is a semiempirical model generally referred to as the vibrational spectroscopic map that is based on a rigorous theory of intermolecular interactions. Recently, genetic algorithm, neural network, and machine learning approaches have been applied to the development of vibrational solvatochromism theory. In this review, we provide comprehensive descriptions of the theoretical foundation and various examples showing its extraordinary successes in the interpretations of experimental observations. In addition, a brief introduction to a newly created repository Web site (http://frequencymap.org) for vibrational spectroscopic maps is presented. We anticipate that a combination of the vibrational frequency map approach and state-of-the-art multidimensional vibrational spectroscopy will be one of the most fruitful ways to study the structure and dynamics of chemical, biological, and functional molecular systems in the future.
Collapse
Affiliation(s)
- Carlos R. Baiz
- Department of Chemistry, University of Texas at Austin, Austin, TX 78712, U.S.A
| | - Bartosz Błasiak
- Department of Physical and Quantum Chemistry, Wroclaw University of Science and Technology, Wybrzeże Wyspiańskiego 27, 50-370 Wrocław, Poland
| | - Jens Bredenbeck
- Johann Wolfgang Goethe-University, Institute of Biophysics, Max-von-Laue-Str. 1, 60438, Frankfurt am Main, Germany
| | - Minhaeng Cho
- Center for Molecular Spectroscopy and Dynamics, Seoul 02841, Republic of Korea
- Department of Chemistry, Korea University, Seoul 02841, Republic of Korea
| | - Jun-Ho Choi
- Department of Chemistry, Gwangju Institute of Science and Technology, Gwangju 61005, Republic of Korea
| | - Steven A. Corcelli
- Department of Chemistry and Biochemistry, University of Notre Dame, Notre Dame, IN 46556, U.S.A
| | - Arend G. Dijkstra
- School of Chemistry and School of Physics and Astronomy, University of Leeds, Leeds LS2 9JT, U.K
| | - Chi-Jui Feng
- Department of Chemistry, James Franck Institute and Institute for Biophysical Dynamics, University of Chicago, Chicago, IL 60637, U.S.A
| | - Sean Garrett-Roe
- Department of Chemistry, University of Pittsburgh, Pittsburgh, PA 15260, U.S.A
| | - Nien-Hui Ge
- Department of Chemistry, University of California at Irvine, Irvine, CA 92697-2025, U.S.A
| | - Magnus W. D. Hanson-Heine
- School of Chemistry, University of Nottingham, Nottingham, University Park, Nottingham, NG7 2RD, U.K
| | - Jonathan D. Hirst
- School of Chemistry, University of Nottingham, Nottingham, University Park, Nottingham, NG7 2RD, U.K
| | - Thomas L. C. Jansen
- University of Groningen, Zernike Institute for Advanced Materials, Nijenborgh 4, 9747 AG Groningen, The Netherlands
| | - Kijeong Kwac
- Center for Molecular Spectroscopy and Dynamics, Seoul 02841, Republic of Korea
| | - Kevin J. Kubarych
- Department of Chemistry, University of Michigan, 930 N. University Ave., Ann Arbor, MI 48109, U.S.A
| | - Casey H. Londergan
- Department of Chemistry, Haverford College, Haverford, Pennsylvania 19041, U.S.A
| | - Hiroaki Maekawa
- Department of Chemistry, University of California at Irvine, Irvine, CA 92697-2025, U.S.A
| | - Mike Reppert
- Chemical Physics Theory Group, Department of Chemistry, University of Toronto, Toronto, Ontario M5S 3H6, Canada
| | - Shinji Saito
- Department of Theoretical and Computational Molecular Science, Institute for Molecular Science, Myodaiji, Okazaki, 444-8585, Japan
| | - Santanu Roy
- Chemical Sciences Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831-6110, U.S.A
| | - James L. Skinner
- Institute for Molecular Engineering, University of Chicago, Chicago, IL 60637, U.S.A
| | - Gerhard Stock
- Biomolecular Dynamics, Institute of Physics, Albert Ludwigs University, 79104 Freiburg, Germany
| | - John E. Straub
- Department of Chemistry, Boston University, Boston, MA 02215, U.S.A
| | - Megan C. Thielges
- Department of Chemistry, Indiana University, 800 East Kirkwood, Bloomington, Indiana 47405, U.S.A
| | - Keisuke Tominaga
- Molecular Photoscience Research Center, Kobe University, Nada, Kobe 657-0013, Japan
| | - Andrei Tokmakoff
- Department of Chemistry, James Franck Institute and Institute for Biophysical Dynamics, University of Chicago, Chicago, IL 60637, U.S.A
| | - Hajime Torii
- Department of Applied Chemistry and Biochemical Engineering, Faculty of Engineering, and Department of Optoelectronics and Nanostructure Science, Graduate School of Science and Technology, Shizuoka University, 3-5-1 Johoku, Naka-Ku, Hamamatsu 432-8561, Japan
| | - Lu Wang
- Department of Chemistry and Chemical Biology, Institute for Quantitative Biomedicine, Rutgers University, 174 Frelinghuysen Road, Piscataway, NJ 08854, U.S.A
| | - Lauren J. Webb
- Department of Chemistry, The University of Texas at Austin, 105 E. 24th Street, STOP A5300, Austin, Texas 78712, U.S.A
| | - Martin T. Zanni
- Department of Chemistry, University of Wisconsin-Madison, Madison, Wisconsin 53706-1396, U.S.A
| |
Collapse
|
25
|
Riera M, Talbot JJ, Steele RP, Paesani F. Infrared signatures of isomer selectivity and symmetry breaking in the Cs+(H2O)3 complex using many-body potential energy functions. J Chem Phys 2020; 153:044306. [DOI: 10.1063/5.0013101] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023] Open
Affiliation(s)
- Marc Riera
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, USA
| | - Justin J. Talbot
- Department of Chemistry, University of Utah, Salt Lake City, Utah 84112, USA
| | - Ryan P. Steele
- Department of Chemistry, University of Utah, Salt Lake City, Utah 84112, USA
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, USA
- Materials Science and Engineering, University of California San Diego, La Jolla, California 92093, USA
- San Diego Supercomputer Center, University of California San Diego, La Jolla, California 92093, USA
| |
Collapse
|
26
|
Mathea T, Rauhut G. Assignment of vibrational states within configuration interaction calculations. J Chem Phys 2020; 152:194112. [DOI: 10.1063/5.0009732] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Tina Mathea
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Guntram Rauhut
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| |
Collapse
|
27
|
Artiukhin DG, Klinting EL, König C, Christiansen O. Adaptive density-guided approach to double incremental potential energy surface construction. J Chem Phys 2020; 152:194105. [PMID: 33687258 DOI: 10.1063/5.0004686] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022] Open
Abstract
We present a combination of the recently developed double incremental expansion of potential energy surfaces with the well-established adaptive density-guided approach to grid construction. This unique methodology is based on the use of an incremental expansion for potential energy surfaces, known as n-mode expansion; an incremental many-body representation of the electronic energy; and an efficient vibrational density-guided approach to automated determination of grid dimensions and granularity. The reliability of the method is validated calculating potential energy surfaces and obtaining fundamental excitation energies for three moderate-size chain-like molecular systems. The use of our methodology leads to considerable computational savings for potential energy surface construction compared to standard approaches while maintaining a high level of accuracy in the resulting potential energy surfaces. Additional investigations indicate that our method can be applied to covalently bound and strongly interacting molecular systems, even though these cases are known to be very unfavorable for fragmentation schemes. We therefore conclude that the presented methodology is a robust and flexible approach to potential energy surface construction, which introduces considerable computational savings without compromising the accuracy of vibrational spectra calculations.
Collapse
Affiliation(s)
- Denis G Artiukhin
- Department of Chemistry, Aarhus Universitet, DK-8000 Aarhus, Denmark
| | | | - Carolin König
- Institute of Physical Chemistry, Christian-Albrechts-University Kiel, Max-Eyth-Straße 1, D-24118 Kiel, Germany
| | - Ove Christiansen
- Department of Chemistry, Aarhus Universitet, DK-8000 Aarhus, Denmark
| |
Collapse
|
28
|
Klinting EL, Lauvergnat D, Christiansen O. Vibrational Coupled Cluster Computations in Polyspherical Coordinates with the Exact Analytical Kinetic Energy Operator. J Chem Theory Comput 2020; 16:4505-4520. [DOI: 10.1021/acs.jctc.0c00261] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Affiliation(s)
| | - David Lauvergnat
- Université Paris-Saclay, CNRS, Institut de Chimie Physique, UMR8000, 91405 Orsay, France
| | - Ove Christiansen
- Department of Chemistry, Aarhus University, DK-8000 Aarhus C, Denmark
| |
Collapse
|
29
|
Kraka E, Zou W, Tao Y. Decoding chemical information from vibrational spectroscopy data: Local vibrational mode theory. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2020. [DOI: 10.1002/wcms.1480] [Citation(s) in RCA: 43] [Impact Index Per Article: 10.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Affiliation(s)
- Elfi Kraka
- Department of Chemistry Southern Methodist University Dallas Texas USA
| | - Wenli Zou
- Institute of Modern Physics Northwest University and Shaanxi Key Laboratory for Theoretical Physics Frontiers, Xi'an Shaanxi PR China
| | - Yunwen Tao
- Department of Chemistry Southern Methodist University Dallas Texas USA
| |
Collapse
|
30
|
Klinting EL, Christiansen O, König C. Toward Accurate Theoretical Vibrational Spectra: A Case Study for Maleimide. J Phys Chem A 2020; 124:2616-2627. [DOI: 10.1021/acs.jpca.9b11915] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Affiliation(s)
| | - Ove Christiansen
- Department of Chemistry, Aarhus University, DK-8000 Aarhus C, Denmark
| | - Carolin König
- Institute of Physical Chemistry, Kiel University, Max-Eyth-Straße 1, D-24118 Kiel, Germany
| |
Collapse
|
31
|
Puzzarini C, Bloino J, Tasinato N, Barone V. Accuracy and Interpretability: The Devil and the Holy Grail. New Routes across Old Boundaries in Computational Spectroscopy. Chem Rev 2019; 119:8131-8191. [DOI: 10.1021/acs.chemrev.9b00007] [Citation(s) in RCA: 114] [Impact Index Per Article: 22.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Affiliation(s)
- Cristina Puzzarini
- Dipartimento di Chimica “Giacomo Ciamician”, Università di Bologna, Via F. Selmi 2, I-40126 Bologna, Italy
| | - Julien Bloino
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| | - Nicola Tasinato
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| | - Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| |
Collapse
|
32
|
Bajaj P, Zhuang D, Paesani F. Specific Ion Effects on Hydrogen-Bond Rearrangements in the Halide-Dihydrate Complexes. J Phys Chem Lett 2019; 10:2823-2828. [PMID: 31082245 DOI: 10.1021/acs.jpclett.9b00899] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
Small aqueous ionic clusters represent ideal systems to investigate the microscopic hydrogen-bonding structure and dynamics in ion hydration shells. In this context, halide-dihydrate complexes are the smallest systems where the interplay between halide-water and water-water interactions can be studied simultaneously. Here, quantum molecular dynamics simulations unravel specific ion effects on the temperature-dependent structural transition in X-(H2O)2 complexes (X = Cl, Br, and I), which is induced by the breaking of the water-water hydrogen bond. A systematic analysis of the hydrogen-bonding rearrangements at low temperature provides fundamental insights into the competition between halide-water and water-water interactions depending on the properties of the halide ion. While the halide-water hydrogen-bond strength decreases going from Cl-(H2O)2 to I-(H2O)2, the opposite trend in observed in the strength of the water-water hydrogen-bond, suggesting that nontrivial many-body effects may also be at play in the hydration shells of halide ions in solution, especially in frustrated systems (e.g., interfaces) where the water molecules can have dangling OH bonds.
Collapse
Affiliation(s)
- Pushp Bajaj
- Department of Chemistry and Biochemistry , University of California, San Diego , La Jolla , California 92093 , United States
| | - Debbie Zhuang
- Department of Chemistry and Biochemistry , University of California, San Diego , La Jolla , California 92093 , United States
| | - Francesco Paesani
- Department of Chemistry and Biochemistry , University of California, San Diego , La Jolla , California 92093 , United States
- Materials Science and Engineering , University of California, San Diego , La Jolla , California 92093 , United States
- San Diego Supercomputer Center , University of California, San Diego , La Jolla , California 92093 , United States
| |
Collapse
|
33
|
Ziegler B, Rauhut G. Localized Normal Coordinates in Accurate Vibrational Structure Calculations: Benchmarks for Small Molecules. J Chem Theory Comput 2019; 15:4187-4196. [DOI: 10.1021/acs.jctc.9b00381] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Benjamin Ziegler
- Institute for Theoretical Chemistry, University of Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Guntram Rauhut
- Institute for Theoretical Chemistry, University of Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| |
Collapse
|
34
|
Hanson-Heine MWD. Excited-State Vibrational Frequencies: Restricted Virtual Space Time-Dependent Density Functional Theory. J Phys Chem A 2019; 123:2949-2956. [DOI: 10.1021/acs.jpca.8b09642] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
35
|
Bajaj P, Riera M, Lin JK, Mendoza Montijo YE, Gazca J, Paesani F. Halide Ion Microhydration: Structure, Energetics, and Spectroscopy of Small Halide–Water Clusters. J Phys Chem A 2019; 123:2843-2852. [DOI: 10.1021/acs.jpca.9b00816] [Citation(s) in RCA: 25] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
|
36
|
Baiardi A, Stein CJ, Barone V, Reiher M. Optimization of highly excited matrix product states with an application to vibrational spectroscopy. J Chem Phys 2019; 150:094113. [DOI: 10.1063/1.5068747] [Citation(s) in RCA: 23] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Alberto Baiardi
- ETH Zürich, Laboratorium für Physikalische Chemie, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Christopher J. Stein
- ETH Zürich, Laboratorium für Physikalische Chemie, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Markus Reiher
- ETH Zürich, Laboratorium für Physikalische Chemie, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| |
Collapse
|
37
|
Sibert EL. Modeling vibrational anharmonicity in infrared spectra of high frequency vibrations of polyatomic molecules. J Chem Phys 2019; 150:090901. [DOI: 10.1063/1.5079626] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Edwin L. Sibert
- Department of Chemistry and Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA
| |
Collapse
|
38
|
Bajaj P, Richardson JO, Paesani F. Ion-mediated hydrogen-bond rearrangement through tunnelling in the iodide–dihydrate complex. Nat Chem 2019; 11:367-374. [DOI: 10.1038/s41557-019-0220-2] [Citation(s) in RCA: 34] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2018] [Accepted: 01/22/2019] [Indexed: 11/09/2022]
|
39
|
Panek PT, Hoeske AA, Jacob CR. On the choice of coordinates in anharmonic theoretical vibrational spectroscopy: Harmonic vs. anharmonic coupling in vibrational configuration interaction. J Chem Phys 2019; 150:054107. [PMID: 30736699 DOI: 10.1063/1.5083186] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023] Open
Abstract
By a suitable choice of coordinates, the computational effort required for calculations of anharmonic vibrational spectra can be reduced significantly. By using suitable localized-mode coordinates obtained from an orthogonal transformation of the conventionally used normal-mode coordinates, anharmonic couplings between modes can be significantly reduced. However, such a transformation introduces harmonic couplings between the localized modes. To elucidate the role of these harmonic couplings, we consider the vibrational self-consistent field (VSCF)/vibrational configuration interaction (VCI) calculations for both few-mode model systems and for ethene as a molecular test case. We show that large harmonic couplings can result in significant errors in localized-mode L-VSCF/L-VCI calculations and study the convergence with respect to the size of the VCI excitation space. To further elucidate the errors introduced by harmonic couplings, we discuss the connection between L-VSCF/L-VCI and vibrational exciton models. With the help of our results, we propose an algorithm for the localization of normal modes in suitable subsets that are chosen to strictly limit the errors introduced by the harmonic couplings while still leading to maximally localized modes.
Collapse
Affiliation(s)
- Paweł T Panek
- Technische Universität Braunschweig, Institute of Physical and Theoretical Chemistry, Gaußstraße 17, 38106 Braunschweig, Germany
| | - Adrian A Hoeske
- Technische Universität Braunschweig, Institute of Physical and Theoretical Chemistry, Gaußstraße 17, 38106 Braunschweig, Germany
| | - Christoph R Jacob
- Technische Universität Braunschweig, Institute of Physical and Theoretical Chemistry, Gaußstraße 17, 38106 Braunschweig, Germany
| |
Collapse
|
40
|
Yagi K, Yamada K, Kobayashi C, Sugita Y. Anharmonic Vibrational Analysis of Biomolecules and Solvated Molecules Using Hybrid QM/MM Computations. J Chem Theory Comput 2019; 15:1924-1938. [PMID: 30730746 PMCID: PMC8864611 DOI: 10.1021/acs.jctc.8b01193] [Citation(s) in RCA: 22] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
![]()
Quantum
mechanics/molecular mechanics (QM/MM) calculations are
applied for anharmonic vibrational analyses of biomolecules and solvated
molecules. The QM/MM method is implemented into a molecular dynamics
(MD) program, GENESIS, by interfacing with external electronic structure
programs. Following the geometry optimization and the harmonic normal-mode
analysis based on a partial Hessian, the anharmonic potential energy
surface (PES) is generated from QM/MM energies and gradients calculated
at grid points. The PES is used for vibrational self-consistent field
(VSCF) and post-VSCF calculations to compute the vibrational spectrum.
The method is first applied to a phosphate ion in solution. With both
the ion and neighboring water molecules taken as a QM region, IR spectra
of representative hydration structures are calculated by the second-order
vibrational quasi-degenerate perturbation theory (VQDPT2) at the level
of B3LYP/cc-pVTZ and TIP3P force field. A weight-average of IR spectra
over the structures reproduces the experimental spectrum with a mean
absolute deviation of 16 cm–1. Then, the method
is applied to an enzyme, P450 nitric oxide reductase (P450nor), with
the NO molecule bound to a ferric (FeIII) heme. Starting
from snapshot structures obtained from MD simulations of P450nor in
solution, QM/MM calculations have been carried out at the level of
B3LYP-D3/def2-SVP(D). The spin state of FeIII(NO) is likely
a closed-shell singlet state based on a ratio of N–O and Fe–NO
stretching frequencies (νN–O and νFe–NO) calculated for closed- and open-shell singlet
states. The calculated νN–O and νFe–NO overestimate the experimental ones by 120 and
75 cm–1, respectively. The electronic structure
and solvation of FeIII(NO) affect the structure around
the heme of P450nor leading to an increase in νN–O and νFe–NO.
Collapse
Affiliation(s)
- Kiyoshi Yagi
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
| | - Kenta Yamada
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
| | - Chigusa Kobayashi
- Computational Biophysics Research Team, RIKEN Center for Computational Science, 7-1-26 Minatojima-Minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
| | - Yuji Sugita
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
- Computational Biophysics Research Team, RIKEN Center for Computational Science, 7-1-26 Minatojima-Minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
- Laboratory for Biomolecular Function Simulation, RIKEN Center for Biosystems Dynamics Research, 1-6-5 Minatojima-Minamimachi,
Chuo-ku, Kobe, Hyogo 650-0047, Japan
| |
Collapse
|
41
|
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
|
42
|
Hunter KM, Shakib FA, Paesani F. Disentangling Coupling Effects in the Infrared Spectra of Liquid Water. J Phys Chem B 2018; 122:10754-10761. [DOI: 10.1021/acs.jpcb.8b09910] [Citation(s) in RCA: 31] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023]
Affiliation(s)
- Kelly M. Hunter
- Department of Chemistry and Biochemistry, University of California, San Diego La Jolla, California 92093, United States
| | - Farnaz A. Shakib
- Department of Chemistry and Biochemistry, University of California, San Diego La Jolla, California 92093, United States
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, University of California, San Diego La Jolla, California 92093, United States
| |
Collapse
|
43
|
Ziegler B, Rauhut G. Rigorous use of symmetry within the construction of multidimensional potential energy surfaces. J Chem Phys 2018; 149:164110. [DOI: 10.1063/1.5047912] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023] Open
Affiliation(s)
- Benjamin Ziegler
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Guntram Rauhut
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| |
Collapse
|
44
|
Kuenzer U, Klotz M, Hofer TS. Probing vibrational coupling via a grid-based quantum approach-an efficient strategy for accurate calculations of localized normal modes in solid-state systems. J Comput Chem 2018; 39:2196-2209. [PMID: 30341952 PMCID: PMC6767160 DOI: 10.1002/jcc.25533] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/16/2018] [Revised: 06/08/2018] [Accepted: 06/29/2018] [Indexed: 11/19/2022]
Abstract
In this work an approach to investigate the properties of strongly localized vibrational modes of functional groups in bulk material and on solid‐state surfaces is presented. The associated normal mode vectors are approximated solely on the basis of structural information and obtained via diagonalization of a reduced Hessian. The grid‐based Numerov procedure in one and two dimensions is then applied to an adequate scan of the respective potential surface yielding the associated vibrational wave functions and energy eigenvalues. This not only provides a detailed description of anharmonic effects but also an accurate inclusion of the coupling between the investigated vibrational states on a quantum mechanical level. All results obtained for the constructed normal modes are benchmarked against their analytical counterparts obtained from the diagonalization of the total Hessian of the entire system. Three increasingly complex systems treated at quantum chemical level of theory have been considered, namely the symmetric and asymmetric stretch vibrations of an isolated water molecule, hydroxyl groups bound to the surface of GeO2 (001), α‐quartz(001) and Rutil (001) as well as crystalline Li2NH serving as an example for a bulk material. While the data obtained for the individual systems verify the applicability of the proposed methodology, comparison to experimental data demonstrates the accuracy of this methodology despite the restriction to limit this methodology to a few selected vibrational modes. The possibility to investigate vibrational phenomena of localized normal modes without the requirement of executing costly harmonic frequency calculations of the entire system enables the application of this method to cases in which the determination of normal modes is prohibitively expensive or not available for a particular level of theory. © 2018 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Ulrich Kuenzer
- University of Innsbruck, Institute of General, Inorganic and Theoretical Chemistry, Innrain 80-82, 6020 Innsbruck, Austria
| | - Martin Klotz
- University of Innsbruck, Institute of General, Inorganic and Theoretical Chemistry, Innrain 80-82, 6020 Innsbruck, Austria
| | - Thomas S Hofer
- University of Innsbruck, Institute of General, Inorganic and Theoretical Chemistry, Innrain 80-82, 6020 Innsbruck, Austria
| |
Collapse
|
45
|
Hanson-Heine MW, Calladine JA, Yang J, Towrie M, Horvath R, Besley NA, George MW. A combined time-resolved infrared and density functional theory study of the lowest excited states of 9-fluorenone and 2-naphthaldehyde. Chem Phys 2018. [DOI: 10.1016/j.chemphys.2018.04.014] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
46
|
Riera M, Brown SE, Paesani F. Isomeric Equilibria, Nuclear Quantum Effects, and Vibrational Spectra of M+(H2O)n=1–3 Clusters, with M = Li, Na, K, Rb, and Cs, through Many-Body Representations. J Phys Chem A 2018; 122:5811-5821. [DOI: 10.1021/acs.jpca.8b04106] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
|
47
|
Avilés-Moreno JR, Berden G, Oomens J, Martínez-Haya B. Intra-cavity proton bonding and anharmonicity in the anionophore cyclen. Phys Chem Chem Phys 2018; 20:8968-8975. [PMID: 29557457 DOI: 10.1039/c8cp00660a] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
Abstract
Proton bonding drives the supramolecular chemistry of a broad range of materials with polar moieties. Proton delocalization and electronic charge redistribution have a profound impact on the structure of proton-bound molecular frameworks, and pose fundamental challenges to quantum chemical modelling. This study provides insights into the structural and spectral signatures of the intramolecular proton bond formed in a benchmark polyazamacrocycle anionophore (cyclen, 1,4,7,10-tetraazacyclododecane). Infrared action spectroscopy is employed to characterize the macrocycle, isolated in protonated form. In its most stable configuration, protonated cyclen adopts an open arrangement of Cs symmetry with a particularly strong NHδ+N bond across the cavity. The quantum chemical analysis of the infrared spectrum reveals intrinsic difficulties for the accurate description of the vibrational modes of the system. The reconciliation of the computational predictions with experiment demands a careful anharmonic treatment of the proton motion, which exposes the limitations of current methods. Best results are obtained with the incorporation of anharmonicity only to the fundamental modes directly related to motions of the proton. However, the full anharmonic treatment of the system fails to describe correctly the vibrations related to the macrocycle backbone. The results should serve as motivation for new developments in the modelling of proton bonded systems.
Collapse
Affiliation(s)
- Juan Ramón Avilés-Moreno
- Department of Physical, Chemical and Natural Systems, Universidad Pablo de Olavide, E-41013 Seville, Spain.
| | - Giel Berden
- Radboud University, Institute for Molecules and Materials, FELIX Laboratory, Toernooiveld 7c, 6525ED Nijmegen, The Netherlands
| | - Jos Oomens
- Radboud University, Institute for Molecules and Materials, FELIX Laboratory, Toernooiveld 7c, 6525ED Nijmegen, The Netherlands
| | - Bruno Martínez-Haya
- Department of Physical, Chemical and Natural Systems, Universidad Pablo de Olavide, E-41013 Seville, Spain.
| |
Collapse
|
48
|
Bajaj P, Wang XG, Carrington T, Paesani F. Vibrational spectra of halide-water dimers: Insights on ion hydration from full-dimensional quantum calculations on many-body potential energy surfaces. J Chem Phys 2018; 148:102321. [DOI: 10.1063/1.5005540] [Citation(s) in RCA: 34] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022] Open
Affiliation(s)
- Pushp Bajaj
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, USA
| | - Xiao-Gang Wang
- Chemistry Department, Queen’s University, Kingston, Ontario K7L3N6, Canada
| | - Tucker Carrington
- Chemistry Department, Queen’s University, Kingston, Ontario K7L3N6, Canada
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, Materials Science and Engineering, and San Diego Supercomputer Center, University of California, San Diego, La Jolla, California 92093, USA
| |
Collapse
|
49
|
Klinting EL, Thomsen B, Godtliebsen IH, Christiansen O. Employing general fit-bases for construction of potential energy surfaces with an adaptive density-guided approach. J Chem Phys 2018; 148:064113. [DOI: 10.1063/1.5016259] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
| | - Bo Thomsen
- Department of Chemistry, Aarhus University, DK-8000 Aarhus C, Denmark
| | | | - Ove Christiansen
- Department of Chemistry, Aarhus University, DK-8000 Aarhus C, Denmark
| |
Collapse
|
50
|
Hanson-Heine MWD. Reduced Basis Set Dependence in Anharmonic Frequency Calculations Involving Localized Coordinates. J Chem Theory Comput 2018; 14:1277-1285. [PMID: 29385338 DOI: 10.1021/acs.jctc.7b01075] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Localized normal coordinates are known to be effective in speeding up anharmonic frequency calculations by reducing the complexity of the nuclear Hamiltonian and wave function. Displacing atoms in localized coordinates can also cause relatively small changes in the electronic structure, which can be exploited for further computational efficiency improvements during ab initio electronic structure calculations of the potential energy surface by reducing the electronic basis set dependence. Three different schemes for reducing the basis set dependence have been investigated in this work. These include combining localized coordinate schemes with general mixed basis sets, distance based force-field reductions, and using coordinate specific basis sets. The importance of accurately describing electronic interactions is found to diminish both for multicoordinate terms involving the displacement of remote atoms and when describing the interactions between more remote atoms within specific coordinates.
Collapse
|