1
|
Rey M, Carrington T. Using nested tensor train contracted basis functions with group theoretical techniques to compute (ro)-vibrational spectra of molecules with non-Abelian groups. J Chem Phys 2024; 161:044102. [PMID: 39037133 DOI: 10.1063/5.0219434] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/17/2024] [Accepted: 06/30/2024] [Indexed: 07/23/2024] Open
Abstract
In this paper, we use nested tensor-train contractions to compute vibrational and ro-vibrational energy levels of molecules with five and six atoms. At each step, we fully exploit symmetry by using symmetry adapted basis functions obtained from an irreducible tensor method. Contracted basis functions are determined by diagonalizing reduced dimensional Hamiltonian matrices. The size of matrices of eigenvectors, used to account for coupling between groups of coordinates, is reduced by discarding rows and columns. The size of the matrices that must be diagonalized is thus substantially reduced, making it possible to use direct eigensolvers, even for molecules with five and six atoms. The symmetry-adapted contracted vibrational basis functions have been used to compute J = 0 energy levels of the CH3CN (C3v) and J > 0 levels of CH4.
Collapse
Affiliation(s)
- Michaël Rey
- Groupe de Spectrométrie Moléculaire et Atmosphérique, UMR CNRS 7331, BP 1039, F-51687 Reims Cedex 2, France
| | - Tucker Carrington
- Chemistry Department, Queen's University, Kingston, Ontario K7L 3N6, Canada
| |
Collapse
|
2
|
Palanisamy N, Banik S. An approximation to the vibrational coupled-cluster method for CH-stretching of large molecules: application to naphthalene and anthracene. Phys Chem Chem Phys 2023; 25:20076-20092. [PMID: 37462438 DOI: 10.1039/d3cp01313h] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/27/2023]
Abstract
We propose an approximation to the vibrational coupled-cluster method (VCCM) to describe the CH-stretching region of the vibrational spectrum of large molecules. The vibrational modes of a molecule are divided into two sets: the target set and the bath set. The target set includes the CH stretches and the modes that are strongly coupled with the CH stretches and/or involve strong Fermi resonances with a CH stretch fundamental. The rest of the modes are in the bath set. First, the effective harmonic oscillator (EHO) approximation is invoked for the whole system to obtain the zeroth-order frequencies and modified potentials. The effects of interaction between the bath set and the target sets are included in the modified potential from the EHO calculation. The VCCM equations are constructed with the modified potential from the EHO calculations and for the target set only. The transition energies and intensities are calculated using such a truncated VCCM approximation. The proposed method is applied to calculate the IR spectra of naphthalene and anthracene. The results with three different criteria for selecting the modes in the target set are compared with the experimental IR spectra.
Collapse
Affiliation(s)
- Nivedhitha Palanisamy
- Department of Chemistry, School of Chemical and Biotechnology, SASTRA Deemed University, Thanjavur, Tamil Nadu, India.
| | - Subrata Banik
- Department of Chemistry, School of Chemical and Biotechnology, SASTRA Deemed University, Thanjavur, Tamil Nadu, India.
| |
Collapse
|
3
|
Harwick OA, Fortenberry RC. Spectroscopic Constants and Anharmonic Vibrational Frequencies of C(O)OC, c-C 2O 2 and Their Silicon-Containing Analogues. Molecules 2023; 28:molecules28114563. [PMID: 37299037 DOI: 10.3390/molecules28114563] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/11/2023] [Revised: 05/31/2023] [Accepted: 06/01/2023] [Indexed: 06/12/2023] Open
Abstract
Comets are likely to contain various carbon oxide molecules potentially including C(O)OC and c-C2O2 on their surfaces and comae, as well as their silicon-substituted analogues possibly playing a role in the formation of interstellar dust grains. In this work, high-level quantum chemical data are provided to support such potential future astrophysical detection through the generation of predicted rovibrational data. Laboratory-based chemistry would also benefit from such aforementioned computational benchmarking considering these molecules' historic computational and experimental elusiveness. Coupled-cluster singles, doubles, and perturbative triples, the F12b formalism, and the cc-pCVTZ-F12 basis set garner the rapid, yet highly trusted F12-TcCR level of theory leveraged presently. This current work points to all four molecules' strong IR activity, coupled with large intensities, thus suggesting the potential for JWST detection. Although Si(O)OSi possesses a permanent dipole moment significantly larger than those of the other molecules of present interest, the significant abundance of the potential precursor carbon monoxide suggests that the dicarbon dioxide molecules may yet be observable in the microwave region of the electromagnetic spectrum. Thus, this present work details the likely existence and detectability of these four cyclic molecules, providing updated implications compared to previous work performed both experimentally and computationally.
Collapse
Affiliation(s)
- Olivia A Harwick
- Department of Chemistry and Biochemistry, University of Mississippi, University, MS 38677, USA
| | - Ryan C Fortenberry
- Department of Chemistry and Biochemistry, University of Mississippi, University, MS 38677, USA
| |
Collapse
|
4
|
Yang Q, Bloino J. An Effective and Automated Processing of Resonances in Vibrational Perturbation Theory Applied to Spectroscopy. J Phys Chem A 2022; 126:9276-9302. [DOI: 10.1021/acs.jpca.2c06460] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/02/2022]
Affiliation(s)
- Qin Yang
- Faculty of Science, Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126Pisa, Italy
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, 16610Prague, Czech Republic
| | - Julien Bloino
- Faculty of Science, Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126Pisa, Italy
| |
Collapse
|
5
|
Rey M. Novel methodology for systematically constructing global effective models from ab initio-based surfaces: A new insight into high-resolution molecular spectra analysis. J Chem Phys 2022; 156:224103. [PMID: 35705402 DOI: 10.1063/5.0089097] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In this paper, a novel methodology is presented for the construction of ab initio effective rotation-vibration spectroscopic models from potential energy and dipole moment surfaces. Non-empirical effective Hamiltonians are obtained via the block-diagonalization of selected variationally computed eigenvector matrices. For the first time, the derivation of an effective dipole moment is carried out in a systematic way. This general approach can be implemented quite easily in most of the variational computer codes and turns out to be a clear alternative to the rather involved Van Vleck perturbation method. Symmetry is exploited at all stages to translate first-principles calculations into a set of spectroscopic parameters to be further refined on experiment. We demonstrate on H2CO, PH3, CH4, C2H4, and SF6 that the proposed effective model can provide crucial information to spectroscopists within a very short time compared to empirical spectroscopic models. This approach brings a new insight into high-resolution spectrum analysis of polyatomic molecules and will be also of great help in the modeling of hot atmospheres where completeness is important.
Collapse
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
|
6
|
Chang X, Dobrolyubov EO, Krasnoshchekov SV. Fundamental studies of vibrational resonance phenomena by multivalued resummation of the divergent Rayleigh-Schrödinger perturbation theory series: deciphering polyad structures of three H 216O isotopologues. Phys Chem Chem Phys 2022; 24:6655-6675. [PMID: 35234755 DOI: 10.1039/d1cp04279c] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023]
Abstract
A fundamental quantitative study of the vibrational resonances of three H216O H,D-isotopologues with a quartic Watson Hamiltonian was carried out using the resummation of the high-order (∼λn, n ≤ 203) divergent Rayleigh-Schrödinger perturbation theory (RSPT) series by quartic Padé-Hermite multivalued diagonal approximants PH[m,m,m,m,m], m ≤ 40. The resonance condition between a pair of states is formulated as the existence of a common complex energy solution branch point inside the unit circle: |E(λj)|, Re(λj)2 + Im(λj)2 ≤ 1. For the matrix formulation of the vibrational problem (VCI), the existence of common branch points is governed by the Katz theorem and they can be found as roots of discriminant polynomials. The main branches of the Padé-Hermite approximants typically reproduce VCI energies with high accuracy while alternative branches often fit nearby resonant states. The resummation of the RSPT series for H2O and D2O (up to the tenth polyad) revealed not only Fermi and Darling-Dennison resonances, but also unusual (0,2,-5) and (5,-2,0) resonance effects matching the (5,2,5) polyad structure, while the (3,2,1) structure was rigorously confirmed for HDO. It is demonstrated that the (5,2,5) polyad structure ensures good organization of high-energy resonating states and breaks down the classic (2,1,2) structure. The advocated methodology of quantitative description of resonance phenomena and revealing polyad structures is suitable for larger molecules and can be adapted to linear molecules and symmetric tops. Its application ensures rigorous classification of vibrational states and can be used in quantitative vibration-rotation spectroscopy.
Collapse
Affiliation(s)
- Xuanhao Chang
- Chemistry Department, Lomonosov Moscow State University, Leninskiye Gory 1-3, Moscow, 119899, Russian Federation.
| | - Egor O Dobrolyubov
- Chemistry Department, Lomonosov Moscow State University, Leninskiye Gory 1-3, Moscow, 119899, Russian Federation.
| | - Sergey V Krasnoshchekov
- Chemistry Department, Lomonosov Moscow State University, Leninskiye Gory 1-3, Moscow, 119899, Russian Federation.
| |
Collapse
|
7
|
Rey M, Chizhmakova IS, Nikitin AV, Tyuterev VG. Towards a complete elucidation of the ro-vibrational band structure in the SF 6 infrared spectrum from full quantum-mechanical calculations. Phys Chem Chem Phys 2021; 23:12115-12126. [PMID: 34032236 DOI: 10.1039/d0cp05727d] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
The first accurate and complete theoretical room-temperature rotationally resolved spectra in the range 300-3000 cm-1 are reported for the three most abundant isotopologues (32SF6, 33SF6 and 34SF6) of the sulfur hexafluoride molecule. The literature reports that SF6 is widely used as a prototype molecule for studying the multi-photon excitation processes with powerful lasers in the infrared range. On the other hand, SF6 is an important greenhouse molecule with a very long lifetime in the atmosphere. Because of relatively low vibrational frequencies, the hot bands of this molecule contribute significantly to the absorption infrared spectra even at room temperature. This makes the calculation of complete rovibrational line lists required for fully converged opacity modeling extremely demanding. In order to reduce the computational costs, symmetry was exploited at all stages of the first global variational nuclear motion calculations by means of irreducible tensor operators. More than 2600 new vibrational band centers were predicted using our empirically refined ab initio potential energy surface. Highly excited rotational states were calculated up to J = 121, resulting in 6 billion transitions computed from an ab initio dipole moment surface and distributed over more than 500 cold and hot bands. The final line lists are made available through the TheoReTS information system (http://theorets.univ-reims.fr, http://theorets.tsu.ru). For the first time, the major (ro)vibrational band structures in the wavenumber range corresponding to the strongest absorption in the infra-red are completely elucidated for a seven-atom molecule, providing excellent agreement with the observed spectral patterns. It is shown that the obtained results are more complete than all available line lists, permitting reliable modelling of the temperature dependence of the molecular opacity.
Collapse
Affiliation(s)
- Michaël Rey
- Groupe de Spectrométrie Moléculaire et Atmosphérique, UMR CNRS 7331, BP 1039, F-51687, Reims Cedex 2, France.
| | - Iana S Chizhmakova
- Laboratory of Quantum Mechanics of Molecules and Radiative Processes, Tomsk State University, 36 Lenin Avenue, 634050 Tomsk, Russia
| | - Andrei V Nikitin
- V.E. Zuev Institute of Atmospheric Optics, Russian Academy of Sciences, 1, Akademichesky Avenue, 634055 Tomsk, Russian Federation
| | - Vladimir G Tyuterev
- Groupe de Spectrométrie Moléculaire et Atmosphérique, UMR CNRS 7331, BP 1039, F-51687, Reims Cedex 2, France. and Laboratory of Quantum Mechanics of Molecules and Radiative Processes, Tomsk State University, 36 Lenin Avenue, 634050 Tomsk, Russia
| |
Collapse
|
8
|
Ferrari BC, Bennett CJ. A computational investigation of the equilibrium geometries, energetics, vibrational frequencies, infrared intensities and Raman activities of C2Oy (y = 3, 4) species. Mol Phys 2020. [DOI: 10.1080/00268976.2020.1837404] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
Affiliation(s)
- Brian C. Ferrari
- Department of Physics, University of Central Florida, Orlando, FL, USA
| | - Chris J. Bennett
- Department of Physics, University of Central Florida, Orlando, FL, USA
| |
Collapse
|
9
|
Nejad A, Crittenden DL. On the separability of large-amplitude motions in anharmonic frequency calculations. Phys Chem Chem Phys 2020; 22:20588-20601. [PMID: 32966420 DOI: 10.1039/d0cp03515g] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Nuclear vibrational theories based upon the Watson Hamiltonian are ubiquitous in quantum chemistry, but are generally unable to model systems in which the wavefunction can delocalise over multiple energy minima, i.e. molecules that have low-energy torsion and inversion barriers. In a 2019 Chemical Reviews article, Puzzarini et al. note that a common workaround is to simply decouple these problematic modes from all other vibrations in the system during anharmonic frequency calculations. They also point out that this approximation can be "ill-suited", but do not quantify the errors introduced. In this work, we present the first systematic investigation into how separating out or constraining torsion and inversion vibrations within potential energy surface (PES) expansions affects the accuracy of computed fundamental wavenumbers for the remaining vibrational modes, using a test set of 19 tetratomic molecules for which high quality analytic potential energy surfaces and fully-coupled anharmonic reference fundamental frequencies are available. We find that the most effective and efficient strategy is to remove the mode in question from the PES expansion entirely. This introduces errors of up to +10 cm-1 in stretching fundamentals that would otherwise couple to the dropped mode, and ±5 cm-1 in all other fundamentals. These errors are approximately commensurate with, but not necessarily additional to, errors due to the choice of electronic structure model used in constructing spectroscopically accurate PES.
Collapse
Affiliation(s)
- Arman Nejad
- Institute of Physical Chemistry, University of Göttingen, Tammannstr. 6, D-37077 Göttingen, Germany.
| | - Deborah L Crittenden
- School of Physical and Chemical Sciences, University of Canterbury, Christchurch, New Zealand
| |
Collapse
|
10
|
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
|
11
|
Sibert EL. Modeling vibrational anharmonicity in infrared spectra of high frequency vibrations of polyatomic molecules. J Chem Phys 2019; 150:090901. [DOI: 10.1063/1.5079626] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Edwin L. Sibert
- Department of Chemistry and Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA
| |
Collapse
|
12
|
Kananenka AA, Skinner JL. Fermi resonance in OH-stretch vibrational spectroscopy of liquid water and the water hexamer. J Chem Phys 2018; 148:244107. [DOI: 10.1063/1.5037113] [Citation(s) in RCA: 40] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Alexei A. Kananenka
- Institute for Molecular Engineering, The University of Chicago, Chicago, Illinois 60637, USA
| | - J. L. Skinner
- Institute for Molecular Engineering, The University of Chicago, Chicago, Illinois 60637, USA
| |
Collapse
|