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
|
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
|
4
|
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
|
5
|
Tran HK, Berkelbach TC. Vibrational heat-bath configuration interaction with semistochastic perturbation theory using harmonic oscillator or VSCF modals. J Chem Phys 2023; 159:194101. [PMID: 37965997 PMCID: PMC10653875 DOI: 10.1063/5.0172702] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/18/2023] [Accepted: 10/23/2023] [Indexed: 11/16/2023] Open
Abstract
Vibrational heat-bath configuration interaction (VHCI)-a selected configuration interaction technique for vibrational structure theory-has recently been developed in two independent works [J. H. Fetherolf and T. C. Berkelbach, J. Chem. Phys. 154, 074104 (2021); A. U. Bhatty and K. R. Brorsen, Mol. Phys. 119, e1936250 (2021)], where it was shown to provide accuracy on par with the most accurate vibrational structure methods with a low computational cost. Here, we eliminate the memory bottleneck of the second-order perturbation theory correction using the same (semi)stochastic approach developed previously for electronic structure theory. This allows us to treat, in an unbiased manner, much larger perturbative spaces, which are necessary for high accuracy in large systems. Stochastic errors are easily controlled to be less than 1 cm-1. We also report two other developments: (i) we propose a new heat-bath criterion and an associated exact implicit sorting algorithm for potential energy surfaces expressible as a sum of products of one-dimensional potentials; (ii) we formulate VHCI to use a vibrational self-consistent field (VSCF) reference, as opposed to the harmonic oscillator reference configuration used in previous reports. Our tests are done with quartic and sextic force fields, for which we find that with VSCF, the minor improvements to accuracy are outweighed by the higher computational cost associated the matrix element evaluations. We expect VSCF-based VHCI to be important for more general potential representations, for which the harmonic oscillator basis function integrals are no longer analytic.
Collapse
Affiliation(s)
- Henry K. Tran
- Department of Chemistry, Columbia University, New York, New York 10027, USA
| | | |
Collapse
|
6
|
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
|
7
|
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
|
8
|
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
|
9
|
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
|
10
|
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
|
11
|
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
|
12
|
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
|
13
|
Yadav S, Banik S, Prasad MD. Understanding of the C-H stretch region of infra-red spectroscopy: an analysis of the final state wavefunctions. Phys Chem Chem Phys 2021; 23:9176-9188. [PMID: 33885051 DOI: 10.1039/d0cp01157f] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The nature of the wavefunctions associated with the final states in the CH stretching region of several medium sized molecules is analysed. The number of optically bright transitions is much larger than the number of CH oscillators present in the molecule, and they are spread over a range of about 300 cm-1. Several of them are clustered together within about 5 cm-1 with near equal intensities. The final states of all these transitions are superpositions of multiple zeroth order states. In almost all of such superpositions, no single zeroth order state has more than 50% weight. Several multiquantum states, with three to four quanta of excitation dominate the final states, with the CH chromophore contributing only a small weightage. Thus the band structure of the CH stretch region is due to several optically bright transitions whose final states are superpositions of low frequency multiquantum states with the CH chromophore contributing only a small weight to make them spectroscopically active.
Collapse
Affiliation(s)
- Swati Yadav
- School of Chemistry, University of Hyderabad, Hyderabad 500046, India.
| | | | | |
Collapse
|
14
|
Huang QR, Endo T, Mishra S, Zhang B, Chen LW, Fujii A, Jiang L, Patwari GN, Matsuda Y, Kuo JL. Understanding Fermi resonances in the complex vibrational spectra of the methyl groups in methylamines. Phys Chem Chem Phys 2021; 23:3739-3747. [PMID: 33533768 DOI: 10.1039/d0cp05745b] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Abstract
Vibrational spectra of the methyl groups in mono-methylamine (MMA), dimethylamine (DMA), and trimethylamine (TMA) monomers and their clusters were measured in three experimental set-ups to capture their complex spectral features as a result of bend/umbrella-stretch Fermi resonance (FR). Multiple bands were observed between 2800 and 3000 cm-1 corresponding to the methyl groups for MMA and DMA. On the other hand, the corresponding spectrum of TMA is relatively simple, exhibiting only four prominent bands in the same frequency window, even though TMA has a larger number of methyl groups. The discrete variable representation (DVR) based ab initio anharmonic algorithm with potential energy surface (PES) at CCSD/aug-cc-pVDZ quality is able to capture all the experimentally observed spectral features across all three amines, and the constructed vibrational Hamiltonian was used to analyze the couplings that give rise to the observed FR patterns. It was observed that the vibrational coupling among CH stretch modes on different methyl groups is weak (less than 2 cm-1) and stronger vibrational coupling is found to localize within a methyl group. In MMA and DMA, the complex feature between 2850 and 2950 cm-1 is a consequence of closely packed overtone states that gain intensities by mixing with the stretching modes. The simplification of the spectral pattern of TMA can be understood by the red-shift of the symmetric CH3 stretching modes by about 80 cm-1 relative to MMA, which causes the symmetric CH3 stretch to shift outside the FR window.
Collapse
Affiliation(s)
- Qian-Rui Huang
- Institute of Atomic and Molecular Sciences, Academia Sinica, Taipei 10617, Taiwan.
| | - Tomoya Endo
- Department of Chemistry, Graduate School of Science, Tohoku University, Aramaki-Aza-Aoba 6-3, Aoba-ku, Sendai 980-8578, Japan.
| | - Saurabh Mishra
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India.
| | - Bingbing Zhang
- State Key Laboratory of Molecular Reaction Dynamics, Collaborative Innovation Center of Chemistry for Energy and Materials (iChEM), Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China.
| | - Li-Wei Chen
- Institute of Atomic and Molecular Sciences, Academia Sinica, Taipei 10617, Taiwan.
| | - Asuka Fujii
- Department of Chemistry, Graduate School of Science, Tohoku University, Aramaki-Aza-Aoba 6-3, Aoba-ku, Sendai 980-8578, Japan.
| | - Ling Jiang
- State Key Laboratory of Molecular Reaction Dynamics, Collaborative Innovation Center of Chemistry for Energy and Materials (iChEM), Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China.
| | - G Naresh Patwari
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India.
| | - Yoshiyuki Matsuda
- Department of Chemistry, Graduate School of Science, Tohoku University, Aramaki-Aza-Aoba 6-3, Aoba-ku, Sendai 980-8578, Japan.
| | - Jer-Lai Kuo
- Institute of Atomic and Molecular Sciences, Academia Sinica, Taipei 10617, Taiwan.
| |
Collapse
|
15
|
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
|
16
|
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
|
17
|
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
|
18
|
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
|
19
|
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
|
20
|
Chaabene M, Agren S, Allouche A, Lahcinie M, Ben Chaâbane R, Baouab MHV. Theoretical and experimental investigations of complexation with BF
3
.Et
2
O effects on electronic structures, energies and photophysical properties of Anil and tetraphenyl (hydroxyl) imidazol. Appl Organomet Chem 2019. [DOI: 10.1002/aoc.5218] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Affiliation(s)
- Marwa Chaabene
- Laboratory of Advanced Materials and Interfaces (LIMA), Faculty of Sciences of MonastirUniversity of Monastir Avenue of the Environment 5000 Monastir Tunisia
| | - Soumaya Agren
- Research Unit Materials and Organic Synthesis (UR17ES31), Preparatory Institute for Engineering Studies of MonastirUniversity of Monastir Avenue of the Environment 5000 Monastir Tunisia
| | - Abdul‐Rahman Allouche
- Institute of Light and MatterUMR5306 University of Lyon 1‐CNRS, University of Lyon 69622 Villeurbanne cedex France
| | - Mohamed Lahcinie
- Laboratory of organometallic and macromolecular chemistry‐composites Materials, Faculty of Sciences and TechnologiesCadi Ayyad University Avenue Abdelhakim Elkhattabi, BP549 40000 Marrakech Morocco
- Mohamed VI Polytechnic University Lot 660, Hay Moulay Rachid 43150 Ben Guerir Morocco
| | - Rafik Ben Chaâbane
- Laboratory of Advanced Materials and Interfaces (LIMA), Faculty of Sciences of MonastirUniversity of Monastir Avenue of the Environment 5000 Monastir Tunisia
| | - Mohamed Hassen V. Baouab
- Research Unit Materials and Organic Synthesis (UR17ES31), Preparatory Institute for Engineering Studies of MonastirUniversity of Monastir Avenue of the Environment 5000 Monastir Tunisia
| |
Collapse
|
21
|
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
|
22
|
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
|
23
|
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]
|
24
|
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
|
25
|
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]
|
26
|
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
|
27
|
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
|
28
|
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
|
29
|
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
|
30
|
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]
|
31
|
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
|
32
|
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
|
33
|
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
|
34
|
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
|
35
|
Biczysko M, Bloino J, Puzzarini C. Computational challenges in Astrochemistry. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2017. [DOI: 10.1002/wcms.1349] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Affiliation(s)
- Malgorzata Biczysko
- International Center for Quantum and Molecular Structures, College of SciencesShanghai University Shanghai China
| | - Julien Bloino
- Consiglio Nazionale delle Ricerche, Istituto di Chimica dei Composti OrganoMetalliciUOS di Pisa, Area della Ricerca CNR Pisa Italy
- Scuola Normale Superiore Classe di Scienze, Pisa Italy
| | - Cristina Puzzarini
- Department of Chemistry “Giacomo Ciamician”University of Bologna Bologna Italy
| |
Collapse
|
36
|
Brown SE, Götz AW, Cheng X, Steele RP, Mandelshtam VA, Paesani F. Monitoring Water Clusters “Melt” Through Vibrational Spectroscopy. J Am Chem Soc 2017; 139:7082-7088. [DOI: 10.1021/jacs.7b03143] [Citation(s) in RCA: 55] [Impact Index Per Article: 7.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Affiliation(s)
| | | | - Xiaolu Cheng
- 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
|
37
|
Yagi K, Thomsen B. Infrared Spectra of Protonated Water Clusters, H+(H2O)4, in Eigen and Zundel Forms Studied by Vibrational Quasi-Degenerate Perturbation Theory. J Phys Chem A 2017; 121:2386-2398. [DOI: 10.1021/acs.jpca.6b11189] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Affiliation(s)
- Kiyoshi Yagi
- Theoretical
Molecular Science Laboratory and ‡iTHES, RIKEN, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
| | - Bo Thomsen
- Theoretical
Molecular Science Laboratory and ‡iTHES, RIKEN, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
| |
Collapse
|