1
|
Wang XG, Carrington T. A two-step quadrature-based variational calculation of ro-vibrational levels and wavefunctions of CO 2 using a bisector- x molecule-fixed frame. Phys Chem Chem Phys 2024; 26:15181-15191. [PMID: 38752328 DOI: 10.1039/d4cp00655k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/30/2024]
Abstract
In this paper, we propose a new two-step strategy for computing ro-vibrational energy levels and wavefunctions of a triatomic molecule and apply it to CO2. A two-step method [J. Tennyson and B. T. Sutcliffe, Mol. Phys., 1986, 58, 1067] uses a basis whose functions are products of K-dependent "vibrational" functions and symmetric top functions. K is the quantum number for the molecule-fixed z component of the angular momentum. For a linear molecule, a two-step method is efficient because the Hamiltonian used to compute the basis functions includes the largest coupling term. The most important distinguishing feature of the two-step method we propose is that it uses an associated Legendre basis and quadrature rather than a K-dependent discrete variable representation. This reduces the cost of the calculation and simplifies the method. We have computed ro-vibrational energy levels with J up to 100 for CO2, on an accurate available potential energy surface which is known as the AMES-2 PES and present a subset of those levels. We have converged most levels up to 20 000 cm-1 to 0.0001 cm-1.
Collapse
Affiliation(s)
- Xiao-Gang Wang
- Chemistry Department, Queen's University, Kingston, Ontario K7L 3N6, Canada.
| | - Tucker Carrington
- Chemistry Department, Queen's University, Kingston, Ontario K7L 3N6, Canada.
| |
Collapse
|
2
|
Martín Santa Daría A, Avila G, Mátyus E. Methane dimer rovibrational states and Raman transition moments. Phys Chem Chem Phys 2024; 26:10254-10264. [PMID: 38497527 DOI: 10.1039/d3cp06222h] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/19/2024]
Abstract
Benchmark-quality rovibrational data are reported for the methane dimer from variational nuclear motion computations using an ab initio intermolecular potential energy surface reported by [M. P. Metz et al., Phys. Chem. Chem. Phys., 2019, 21, 13504-13525]. A simple polarizability model is used to compute Raman transition moments that may be relevant for future direct observation of the intermolecular dynamics. Non-negligible ΔK ≠ 0 transition moments arise in this symmetric top system due to strong rovibrational couplings.
Collapse
Affiliation(s)
- Alberto Martín Santa Daría
- Departamento de Química Física, University of Salamanca, 37008 Salamanca, Spain
- 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.
| | - 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.
| |
Collapse
|
3
|
Wang XG, Carrington T. Computing excited OH stretch states of water dimer in 12D using contracted intermolecular and intramolecular basis functions. J Chem Phys 2023; 158:084107. [PMID: 36859104 DOI: 10.1063/5.0139586] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023] Open
Abstract
Due to the ubiquity and importance of water, water dimer has been intensively studied. Computing the (ro-)vibrational spectrum of water dimer is challenging. The potential has eight wells separated by low barriers, which makes harmonic approximations of limited utility. A variational approach is imperative, but difficult because there are 12 coupled vibrational coordinates. In this paper, we use a product contracted basis whose functions are products of intramolecular and intermolecular functions computed using an iterative eigensolver. An intermediate matrix F facilitates calculating matrix elements. Using F, it is possible to do calculations on a general potential without storing the potential on the full quadrature grid. We find that surprisingly many intermolecular functions are required. This is due to the importance of coupling between inter- and intra-molecular coordinates. The full G16 symmetry of water dimer is exploited. We calculate, for the first time, monomer excited stretch states and compare P(1) transition frequencies with their experimental counterparts. We also compare with experimental vibrational shifts and tunneling splittings. Surprisingly, we find that the largest tunneling splitting, which does not involve the interchange of the two monomers, is smaller in the asymmetric stretch excited state than in the ground state. Differences between levels we compute and those obtained with a [6+6]D adiabatic approximation [Leforestier et al. J. Chem. Phys. 137 014305 (2012)] are ∼0.6 cm-1 for states without monomer excitation, ∼4 cm-1 for monomer excited bend states, and as large as ∼10 cm-1 for monomer excited stretch states.
Collapse
Affiliation(s)
- Xiao-Gang Wang
- Chemistry Department, Queen's University, Kingston, Ontario K7L 3N6, Canada
| | - Tucker Carrington
- Chemistry Department, Queen's University, Kingston, Ontario K7L 3N6, Canada
| |
Collapse
|
4
|
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
|
5
|
Wang XG, Carrington T. Using nondirect product Wigner D basis functions and the symmetry-adapted Lanczos algorithm to compute the ro-vibrational spectrum of CH 4-H 2O. J Chem Phys 2021; 154:124112. [PMID: 33810654 DOI: 10.1063/5.0044010] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
By doing calculations on the methane-water van der Waals complex, we demonstrate that highly converged energy levels and wavefunctions can be obtained using Wigner D basis functions and the Symmetry-Adapted Lanczos (SAL) method. The Wigner D basis is a nondirect product basis and, therefore, efficient when the kinetic energy operator has accessible singularities. The SAL method makes it possible to exploit symmetry to label energy levels and reduce the cost of the calculation, without explicitly using symmetry-adapted basis functions. Line strengths are computed, and new bands are identified. In particular, we find unusually strong transitions between states associated with the isomers of the global minimum and the secondary minimum.
Collapse
Affiliation(s)
- Xiao-Gang Wang
- Chemistry Department, Queen's University, Kingston, Ontario K7L 3N6, Canada
| | - Tucker Carrington
- Chemistry Department, Queen's University, Kingston, Ontario K7L 3N6, Canada
| |
Collapse
|
6
|
Zhao B, Manthe U. Direct product-type grid representations for angular coordinates in extended space and their application in the MCTDH approach. J Chem Phys 2021; 154:104115. [PMID: 33722051 DOI: 10.1063/5.0045054] [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/14/2022] Open
Abstract
Multi-configurational time-dependent Hartree (MCTDH) calculations using time-dependent grid representations can be used to accurately simulate high-dimensional quantum dynamics on general ab initio potential energy surfaces. Employing the correlation discrete variable representation, sets of direct product type grids are employed in the calculation of the required potential energy matrix elements. This direct product structure can be a problem if the coordinate system includes polar and azimuthal angles that result in singularities in the kinetic energy operator. In the present work, a new direct product-type discrete variable representation (DVR) for arbitrary sets of polar and azimuthal angles is introduced. It employs an extended coordinate space where the range of the polar angles is taken to be [-π, π]. The resulting extended space DVR resolves problems caused by the singularities in the kinetic energy operator without generating a very large spectral width. MCTDH calculations studying the F·CH4 complex are used to investigate important properties of the new scheme. The scheme is found to allow for more efficient integration of the equations of motion compared to the previously employed cot-DVR approach [G. Schiffel and U. Manthe, Chem. Phys. 374, 118 (2010)] and decreases the required central processing unit times by about an order of magnitude.
Collapse
Affiliation(s)
- Bin Zhao
- Theoretische Chemie, Fakultät für Chemie, Universität Bielefeld, Universitätsstr. 25, D-33615 Bielefeld, Germany
| | - Uwe Manthe
- Theoretische Chemie, Fakultät für Chemie, Universität Bielefeld, Universitätsstr. 25, D-33615 Bielefeld, Germany
| |
Collapse
|
7
|
Viglaska D, Rey M, Nikitin AV, Tyuterev VG. Derivation of ρ-dependent coordinate transformations for nonrigid molecules in the Hougen-Bunker-Johns formalism. J Chem Phys 2020; 153:084102. [PMID: 32872870 DOI: 10.1063/5.0016365] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
In this paper, we report a series of transformations for the construction of a Hamiltonian model for nonrigid polyatomic molecules in the framework of the Hougen-Bunker-Johns formalism (HBJ). This model is expressed in normal mode coordinates for small vibrations and in a specific coordinate ρ to describe the large amplitude motion. For the first time, a general procedure linking the "true" curvilinear coordinates to ρ is proposed, allowing the expression of the potential energy part in the same coordinate representation as the kinetic energy operator, whatever the number of atoms. A Lie group-based method is also proposed for the derivation of the reference configuration in the internal axis system. This work opens new perspectives for future high-resolution spectroscopy studies of nonrigid, medium-sized molecules using HBJ-type Hamiltonians. Illustrative examples and computation of vibrational energy levels on semirigid and nonrigid molecules are given to validate this method.
Collapse
Affiliation(s)
- Dominika Viglaska
- Chemistry Department, Queen's University, Kingston, Ontario K7L 3N6, Canada
| | - Michael Rey
- Groupe de Spectrométrie Moléculaire et Atmosphérique, UMR CNRS 7331, BP 1039, F-51687 Reims Cedex 2, France
| | - Andrei V Nikitin
- Laboratory of Theoretical Spectroscopy, Institute of Atmospheric Optics, SB RAS, 634055 Tomsk, Russia
| | - Vladimir G Tyuterev
- Groupe de Spectrométrie Moléculaire et Atmosphérique, UMR CNRS 7331, BP 1039, F-51687 Reims Cedex 2, France
| |
Collapse
|
8
|
Wang XG, Carrington T. Using quadrature and an iterative eigensolver to compute fine-structure ro-vibrational levels of Van der Waals complexes: NH(Σ−3)–He, O 2(Σg−3)–Ar, and O 2(Σg−3)–He. J Chem Phys 2019. [DOI: 10.1063/1.5110873] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023] Open
Affiliation(s)
- Xiao-Gang Wang
- Chemistry Department, Queen’s University, Kingston, Ontario K7L 3N6, Canada
| | - Tucker Carrington
- Chemistry Department, Queen’s University, Kingston, Ontario K7L 3N6, Canada
| |
Collapse
|
9
|
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
|
10
|
An efficient pseudo-spectral method for the description of atomic electronic wave functions – Application to the hydrogen atom in a uniform magnetic field. Chem Phys 2018. [DOI: 10.1016/j.chemphys.2018.09.025] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
|
11
|
Wang XG, Carrington T. Using monomer vibrational wavefunctions as contracted basis functions to compute rovibrational levels of an H2O-atom complex in full dimensionality. J Chem Phys 2017; 146:104105. [DOI: 10.1063/1.4977179] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Xiao-Gang Wang
- Chemistry Department, Queen’s University, Kingston, Ontario K7L 3N6, Canada
| | - Tucker Carrington
- Chemistry Department, Queen’s University, Kingston, Ontario K7L 3N6, Canada
| |
Collapse
|
12
|
Ndengué S, Dawes R, Wang XG, Carrington T, Sun Z, Guo H. Calculated vibrational states of ozone up to dissociation. J Chem Phys 2016; 144:074302. [DOI: 10.1063/1.4941559] [Citation(s) in RCA: 36] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Steve Ndengué
- Department of Chemistry, Missouri University of Science and Technology, Rolla, Missouri 65409, USA
| | - Richard Dawes
- Department of Chemistry, Missouri University of Science and Technology, Rolla, Missouri 65409, USA
| | - Xiao-Gang Wang
- Chemistry Department, Queen’s University, Kingston, Ontario K7L 3N6, Canada
| | - Tucker Carrington
- Chemistry Department, Queen’s University, Kingston, Ontario K7L 3N6, Canada
| | - Zhigang Sun
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China and Center for Advanced Chemical Physics and 2011 Frontier Center for Quantum Science and Technology, University of Science and Technology of China, 96 Jinzhai Road, Hefei 230026, China
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, USA
| |
Collapse
|
13
|
Wang XG, Carrington T. Computing rovibrational levels of polyatomic molecules with polyspherical coordinates and a contracted basis built with a K-independent vibrational primitive basis. Mol Phys 2013. [DOI: 10.1080/00268976.2013.808387] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
Affiliation(s)
- Xiao-Gang Wang
- Chemistry Department, Queen’s University , Kingston, Ontario, Canada
| | - Tucker Carrington
- Chemistry Department, Queen’s University , Kingston, Ontario, Canada
| |
Collapse
|
14
|
Jaquet R, Carrington T. Using a Nondirect Product Basis to Compute J > 0 Rovibrational States of H3+. J Phys Chem A 2013; 117:9493-500. [DOI: 10.1021/jp312027s] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Ralph Jaquet
- Theoretische Chemie, Universität Siegen, D-57068 Siegen, Germany
| | - Tucker Carrington
- Chemistry Department, Queens’s University, Kingston, Ontario, K7L
3N6, Canada
| |
Collapse
|
15
|
Wang XG, Carrington T. Computing rovibrational levels of methane with curvilinear internal vibrational coordinates and an Eckart frame. J Chem Phys 2013; 138:104106. [DOI: 10.1063/1.4793474] [Citation(s) in RCA: 66] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|