1
|
Schneider M, Rauhut G. Comparison of curvilinear coordinates within vibrational structure calculations based on automatically generated potential energy surfaces. J Chem Phys 2024; 161:094102. [PMID: 39225517 DOI: 10.1063/5.0225991] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2024] [Accepted: 08/12/2024] [Indexed: 09/04/2024] Open
Abstract
For floppy molecules showing internal rotations and/or large amplitude motions, curvilinear internal coordinates are known to be superior to rectilinear normal coordinates within vibrational structure calculations. Due to the myriad definitions of internal coordinates, automated and efficient potential energy surface generators necessitate a high degree of flexibility, supporting the properties arising from these coordinates. Within this work, an approach to deal with these challenges is presented, including key elements, such as the selection of appropriate fit functions, the exploitation of symmetry, the positioning of grid points, or elongation limits for different coordinates. These elements are tested for five definitions of curvilinear coordinates, with three of them being generated in an automated manner. Calculations for semi-rigid molecules, namely H2O, H2CO, CH2F2, and H2CNH, demonstrate the general functionality of the implemented algorithms. Additional calculations for the HOPO molecule highlight the benefits of these curvilinear coordinates for systems with large amplitude motions. This new implementation allowed us to compare the performance of these different coordinate systems with respect to the convergence of the underlying expansion of the potential energy surface and subsequent vibrational configuration interaction calculations.
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
|
Mendolicchio M, Barone V. Accurate Vibrational and Ro-Vibrational Contributions to the Properties of Large Molecules by a New Engine Employing Curvilinear Internal Coordinates and Vibrational Perturbation Theory to Second Order. J Chem Theory Comput 2024. [PMID: 39215708 DOI: 10.1021/acs.jctc.4c00857] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/04/2024]
Abstract
The unbiased comparison between theory and experiment requires approaches more sophisticated than the basic harmonic-oscillator rigid-rotor model, for taking into account vibrational averaging effects and ro-vibrational couplings in molecules of increasing size. Second-order vibrational perturbation theory based on curvilinear internal coordinates (ICs) offers a remarkable compromise between accuracy and computational cost, thanks to the reduction of mode-mode couplings with respect to their counterparts based on Cartesian coordinates. Therefore, we have developed, implemented, and validated a general engine employing ICs, which allows the accurate evaluation of vibrational averages and ro-vibrational couplings for molecules containing up to about 50 atoms beyond the harmonic approximation. After validation of the new tool for relatively small molecules, the effectiveness of ICs has been demonstrated for some flexible and/or quite large molecular bricks of life.
Collapse
|
3
|
Sunaga A, Avila G, Mátyus E. Variational Vibrational States of Methanol (12D). J Chem Theory Comput 2024. [PMID: 39213609 DOI: 10.1021/acs.jctc.4c00647] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/04/2024]
Abstract
Full-dimensional (12D) vibrational states of the methanol molecule (CH3OH) have been computed using the GENIUSH-Smolyak approach and the potential energy surface from Qu and Bowman (2013). All vibrational energies are converged better than 0.5 cm-1 with respect to the basis and grid size up to the first overtone of the CO stretch, ca. 2000 cm-1 beyond the zero-point vibrational energy. About 70 torsion-vibration states are reported and assigned. The computed vibrational energies agree with the available experimental data within less than a few cm-1 in most cases, which confirms the good accuracy of the potential energy surface. The computations are carried out using curvilinear normal coordinates with the option of path-following coefficients, which minimize the coupling of the small- and large-amplitude motions. It is important to ensure tight numerical fulfillment of the C3v(M) molecular symmetry for every geometry and coefficient set used to define the curvilinear normal coordinates along the torsional coordinate to obtain a faithful description of degeneracy in this floppy system. The reported values may provide a computational reference for fundamental spectroscopy, astrochemistry, and for the search of the proton-to-electron mass ratio variation using the methanol molecule.
Collapse
Affiliation(s)
- Ayaki Sunaga
- ELTE, Eötvös Loránd University, Institute of Chemistry, Budapest 1117, Hungary
| | - Gustavo Avila
- ELTE, Eötvös Loránd University, Institute of Chemistry, Budapest 1117, Hungary
| | - Edit Mátyus
- ELTE, Eötvös Loránd University, Institute of Chemistry, Budapest 1117, Hungary
| |
Collapse
|
4
|
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
|
5
|
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
|
6
|
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
|
7
|
Bader F, Lauvergnat D, Christiansen O. Vibrationally correlated calculations in polyspherical coordinates: Taylor expansion-based kinetic energy operators. J Chem Phys 2023; 159:214107. [PMID: 38047511 DOI: 10.1063/5.0171912] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/11/2023] [Accepted: 11/06/2023] [Indexed: 12/05/2023] Open
Abstract
The efficiency of quantum chemical simulations of nuclear motion can in many cases greatly benefit from the application of curvilinear coordinate systems. This is rooted in the fact that a set of smartly selected curvilinear coordinates may represent the motion naturally well, thus decreasing the couplings between motions in these coordinates. In this study, we assess the validity of different Taylor expansion-based approximations of kinetic energy operators in a (curvilinear) polyspherical parametrization. To this end, we investigate the accuracy as well as the numerical performance of the approximations in time-independent vibrational coupled cluster and full vibrational interaction calculations for several test cases ranging from tri- to penta-atomic molecules. We find that several of the proposed schemes reproduce the vibrational ground state and excitation energies to a decent accuracy, justifying their application in future investigations. Furthermore, due to the restricted mode coupling and their inherent sum-of-products form, the new approximations open up the possibility of treating large molecular systems with efficient vibrational coupled cluster schemes in general coordinates.
Collapse
Affiliation(s)
- F Bader
- Department of Chemistry, Aarhus University, DK-8000 Aarhus C, Denmark
| | - D Lauvergnat
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, Orsay 91405, France
| | - O Christiansen
- Department of Chemistry, Aarhus University, DK-8000 Aarhus C, Denmark
| |
Collapse
|
8
|
Lauvergnat D, Nauts A. Smolyak Scheme for solving the Schrödinger equation: Application to Malonaldehyde in Full Dimensionality. Chemphyschem 2023; 24:e202300501. [PMID: 37555577 DOI: 10.1002/cphc.202300501] [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: 07/14/2023] [Revised: 08/08/2023] [Accepted: 08/09/2023] [Indexed: 08/10/2023]
Abstract
In 1963 Smolyak introduced an approach to overcome the exponential scaling with respect to the number of variables of the direct product size [S. A. Smolyak Soviet Mathematics Doklady, 4, 240 (1963)]. The main idea is to replace a single large direct product by a sum of selected small direct products. It was first used in quantum dynamics in 2009 by Avila and Carrington [G. Avila and T. Carrington, J. Chem. Phys., 131, 174103 (2009)]. Since then, several calculations have been published by Avila and Carrington and by other groups. In the present study, and to push the limit to larger and more complex systems, this scheme is combined with the use of an on-the-fly calculation of the kinetic energy operator and a Block-Davidson procedure to obtain eigenstates in our home-made Fortran codes, ElVibRot and Tnum-Tana. This was applied to compute the tunneling splitting of malonaldehyde in full dimensionality (21D) using the potential of Mizukami et al. [W. Mizukami, S. Habershon, and D.P. Tew, J. Chem. Phys. 141, 1443-10 (2014)]. Our tunneling splitting calculations, 21.7±0.3 cm-1 and 2.9±0.1 cm-1 , show excellent agreement with the experimental values, 21.6 cm-1 and 2.9 cm-1 for the normal isotopologue and the mono-deuterated one, respectively.
Collapse
Affiliation(s)
- David Lauvergnat
- Institut de Chimie Physique, CNRS, Univ. Paris-Sud, Université Paris-Saclay, 91405, Orsay, France
| | - André Nauts
- Institut de Chimie Physique, CNRS, Univ. Paris-Sud, Université Paris-Saclay, 91405, Orsay, France
- Institute of Condensed Matter and Nanosciences (NAPS), Université Catholique de Louvain, 2 Chemin du Cyclotron, 1348, Louvain-la-Neuve, Belgium
| |
Collapse
|
9
|
Manzhos S, Ihara M, Carrington T. Using Collocation to Solve the Schrödinger Equation. J Chem Theory Comput 2023; 19:1641-1656. [PMID: 36974479 DOI: 10.1021/acs.jctc.2c01232] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/10/2023]
Abstract
We review the collocation approach to the solution of the Schrödinger equation and its uses in applications. Interrelations between collocation and other methods are highlighted. We also stress advantages and disadvantages of the rectangular collocation formulation. Using collocation makes it possible to use any, e.g. optimized, coordinates and basis functions, including nonintegrable basis functions, and provides a straightforward way of dealing with singularities in the potential. In addition, we stress that using collocation facilitates tuning the shape of basis functions and the placement of points, both of which can be done with machine-learning methods. Applications to electronic and vibrational problems are reviewed focusing on calculations for molecules on surfaces for which there are few variational calculations. Collocation has advantages when potential energy surfaces are unavailable, in particular, for molecule-surface systems, and for systems for which standard direct product quadrature grids, often used with variational methods, are costly.
Collapse
Affiliation(s)
- Sergei Manzhos
- School of Materials and Chemical Technology, Tokyo Institute of Technology, Ookayama 2-12-1, Meguro-ku, Tokyo 152-8552, Japan
| | - Manabu Ihara
- School of Materials and Chemical Technology, Tokyo Institute of Technology, Ookayama 2-12-1, Meguro-ku, Tokyo 152-8552, Japan
| | - Tucker Carrington
- Department of Chemistry, Queen’s University, 90 Bader Lane, Kingston, Ontario K7L 3N6, Canada
| |
Collapse
|
10
|
Mátyus E, Martín Santa Daría A, Avila G. Exact quantum dynamics developments for floppy molecular systems and complexes. Chem Commun (Camb) 2023; 59:366-381. [PMID: 36519578 DOI: 10.1039/d2cc05123k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
Abstract
Molecular rotation, vibration, internal rotation, isomerization, tunneling, intermolecular dynamics of weakly and strongly interacting systems, intra-to-inter-molecular energy transfer, hindered rotation and hindered translation over surfaces are important types of molecular motions. Their fundamentally correct and detailed description can be obtained by solving the nuclear Schrödinger equation on a potential energy surface. Many of the chemically interesting processes involve quantum nuclear motions which are 'delocalized' over multiple potential energy wells. These 'large-amplitude' motions in addition to the high dimensionality of the vibrational problem represent challenges to the current (ro)vibrational methodology. A review of the quantum nuclear motion methodology is provided, current bottlenecks of solving the nuclear Schrödinger equation are identified, and solution strategies are reviewed. Technical details, computational results, and analysis of these results in terms of limiting models and spectroscopically relevant concepts are highlighted for selected numerical examples.
Collapse
Affiliation(s)
- Edit Mátyus
- ELTE, Eötvös Loránd University, Institute of Chemistry, Pázmány Péter sétány 1/A, 1117 Budapest, Hungary.
| | - Alberto Martín Santa Daría
- ELTE, Eötvös Loránd University, Institute of Chemistry, Pázmány Péter sétány 1/A, 1117 Budapest, Hungary.
| | - Gustavo Avila
- ELTE, Eötvös Loránd University, Institute of Chemistry, Pázmány Péter sétány 1/A, 1117 Budapest, Hungary.
| |
Collapse
|
11
|
Mendolicchio M, Bloino J, Barone V. Perturb-Then-Diagonalize Vibrational Engine Exploiting Curvilinear Internal Coordinates. J Chem Theory Comput 2022; 18:7603-7619. [PMID: 36322968 DOI: 10.1021/acs.jctc.2c00773] [Citation(s) in RCA: 15] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
The present paper is devoted to the implementation and validation of a second-order perturbative approach to anharmonic vibrations, followed by variational treatment of strong couplings (GVPT2) based on curvilinear internal coordinates. The main difference with respect to the customary Cartesian-based formulation is that the kinetic energy operator is no longer diagonal, and has to be expanded as well, leading to additional terms which have to be taken into proper account. It is, however, possible to recast all the equations as well-defined generalizations of the corresponding Cartesian-based counterparts, thus achieving a remarkable simplification of the new implementation. Particular attention is paid to the treatment of Fermi resonances with significant number of test cases analyzed fully, validating the new implementation. The results obtained in this work confirm that curvilinear coordinates strongly reduce the strength of inter-mode couplings compared to their Cartesian counterparts. This increases the reliability of low-order perturbative treatments for semi-rigid molecules and paves the way toward the reliable representation of more flexible molecules where small- and large-amplitude motions can be safely decoupled and treated at different levels of theory.
Collapse
Affiliation(s)
- Marco Mendolicchio
- Scuola Superiore Meridionale, Largo S. Marcellino 10, Napoli I-80138, Italy
| | - Julien Bloino
- Scuola Normale Superiore, Piazza dei Cavalieri 7, Pisa I-56126, Italy
| | - Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, Pisa I-56126, Italy
| |
Collapse
|
12
|
Papp D, Tajti V, Avila G, Mátyus E, Czakó G. CH 4·F − revisited: full-dimensional ab initio potential energy surface and variational vibrational states. Mol Phys 2022. [DOI: 10.1080/00268976.2022.2113565] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/15/2022]
Affiliation(s)
- Dóra Papp
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group, Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged, Szeged, Hungary
| | - Viktor Tajti
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group, Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged, Szeged, Hungary
| | - Gustavo Avila
- ELTE, Eötvös Loránd University, Institute of Chemistry, Budapest, Hungary
| | - Edit Mátyus
- ELTE, Eötvös Loránd University, Institute of Chemistry, Budapest, Hungary
| | - Gábor Czakó
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group, Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged, Szeged, Hungary
| |
Collapse
|
13
|
Aerts A, Brown A, Gatti F. Intramolecular Vibrational Redistribution in Formic Acid and its Deuterated Forms. J Chem Phys 2022; 157:014306. [DOI: 10.1063/5.0098819] [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 intramolecular vibrational relaxation dynamics of formic acid and its deuterated isotopologues is simulated on the full-dimensional potential energy surface of Richter and Carbonnière [F. Richter and P. Carbonnière, J. Chem. Phys. 148, 064303 (2018)] using the Heidelberg MCTDH package. We focus on couplings with the torsion vibrational modes close to the trans- cis isomerisation coordinate from the dynamics of artificially excited vibrational mode overtones. The C-O stretch bright vibrational mode is coupled to the out-of-the plane torsion mode in HCOOH, where this coupling could be exploited for laser-induced trans-to- cis isomerisation. Strong isotopic effects are observed: deuteration of the hydroxyl group, i.e., in HCOOD and DCOOD, destroys the C-O stretch to torsion mode coupling whereas in DCOOH, little to no effect is observed.
Collapse
Affiliation(s)
| | - Alex Brown
- Department of Chemistry, University of Alberta, Canada
| | - Fabien Gatti
- Institut des Sciences Moléculaires d'Orsay, Paris-Saclay University, France
| |
Collapse
|
14
|
Marsili E, Agostini F, Nauts A, Lauvergnat D. Quantum dynamics with curvilinear coordinates: models and kinetic energy operator. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2022; 380:20200388. [PMID: 35341305 DOI: 10.1098/rsta.2020.0388] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/29/2021] [Accepted: 11/01/2021] [Indexed: 06/14/2023]
Abstract
In order to simplify the numerical solution of the time-dependent or time-independent Schrödinger equations associated with atomic and molecular motions, the use of well-adapted coordinates is essential. Usually, this set of curvilinear coordinates leads to a Hamiltonian operator that is as separable as possible. Although their corresponding kinetic energy operator (KEO) expressions can be derived analytically for small systems or special kinds of coordinates, a numerical and exact approach allows one to compute them in terms of sophisticated curvilinear coordinates. Furthermore, the numerical approach enables one to easily define reduced-dimensionality or constrained models. We present here a recent implementation of this numerical approach that allows nested coordinate transformations, therefore leading to great flexibility in the definition of the curvilinear coordinates. Furthermore, this implementation has no limitations in terms of numbers of atoms or coordinate transformations. The quantum dynamics of the cis-trans photoisomerization of part of the retinal chromophore illustrates the construction of the coordinates and KEO part of a three-dimensional model. This article is part of the theme issue 'Chemistry without the Born-Oppenheimer approximation'.
Collapse
Affiliation(s)
- Emanuele Marsili
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, Orsay 91405, France
- Department of Chemistry, Durham University, Durham DH1 3LE, UK
| | - Federica Agostini
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, Orsay 91405, France
| | - André Nauts
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, Orsay 91405, France
- Institute of Condensed Matter and Nanosciences (NAPS), Université Catholique de Louvain, BE-1348 Louvain-la-Neuve, Belgium
| | - David Lauvergnat
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, Orsay 91405, France
| |
Collapse
|
15
|
Martín Santa Daría A, Avila G, Mátyus E. Performance of a black-box-type rovibrational method in comparison with a tailor-made approach: Case study for the methane-water dimer. J Chem Phys 2021; 154:224302. [PMID: 34241197 DOI: 10.1063/5.0054512] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The present work intends to join and respond to the excellent and thoroughly documented rovibrational study of X. G. Wang and T. Carrington, Jr. [J. Chem. Phys. 154, 124112 (2021)] that used an approach tailored for floppy dimers with an analytic dimer Hamiltonian and a non-product basis set including Wigner D functions. It is shown in the present work that the GENIUSH black-box-type rovibrational method can approach the performance of the tailor-made computation for the example of the floppy methane-water dimer. Rovibrational transition energies and intensities are obtained in the black-box-type computation with a twice as large basis set and in excellent numerical agreement in comparison with the more efficient tailor-made approach.
Collapse
Affiliation(s)
| | - Gustavo Avila
- Institute of Chemistry, ELTE Eötvös Loránd University, Pázmány Péter sétány 1/A, 1117 Budapest, Hungary
| | - Edit Mátyus
- Institute of Chemistry, ELTE Eötvös Loránd University, Pázmány Péter sétány 1/A, 1117 Budapest, Hungary
| |
Collapse
|
16
|
Paoloni L, Maris A. Interplay of Rotational and Pseudorotational Motions in Flexible Cyclic Molecules. J Phys Chem A 2021; 125:4098-4113. [PMID: 33973473 PMCID: PMC8279653 DOI: 10.1021/acs.jpca.1c01472] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/17/2021] [Revised: 04/23/2021] [Indexed: 11/29/2022]
Abstract
Solutions to the time-independent nuclear Schrödinger equation associated with the pseudorotational motion of three flexible cyclic molecules are presented and discussed. Structural relaxations related to the pseudorotational motion are described as functions of a pseudorotation angle ϕ which is formulated according to the definition of ring-puckering coordinates originally proposed by Cremer and Pople ( J. Am. Chem. Soc. 1975, 97 (6), 1354-1358). In order to take into account the interplay between pseudorotational and rotational motions, the rovibrational Hamiltonian matrices are formulated for the rotational quantum numbers J = 0 and J = 1. The rovibrational Hamiltonian matrices are constructed and diagonalized using a Python program developed by the authors. Suitable algorithms for (i) the construction of one-dimensional cuts of potential energy surfaces along the pseudorotation angle ϕ and (ii) the assignment of the vibrorotational wave functions (which are needed for the automatic calculation of rotational transition energies J = 0 → J = 1) are described and discussed.
Collapse
Affiliation(s)
- Lorenzo Paoloni
- Dipartimento
di Fisica e Astronomia, Università
di Padova, via Marzolo 8, I-35131 Padova, Italy
| | - Assimo Maris
- Dipartimento
di Chimica G. Ciamician, Università
di Bologna, via Selmi 2, I-40126 Bologna, Italy
| |
Collapse
|
17
|
Sarka J, Poirier B, Szalay V, Császár AG. On neglecting Coriolis and related couplings in first-principles rovibrational spectroscopy: Considerations of symmetry, accuracy, and simplicity. II. Case studies for H 2O isotopologues, H 3+, O 3, and NH 3. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2021; 250:119164. [PMID: 33493950 DOI: 10.1016/j.saa.2020.119164] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/05/2020] [Accepted: 10/26/2020] [Indexed: 06/12/2023]
Abstract
For centuries, it has been known that vibrational and rotational degrees of freedom are in general not separable. Nevertheless, surprisingly little is known about the best strategies for approximately separating these degrees of freedom in practice-even in the case of semirigid molecules, where the separation is most meaningful. There is also some confusion in the literature about the proper way to quantify the magnitude of the Coriolis (i.e., rotation-vibration) coupling in rovibrational Hamiltonians or its effect on the rovibrational eigenenergies. In this study, a vibrational-coordinate-independent metric is proposed to quantify the magnitude of the Coriolis contribution to the rovibrational Hamiltonian. The impact of Coriolis coupling on the rovibrational eigenenergies is computed numerically exactly, using both full and various truncated Hamiltonians. The role played by the choice of the vibrational coordinate system-and especially by the choice of "embedding" or body-fixed frame-is examined extensively, both numerically and analytically. This investigation targets several molecular prototypes, all of which serve as important benchmarks for the high-resolution spectroscopic community. Most of these are triatomic molecules, including water (H216O), its deuterated isotopologues (D216O and HD16O), H3+, and ozone (16O3), but the tetratomic ammonia molecule (14NH3) is also investigated. These studies provide important insight into the nature of Coriolis coupling under various circumstances. The findings of this study also have significant practical ramifications, vis-à-vis the use of simplifying numerical approximation techniques in nuclear-motion computations.
Collapse
Affiliation(s)
- János Sarka
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, TX 79409, USA.
| | - Bill Poirier
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, TX 79409, USA.
| | - Viktor Szalay
- Institute for Solid State Physics and Optics, Wigner Research Centre for Physics, Hungarian Academy of Sciences, P.O. Box 49, H-1525 Budapest, Hungary
| | - Attila G Császár
- MTA-ELTE Complex Chemical Systems Research Group and Laboratory of Molecular Structure and Dynamics, Institute of Chemistry, ELTE Eötvös Loránd University, H-1117 Budapest, Pázmány Péter sétány 1/A, Hungary.
| |
Collapse
|
18
|
Martín Santa Daría A, Avila G, Mátyus E. Fingerprint region of the formic acid dimer: variational vibrational computations in curvilinear coordinates. Phys Chem Chem Phys 2021; 23:6526-6535. [PMID: 33690754 DOI: 10.1039/d0cp06289h] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Curvilinear kinetic energy models are developed for variational nuclear motion computations including the inter- and the low-frequency intra-molecular degrees of freedom of the formic acid dimer. The coupling of the inter- and intra-molecular modes is studied by solving the vibrational Schrödinger equation for a series of vibrational models, from two up to ten active vibrational degrees of freedom by selecting various combinations of active modes and constrained coordinate values. Vibrational states, nodal assignment, and infrared vibrational intensity information is computed using the full-dimensional potential energy surface (PES) and electric dipole moment surface developed by Qu and Bowman [Phys. Chem. Chem. Phys., 2016, 18, 24835; J. Chem. Phys., 2018, 148, 241713]. Good results are obtained for several fundamental and combination bands in comparison with jet-cooled vibrational spectroscopy experiments, but the description of the ν8 and ν9 fundamental vibrations, which are close in energy and have the same symmetry, appears to be problematic. For further progress in comparison with experiment, the potential energy surface, and in particular, its multi-dimensional couplings representation, requires further improvement.
Collapse
Affiliation(s)
- Alberto Martín Santa Daría
- Institute of Chemistry, ELTE Eötvös Loránd University, Pázmány Péter sétány 1/A, 1117 Budapest, Hungary.
| | | | | |
Collapse
|
19
|
Carrington T. Using collocation to study the vibrational dynamics of molecules. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2021; 248:119158. [PMID: 33218875 DOI: 10.1016/j.saa.2020.119158] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/23/2020] [Accepted: 10/26/2020] [Indexed: 06/11/2023]
Abstract
In this paper, I review collocation methods for solving the time-independent and the time-dependent Schroedinger equation. Unlike traditional variational methods, collocation methods do not require integrals and quadrature. Either collocation or quadrature is necessary if the potential does not have a special form. If the basis is a direct product of univariate bases and the quadrature grid is also a direct product, there exist variational methods that do not require quadrature approximations for potential energy matrix elements. These methods, however, do require storing, in computer memory, vectors with as many components as there are quadrature points. For this reason direct-product variational methods are poor for problems with more than five atoms. There are well established ideas for reducing the size of the basis in a variational calculation. Three such ideas are: 1) prune the direct product basis; 2) use basis functions that are products of multivariate functions; 3) optimise the basis functions (e.g. Multiconfiguration time-dependent Hartree). Reducing the basis size, however, is not enough to the make variational methods tractable because, for all three of these ideas, quadrature rears its ugly head. Collocation is an attractive alternative to variational methods.
Collapse
Affiliation(s)
- Tucker Carrington
- Department of Chemistry, Queen's University, Kingston, Ontario K7L 3N6, Canada.
| |
Collapse
|
20
|
Panadés-Barrueta RL, Peláez D. Low-rank sum-of-products finite-basis-representation (SOP-FBR) of potential energy surfaces. J Chem Phys 2020; 153:234110. [PMID: 33353311 DOI: 10.1063/5.0027143] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023] Open
Abstract
The sum-of-products finite-basis-representation (SOP-FBR) approach for the automated multidimensional fit of potential energy surfaces (PESs) is presented. In its current implementation, the method yields a PES in the so-called Tucker sum-of-products form, but it is not restricted to this specific ansatz. The novelty of our algorithm lies in the fact that the fit is performed in terms of a direct product of a Schmidt basis, also known as natural potentials. These encode in a non-trivial way all the physics of the problem and, hence, circumvent the usual extra ad hoc and a posteriori adjustments (e.g., damping functions) of the fitted PES. Moreover, we avoid the intermediate refitting stage common to other tensor-decomposition methods, typically used in the context of nuclear quantum dynamics. The resulting SOP-FBR PES is analytical and differentiable ad infinitum. Our ansatz is fully general and can be used in combination with most (molecular) dynamics codes. In particular, it has been interfaced and extensively tested with the Heidelberg implementation of the multiconfiguration time-dependent Hartree quantum dynamical software package.
Collapse
Affiliation(s)
- Ramón L Panadés-Barrueta
- Laboratoire de Physique des Lasers, Atomes et Molécules (PhLAM), Université Lille 1, Villeneuve d'Ascq Cedex, France
| | - Daniel Peláez
- Institut des Sciences Moléculaires d'Orsay (ISMO) - UMR 8214, Bât. 520, Université Paris-Saclay, 91405 Orsay Cedex, France
| |
Collapse
|
21
|
Di Maiolo F, Brey D, Binder R, Burghardt I. Quantum dynamical simulations of intra-chain exciton diffusion in an oligo (para-phenylene vinylene) chain at finite temperature. J Chem Phys 2020; 153:184107. [PMID: 33187420 DOI: 10.1063/5.0027588] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022] Open
Abstract
We report on quantum dynamical simulations of exciton diffusion in an oligo(para-phenylene vinylene) chain segment with 20 repeat units (OPV-20) at finite temperature, complementary to our recent study of the same system at T = 0 K [R. Binder and I. Burghardt, J. Chem. Phys. 152, 204120 (2020)]. Accurate quantum dynamical simulations are performed using the multi-layer multi-configuration time-dependent Hartree method as applied to a site-based Hamiltonian comprising 20 electronic states of Frenkel type and 460 vibrational modes, including site-local quinoid-distortion modes along with site-correlated bond-length alternation (BLA) modes, ring torsional modes, and an explicit harmonic-oscillator bath. A first-principles parameterized Frenkel-Holstein type Hamiltonian is employed, which accounts for correlations between the ring torsional modes and the anharmonically coupled BLA coordinates located at the same junction. Thermally induced fluctuations of the torsional modes are described by a stochastic mean-field approach, and their impact on the excitonic motion is characterized in terms of the exciton mean-squared displacement. A normal diffusion regime is observed under periodic boundary conditions, apart from transient localization features. Even though the polaronic exciton species are comparatively weakly bound, exciton diffusion is found to be a coherent-rather than hopping type-process, driven by the fluctuations of the soft torsional modes. Similar to the previous observations for oligothiophenes, the evolution for the most part exhibits a near-adiabatic dynamics of local exciton ground states (LEGSs) that adjust to the local conformational dynamics. However, a second mechanism, involving resonant transitions between neighboring LEGSs, gains importance at higher temperatures.
Collapse
Affiliation(s)
- Francesco Di Maiolo
- Institute for Physical and Theoretical Chemistry, Goethe University, Max-von-Laue-Str. 7, 60438 Frankfurt/Main, Germany
| | - Dominik Brey
- Institute for Physical and Theoretical Chemistry, Goethe University, Max-von-Laue-Str. 7, 60438 Frankfurt/Main, Germany
| | - Robert Binder
- Institute for Physical and Theoretical Chemistry, Goethe University, Max-von-Laue-Str. 7, 60438 Frankfurt/Main, Germany
| | - Irene Burghardt
- Institute for Physical and Theoretical Chemistry, Goethe University, Max-von-Laue-Str. 7, 60438 Frankfurt/Main, Germany
| |
Collapse
|
22
|
Marsili E, Olivucci M, Lauvergnat D, Agostini F. Quantum and Quantum-Classical Studies of the Photoisomerization of a Retinal Chromophore Model. J Chem Theory Comput 2020; 16:6032-6048. [PMID: 32931266 DOI: 10.1021/acs.jctc.0c00679] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
We report an in-depth analysis of the photo-induced isomerization of the 2-cis-penta-2,4-dieniminium cation: a minimal model of the 11-cis retinal protonated Schiff base chromophore of the dim-light photoreceptor rhodopsin. Based on recently developed three-dimensional potentials parametrized on ab initio multi-state multi-configurational second-order perturbation theory data, we perform quantum-dynamical studies. In addition, simulations based on various quantum-classical methods, among which Tully surface hopping and the coupled-trajectory approach derived from the exact factorization, allow us to validate their performance against vibronic wavepacket propagation and, therefore, a purely quantum treatment. Quantum-dynamics results uncover qualitative differences with respect to the two-dimensional Hahn-Stock potentials, widely used as model potentials for the isomerization of the same chromophore, due to the increased dimensionality and three-mode correlation. Quantum-classical simulations show, instead, that three-dimensional model potentials are capable of capturing a number of features revealed by atomistic simulations and experimental observations. In particular, a recently reported vibrational phase relationship between double-bond torsion and hydrogen-out-of-plane modes critical for rhodopsin isomerization efficiency is correctly reproduced.
Collapse
Affiliation(s)
- Emanuele Marsili
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, Orsay 91405, France.,Department of Chemistry, Durham University, South Road, Durham DH1 3LE, United Kingdom
| | - Massimo Olivucci
- Department of Biotechnology, Chemistry and Pharmacy, Università degli Studi di Siena, Via A. Moro 2, I-53100 Siena, Italy.,Department of Chemistry, Bowling Green State University, Bowling Green, Ohio 43403, United States
| | - David Lauvergnat
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, Orsay 91405, France
| | - Federica Agostini
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, Orsay 91405, France
| |
Collapse
|
23
|
Hegger R, Binder R, Burghardt I. First-Principles Quantum and Quantum-Classical Simulations of Exciton Diffusion in Semiconducting Polymer Chains at Finite Temperature. J Chem Theory Comput 2020; 16:5441-5455. [PMID: 32786907 DOI: 10.1021/acs.jctc.0c00351] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
We report on first-principles quantum-dynamical and quantum-classical simulations of photoinduced exciton dynamics in oligothiophene chain segments, representative of intrachain exciton migration in the poly(3-hexylthiophene) (P3HT) polymer. Following up on our recent study (Binder R.; Burghardt, I. Faraday Discuss. 2020, 221, 406), multilayer multiconfiguration time-dependent Hartree calculations for a short oligothiophene segment comprising 20 monomer units (OT-20) are carried out to obtain full quantum-dynamical simulations at finite temperature. These are employed to benchmark mean-field Ehrenfest calculations, which are shown to give qualitatively correct results for the present system. Periodic boundary conditions turn out to significantly improve earlier estimates of diffusion coefficients. Using the Ehrenfest approach, a series of calculations are subsequently carried out for larger lattices (OT-40 to OT-80), leading to estimates for temperature-dependent mean-squared displacements, which are found to exhibit a near-linear dependence as a function of time. The resulting diffusion coefficient estimates are an increasing function of temperature, whose detailed functional form depends on the degree of static disorder. With a realistic static disorder parameter (σs ≃ 0.06 eV), the diffusion coefficients decrease from D ∼ 1 × 10-2 cm2 s-1 to D ∼ 1 × 10-3 cm2 s-1, in qualitative agreement with experimental data for P3HT. The dynamical scenario obtained from our simulations shows that exciton migration in P3HT-type chains is a largely adiabatic process throughout the temperature regime we investigated (i.e., T = 50-300 K). The resulting picture of exciton migration is a coherent, but not bandlike, motion of an exciton-polaron driven by fluctuations induced by low-frequency modes. This process acquires partial hopping character if static disorder becomes prominent and Anderson localization sets in.
Collapse
Affiliation(s)
- Rainer Hegger
- Institute for Physical and Theoretical Chemistry, Goethe University, Max-von-Laue-Strasse 7, 60438 Frankfurt/Main, Germany
| | - Robert Binder
- Institute for Physical and Theoretical Chemistry, Goethe University, Max-von-Laue-Strasse 7, 60438 Frankfurt/Main, Germany
| | - Irene Burghardt
- Institute for Physical and Theoretical Chemistry, Goethe University, Max-von-Laue-Strasse 7, 60438 Frankfurt/Main, Germany
| |
Collapse
|
24
|
Binder R, Bonfanti M, Lauvergnat D, Burghardt I. First-principles description of intra-chain exciton migration in an oligo(para-phenylene vinylene) chain. I. Generalized Frenkel-Holstein Hamiltonian. J Chem Phys 2020; 152:204119. [PMID: 32486686 DOI: 10.1063/5.0004510] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Abstract
A generalized Frenkel-Holstein Hamiltonian is constructed to describe exciton migration in oligo(para-phenylene vinylene) chains, based on excited state electronic structure data for an oligomer comprising 20 monomer units (OPV-20). Time-dependent density functional theory calculations using the ωB97XD hybrid functional are employed in conjunction with a transition density analysis to study the low-lying singlet excitations and demonstrate that these can be characterized to a good approximation as a Frenkel exciton manifold. Based on these findings, we employ the analytic mapping procedure of Binder et al. [J. Chem. Phys. 141, 014101 (2014)] to translate one-dimensional (1D) and two-dimensional (2D) potential energy surface (PES) scans to a fully anharmonic, generalized Frenkel-Holstein (FH) Hamiltonian. A 1D PES scan is carried out for intra-ring quinoid distortion modes, while 2D PES scans are performed for the anharmonically coupled inter-monomer torsional and vinylene bridge bond length alternation modes. The kinetic energy is constructed in curvilinear coordinates by an exact numerical procedure, using the TNUM Fortran code. As a result, a fully molecular-based, generalized FH Hamiltonian is obtained, which is subsequently employed for quantum exciton dynamics simulations, as shown in Paper II [R. Binder and I. Burghardt, J. Chem. Phys. 152, 204120 (2020)].
Collapse
Affiliation(s)
- Robert Binder
- Institute for Physical and Theoretical Chemistry, Goethe University, Max-von-Laue-Str. 7, 60438 Frankfurt/Main, Germany
| | - Matteo Bonfanti
- Institute for Physical and Theoretical Chemistry, Goethe University, Max-von-Laue-Str. 7, 60438 Frankfurt/Main, Germany
| | - David Lauvergnat
- Université Paris-Saclay, CNRS, Institut de Chimie Physique, UMR 8000, 91405 Orsay, France
| | - Irene Burghardt
- Institute for Physical and Theoretical Chemistry, Goethe University, Max-von-Laue-Str. 7, 60438 Frankfurt/Main, Germany
| |
Collapse
|
25
|
Binder R, Burghardt I. First-principles description of intra-chain exciton migration in an oligo(para-phenylene vinylene) chain. II. ML-MCTDH simulations of exciton dynamics at a torsional defect. J Chem Phys 2020; 152:204120. [DOI: 10.1063/5.0004511] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023] Open
Affiliation(s)
- Robert Binder
- Institute for Physical and Theoretical Chemistry, Goethe University, Max-von-Laue-Str. 7, 60438 Frankfurt/Main, Germany
| | - Irene Burghardt
- Institute for Physical and Theoretical Chemistry, Goethe University, Max-von-Laue-Str. 7, 60438 Frankfurt/Main, Germany
| |
Collapse
|
26
|
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
|
27
|
Sarka J, Poirier B, Szalay V, Császár AG. On neglecting Coriolis and related couplings in first-principles rovibrational spectroscopy: considerations of symmetry, accuracy, and simplicity. Sci Rep 2020; 10:4872. [PMID: 32184431 PMCID: PMC7078231 DOI: 10.1038/s41598-020-60971-x] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/09/2019] [Accepted: 02/19/2020] [Indexed: 11/20/2022] Open
Abstract
The rotation-vibration (Coriolis) coupling contribution to variationally computed rovibrational energy levels is investigated, employing triatomic AB[Formula: see text] molecules as models. In particular, calculations are performed for H[Formula: see text][Formula: see text]O, across a range of vibrational and rotational excitations, both with and without the Coriolis contribution. A variety of different embedding choices are considered, together with a hierarchy of increasingly severe approximations culminating in a generalized version of the so-called "centrifugal sudden" method. Several surprising and remarkable conclusions are found, including that the Eckart embedding is not the best embedding choice.
Collapse
Affiliation(s)
- János Sarka
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas, 79409, USA
| | - Bill Poirier
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas, 79409, USA.
| | - Viktor Szalay
- Institute for Solid State Physics and Optics, Wigner Research Centre for Physics, P.O. Box 49, H-1525, Budapest, Hungary
| | - Attila G Császár
- MTA-ELTE Complex Chemical Systems Research Group and Laboratory of Molecular Structure and Dynamics, Institute of Chemistry, ELTE Eötvös Loránd University, H-1117, Budapest, Pázmány Péter sétány 1/A, Hungary.
| |
Collapse
|
28
|
Binder R, Burghardt I. First-principles quantum simulations of exciton diffusion on a minimal oligothiophene chain at finite temperature. Faraday Discuss 2020; 221:406-427. [PMID: 31596291 DOI: 10.1039/c9fd00066f] [Citation(s) in RCA: 24] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023]
Abstract
High-dimensional multiconfigurational quantum dynamics simulations are carried out at finite temperature to simulate exciton diffusion on an oligothiophene chain, representative of a segment of the poly(3-hexylthiophene) (P3HT) polymer. The ab initio parametrized site-based Hamiltonian of Binder et al. [Phys. Rev. Lett., 2018, 120, 227401] is employed to model a 20-site system, including intra-ring and inter-ring high-frequency modes as well as torsional modes which undergo thermal fluctuations induced by an explicit harmonic oscillator bath. The system-bath dynamics is treated within the setting of a stochastic mean-field Schrödinger equation. For the 20-site excitonic system, a total of 20 Frenkel states and 248 modes are propagated using the multi-layer multi-configuration time-dependent Hartree (ML-MCTDH) method. The resulting dynamics can be interpreted in terms of the coherent motion of an exciton-polaron quasi-particle stochastically driven by torsional fluctuations. This dynamics yields a near-linear mean squared displacement (MSD) as a function of time, from which a diffusion coefficient can be deduced which increases with temperature, up to 5.7 × 10-3 cm2 s-1 at T = 300 K.
Collapse
Affiliation(s)
- Robert Binder
- Institute of Physical and Theoretical Chemistry, Goethe University Frankfurt, Max-von-Laue-Str. 7, 60438 Frankfurt am Main, Germany.
| | | |
Collapse
|
29
|
Rey M. Group-theoretical formulation of an Eckart-frame kinetic energy operator in curvilinear coordinates for polyatomic molecules. J Chem Phys 2019; 151:024101. [DOI: 10.1063/1.5109482] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Michael Rey
- Groupe de Spectrométrie Moléculaire et Atmosphérique, UMR CNRS 7331, BP 1039, F-51687 Reims Cedex 2, France
| |
Collapse
|
30
|
Puzzarini C, Bloino J, Tasinato N, Barone V. Accuracy and Interpretability: The Devil and the Holy Grail. New Routes across Old Boundaries in Computational Spectroscopy. Chem Rev 2019; 119:8131-8191. [DOI: 10.1021/acs.chemrev.9b00007] [Citation(s) in RCA: 114] [Impact Index Per Article: 22.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Affiliation(s)
- Cristina Puzzarini
- Dipartimento di Chimica “Giacomo Ciamician”, Università di Bologna, Via F. Selmi 2, I-40126 Bologna, Italy
| | - Julien Bloino
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| | - Nicola Tasinato
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| | - Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| |
Collapse
|
31
|
Polimeno A, Zerbetto M, Abergel D. Stochastic modeling of macromolecules in solution. I. Relaxation processes. J Chem Phys 2019; 150:184107. [PMID: 31091939 DOI: 10.1063/1.5077065] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Abstract
A framework for the stochastic description of relaxation processes in flexible macromolecules, including dissipative effects, is introduced from an atomistic point of view. Projection-operator techniques are employed to obtain multidimensional Fokker-Planck operators governing the relaxation of internal coordinates and global degrees of freedom and depending upon parameters fully recoverable from classic force fields (energetics) and continuum models (friction tensors). A hierarchy of approaches of different complexity is proposed in this unified context, aimed primarily at the interpretation of magnetic resonance relaxation experiments. In particular, a model based on a harmonic internal Hamiltonian is discussed as a test case.
Collapse
Affiliation(s)
- Antonino Polimeno
- Dipartimento di Scienze Chimiche, Università degli Studi di Padova, I-35131 Padova, Italy
| | - Mirco Zerbetto
- Dipartimento di Scienze Chimiche, Università degli Studi di Padova, I-35131 Padova, Italy
| | - Daniel Abergel
- Laboratoire des Biomolécules, LBM, Département de Chimie, Ecole Normale Supérieure, PSL University, Sorbonne Université, CNRS, 75005 Paris, France
| |
Collapse
|
32
|
Avila G, Mátyus E. Toward breaking the curse of dimensionality in (ro)vibrational computations of molecular systems with multiple large-amplitude motions. J Chem Phys 2019; 150:174107. [PMID: 31067897 DOI: 10.1063/1.5090846] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Methodological progress is reported in the challenging direction of a black-box-type variational solution of the (ro)vibrational Schrödinger equation applicable to floppy, polyatomic systems with multiple large-amplitude motions. This progress is achieved through the combination of (i) the numerical kinetic-energy operator (KEO) approach of Mátyus et al. [J. Chem. Phys. 130, 134112 (2009)] and (ii) the Smolyak nonproduct grid method of Avila and Carrington, Jr. [J. Chem. Phys. 131, 174103 (2009)]. The numerical representation of the KEO makes it possible to choose internal coordinates and a body-fixed frame best suited for the molecular system. The Smolyak scheme reduces the size of the direct-product grid representation by orders of magnitude, while retaining some of the useful features of it. As a result, multidimensional (ro)vibrational states are computed with system-adapted coordinates, a compact basis- and grid-representation, and an iterative eigensolver. Details of the methodological developments and the first numerical applications are presented for the CH4·Ar complex treated in full (12D) vibrational dimensionality.
Collapse
Affiliation(s)
- Gustavo Avila
- Institute of Chemistry, ELTE, Eötvös Loránd University, Pázmány Péter sétány 1/A, Budapest 1117, Hungary
| | - Edit Mátyus
- Institute of Chemistry, ELTE, Eötvös Loránd University, Pázmány Péter sétány 1/A, Budapest 1117, Hungary
| |
Collapse
|
33
|
Mátyus E. Non-adiabatic mass correction to the rovibrational states of molecules: Numerical application for the H 2 + molecular ion. J Chem Phys 2018; 149:194111. [PMID: 30466265 DOI: 10.1063/1.5050401] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
General transformation expressions of the second-order non-adiabatic Hamiltonian of the atomic nuclei, including the kinetic-energy correction terms, are derived upon the change from laboratory-fixed Cartesian coordinates to general curvilinear coordinate systems commonly used in rovibrational computations. The kinetic-energy or so-called "mass-correction" tensor elements are computed with the stochastic variational method and floating explicitly correlated Gaussian functions for the H 2 + molecular ion in its ground electronic state. {Further numerical applications for the 4 He 2 + molecular ion are presented in the forthcoming paper, Paper II [E. Mátyus, J. Chem. Phys. 149, 194112 (2018)]}. The general, curvilinear non-adiabatic kinetic energy operator expressions are used in the examples, and non-adiabatic rovibrational energies and corrections are determined by solving the rovibrational Schrödinger equation including the diagonal Born-Oppenheimer as well as the mass-tensor corrections.
Collapse
Affiliation(s)
- Edit Mátyus
- Institute of Chemistry, Eötvös Loránd University, Pázmány Péter sétány 1/A, Budapest H-1117, Hungary
| |
Collapse
|
34
|
Affiliation(s)
- Edit Mátyus
- Institute of Chemistry, ELTE Eötvös Loránd University, Budapest, Hungary
| |
Collapse
|
35
|
Szekeres B, Pártay LB, Mátyus E. Direct Computation of the Quantum Partition Function by Path-Integral Nested Sampling. J Chem Theory Comput 2018; 14:4353-4359. [PMID: 29944376 DOI: 10.1021/acs.jctc.8b00368] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
In the present work we introduce a computational approach to the absolute rovibrational quantum partition function using the path-integral formalism of quantum mechanics in combination with the nested sampling technique. The numerical applicability of path-integral nested sampling is demonstrated for small molecules of spectroscopic interest. The computational cost of the method is determined by the evaluation time of a point on the potential energy surface (PES). For efficient PES implementations, the path-integral nested sampling method can be a viable alternative to the direct-Boltzmann-summation technique of variationally computed rovibrational energies, especially for medium-sized molecules and at elevated temperatures.
Collapse
Affiliation(s)
- Béla Szekeres
- Institute of Chemistry , Eötvös Loránd University , Pázmány Péter sétány 1/A , Budapest H-1117 , Hungary.,Department of Numerical Analysis, Faculty of Informatics , Eötvös Loránd University , Pázmány Péter sétány 1/C , Budapest H-1117 , Hungary
| | - Lívia B Pártay
- Department of Chemistry , University of Reading , Whiteknights, Reading RG6 6AD , United Kingdom
| | - Edit Mátyus
- Institute of Chemistry , Eötvös Loránd University , Pázmány Péter sétány 1/A , Budapest H-1117 , Hungary
| |
Collapse
|
36
|
Liang R, Cotton SJ, Binder R, Hegger R, Burghardt I, Miller WH. The symmetrical quasi-classical approach to electronically nonadiabatic dynamics applied to ultrafast exciton migration processes in semiconducting polymers. J Chem Phys 2018; 149:044101. [DOI: 10.1063/1.5037815] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
Affiliation(s)
- Ruibin Liang
- Department of Chemistry and Kenneth S. Pitzer Center for Theoretical Chemistry, University of California, Berkeley,
California 94720, USA and Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720,
USA
| | - Stephen J. Cotton
- Department of Chemistry and Kenneth S. Pitzer Center for Theoretical Chemistry, University of California, Berkeley,
California 94720, USA and Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720,
USA
| | - Robert Binder
- Institute of Physical and Theoretical Chemistry, Goethe University Frankfurt, Max-von-Laue-Str. 7, 60438 Frankfurt/Main,
Germany
| | - Rainer Hegger
- Institute of Physical and Theoretical Chemistry, Goethe University Frankfurt, Max-von-Laue-Str. 7, 60438 Frankfurt/Main,
Germany
| | - Irene Burghardt
- Institute of Physical and Theoretical Chemistry, Goethe University Frankfurt, Max-von-Laue-Str. 7, 60438 Frankfurt/Main,
Germany
| | - William H. Miller
- Department of Chemistry and Kenneth S. Pitzer Center for Theoretical Chemistry, University of California, Berkeley,
California 94720, USA and Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720,
USA
| |
Collapse
|
37
|
Khoma M, Jaquet R. The kinetic energy operator for distance-dependent effective nuclear masses: Derivation for a triatomic molecule. J Chem Phys 2018; 147:114106. [PMID: 28938805 DOI: 10.1063/1.5000267] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The kinetic energy operator for triatomic molecules with coordinate or distance-dependent nuclear masses has been derived. By combination of the chain rule method and the analysis of infinitesimal variations of molecular coordinates, a simple and general technique for the construction of the kinetic energy operator has been proposed. The asymptotic properties of the Hamiltonian have been investigated with respect to the ratio of the electron and proton mass. We have demonstrated that an ad hoc introduction of distance (and direction) dependent nuclear masses in Cartesian coordinates preserves the total rotational invariance of the problem. With the help of Wigner rotation functions, an effective Hamiltonian for nuclear motion can be derived. In the derivation, we have focused on the effective trinuclear Hamiltonian. All necessary matrix elements are given in closed analytical form. Preliminary results for the influence of non-adiabaticity on vibrational band origins are presented for H3+.
Collapse
Affiliation(s)
- Mykhaylo Khoma
- Theoretical Chemistry, Siegen University, Siegen, Germany
| | - Ralph Jaquet
- Theoretical Chemistry, Siegen University, Siegen, Germany
| |
Collapse
|
38
|
|
39
|
Binder R, Lauvergnat D, Burghardt I. Conformational Dynamics Guides Coherent Exciton Migration in Conjugated Polymer Materials: First-Principles Quantum Dynamical Study. PHYSICAL REVIEW LETTERS 2018; 120:227401. [PMID: 29906150 DOI: 10.1103/physrevlett.120.227401] [Citation(s) in RCA: 35] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/06/2017] [Indexed: 06/08/2023]
Abstract
We report on high-dimensional quantum dynamical simulations of photoinduced exciton migration in a single-chain oligothiophene segment, in view of elucidating the controversial nature of the elementary exciton transport steps in semiconducting polymers. A novel first-principles parametrized Frenkel J aggregate Hamiltonian is employed that goes significantly beyond the standard Frenkel-Holstein Hamiltonian. Departing from a nonequilibrium state created by photoexcitation, these simulations provide evidence of an ultrafast two-timescale process at low temperatures, involving exciton-polaron formation within tens of femtoseconds (fs), followed by torsional relaxation on an ∼400 fs timescale. The second step is the driving force for exciton migration, as initial conjugation breaks are removed by dynamical planarization. The quantum coherent nature of the elementary exciton migration step is consistent with experimental observations highlighting the correlated and vibrationally coherent nature of the dynamics on ultrafast timescales.
Collapse
Affiliation(s)
- Robert Binder
- Institute of Physical and Theoretical Chemistry, Goethe University Frankfurt, Max-von-Laue-Straße 7, 60438 Frankfurt, Germany
| | - David Lauvergnat
- Laboratoire de Chimie Physique, Université Paris-Sud, UMR 8000, 91405 Orsay, France
| | - Irene Burghardt
- Institute of Physical and Theoretical Chemistry, Goethe University Frankfurt, Max-von-Laue-Straße 7, 60438 Frankfurt, Germany
| |
Collapse
|
40
|
Manzhos S, Wang X, Carrington T. A multimode-like scheme for selecting the centers of Gaussian basis functions when computing vibrational spectra. Chem Phys 2018. [DOI: 10.1016/j.chemphys.2017.10.006] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
41
|
Nauts A, Lauvergnat D. Numerical on-the-fly implementation of the action of the kinetic energy operator on a vibrational wave function: application to methanol. Mol Phys 2018. [DOI: 10.1080/00268976.2018.1473652] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/16/2022]
Affiliation(s)
- André Nauts
- Institute of Condensed Matter and Nanosciences (NAPS), Université Catholique de Louvain, Louvain-la-Neuve, Belgium
- Laboratoire de Chimie Physique, CNRS, Univ. Paris-Sud, Université Paris-Saclay, Orsay, France
| | - David Lauvergnat
- Laboratoire de Chimie Physique, CNRS, Univ. Paris-Sud, Université Paris-Saclay, Orsay, France
| |
Collapse
|
42
|
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
|
43
|
Perveaux A, Lorphelin M, Lasorne B, Lauvergnat D. Fast and slow excited-state intramolecular proton transfer in 3-hydroxychromone: a two-state story? Phys Chem Chem Phys 2018; 19:6579-6593. [PMID: 28203670 DOI: 10.1039/c6cp06603h] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
The photodynamics of 3-hydroxychromone in its first-excited singlet electronic state (bright state of ππ* character) is investigated with special emphasis given to two types of reaction pathways: the excited-state intramolecular-proton-transfer coordinate and the hydrogen-torsion coordinate linking the excited cis and trans isomers. A newly-found conical intersection with the second-excited singlet electronic state (dark state of nπ* character) is suspected to be, to some extent, the reason for the slower rate constant. This hypothesis based on quantum-chemistry calculations is supported by quantum-dynamics simulations in full dimensionality. They show significant transfer of electronic population and provide consistently a vibronic interpretation for the forbidden band in the UV absorption spectrum.
Collapse
Affiliation(s)
- Aurelie Perveaux
- Laboratoire de Chimie Physique (UMR 8000), CNRS, Université Paris-Sud, F-91405 Orsay, France. and Institut Charles Gerhardt (UMR 5253), CNRS, Université de Montpellier, F-34095 Montpellier, France
| | - Maxime Lorphelin
- Laboratoire de Chimie Physique (UMR 8000), CNRS, Université Paris-Sud, F-91405 Orsay, France.
| | - Benjamin Lasorne
- Institut Charles Gerhardt (UMR 5253), CNRS, Université de Montpellier, F-34095 Montpellier, France
| | - David Lauvergnat
- Laboratoire de Chimie Physique (UMR 8000), CNRS, Université Paris-Sud, F-91405 Orsay, France.
| |
Collapse
|
44
|
Yachmenev A, Küpper J. Communication: General variational approach to nuclear-quadrupole coupling in rovibrational spectra of polyatomic molecules. J Chem Phys 2017; 147:141101. [PMID: 29031262 DOI: 10.1063/1.5002533] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A general algorithm for computing the quadrupole-hyperfine effects in the rovibrational spectra of polyatomic molecules is presented for the case of ammonia (NH3). The method extends the general variational approach TROVE [J. Mol. Spectrosc. 245, 126-140 (2007)] by adding the extra term in the Hamiltonian that describes the nuclear quadrupole coupling, with no inherent limitation on the number of quadrupolar nuclei in a molecule. We applied the new approach to compute the nitrogen-nuclear-quadrupole hyperfine structure in the rovibrational spectrum of NH314. These results agree very well with recent experimental spectroscopic data for the pure rotational transitions in the ground vibrational and ν2 states and the rovibrational transitions in the ν1, ν3, 2ν4, and ν1 + ν3 bands. The computed hyperfine-resolved rovibrational spectrum of ammonia will be beneficial for the assignment of experimental rovibrational spectra, further detection of ammonia in interstellar space, and studies of the proton-to-electron mass variation.
Collapse
Affiliation(s)
- Andrey Yachmenev
- Center for Free-Electron Laser Science, Deutsches Elektronen-Synchrotron DESY, Notkestrasse 85, 22607 Hamburg, Germany
| | - Jochen Küpper
- Center for Free-Electron Laser Science, Deutsches Elektronen-Synchrotron DESY, Notkestrasse 85, 22607 Hamburg, Germany
| |
Collapse
|
45
|
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
|
46
|
Bulik IW, Frisch MJ, Vaccaro PH. Vibrational self-consistent field theory using optimized curvilinear coordinates. J Chem Phys 2017; 147:044110. [DOI: 10.1063/1.4995440] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Ireneusz W. Bulik
- Department of Chemistry, Yale University, New Haven, Connecticut 06511, USA
| | - Michael J. Frisch
- Gaussian, Inc., 340 Quinnipiac St. Bldg. 40, Wallingford, Connecticut 06492, USA
| | - Patrick H. Vaccaro
- Department of Chemistry, Yale University, New Haven, Connecticut 06511, USA
| |
Collapse
|
47
|
Carrington T. Perspective: Computing (ro-)vibrational spectra of molecules with more than four atoms. J Chem Phys 2017; 146:120902. [DOI: 10.1063/1.4979117] [Citation(s) in RCA: 67] [Impact Index Per Article: 9.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Tucker Carrington
- Chemistry Department, Queen’s University, Kingston, Ontario
K7L 3N6, Canada
| |
Collapse
|
48
|
Peláez D, Meyer HD. On the infrared absorption spectrum of the hydrated hydroxide (H3O2-) cluster anion. Chem Phys 2017. [DOI: 10.1016/j.chemphys.2016.08.025] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
49
|
Kopp WA, Leonhard K. General formulation of rovibrational kinetic energy operators and matrix elements in internal bond-angle coordinates using factorized Jacobians. J Chem Phys 2016; 145:234102. [DOI: 10.1063/1.4972224] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Wassja A. Kopp
- Lehrstuhl für Technische Thermodynamik, RWTH Aachen University, Aachen 52062, Germany
| | - Kai Leonhard
- Lehrstuhl für Technische Thermodynamik, RWTH Aachen University, Aachen 52062, Germany
| |
Collapse
|
50
|
Manzhos S, Carrington T. Using an internal coordinate Gaussian basis and a space-fixed Cartesian coordinate kinetic energy operator to compute a vibrational spectrum with rectangular collocation. J Chem Phys 2016; 145:224110. [DOI: 10.1063/1.4971295] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Sergei Manzhos
- Department of Mechanical Engineering, National University of Singapore, Block EA #07-08, 9 Engineering Drive 1, 117576 Singapore
| | - Tucker Carrington
- Chemistry Department, Queen’s University, Kingston, Ontario K7L 3N6, Canada
| |
Collapse
|