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
|
Voronin BA, Tennyson J, Yurchenko SN, Chesnokova TY, Chentsov AV, Bykov AD, Makarova MV, Voronina SS, Cruz FC. The infrared absorption spectrum of radioactive water isotopologue H 215O. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2024; 311:124007. [PMID: 38354674 DOI: 10.1016/j.saa.2024.124007] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/03/2023] [Revised: 01/16/2024] [Accepted: 02/05/2024] [Indexed: 02/16/2024]
Abstract
A room temperature line list for the H215O radioactive isotopologue of the water molecule is computed using the variational nuclear-motion DVR3D program suite and an empirical high-precision potential energy function. The line list consists of rotation-vibrational energies and Einstein-A coefficients, covering a wide spectral range from 0 to 25000 cm-1 and the total angular momenta J up to 30. Estimates of air-broadening coefficients are provided. Experimentally derived energies of H216O, H217O and H218O from the literature are used to provide improved energies for important states with uncertainty estimates for the H215O. A number of the wmost promising spectroscopic ranges for the detection of H215O are proposed. The calculated absorption spectrum should be useful for the study gaseous radioactive water at IR region, determining concentration, etc.
Collapse
Affiliation(s)
- Boris A Voronin
- Instituto de Fisica Gleb Wataghin, Universidade Estadual de Campinas, Campinas, SP, 13083-859, Brazil; V.E. Zuev Institute of Atmospheric Optics SB RAS, sq. Ak. Zuev 1, 643021 Tomsk, Russia
| | - Jonathan Tennyson
- Department of Physics and Astronomy, University College London, Gower Street, London, WC1E 6BT, UK.
| | - Sergey N Yurchenko
- Department of Physics and Astronomy, University College London, Gower Street, London, WC1E 6BT, UK
| | - Tatyana Yu Chesnokova
- V.E. Zuev Institute of Atmospheric Optics SB RAS, sq. Ak. Zuev 1, 643021 Tomsk, Russia
| | - Aleksei V Chentsov
- V.E. Zuev Institute of Atmospheric Optics SB RAS, sq. Ak. Zuev 1, 643021 Tomsk, Russia
| | - Aleksandr D Bykov
- V.E. Zuev Institute of Atmospheric Optics SB RAS, sq. Ak. Zuev 1, 643021 Tomsk, Russia
| | - Maria V Makarova
- Faculty of Physics, St. Petersburg State University, 198504 St. Petersburg, Russia
| | - Svetlana S Voronina
- V.E. Zuev Institute of Atmospheric Optics SB RAS, sq. Ak. Zuev 1, 643021 Tomsk, Russia
| | - Flávio C Cruz
- Instituto de Fisica Gleb Wataghin, Universidade Estadual de Campinas, Campinas, SP, 13083-859, Brazil
| |
Collapse
|
3
|
Egorov O, Rey M, Viglaska D, Nikitin AV. Accurate ab initio potential energy surface, rovibrational energy levels and resonance interactions of triplet ( X ~ 3 B 1 ) methylene. J Comput Chem 2024; 45:83-100. [PMID: 37681560 DOI: 10.1002/jcc.27220] [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: 07/14/2023] [Revised: 08/15/2023] [Accepted: 08/17/2023] [Indexed: 09/09/2023]
Abstract
In this work, we report rovibrational energy levels for four isotopologues of methylene (CH2 , CHD, CD2 , and 13 CH2 ) in their ground triplet electronic state (X ~ 3 B1 ) from variational calculation up to ~10,000 cm-1 and using a new accurate ab initio potential energy surface (PES). Triplet methylene exhibits a large-amplitude bending vibration and can reach a quasilinear configuration due to its low barrier (~2000 cm-1 ). To construct the ab initio PES, the Dunning's augmented correlation-consistent core-valence orbital basis sets were employed up to the sextuple-ζ quality [aug-cc-pCVXZ, X = T, Q, 5, and 6] combined with the single- and double-excitation unrestricted coupled cluster approach with a perturbative treatment of triple excitations [RHF-UCCSD(T)]. We have shown that the accuracy of the ab initio energies is further improved by including the corrections due to the scalar relativistic effects, DBOC and high-order electronic correlations. For the first time, all the available experimental rovibrational transitions were reproduced with errors less than 0.12 cm-1 , without any empirical corrections. Unlike more "traditional" nonlinear triatomic molecules, we have shown that even the energies of the ground vibrational state (000) with rather small rotational quantum numbers are strongly affected by the very pronounced rovibrational resonance interactions. Accordingly, the polyad structure of the vibrational levels of CH2 and CD2 was analyzed and discussed. The comparison between the energy levels obtained from the effective Watson A-reduced Hamiltonian, from the generating-function approach and from a variational calculation was given.
Collapse
Affiliation(s)
- Oleg Egorov
- Laboratory of Theoretical Spectroscopy, V.E. Zuev Institute of Atmospheric Optics SB RAS, Tomsk, Russia
- Faculty of Physics, Tomsk State University, Tomsk, Russia
| | - Michaël Rey
- Groupe de Spectrométrie Moléculaire et Atmosphérique, UMR CNRS 7331, UFR Sciences, Reims, France
| | - Dominika Viglaska
- Groupe de Spectrométrie Moléculaire et Atmosphérique, UMR CNRS 7331, UFR Sciences, Reims, France
| | - Andrei V Nikitin
- Laboratory of Theoretical Spectroscopy, V.E. Zuev Institute of Atmospheric Optics SB RAS, Tomsk, Russia
| |
Collapse
|
4
|
High-level ab initio study of disulfur monoxide: ground state potential energy surface and band origins for six isotopic species. Chem Phys Lett 2022. [DOI: 10.1016/j.cplett.2022.140216] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|
5
|
Egorov O, Rey M, Nikitin AV, Viglaska D. New Theoretical Infrared Line List for the Methyl Radical with Accurate Vibrational Band Origins from High-Level Ab Initio Calculations. J Phys Chem A 2022; 126:6429-6442. [PMID: 36073030 DOI: 10.1021/acs.jpca.2c04822] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
In the present work, high-level ab initio calculations were carried out for the ground electronic state of the methyl radical (CH3). Dunning's augmented correlation-consistent orbital basis sets were employed up to the quintuple-ζ valence quality with the core-valence electron correlation [aug-cc-pCV5Z] combined with the single- and double-excitation unrestricted coupled-cluster approach with a perturbative treatment of triple excitations [RHF-UCCSD(T)]. The explicitly correlated version of the coupled-cluster approach [RHF-UCCSD(T)-F12x{x = a, b}] was additionally applied with the core-valence cc-pCVQZ-F12 basis set in order to study convergence with respect to the basis set size. The contributions beyond the coupled-cluster level of the theory like Douglas-Kroll-Hess scalar relativistic Hamiltonian, diabatic Born-Oppenheimer corrections, and high-order electronic correlations have been included into the ab initio potential energy surfaces (PESs). It is shown that the theoretical band origins of CH3 converge gradually to the experimental values when applying the ab initio PESs using the aug-cc-pCVXZ [X = T, Q, and 5] basis sets. For the first time, all available experimental band origins of the gaseous CH3 are reproduced within an accuracy of 0.2 cm-1 using a newly developed PES extrapolated to the complete basis set limit [CBS(TQ5Z)]. The reached accuracy is one order of magnitude better than that of the best available calculations. A new theoretical infrared line list was generated for astrophysical applications using an ab initio dipole moment surface computed at the RHF-UCCSD(T)/aug-cc-pCVQZ level of the theory. The manifestation of a large-amplitude motion in CH3 is also discussed.
Collapse
Affiliation(s)
- Oleg Egorov
- Laboratory of Theoretical Spectroscopy, V.E. Zuev Institute of Atmospheric Optics SB RAS 1, Akademician Zuev Sq., Tomsk 634055 Russia
| | - Michaël Rey
- Groupe de Spectrométrie Moléculaire et Atmosphérique UMR CNRS 7331, UFR Sciences BP 1039, 51687 Reims Cedex 2, France
| | - Andrei V Nikitin
- Laboratory of Theoretical Spectroscopy, V.E. Zuev Institute of Atmospheric Optics SB RAS 1, Akademician Zuev Sq., Tomsk 634055 Russia
| | - Dominika Viglaska
- Groupe de Spectrométrie Moléculaire et Atmosphérique UMR CNRS 7331, UFR Sciences BP 1039, 51687 Reims Cedex 2, France
| |
Collapse
|
6
|
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
|
7
|
Egorov O, Rey M, Nikitin AV, Viglaska D. New Ab Initio Potential Energy Surfaces for NH 3 Constructed from Explicitly Correlated Coupled-Cluster Methods. J Phys Chem A 2021; 125:10568-10579. [PMID: 34874738 DOI: 10.1021/acs.jpca.1c08717] [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/28/2022]
Abstract
Accurate calculation of ab initio potential energy surfaces (PESs) for the NH3 molecule is a difficult task because of the poor convergence of the standard CCSD(T) method with respect to the basis size. Currently, the best available NH3 PESs contain empirically refined parameters. In this paper, we show that CCSD(T)/aug-cc-pCV6Z calculations are not sufficient to properly describe the PES over a large range of nuclear configurations. However, the PES obtained by the extrapolation of the CCSD(T)/aug-cc-pCVXZ (X = T, Q, 5, and 6) energies to the complete basis set limit is closer to that based on the explicitly correlated CCSD(T)-F12a method using the cc-pCVQZ-F12 orbital basis set. All of the ab initio PESs constructed in this work include the following corrections: one electron relativistic effects, diagonal Born-Oppenheimer correction, and high-order electronic correlations (CCSDT, CCSDTQ, and CCSDTQP). Finally, the root-mean-square deviation between the predicted band centers obtained from our final "pure" ab initio PES and the experimental ones in the spectral region [0-7000] cm-1 is divided by two compared to the most accurate ab initio PES available in the literature.
Collapse
Affiliation(s)
- Oleg Egorov
- Laboratory of Quantum Mechanics of Molecules and Radiative Processes, Tomsk State University 36, Lenin Avenue, Tomsk 634050, Russia.,Laboratory of Theoretical Spectroscopy, V.E. Zuev Institute of Atmospheric Optics SB RAS, 1, Akademician Zuev Square, Tomsk 634055, Russia
| | - Michaël Rey
- Groupe de Spectrométrie Moléculaire et Atmosphérique, UMR CNRS 7331, UFR Sciences BP 1039, 51687 Reims Cedex 2, France
| | - Andrei V Nikitin
- Laboratory of Quantum Mechanics of Molecules and Radiative Processes, Tomsk State University 36, Lenin Avenue, Tomsk 634050, Russia.,Laboratory of Theoretical Spectroscopy, V.E. Zuev Institute of Atmospheric Optics SB RAS, 1, Akademician Zuev Square, Tomsk 634055, Russia
| | - Dominika Viglaska
- Groupe de Spectrométrie Moléculaire et Atmosphérique, UMR CNRS 7331, UFR Sciences BP 1039, 51687 Reims Cedex 2, France
| |
Collapse
|
8
|
Copeland RA, Zhou S, Schaperdoth I, Shoda TKC, Bollinger JM, Krebs C. Hybrid radical-polar pathway for excision of ethylene from 2-oxoglutarate by an iron oxygenase. Science 2021; 373:1489-1493. [PMID: 34385355 DOI: 10.1126/science.abj4290] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/02/2022]
Abstract
[Figure: see text].
Collapse
Affiliation(s)
- Rachelle A Copeland
- Department of Chemistry, The Pennsylvania State University, University Park, PA 16802, USA
| | - Shengbin Zhou
- Department of Chemistry, The Pennsylvania State University, University Park, PA 16802, USA
| | - Irene Schaperdoth
- Department of Chemistry, The Pennsylvania State University, University Park, PA 16802, USA
| | - Tokufu Kent C Shoda
- Department of Chemistry, The Pennsylvania State University, University Park, PA 16802, USA
| | - J Martin Bollinger
- Department of Chemistry, The Pennsylvania State University, University Park, PA 16802, USA.,Department of Biochemistry and Molecular Biology, The Pennsylvania State University, University Park, PA 16802, USA
| | - Carsten Krebs
- Department of Chemistry, The Pennsylvania State University, University Park, PA 16802, USA.,Department of Biochemistry and Molecular Biology, The Pennsylvania State University, University Park, PA 16802, USA
| |
Collapse
|
9
|
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
|
10
|
Furuya K, Sakamoto A, Tasumi M. Anharmonic Vibrational Analysis of Ethylene, Butadiene, and Longer Conjugated Alkenes by Double-Hybrid Density-Functional Theory Calculations. BULLETIN OF THE CHEMICAL SOCIETY OF JAPAN 2020. [DOI: 10.1246/bcsj.20200032] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Affiliation(s)
- Kazuhiko Furuya
- HPC SYSTEMS Inc., 8F LOOP-X, 3-9-15 Kaigan, Minato-ku, Tokyo 108-0022, Japan
| | - Akira Sakamoto
- College of Science and Engineering, Aoyama Gakuin University, 5-10-1 Fuchinobe, Chuo-ku, Sagamihara, Kanagawa 252-5258, Japan
| | - Mitsuo Tasumi
- Sapiarc Research Laboratory, 5-14-19-101 Roppongi, Minato-ku, Tokyo 106-0032, Japan
| |
Collapse
|