1
|
Schneider M, Rauhut G. VSCF/VCI theory based on the Podolsky Hamiltonian. J Chem Phys 2024; 160:214118. [PMID: 38832735 DOI: 10.1063/5.0213401] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/10/2024] [Accepted: 05/14/2024] [Indexed: 06/05/2024] Open
Abstract
While the vibrational spectra of semi-rigid molecules can be computed on approaches relying on the Watson Hamiltonian, floppy molecules or molecular clusters are better described by Hamiltonians, which are capable of dealing with any curvilinear coordinates. It is the kinetic energy operator (KEO) of these Hamiltonians, which render the correlated calculations relying on them rather costly. Novel implementation of vibrational self-consistent field theory and vibrational configuration interaction theory on the basis of the Podolsky Hamiltonian are reported, in which the inverse of the metric tensor, i.e., the G matrix, is represented by an n-mode expansion expressed in terms of polynomials. An analysis of the importance of the individual terms of the KEO with respect to the truncation orders of the n-mode expansion is provided. Benchmark calculations have been performed for the cis-HOPO and methanimine, H2CNH, molecules and are compared to experimental data and to calculations based on the Watson Hamiltonian and the internal coordinate path Hamiltonian.
Collapse
Affiliation(s)
- Moritz Schneider
- 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
|
2
|
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
|
3
|
Bader F, Lauvergnat D, Christiansen O. Efficient vibrationally correlated calculations using n-mode expansion-based kinetic energy operators. Phys Chem Chem Phys 2024; 26:11469-11481. [PMID: 38546727 DOI: 10.1039/d4cp00423j] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/18/2024]
Abstract
Due to its efficiency and flexibility, the n-mode expansion is a frequently used tool for representing molecular potential energy surfaces in quantum chemical simulations. In this work, we investigate the performance of n-mode expansion-based models of kinetic energy operators in general polyspherical coordinate systems. In particular, we assess the operators with respect to accuracy in vibrationally correlated calculations and their effect on potential energy surface construction with the adaptive density guided approach. Our results show that the n-mode expansion-based operator variants are reliable and systematically improvable approximations of the full kinetic energy operator. Moreover, we introduce a workflow to generate the n-mode expanded kinetic energy operators on-the-fly within the adaptive density guided approach. This scheme can be applied in studies of species and coordinate systems, for which an analytical form of the kinetic energy operator is not available.
Collapse
Affiliation(s)
- Frederik Bader
- Department of Chemistry, Aarhus University, DK-8000 Aarhus C, Denmark.
| | - David Lauvergnat
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, Orsay 91405, France.
| | - Ove Christiansen
- Department of Chemistry, Aarhus University, DK-8000 Aarhus C, Denmark.
| |
Collapse
|
4
|
Schröder B, Rauhut G. From the Automated Calculation of Potential Energy Surfaces to Accurate Infrared Spectra. J Phys Chem Lett 2024; 15:3159-3169. [PMID: 38478898 PMCID: PMC10961845 DOI: 10.1021/acs.jpclett.4c00186] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/19/2024] [Revised: 02/20/2024] [Accepted: 02/28/2024] [Indexed: 03/22/2024]
Abstract
Advances in the development of quantum chemical methods and progress in multicore architectures in computer science made the simulation of infrared spectra of isolated molecules competitive with respect to established experimental methods. Although it is mainly the multidimensional potential energy surface that controls the accuracy of these calculations, the subsequent vibrational structure calculations need to be carefully converged in order to yield accurate results. As both aspects need to be considered in a balanced way, we focus on approaches for molecules of up to 12-15 atoms with respect to both parts, which have been automated to some extent so that they can be employed in routine applications. Alternatives to machine learning will be discussed, which appear to be attractive, as long as local regions of the potential energy surface are sufficient. The automatization of these methods is still in its infancy, and the generalization to molecules with large amplitude motions or molecular clusters is far from trivial, but many systems relevant for astrophysical studies are already in reach.
Collapse
Affiliation(s)
- Benjamin Schröder
- Institute
of Physical Chemistry, University of Goettingen, Tammannstrasse 6, Göttingen 37077, Germany
| | - Guntram Rauhut
- Institute
for Theoretical Chemistry, University of
Stuttgart, Pfaffenwaldring 55, Stuttgart 70569, Germany
| |
Collapse
|
5
|
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
|
6
|
Qin X, Hirata S. Finite-temperature many-body perturbation theory for anharmonic vibrations: Recursions, algebraic reduction, second-quantized reduction, diagrammatic rules, linked-diagram theorem, finite-temperature self-consistent field, and general-order algorithm. J Chem Phys 2023; 159:084114. [PMID: 37638629 DOI: 10.1063/5.0164326] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/22/2023] [Accepted: 07/24/2023] [Indexed: 08/29/2023] Open
Abstract
A unified theory is presented for finite-temperature many-body perturbation expansions of the anharmonic vibrational contributions to thermodynamic functions, i.e., the free energy, internal energy, and entropy. The theory is diagrammatically size-consistent at any order, as ensured by the linked-diagram theorem proved in this study, and, thus, applicable to molecular gases and solids on an equal footing. It is also a basis-set-free formalism, just like its underlying Bose-Einstein theory, capable of summing anharmonic effects over an infinite number of states analytically. It is formulated by the Rayleigh-Schrödinger-style recursions, generating sum-over-states formulas for the perturbation series, which unambiguously converges at the finite-temperature vibrational full-configuration-interaction limits. Two strategies are introduced to reduce these sum-over-states formulas into compact sum-over-modes analytical formulas. One is a purely algebraic method that factorizes each many-mode thermal average into a product of one-mode thermal averages, which are then evaluated by the thermal Born-Huang rules. Canonical forms of these rules are proposed, dramatically expediting the reduction process. The other is finite-temperature normal-ordered second quantization, which is fully developed in this study, including a proof of thermal Wick's theorem and the derivation of a normal-ordered vibrational Hamiltonian at finite temperature. The latter naturally defines a finite-temperature extension of size-extensive vibrational self-consistent field theory. These reduced formulas can be represented graphically as Feynman diagrams with resolvent lines, which include anomalous and renormalization diagrams. Two order-by-order and one general-order algorithms of computing these perturbation corrections are implemented and applied up to the eighth order. The results show no signs of Kohn-Luttinger-type nonconvergence.
Collapse
Affiliation(s)
- Xiuyi Qin
- Department of Chemistry, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| | - So Hirata
- Department of Chemistry, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| |
Collapse
|
7
|
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
|
8
|
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
|
9
|
Tan JA, Barbarona RF, Kuo JL. Approaching the "Zundel" Limit: Tuning the Vibrational Coupling in N 2H +Ng, Ng = {He, Ne, Ar, Kr, Xe, and Rn}. J Phys Chem A 2023. [PMID: 37418837 DOI: 10.1021/acs.jpca.3c01956] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/09/2023]
Abstract
The diazenylium ion (N2H+) is a ubiquitous ion in dense molecular clouds. This ion is often used as a dense gas tracer in outer space. Most of the previous works on diazenylium ion have focused on the shared-proton stretch band, νH+. In this work, we have performed reduced-dimensional calculations to investigate the vibrational structure of N2H+Ng, Ng = {He, Ne, Ar, Kr, Xe, and Rn}. We demonstrate a few interesting things about this system. First, the vibrational coupling in N2H+ can be tuned to switch on interesting anharmonic effects such as Fermi resonance or combination bands by tagging it with different noble gases. Second, a comparison of the vibrational spectrum from N2H+He to N2H+Rn shows that the νH+ can be swept from an "Eigen-like" to a "Zundel-like" limiting case. Anharmonic calculations were performed using a multilevel approach, which utilized the MP2 and CCSD(T) levels of theories. Binding energies for the elimination of Ng in N2H+Ng are also reported.
Collapse
Affiliation(s)
- Jake A Tan
- Institute of Atomic and Molecular Science, Academia Sinica, Taipei 10617, Taiwan
| | - Rona F Barbarona
- Institute of Atomic and Molecular Science, Academia Sinica, Taipei 10617, Taiwan
- Molecular Science and Technology Program, Taiwan International Graduate Program, Academia Sinica, Taipei 11529, Taiwan
- Department of Physics, National Central University, Taoyuan 32001, Taiwan
| | - Jer-Lai Kuo
- Institute of Atomic and Molecular Science, Academia Sinica, Taipei 10617, Taiwan
- Molecular Science and Technology Program, Taiwan International Graduate Program, Academia Sinica, Taipei 11529, Taiwan
| |
Collapse
|
10
|
Schiltz C, Rappoport D, Mandelshtam VA. Implementation of the self-consistent phonons method with ab initio potentials (AI-SCP). J Chem Phys 2023; 158:2890485. [PMID: 37184023 DOI: 10.1063/5.0146682] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2023] [Accepted: 04/24/2023] [Indexed: 05/16/2023] Open
Abstract
The self-consistent phonon (SCP) method allows one to include anharmonic effects when treating a many-body quantum system at thermal equilibrium. The system is then described by an effective temperature-dependent harmonic Hamiltonian, which can be used to estimate its various dynamic and static properties. In this paper, we combine SCP with ab initio (AI) potential energy evaluation in which case the numerical bottleneck of AI-SCP is the evaluation of Gaussian averages of the AI potential energy and its derivatives. These averages are computed efficiently by the quasi-Monte Carlo method utilizing low-discrepancy sequences leading to a fast convergence with respect to the number, S, of the AI energy evaluations. Moreover, a further substantial (an-order-of-magnitude) improvement in efficiency is achieved once a numerically cheap approximation of the AI potential is available. This is based on using a perturbation theory-like (the two-grid) approach in which it is the average of the difference between the AI and the approximate potential that is computed. The corresponding codes and scripts are provided.
Collapse
Affiliation(s)
- Colin Schiltz
- Department of Chemistry, University of California, Irvine, California 92697, USA
| | - Dmitrij Rappoport
- Department of Chemistry, University of California, Irvine, California 92697, USA
| | | |
Collapse
|
11
|
Schneider M, Born D, Kästner J, Rauhut G. Positioning of grid points for spanning potential energy surfaces-How much effort is really needed? J Chem Phys 2023; 158:144118. [PMID: 37061506 DOI: 10.1063/5.0146020] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/17/2023] Open
Abstract
The positions of grid points for representing a multidimensional potential energy surface (PES) have a non-negligible impact on its accuracy and the associated computational effort for its generation. Six different positioning schemes were studied for PESs represented by n-mode expansions as needed for the accurate calculation of anharmonic vibrational frequencies by means of vibrational configuration interaction theory. A static approach, which has successfully been used in many applications, and five adaptive schemes based on Gaussian process regression have been investigated with respect to the number of necessary grid points and the accuracy of the fundamental modes for a small set of test molecules. A comparison with a related, more sophisticated, and consistent approach by Christiansen et al. is provided. The impact of the positions of the ab initio grid points is discussed for multilevel PESs, for which the computational effort of the individual electronic structure calculations decreases for increasing orders of the n-mode expansion. As a result of that, the ultimate goal is not the maximal reduction of grid points but rather the computational cost, which is not directly related.
Collapse
Affiliation(s)
- Moritz Schneider
- Institute for Theoretical Chemistry, University of Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Daniel Born
- Institute for Theoretical Chemistry, University of Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Johannes Kästner
- 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
|
12
|
Schneider M, Rauhut G. Quantum chemical rovibrational analysis of aminoborane and its isotopologues. J Comput Chem 2023; 44:298-306. [PMID: 35582830 DOI: 10.1002/jcc.26893] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/25/2022] [Revised: 04/22/2022] [Accepted: 04/28/2022] [Indexed: 12/31/2022]
Abstract
Aminoborane, H2 NBH2 and its isotopologues, H2 N10 BH2 , D2 NBD2 , and D2 N10 BD2 , have been studied by high-level ab initio methods. All calculations rely on multidimensional potential energy surfaces and dipole moment surfaces including high-order mode coupling terms, which have been obtained from electronic structure calculations at the level of explicitly correlated coupled-cluster theory, CCSD(T)-F12, or the distinguishable cluster approximation, DCSD. Subsequent vibrational structure calculations based on second-order vibrational perturbation theory, VPT2, and vibrational configuration interaction theory, VCI, were used to determine rotational constants, centrifugal distortion constants, vibrationally averaged geometrical parameters and (an)harmonic vibrational frequencies. The impact of core-correlation effects is discussed in detail. Rovibrational VCI calculations were used to simulate the gas phase spectra of these species and an in-depth analysis of the ν7 band of aminoborane is provided. Color-coding is used to reveal the identity of the individual progressions of the rovibrational transitions for this particular mode.
Collapse
Affiliation(s)
- Moritz Schneider
- Institute for Theoretical Chemistry, University of Stuttgart, Stuttgart, Germany
| | - Guntram Rauhut
- Institute for Theoretical Chemistry, University of Stuttgart, Stuttgart, Germany
| |
Collapse
|
13
|
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
|
14
|
Tan JA, Kuo JL. Spectral Signatures of Protonated Noble Gas Clusters of Ne, Ar, Kr, and Xe: From Monomers to Trimers. MOLECULES (BASEL, SWITZERLAND) 2022; 27:molecules27103198. [PMID: 35630674 PMCID: PMC9143425 DOI: 10.3390/molecules27103198] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 04/01/2022] [Revised: 05/08/2022] [Accepted: 05/13/2022] [Indexed: 11/25/2022]
Abstract
The structures and spectral features of protonated noble gas clusters are examined using a first principles approach. Protonated noble gas monomers (NgH+) and dimers (NgH+Ng) have a linear structure, while the protonated noble gas trimers (Ng3H+) can have a T-shaped or linear structure. Successive binding energies for these complexes are calculated at the CCSD(T)/CBS level of theory. Anharmonic simulations for the dimers and trimers unveil interesting spectral features. The symmetric NgH+Ng are charactized by a set of progression bands, which involves one quantum of the asymmetric Ng-H+ stretch with multiple quanta of the symmetric Ng-H+ stretch. Such a spectral signature is very robust and is predicted to be observed in both T-shaped and linear isomers of Ng3H+. Meanwhile, for selected asymmetric NgH+Ng’, a Fermi resonance interaction involving the first overtone of the proton bend with the proton stretch is predicted to occur in ArH+Kr and XeH+Kr.
Collapse
|
15
|
Schröder B, Rauhut G. Comparison of body definitions for incremental vibrational configuration interaction theory (iVCI). J Chem Phys 2022; 156:174103. [DOI: 10.1063/5.0085082] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Within incremental vibrational configuration interaction theory (iVCI), the vibrational state energy is determined by means of a many-body expansion, i.e., it is a sum of terms of increasing order, which allow for an embarrassingly parallel evaluation. The convergence of this expansion depends strongly on the definition of the underlying bodies, which essentially decompose the correlation space into fragments. The different definitions considered here comprise mode-based bodies, excitation level-based bodies, and energy-based bodies. An analysis of the convergence behavior revealed that accounting for resonances within these definitions is mandatory and leads to a substantial improvement of the convergence, that is, the expansions can be truncated at lower orders. Benchmark calculations and systematic comparisons of the different body definitions for a small set of molecules, i.e., ketene, ethene, and diborane, have been conducted to study the overall performance of these iVCI implementations with respect to accuracy and central processing unit time.
Collapse
Affiliation(s)
- Benjamin Schröder
- 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
|
Erfort S, Tschoepe M, Rauhut G. Efficient and Automated Quantum Chemical Calculation of Rovibrational Nonresonant Raman Spectra. J Chem Phys 2022; 156:124102. [DOI: 10.1063/5.0087359] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
An outline of a newly developed program for the simulation of rovibrational nonresonant Raman spectra is presented. This program is an extension of our recently developed code for rovibrational infrared spectra [J. Chem Phys. 152 (2020) 244104] and relies on vibrational wavefunctions from variational configuration interaction theory to allow for an almost fully automated calculation of such spectra in pure ab initio fashion. Due to efficient contraction schemes this program requires modest computational resources and it can be controlled by only a few lines of input. As the required polarizability surfaces are also computed in an automated fashion, this implementation enables the routine application to small molecules. For demonstrating its capabilities, benchmark calculations for water H216O are compared to reference data and spectra for the beryllium dihydride dimer, Be2H4 (D2h), are predicted. The inversion symmetry of the D2h systems lead to complementary infrared and Raman spectra, which are needed both for a comprehensive investigation of this system.
Collapse
Affiliation(s)
- Sebastian Erfort
- Institute for Theoretical Chemistry, University of Stuttgart Faculty of Chemistry, Germany
| | | | - Guntram Rauhut
- Institut fuer Theoretische Chemie, University of Stuttgart Faculty of Chemistry, Germany
| |
Collapse
|
17
|
Mathea T, Rauhut G. Advances in vibrational configuration interaction theory - part 1: Efficient calculation of vibrational angular momentum terms. J Comput Chem 2021; 42:2321-2333. [PMID: 34651703 DOI: 10.1002/jcc.26762] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/10/2021] [Revised: 09/20/2021] [Accepted: 09/20/2021] [Indexed: 02/02/2023]
Abstract
Finite basis vibrational configuration interaction theory (VCI) is a highly accurate method for the variational calculation of state energies and related properties, but suffers from fast growing computational costs in dependence of the size of the correlation space. In this series of papers, concepts and techniques will be presented, which diminish the computational demands and thus broaden the applicability of this method to larger molecules or more complex situations. This first part focuses on a highly efficient implementation of the vibrational angular momentum (VAM) terms as occurring in the Watson Hamiltonian and the prediagonalization of initial subspaces within an iterative configuration selective VCI implementation. Working equations and benchmark calculations are provided, the latter demonstrating the increased performance of the new algorithm.
Collapse
Affiliation(s)
- Tina Mathea
- Institute for Theoretical Chemistry, University of Stuttgart, Stuttgart, Germany
| | - Guntram Rauhut
- Institute for Theoretical Chemistry, University of Stuttgart, Stuttgart, Germany
| |
Collapse
|
18
|
Tan JA, Teh S, Kuo JL. Structural and vibrational characterization of HCO + and Rg-HCO +, Rg = {He, Ne, Ar, Kr, and Xe}. J Chem Phys 2021; 155:174306. [PMID: 34742194 DOI: 10.1063/5.0069726] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The structures of the formyl ion (HCO+) and its rare gas tagged counterparts (Rg-HCO+, Rg = He, Ne, Ar, Kr, and Xe) were studied at the coupled-cluster singles, doubles, and perturbative triples [CCSD(T)]/aug-cc-pVTZ level of theory and basis set. A linear structure for these tagged complexes was predicted. The Rg binding energies for Rg-HCO+ are also examined at the CCSD(T) level. It was found that the binding interaction increases from He-HCO+ to Xe-HCO+. A multilevel potential energy surface built at the CCSD(T) and second-order Møller-Plesset perturbation levels of theory were used to study these species' vibrational spectra. By changing the Rg in the first-solvation shell for HCO+, the Fermi resonance interaction between the first H+ bend overtone and the asymmetric and symmetric H-C-O stretches can be modulated. This Fermi resonance modulation is demonstrated by examining a series of rare gas solvated HCO+.
Collapse
Affiliation(s)
- Jake A Tan
- Institute of Atomic and Molecular Sciences, Academia Sinica, No. 1 Roosevelt Road, Section 4, Da-an District, Taipei City 10617, Taiwan
| | - Soon Teh
- Institute of Atomic and Molecular Sciences, Academia Sinica, No. 1 Roosevelt Road, Section 4, Da-an District, Taipei City 10617, Taiwan
| | - Jer-Lai Kuo
- Institute of Atomic and Molecular Sciences, Academia Sinica, No. 1 Roosevelt Road, Section 4, Da-an District, Taipei City 10617, Taiwan
| |
Collapse
|
19
|
Towards complete assignment of the infrared spectrum of the protonated water cluster H +(H 2O) 21. Nat Commun 2021; 12:6141. [PMID: 34686665 PMCID: PMC8536673 DOI: 10.1038/s41467-021-26284-x] [Citation(s) in RCA: 27] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/24/2021] [Accepted: 09/22/2021] [Indexed: 11/08/2022] Open
Abstract
The spectroscopic features of protonated water species in dilute acid solutions have been long sought after for understanding the microscopic behavior of the proton in water with gas-phase water clusters H+(H2O)n extensively studied as bottom-up model systems. We present a new protocol for the calculation of the infrared (IR) spectra of complex systems, which combines the fragment-based Coupled Cluster method and anharmonic vibrational quasi-degenerate perturbation theory, and demonstrate its accuracy towards the complete and accurate assignment of the IR spectrum of the H+(H2O)21 cluster. The site-specific IR spectral signatures reveal two distinct structures for the internal and surface four-coordinated water molecules, which are ice-like and liquid-like, respectively. The effect of inter-molecular interaction between water molecules is addressed, and the vibrational resonance is found between the O-H stretching fundamental and the bending overtone of the nearest neighboring water molecule. The revelation of the spectral signature of the excess proton offers deeper insight into the nature of charge accommodation in the extended hydrogen-bonding network underpinning this aqueous cluster. Protonated water species have been the subject of numerous experimental and computational studies. Here the authors provide a nearly complete assignment of the experimental IR spectrum of the H+(H2O)21 water cluster based on high-level wavefunction theory and anharmonic vibrational quasi-degenerate perturbation theory.
Collapse
|
20
|
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: 8] [Impact Index Per Article: 2.7] [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
|
21
|
Tan JA, Kuo JL. Fermi resonance switching in KrH +Rg and XeH +Rg (Rg = Ne, Ar, Kr, and Xe). J Chem Phys 2021; 154:134302. [PMID: 33832263 DOI: 10.1063/5.0044703] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Abstract
Matrix isolation experiments have been successfully employed to extensively study the infrared spectrum of several proton-bound rare gas complexes. Most of these studies have focused on the spectral signature for the H+ stretch (ν3) and its combination bands with the intermolecular stretch coordinate (ν1). However, little attention has been paid to the Fermi resonance interaction between the H+ stretch (ν3) and H+ bend overtone (2ν2) in the asymmetric proton-bound rare gas dimers, RgH+Rg'. In this work, we have investigated this interaction on KrH+Rg and XeH+Rg with Rg = (Ne, Ar, Kr, and Xe). A multilevel potential energy surface (PES) was used to simulate the vibrational structure of these complexes. This PES is a dual-level comprising of second-order Møller-Plesset perturbation theory and coupled-cluster singles doubles with perturbative triples [CCSD(T)] levels of ab initio theories. We found that when both the combination bands (nν1 + ν3) and bend overtone 2ν2 compete to borrow intensity from the ν3 band, the latter wins over the former, which then results in the suppression of the nν1 + ν3 bands. The current simulations offer new assignments for the ArH+Xe and KrH+Xe spectra. Complete basis set (CBS) binding energies for these complexes were also calculated at the CCSD(T)/CBS level.
Collapse
Affiliation(s)
- Jake A Tan
- Institute of Atomic and Molecular Sciences, Academia Sinica, No. 1 Roosevelt Road, Section 4, Da-an District, Taipei City 10617, Taiwan
| | - Jer-Lai Kuo
- Institute of Atomic and Molecular Sciences, Academia Sinica, No. 1 Roosevelt Road, Section 4, Da-an District, Taipei City 10617, Taiwan
| |
Collapse
|
22
|
Tschöpe M, Schröder B, Erfort S, Rauhut G. High-Level Rovibrational Calculations on Ketenimine. Front Chem 2021; 8:623641. [PMID: 33585403 PMCID: PMC7873934 DOI: 10.3389/fchem.2020.623641] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2020] [Accepted: 12/03/2020] [Indexed: 11/13/2022] Open
Abstract
From an astrochemical point of view ketenimine (CH2CNH) is a complex organic molecule (COM) and therefore likely to be a building block for biologically relevant molecules. Since it has been detected in the star-forming region Sagittarius B2(N), it is of high relevance in this field. Although experimental data are available for certain bands, for some energy ranges such as above 1200 cm-1 reliable data virtually do not exist. In addition, high-level ab initio calculations are neither reported for ketenimine nor for one of its deuterated isotopologues. In this paper, we provide for the first time data from accurate quantum chemical calculations and a thorough analysis of the full rovibrational spectrum. Based on high-level potential energy surfaces obtained from explicitly correlated coupled-cluster calculations including up to 4-mode coupling terms, the (ro)vibrational spectrum of ketenimine has been studied in detail by variational calculations relying on rovibrational configuration interaction (RVCI) theory. Strong Fermi resonances were found for all isotopologues. Rovibrational infrared intensities have been obtained from dipole moment surfaces determined from the distinguishable cluster approximation. A comparison of the spectra of the CH2CNH molecule with experimental data validates our results, but also reveals new insight about the system, which shows very strong Coriolis coupling effects.
Collapse
Affiliation(s)
- Martin Tschöpe
- Institute for Theoretical Chemistry, University of Stuttgart, Stuttgart, Germany
| | - Benjamin Schröder
- Institute for Theoretical Chemistry, University of Stuttgart, Stuttgart, Germany
| | - Sebastian Erfort
- Institute for Theoretical Chemistry, University of Stuttgart, Stuttgart, Germany
| | - Guntram Rauhut
- Institute for Theoretical Chemistry, University of Stuttgart, Stuttgart, Germany
| |
Collapse
|
23
|
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
|
24
|
Nejad A, Crittenden DL. On the separability of large-amplitude motions in anharmonic frequency calculations. Phys Chem Chem Phys 2020; 22:20588-20601. [PMID: 32966420 DOI: 10.1039/d0cp03515g] [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
Nuclear vibrational theories based upon the Watson Hamiltonian are ubiquitous in quantum chemistry, but are generally unable to model systems in which the wavefunction can delocalise over multiple energy minima, i.e. molecules that have low-energy torsion and inversion barriers. In a 2019 Chemical Reviews article, Puzzarini et al. note that a common workaround is to simply decouple these problematic modes from all other vibrations in the system during anharmonic frequency calculations. They also point out that this approximation can be "ill-suited", but do not quantify the errors introduced. In this work, we present the first systematic investigation into how separating out or constraining torsion and inversion vibrations within potential energy surface (PES) expansions affects the accuracy of computed fundamental wavenumbers for the remaining vibrational modes, using a test set of 19 tetratomic molecules for which high quality analytic potential energy surfaces and fully-coupled anharmonic reference fundamental frequencies are available. We find that the most effective and efficient strategy is to remove the mode in question from the PES expansion entirely. This introduces errors of up to +10 cm-1 in stretching fundamentals that would otherwise couple to the dropped mode, and ±5 cm-1 in all other fundamentals. These errors are approximately commensurate with, but not necessarily additional to, errors due to the choice of electronic structure model used in constructing spectroscopically accurate PES.
Collapse
Affiliation(s)
- Arman Nejad
- Institute of Physical Chemistry, University of Göttingen, Tammannstr. 6, D-37077 Göttingen, Germany.
| | - Deborah L Crittenden
- School of Physical and Chemical Sciences, University of Canterbury, Christchurch, New Zealand
| |
Collapse
|
25
|
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
|
26
|
Werner HJ, Knowles PJ, Manby FR, Black JA, Doll K, Heßelmann A, Kats D, Köhn A, Korona T, Kreplin DA, Ma Q, Miller TF, Mitrushchenkov A, Peterson KA, Polyak I, Rauhut G, Sibaev M. The Molpro quantum chemistry package. J Chem Phys 2020; 152:144107. [PMID: 32295355 DOI: 10.1063/5.0005081] [Citation(s) in RCA: 475] [Impact Index Per Article: 118.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Abstract
Molpro is a general purpose quantum chemistry software package with a long development history. It was originally focused on accurate wavefunction calculations for small molecules but now has many additional distinctive capabilities that include, inter alia, local correlation approximations combined with explicit correlation, highly efficient implementations of single-reference correlation methods, robust and efficient multireference methods for large molecules, projection embedding, and anharmonic vibrational spectra. In addition to conventional input-file specification of calculations, Molpro calculations can now be specified and analyzed via a new graphical user interface and through a Python framework.
Collapse
Affiliation(s)
- Hans-Joachim Werner
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Peter J Knowles
- School of Chemistry, Cardiff University, Main Building, Park Place, Cardiff CF10 3AT, United Kingdom
| | - Frederick R Manby
- School of Chemistry, University of Bristol, Cantock's Close, Bristol BS8 1TS, United Kingdom
| | - Joshua A Black
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Klaus Doll
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Andreas Heßelmann
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Daniel Kats
- Max-Planck Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | - Andreas Köhn
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Tatiana Korona
- Faculty of Chemistry, University of Warsaw, L. Pasteura 1 St., 02-093 Warsaw, Poland
| | - David A Kreplin
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Qianli Ma
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Thomas F Miller
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | | | - Kirk A Peterson
- Department of Chemistry, Washington State University, Pullman, Washington 99164-4630, USA
| | - Iakov Polyak
- School of Chemistry, Cardiff University, Main Building, Park Place, Cardiff CF10 3AT, United Kingdom
| | - Guntram Rauhut
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Marat Sibaev
- School of Chemistry, Cardiff University, Main Building, Park Place, Cardiff CF10 3AT, United Kingdom
| |
Collapse
|
27
|
Petrenko T, Rauhut G. Account of non-Condon effects in time-independent Raman wavefunction theory: Calculation of the S 1 ← S 0 vibronic absorption spectrum of formaldehyde. J Chem Phys 2020; 152:114109. [DOI: 10.1063/5.0003272] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Taras Petrenko
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Guntram Rauhut
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| |
Collapse
|
28
|
|
29
|
Schmitz G, Godtliebsen IH, Christiansen O. Machine learning for potential energy surfaces: An extensive database and assessment of methods. J Chem Phys 2019; 150:244113. [DOI: 10.1063/1.5100141] [Citation(s) in RCA: 32] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022] Open
Affiliation(s)
- Gunnar Schmitz
- Department of Chemistry, Aarhus Universitet, DK-8000 Aarhus, Denmark
| | | | - Ove Christiansen
- Department of Chemistry, Aarhus Universitet, DK-8000 Aarhus, Denmark
| |
Collapse
|
30
|
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
|
31
|
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
|
32
|
Tan JA, Kuo JL. Multilevel Approach for Direct VSCF/VCI MULTIMODE Calculations with Applications to Large “Zundel” Cations. J Chem Theory Comput 2018; 14:6405-6416. [DOI: 10.1021/acs.jctc.8b00679] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Jake A. Tan
- Institute of Atomic and Molecular Sciences, Academia Sinica, No. 1 Roosevelt Road, Section 4, Taipei 10617, Taiwan (ROC)
| | - Jer-Lai Kuo
- Institute of Atomic and Molecular Sciences, Academia Sinica, No. 1 Roosevelt Road, Section 4, Taipei 10617, Taiwan (ROC)
| |
Collapse
|
33
|
Morisawa Y, Suga A. Effects of intermolecular interactions on absorption intensities of the fundamental and the first, second, and third overtones of OH stretching vibrations of methanol and t-butanol‑d 9 in n-hexane studied by visible/near-infrared/infrared spectroscopy. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2018; 197:121-125. [PMID: 29221935 DOI: 10.1016/j.saa.2017.11.059] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/19/2017] [Revised: 11/24/2017] [Accepted: 11/27/2017] [Indexed: 06/07/2023]
Abstract
Visible (Vis), near-infrared (NIR) and IR spectra in the 15,600-2500cm-1 region were measured for methanol, methanol-d3, and t-butanol-d9 in n-hexane to investigate effects of intermolecular interaction on absorption intensities of the fundamental and the first, second, and third overtones of their OH stretching vibrations. The relative area intensities of OH stretching bands of free and hydrogen-bonded species were plotted versus the vibrational quantum number using logarithm plots (V=1-4) for 0.5M methanol, 0.5M methanol‑d3, and 0.5M t-butanol-d9 in n-hexane. In the logarithm plots the relative intensities of free species yield a linear dependence irrespective of the solutes while those of hydrogen-bonded species deviate significantly from the linearity. The observed results suggest that the modifications in dipole moment functions of the OH bond induced by the formation of the hydrogen bondings change transient dipole moment, leading to the deviations of the dependences of relative absorption intensities on the vibrational quantum number from the linearity.
Collapse
Affiliation(s)
- Yusuke Morisawa
- Department of Chemistry, School of Science and Engineering, Kindai University, Higashi-Osaka, Osaka, Japan.
| | - Arisa Suga
- Department of Chemistry, School of Science and Technology, Kwansei Gakuin University, Gakuen, Sanda, Hyogo 669-1337, Japan
| |
Collapse
|
34
|
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
|
35
|
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
|
36
|
Petrenko T, Rauhut G. A General Approach for Calculating Strongly Anharmonic Vibronic Spectra with a High Density of States: The X̃2B1 ← X̃1A1 Photoelectron Spectrum of Difluoromethane. J Chem Theory Comput 2017; 13:5515-5527. [DOI: 10.1021/acs.jctc.7b00468] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Taras Petrenko
- Institute for Theoretical Chemistry, Pfaffenwaldring
55, 70569 Stuttgart, Germany
| | - Guntram Rauhut
- Institute for Theoretical Chemistry, Pfaffenwaldring
55, 70569 Stuttgart, Germany
| |
Collapse
|
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
|
38
|
Rai P, Sargsyan K, Najm H, Hermes MR, Hirata S. Low-rank canonical-tensor decomposition of potential energy surfaces: application to grid-based diagrammatic vibrational Green's function theory. Mol Phys 2017. [DOI: 10.1080/00268976.2017.1288937] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
Affiliation(s)
- Prashant Rai
- Sandia National Laboratories, Livermore, CA, USA
| | | | - Habib Najm
- Sandia National Laboratories, Livermore, CA, USA
| | - Matthew R. Hermes
- Department of Chemistry, University of Illinois at Urbana-Champaign, Urbana, IL, USA
| | - So Hirata
- Department of Chemistry, University of Illinois at Urbana-Champaign, Urbana, IL, USA
| |
Collapse
|
39
|
Panek PT, Jacob CR. On the benefits of localized modes in anharmonic vibrational calculations for small molecules. J Chem Phys 2017; 144:164111. [PMID: 27131535 DOI: 10.1063/1.4947213] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Anharmonic vibrational calculations can already be computationally demanding for relatively small molecules. The main bottlenecks lie in the construction of the potential energy surface and in the size of the excitation space in the vibrational configuration interaction (VCI) calculations. To address these challenges, we use localized-mode coordinates to construct potential energy surfaces and perform vibrational self-consistent field and L-VCI calculations [P. T. Panek and C. R. Jacob, ChemPhysChem 15, 3365 (2014)] for all vibrational modes of two prototypical test cases, the ethene and furan molecules. We find that the mutual coupling between modes is reduced when switching from normal-mode coordinates to localized-mode coordinates. When using such localized-mode coordinates, we observe a faster convergence of the n-mode expansion of the potential energy surface. This makes it possible to neglect higher-order contributions in the n-mode expansion of the potential energy surface or to approximate higher-order contributions in hybrid potential energy surfaces, which reduced the computational effort for the construction of the anharmonic potential energy surface significantly. Moreover, we find that when using localized-mode coordinates, the convergence with respect to the VCI excitation space proceeds more smoothly and that the error at low orders is reduced significantly. This makes it possible to devise low-cost models for obtaining a first approximation of anharmonic corrections. This demonstrates that the use of localized-mode coordinates can be beneficial already in anharmonic vibrational calculations of small molecules and provides a possible avenue for enabling such accurate calculations also for larger molecules.
Collapse
Affiliation(s)
- Paweł T Panek
- TU Braunschweig, Institute of Physical and Theoretical Chemistry, Hans-Sommer-Str. 10, 38106 Braunschweig, Germany
| | - Christoph R Jacob
- TU Braunschweig, Institute of Physical and Theoretical Chemistry, Hans-Sommer-Str. 10, 38106 Braunschweig, Germany
| |
Collapse
|
40
|
Otaki H, Yagi K, Ishiuchi SI, Fujii M, Sugita Y. Anharmonic Vibrational Analyses of Pentapeptide Conformations Explored with Enhanced Sampling Simulations. J Phys Chem B 2016; 120:10199-10213. [DOI: 10.1021/acs.jpcb.6b06672] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
| | | | - Shun-ichi Ishiuchi
- Laboratory
for Chemistry and Life Science, Institute for Innovative Research, Tokyo Institute of Technology, 4259 Nagatsuta-cho, Midori-ku, Yokohama 226-8503, Japan
| | - Masaaki Fujii
- Laboratory
for Chemistry and Life Science, Institute for Innovative Research, Tokyo Institute of Technology, 4259 Nagatsuta-cho, Midori-ku, Yokohama 226-8503, Japan
| | - Yuji Sugita
- RIKEN Advanced Institute for Computational Science, 7-1-26 Minatojima-Minamimachi,
Chuo-ku, Kobe, Hyogo 650-0047, Japan
- RIKEN Quantitative Biology Center, 1-6-5 Minatojima-Minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
| |
Collapse
|
41
|
König C, Christiansen O. Linear-scaling generation of potential energy surfaces using a double incremental expansion. J Chem Phys 2016. [DOI: 10.1063/1.4960189] [Citation(s) in RCA: 35] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Carolin König
- Department of Chemistry, Aarhus University, DK-8000 Aarhus C, Denmark
| | - Ove Christiansen
- Department of Chemistry, Aarhus University, DK-8000 Aarhus C, Denmark
| |
Collapse
|
42
|
Sibaev M, Crittenden DL. An efficient and numerically stable procedure for generating sextic force fields in normal mode coordinates. J Chem Phys 2016; 144:214107. [DOI: 10.1063/1.4953080] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- M. Sibaev
- Department of Chemistry, University of Canterbury, Christchurch, New Zealand
| | - D. L. Crittenden
- Department of Chemistry, University of Canterbury, Christchurch, New Zealand
| |
Collapse
|
43
|
Asami H, Tokugawa M, Masaki Y, Ishiuchi SI, Gloaguen E, Seio K, Saigusa H, Fujii M, Sekine M, Mons M. Effective Strategy for Conformer-Selective Detection of Short-Lived Excited State Species: Application to the IR Spectroscopy of the N1H Keto Tautomer of Guanine. J Phys Chem A 2016; 120:2179-84. [DOI: 10.1021/acs.jpca.6b01194] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Hiroya Asami
- Department
of Life Science, Graduate School of Bioscience and Biotechnology, Tokyo Institute of Technology, 4259-J2-12, Nagatsuta-cho Midori, Yokohama 226-8501, Japan
- LIDYL,
CEA, CNRS, Université Paris-Saclay, CEA Saclay, 91191 Gif-sur-Yvette, France
| | - Munefumi Tokugawa
- Department
of Life Science, Graduate School of Bioscience and Biotechnology, Tokyo Institute of Technology, 4259-J2-12, Nagatsuta-cho Midori, Yokohama 226-8501, Japan
| | - Yoshiaki Masaki
- Department
of Life Science, Graduate School of Bioscience and Biotechnology, Tokyo Institute of Technology, 4259-J2-12, Nagatsuta-cho Midori, Yokohama 226-8501, Japan
| | - Shun-ichi Ishiuchi
- Laboratory
for Chemistry and Life Science, Tokyo Institute of Technology, 4259-R1-15,
Nagatsuta-cho Midori, Yokohama 226-8503, Japan
| | - Eric Gloaguen
- LIDYL,
CEA, CNRS, Université Paris-Saclay, CEA Saclay, 91191 Gif-sur-Yvette, France
| | - Kohji Seio
- Department
of Life Science, Graduate School of Bioscience and Biotechnology, Tokyo Institute of Technology, 4259-J2-12, Nagatsuta-cho Midori, Yokohama 226-8501, Japan
| | - Hiroyuki Saigusa
- Graduate
School of Bio- and Nanosystem Science, Yokohama City University, Yokohama 236-0027, Japan
| | - Masaaki Fujii
- Laboratory
for Chemistry and Life Science, Tokyo Institute of Technology, 4259-R1-15,
Nagatsuta-cho Midori, Yokohama 226-8503, Japan
| | - Mitsuo Sekine
- Department
of Life Science, Graduate School of Bioscience and Biotechnology, Tokyo Institute of Technology, 4259-J2-12, Nagatsuta-cho Midori, Yokohama 226-8501, Japan
| | - Michel Mons
- LIDYL,
CEA, CNRS, Université Paris-Saclay, CEA Saclay, 91191 Gif-sur-Yvette, France
| |
Collapse
|
44
|
Ziegler B, Rauhut G. Efficient generation of sum-of-products representations of high-dimensional potential energy surfaces based on multimode expansions. J Chem Phys 2016; 144:114114. [DOI: 10.1063/1.4943985] [Citation(s) in RCA: 47] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] 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
|
45
|
Yagi K. Development of Molecular Vibrational Structure Theory with an Explicit Account of Anharmonicity. ACTA ACUST UNITED AC 2016. [DOI: 10.3175/molsci.10.a0085] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
46
|
Meier P, Oschetzki D, Pfeiffer F, Rauhut G. Towards an automated and efficient calculation of resonating vibrational states based on state-averaged multiconfigurational approaches. J Chem Phys 2015; 143:244111. [DOI: 10.1063/1.4938280] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Patrick Meier
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Dominik Oschetzki
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Florian Pfeiffer
- 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
|
47
|
Czarnecki MA, Morisawa Y, Futami Y, Ozaki Y. Advances in Molecular Structure and Interaction Studies Using Near-Infrared Spectroscopy. Chem Rev 2015; 115:9707-44. [DOI: 10.1021/cr500013u] [Citation(s) in RCA: 146] [Impact Index Per Article: 16.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
| | - Yusuke Morisawa
- Department
of Chemistry, School of Science and Engineering, Kinki University, Higashi-Osaka, Osaka 577-8502, Japan
| | - Yoshisuke Futami
- Department
of Biological and Chemical Systems Engineering, National Institute of Technology, Kumamoto College, Yatsushiro, Kumamoto 866-8501, Japan
| | - Yukihiro Ozaki
- Department
of Chemistry, School of Science and Technology, Kwansei Gakuin University, Sanda, Hyogo 669-1337, Japan
| |
Collapse
|
48
|
Oschetzki D, Rauhut G. Pushing the limits in accurate vibrational structure calculations: anharmonic frequencies of lithium fluoride clusters (LiF)n, n = 2-10. Phys Chem Chem Phys 2015; 16:16426-35. [PMID: 24981078 DOI: 10.1039/c4cp02264e] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The vibrational spectra of a series of small lithium fluoride clusters, i.e. (LiF)n, n = 2-10, were studied by vibrational configuration interaction (VCI) calculations relying on potential energy surfaces including three-mode coupling terms and being obtained from explicitly correlated local coupled cluster calculations. Due to the account for anharmonicity effects, the simulated spectra allow for a direct comparison with experimental data and may thus help to identify clusters in the experiments. Even structurally closely related clusters can clearly be distinguished by infrared spectroscopy. The largest system in this study required more than 1000 basis functions in the electronic structure calculations and more than 10(7) configurations in the vibrational structure calculations and became computationally feasible only due to a combination of different approximations and highly parallelized algorithms.
Collapse
Affiliation(s)
- Dominik Oschetzki
- Institute of Theoretical Chemistry, University of Stuttgart, Pfaffenwaldring 55, Stuttgart D-70569, Germany.
| | | |
Collapse
|
49
|
Ramakrishnan R, Rauhut G. Semi-quartic force fields retrieved from multi-mode expansions: Accuracy, scaling behavior, and approximations. J Chem Phys 2015; 142:154118. [DOI: 10.1063/1.4918587] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Raghunathan Ramakrishnan
- Institute of Physical Chemistry and National Center for Computational Design and Discovery of Novel Materials (MARVEL), Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| | - Guntram Rauhut
- Institute for Theoretical Chemistry, University of Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
| |
Collapse
|
50
|
König C, Christiansen O. Automatic determination of important mode–mode correlations in many-mode vibrational wave functions. J Chem Phys 2015; 142:144115. [DOI: 10.1063/1.4916518] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Carolin König
- Department of Chemistry, Aarhus University, DK-8000 Aarhus C, Denmark
| | - Ove Christiansen
- Department of Chemistry, Aarhus University, DK-8000 Aarhus C, Denmark
| |
Collapse
|