1
|
Faraji S, Liu M. Transferable machine learning interatomic potential for carbon hydrogen systems. Phys Chem Chem Phys 2024; 26:22346-22358. [PMID: 39140158 DOI: 10.1039/d4cp02300e] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/15/2024]
Abstract
In this study, we developed a machine learning interatomic potential based on artificial neural networks (ANN) to model carbon-hydrogen (C-H) systems. The ANN potential was trained on a dataset of C-H clusters obtained through density functional theory (DFT) calculations. Through comprehensive evaluations against DFT results, including predictions of geometries and formation energies across 0D-3D systems comprising C and C-H, as well as modeling various chemical processes, the ANN potential demonstrated exceptional accuracy and transferability. Its capability to accurately predict lattice dynamics, crucial for stability assessment in crystal structure prediction, was also verified through phonon dispersion analysis. Notably, its accuracy and computational efficiency in calculating force constants facilitated the exploration of complex energy landscapes, leading to the discovery of a novel C polymorph. These results underscore the robustness and versatility of the ANN potential, highlighting its efficacy in advancing computational materials science by conducting precise atomistic simulations on a wide range of C-H materials.
Collapse
Affiliation(s)
- Somayeh Faraji
- Department of Chemistry, University of Florida, Gainesville, FL 32611, USA.
| | - Mingjie Liu
- Department of Chemistry, University of Florida, Gainesville, FL 32611, USA.
| |
Collapse
|
2
|
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
|
3
|
Owens A. A highly accurate potential energy surface for carbonyl sulphide (OCS): how important are the ab initio calculations? Phys Chem Chem Phys 2024; 26:17684-17694. [PMID: 38869020 DOI: 10.1039/d4cp01205d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/14/2024]
Abstract
Ab initio quantum chemical methods can produce accurate molecular potential energy surfaces (PESs) capable of predicting the fundamental vibrational wavenumbers to within 1 cm-1. However, for high-resolution applications this is simply not good enough and empirical refinement is necessary, i.e. adjusting the PES to better match laboratory spectroscopic data. Here, the impact of the underlying ab initio calculations is rigorously investigated within the context of empirical refinement. For carbonyl sulphide (OCS), state-of-the-art electronic structure calculations are employed to construct higher- and lower-level ab initio PESs, which are then empirically refined in near-identical procedures. The initial ab initio calculations are shown to considerably affect the accuracy of the final refined PES, with an order-of-magnitude improvement in computed rotation-vibration energy levels achieved for OCS. In demonstrating this, the most accurate PES of the electronic ground state of OCS is produced, reproducing the fundamentals with a root-mean-square error (RMSE) of 0.004 cm-1, and 884 rovibrational energy levels below 14 000 cm-1 with an RMSE of 0.060 cm-1.
Collapse
Affiliation(s)
- Alec Owens
- Department of Physics and Astronomy, University College London, Gower Street, WC1E 6BT, London, UK.
| |
Collapse
|
4
|
Rey M, Viglaska D, Egorov O, Nikitin AV. A numerical-tensorial "hybrid" nuclear motion Hamiltonian and dipole moment operator for spectra calculation of polyatomic nonrigid molecules. J Chem Phys 2023; 159:114103. [PMID: 37712781 DOI: 10.1063/5.0166657] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/06/2023] [Accepted: 08/24/2023] [Indexed: 09/16/2023] Open
Abstract
The analysis and modeling of high-resolution spectra of nonrigid molecules require a specific Hamiltonian and group-theoretical formulation that differs significantly from that of more familiar rigid systems. Within the framework of Hougen-Bunker-Johns (HBJ) theory, this paper is devoted to the construction of a nonrigid Hamiltonian based on a suitable combination of numerical calculations for the nonrigid part in conjunction with the irreducible tensor operator method for the rigid part. For the first time, a variational calculation from ab initio potential energy surfaces is performed using the HBJ kinetic energy operator built from vibrational, large-amplitude motion, and rotational tensor operators expressed in terms of curvilinear and normal coordinates. Group theory for nonrigid molecules plays a central role in the characterization of the overall tunneling splittings and is discussed in the present approach. The construction of the dipole moment operator is also examined. Validation tests consisting of a careful convergence study of the energy levels as well as a comparison of results obtained from independent computer codes are given for the nonrigid molecules CH2, CH3, NH3, and H2O2. This work paves the way for the modeling of high-resolution spectra of larger nonrigid systems.
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
| | - Dominika Viglaska
- Groupe de Spectrométrie Moléculaire et Atmosphérique, UMR CNRS 7331, BP 1039, F-51687 Reims Cedex 2, France
| | - Oleg Egorov
- Laboratory of Theoretical Spectroscopy, V.E. Zuev Institute of Atmospheric Optics SB RAS, 1, Akademician Zuev Sq., Tomsk 634055, Russia
- Tomsk State University, 36, Lenin Ave., Tomsk 634050, Russia
| | - Andrei V Nikitin
- Laboratory of Theoretical Spectroscopy, V.E. Zuev Institute of Atmospheric Optics SB RAS, 1, Akademician Zuev Sq., Tomsk 634055, Russia
| |
Collapse
|
5
|
Tyuterev V, Tashkun S, Rey M, Nikitin A. High-order contact transformations of molecular Hamiltonians: general approach, fast computational algorithm and convergence of ro-vibrational polyad models. Mol Phys 2022. [DOI: 10.1080/00268976.2022.2096140] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/16/2022]
Affiliation(s)
- Vladimir Tyuterev
- Groupe de Spectrométrie Moléculaire et Atmosphérique, Université de Reims, Reims, France
- Laboratory of Molecular Quantum Mechanics and Radiative Transfer, Tomsk State University, Tomsk, Russia
| | - Sergey Tashkun
- Laboratory of Theoretical Spectroscopy, V.E. Zuev Institute of Atmospheric Optics, Siberian Branch, Russian Academy of Sciences, Tomsk, Russia
| | - Michael Rey
- Groupe de Spectrométrie Moléculaire et Atmosphérique, Université de Reims, Reims, France
| | - Andrei Nikitin
- Laboratory of Theoretical Spectroscopy, V.E. Zuev Institute of Atmospheric Optics, Siberian Branch, Russian Academy of Sciences, Tomsk, Russia
| |
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
|
Zhang L, Zhang S, Owens A, Yurchenko SN, Dral PO. VIB5 database with accurate ab initio quantum chemical molecular potential energy surfaces. Sci Data 2022; 9:84. [PMID: 35277513 PMCID: PMC8917215 DOI: 10.1038/s41597-022-01185-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/02/2021] [Accepted: 01/19/2022] [Indexed: 11/09/2022] Open
Abstract
High-level ab initio quantum chemical (QC) molecular potential energy surfaces (PESs) are crucial for accurately simulating molecular rotation-vibration spectra. Machine learning (ML) can help alleviate the cost of constructing such PESs, but requires access to the original ab initio PES data, namely potential energies computed on high-density grids of nuclear geometries. In this work, we present a new structured PES database called VIB5, which contains high-quality ab initio data on 5 small polyatomic molecules of astrophysical significance (CH3Cl, CH4, SiH4, CH3F, and NaOH). The VIB5 database is based on previously used PESs, which, however, are either publicly unavailable or lacking key information to make them suitable for ML applications. The VIB5 database provides tens of thousands of grid points for each molecule with theoretical best estimates of potential energies along with their constituent energy correction terms and a data-extraction script. In addition, new complementary QC calculations of energies and energy gradients have been performed to provide a consistent database, which, e.g., can be used for gradient-based ML methods. Measurement(s) | potential energy surfaces | Technology Type(s) | quantum chemistry computational methods |
Collapse
Affiliation(s)
- Lina Zhang
- State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, Department of Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, 361005, China
| | - Shuang Zhang
- State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, Department of Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, 361005, China
| | - Alec Owens
- Department of Physics and Astronomy, University College London, Gower Street, WC1E 6BT, London, United Kingdom.
| | - Sergei N Yurchenko
- Department of Physics and Astronomy, University College London, Gower Street, WC1E 6BT, London, United Kingdom
| | - Pavlo O Dral
- State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, Department of Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, 361005, China.
| |
Collapse
|
8
|
Abstract
Linear molecules usually represent a special case in rotational-vibrational calculations due to a singularity of the kinetic energy operator that arises from the rotation about the a (the principal axis of least moment of inertia, becoming the molecular axis at the linear equilibrium geometry) being undefined. Assuming the standard ro-vibrational basis functions, in the 3N−6 approach, of the form ∣ν1,ν2,ν3ℓ3;J,k,m⟩, tackling the unique difficulties of linear molecules involves constraining the vibrational and rotational functions with k=ℓ3, which are the projections, in units of ℏ, of the corresponding angular momenta onto the molecular axis. These basis functions are assigned to irreducible representations (irreps) of the C2v(M) molecular symmetry group. This, in turn, necessitates purpose-built codes that specifically deal with linear molecules. In the present work, we describe an alternative scheme and introduce an (artificial) group that ensures that the condition ℓ3=k is automatically applied solely through symmetry group algebra. The advantage of such an approach is that the application of symmetry group algebra in ro-vibrational calculations is ubiquitous, and so this method can be used to enable ro-vibrational calculations of linear molecules in polyatomic codes with fairly minimal modifications. To this end, we construct a—formally infinite—artificial molecular symmetry group D∞h(AEM), which consists of one-dimensional (non-degenerate) irreducible representations and use it to classify vibrational and rotational basis functions according to ℓ and k. This extension to non-rigorous, artificial symmetry groups is based on cyclic groups of prime-order. Opposite to the usual scenario, where the form of symmetry adapted basis sets is dictated by the symmetry group the molecule belongs to, here the symmetry group D∞h(AEM) is built to satisfy properties for the convenience of the basis set construction and matrix elements calculations. We believe that the idea of purpose-built artificial symmetry groups can be useful in other applications.
Collapse
|
9
|
Tajti A, Szalay PG, Kochanov R, Tyuterev VG. Diagonal Born-Oppenheimer corrections to the ground electronic state potential energy surfaces of ozone: improvement of ab initio vibrational band centers for the 16O 3, 17O 3 and 18O 3 isotopologues. Phys Chem Chem Phys 2020; 22:24257-24269. [PMID: 33089270 DOI: 10.1039/d0cp02457k] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Mass-dependent diagonal Born-Oppenheimer corrections (DBOCs) to the ab initio electronic ground state potential energy surface for the main 16O3 isotopologue and for homogeneous isotopic substitutions 17O3 and 18O3 of the ozone molecule are reported for the first time. The system being of strongly multiconfigurational character, multireference configuration interaction wave function ansatz with different complete active spaces was used. The reliable DBOC calculations with the targeted accuracy were possible to carry out up to about half of the dissociation threshold D0. The comparison with the experimental band centers shows a significant improvement of the accuracy with respect to the best Born-Oppenheimer (BO) ab initio calculations reducing the total root-mean-squares (calculated-observed) deviations by about a factor of two. For the set of 16O3 vibrations up to five bending and four stretching quanta, the mean (calculated-observed) deviations drop down from 0.7 cm-1 (BO) to about 0.1 cm-1, with the most pronounced improvement seen for bending states and for mixed bending-stretching polyads. In the case of bending band centers directly observed under high spectral resolutions, the errors are reduced by more than an order of magnitude down to 0.02 cm-1 from the observed levels, approaching nearly experimental accuracy. A similar improvement for heavy isotopologues shows that the reported DBOC corrections almost remove the systematic BO errors in vibrational levels below D0/2, though the scatter increases towards higher energies. The possible reasons for this finding, as well as remaining issues are discussed in detail. The reported results provide an encouraging accuracy validation for the multireference methods of the ab initio theory. New sets of ab initio vibrational states can be used for improving effective spectroscopic models for analyses of the observed high-resolution spectra, particularly in the cases of accidental resonances with "dark" states requiring accurate theoretical predictions.
Collapse
Affiliation(s)
- Attila Tajti
- ELTE Eötvös Loránd University, Institute of Chemistry, Laboratory of Theoretical Chemistry, P. O. Box 32, H-1518, Budapest 112, Hungary.
| | | | | | | |
Collapse
|
10
|
Nikitin AV, Rey M, Chizhmakova IS, Tyuterev VG. First Full-Dimensional Potential Energy and Dipole Moment Surfaces of SF6. J Phys Chem A 2020; 124:7014-7023. [DOI: 10.1021/acs.jpca.0c02733] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Andrei V. Nikitin
- V.E. Zuev Institute of Atmospheric Optics, Russian Academy of Sciences, 1, Akademichesky Avenue, 634055 Tomsk, Russian Federation
| | - Michael Rey
- Groupe de Spectrométrie Moléculaire et Atmosphérique, UMR CNRS 6089, Université de Reims, U.F.R. Sciences, B.P. 1039, 51687 Reims Cedex 2, France
| | - Iana S. Chizhmakova
- Institute of Monitoring of Climatic and Ecological Systems, Russian Academy of Sciences, 10/3, Academichesky Avenue, 634055 Tomsk, Russian Federation
- QUAMER Laboratory, Tomsk State University, 36 Lenin Avenue, 634050 Tomsk, Russian Federation
| | - Vladimir G. Tyuterev
- Groupe de Spectrométrie Moléculaire et Atmosphérique, UMR CNRS 6089, Université de Reims, U.F.R. Sciences, B.P. 1039, 51687 Reims Cedex 2, France
- QUAMER Laboratory, Tomsk State University, 36 Lenin Avenue, 634050 Tomsk, Russian Federation
| |
Collapse
|
11
|
Dral PO, Owens A, Dral A, Csányi G. Hierarchical machine learning of potential energy surfaces. J Chem Phys 2020; 152:204110. [DOI: 10.1063/5.0006498] [Citation(s) in RCA: 31] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022] Open
Affiliation(s)
- Pavlo O. Dral
- State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, Department of Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, China
| | - Alec Owens
- Department of Physics and Astronomy, University College London, Gower Street, WC1E 6BT London, United Kingdom
| | - Alexey Dral
- BigData Team, 1A Tormoznoye Shosse Off 17, Yaroslavl, Yaroslavl 150022, Russian Federation
| | - Gábor Csányi
- Department of Engineering, University of Cambridge, Cambridge CB2 1PZ, United Kingdom
| |
Collapse
|
12
|
Affiliation(s)
- James H. Thorpe
- The Quantum Theory Project, Department of Chemistry, The University of Florida, Gainesville, FL, USA
| | - John F. Stanton
- The Quantum Theory Project, Department of Chemistry, The University of Florida, Gainesville, FL, USA
| |
Collapse
|
13
|
Sahnoun E, Wiesenfeld L, Hammami K. Potential Energy Surface for the CH4–H2 van der Waals Interaction. J Phys Chem A 2020; 124:3242-3248. [DOI: 10.1021/acs.jpca.9b10499] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Emna Sahnoun
- Laboratoire Aimé-Cotton, , CNRS, Université Paris-Saclay91405 Orsay, France
- LSAMA, Université de Tunis El Manar, Tunis, Tunisia
| | - Laurent Wiesenfeld
- Laboratoire Aimé-Cotton, , CNRS, Université Paris-Saclay91405 Orsay, France
| | | |
Collapse
|
14
|
Transformation Properties under the Operations of
the Molecular Symmetry Groups G36 and G36(EM) of
Ethane H3CCH3. Symmetry (Basel) 2019. [DOI: 10.3390/sym11070862] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022] Open
Abstract
In the present work, we report a detailed description of the symmetry propertiesof the eight-atomic molecule ethane, with the aim of facilitating the variational calculations ofrotation-vibration spectra of ethane and related molecules. Ethane consists of two methyl groupsCH3 where the internal rotation (torsion) of one CH3 group relative to the other is of large amplitudeand involves tunnelling between multiple minima of the potential energy function. The molecularsymmetry group of ethane is the 36-element group G36, but the construction of symmetrised basisfunctions is most conveniently done in terms of the 72-element extended molecular symmetrygroup G36(EM). This group can subsequently be used in the construction of block-diagonal matrixrepresentations of the ro-vibrational Hamiltonian for ethane. The derived transformation matricesassociated with G36(EM) have been implemented in the variational nuclear motion program TROVE(Theoretical ROVibrational Energies). TROVE variational calculations are used as a practical exampleof a G36(EM) symmetry adaptation for large systems with a non-rigid, torsional degree of freedom.We present the derivation of irreducible transformation matrices for all 36 (72) operations of G36(M)(G36(EM)) and also describe algorithms for a numerical construction of these matrices based on aset of four (five) generators. The methodology presented is illustrated on the construction of thesymmetry-adapted representations both of the potential energy function of ethane and of the rotation,torsion and vibration basis set functions.
Collapse
|
15
|
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
|
16
|
Owens A, Yurchenko SN. Theoretical rotation-vibration spectroscopy of cis- and trans-diphosphene (P 2H 2) and the deuterated species P 2HD. J Chem Phys 2019; 150:194308. [PMID: 31117776 DOI: 10.1063/1.5092767] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Growing astronomical interest in phosphorous (P) chemistry is stimulating the search for new interstellar P-bearing molecules, a task requiring detailed knowledge of the microwave and infrared molecular spectrum. In this work, we present comprehensive rotation-vibration line lists of the cis- and trans-isomers of diphosphene (P2H2). The line lists have been generated using robust, first-principles methodologies based on newly computed, high-level ab initio potential energy and dipole moment surfaces. Transitions are considered between states with energies up to 8000 cm-1 and total angular momentum J ≤ 25. These are the first-ever line lists to be reported for P2H2, and they should significantly facilitate future spectroscopic characterization of this system. The deuterated species trans-P2HD and the effect of its dynamic dipole moment on the rovibrational spectrum are also discussed.
Collapse
Affiliation(s)
- Alec Owens
- Department of Physics and Astronomy, University College London, Gower Street, WC1E 6BT London, United Kingdom
| | - Sergei N Yurchenko
- Department of Physics and Astronomy, University College London, Gower Street, WC1E 6BT London, United Kingdom
| |
Collapse
|
17
|
Metz MP, Szalewicz K, Sarka J, Tóbiás R, Császár AG, Mátyus E. Molecular dimers of methane clathrates: ab initio potential energy surfaces and variational vibrational states. Phys Chem Chem Phys 2019; 21:13504-13525. [DOI: 10.1039/c9cp00993k] [Citation(s) in RCA: 23] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Motivated by the energetic and environmental relevance of methane clathrates, highly accurate ab initio potential energy surfaces (PESs) have been developed for the three possible dimers of the methane and water molecules: (H2O)2, CH4·H2O, and (CH4)2.
Collapse
Affiliation(s)
- Michael P. Metz
- Department of Physics and Astronomy
- University of Delaware
- Newark
- USA
| | | | - János Sarka
- Department of Chemistry and Biochemistry
- Texas Tech University
- Lubbock
- USA
| | - Roland Tóbiás
- Institute of Chemistry
- ELTE Eötvös Loránd University
- Budapest
- Hungary
- MTA-ELTE Complex Chemical Systems Research Group
| | - Attila G. Császár
- Institute of Chemistry
- ELTE Eötvös Loránd University
- Budapest
- Hungary
- MTA-ELTE Complex Chemical Systems Research Group
| | - Edit Mátyus
- Institute of Chemistry
- ELTE Eötvös Loránd University
- Budapest
- Hungary
| |
Collapse
|
18
|
Owens A, Yachmenev A, Küpper J, Yurchenko SN, Thiel W. The rotation–vibration spectrum of methyl fluoride from first principles. Phys Chem Chem Phys 2019; 21:3496-3505. [DOI: 10.1039/c8cp01721b] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
State-of-the-art theoretical methods are used to investigate the rotation–vibration spectrum of methyl fluoride.
Collapse
Affiliation(s)
- Alec Owens
- Center for Free-Electron Laser Science (CFEL)
- Deutsches Elektronen-Synchrotron DESY
- 22607 Hamburg
- Germany
- The Hamburg Center for Ultrafast Imaging
| | - Andrey Yachmenev
- Center for Free-Electron Laser Science (CFEL)
- Deutsches Elektronen-Synchrotron DESY
- 22607 Hamburg
- Germany
- The Hamburg Center for Ultrafast Imaging
| | - Jochen Küpper
- Center for Free-Electron Laser Science (CFEL)
- Deutsches Elektronen-Synchrotron DESY
- 22607 Hamburg
- Germany
- The Hamburg Center for Ultrafast Imaging
| | - Sergei N. Yurchenko
- Department of Physics and Astronomy
- University College London
- WC1E 6BT London
- UK
| | - Walter Thiel
- Max-Planck-Institut für Kohlenforschung
- 45470 Mülheim an der Ruhr
- Germany
| |
Collapse
|
19
|
Nikitin AV, Protasevich AE, Rey M, Tyuterev VG. Highly excited vibrational levels of methane up to 10 300 cm -1: Comparative study of variational methods. J Chem Phys 2018; 149:124305. [PMID: 30278662 DOI: 10.1063/1.5042154] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023] Open
Abstract
In this work, we report calculated vibrational energy levels of the methane molecule up to 10 300 cm-1. Two potential energy surfaces constructed in quite different coordinate systems with different analytical representations are employed in order to evaluate the uncertainty of vibrational predictions. To calculate methane energy levels, we used two independent techniques of the variational method. One method uses an exact kinetic energy operator in internal curvilinear coordinates. Another one uses an expansion of Eckart-Watson nuclear motion Hamiltonian in rectilinear normal coordinates. In the Icosad range (up to five vibrational quanta bands-below 7800 cm-1), the RMS standard deviations between calculated and observed energy levels were 0.22 cm-1 and 0.41 cm-1 for these two quite different approaches. For experimentally well-known 3v3 sub-levels, the calculation accuracy is estimated to be ∼1 cm-1. In the Triacontad range (7660-9188 cm-1), the average error of the calculation is about 0.5 cm-1. The accuracy and convergence issues for higher energy ranges are discussed.
Collapse
Affiliation(s)
- Andrei V Nikitin
- V.E. Zuev Institute of Atmospheric Optics, Russian Academy of Sciences, 1, Akademichesky Avenue, 634055 Tomsk, Russian Federation
| | - Alexander E Protasevich
- V.E. Zuev Institute of Atmospheric Optics, Russian Academy of Sciences, 1, Akademichesky Avenue, 634055 Tomsk, Russian Federation
| | - Michael Rey
- Groupe de Spectrométrie Moléculaire et Atmosphérique, UMR CNRS 7331, Université de Reims, U.F.R. Sciences, B.P. 1039, 51687 Reims Cedex 2, France
| | - Vladimir G Tyuterev
- Groupe de Spectrométrie Moléculaire et Atmosphérique, UMR CNRS 7331, Université de Reims, U.F.R. Sciences, B.P. 1039, 51687 Reims Cedex 2, France
| |
Collapse
|
20
|
Rey M, Chizhmakova IS, Nikitin AV, Tyuterev VG. Understanding global infrared opacity and hot bands of greenhouse molecules with low vibrational modes from first-principles calculations: the case of CF 4. Phys Chem Chem Phys 2018; 20:21008-21033. [PMID: 30070661 DOI: 10.1039/c8cp03252a] [Citation(s) in RCA: 29] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Fluorine containing molecules have a particularly long atmospheric lifetime and their very big estimated global warming potentials are expected to rapidly increase in the future. This work is focused on the global theoretical prediction of infrared spectra of the tetrafluoromethane molecule that is considered as a potentially powerful greenhouse gas having the largest estimated lifetime of over 50 000 years in the atmosphere. The presence of relatively low vibrational frequencies makes the Boltzmann population of the excited levels important. Consequently, the "hot bands" corresponding to transitions among excited rovibrational states contribute significantly to the CF4 opacity in the infrared even at room temperature conditions but the existing laboratory data analyses are not sufficiently complete. In this work, we construct the first accurate and complete ab initio based line lists for CF4 in the range 0-4000 cm-1, containing rovibrational bands that are the most active in absorption. An efficient basis set compression method was applied to predict more than 700 new bands and subbands via variational nuclear motion calculations. We show that already at room temperature a quasi-continuum of overlapping weak lines appears in the CF4 infrared spectra due to the increasing density of bands and transitions. In order to converge the infrared opacity at room temperature, it was necessary to include a high rotational quantum number up to J = 80 resulting in 2 billion rovibrational transitions. In order to make the cross-section simulation faster, we have partitioned our data into two parts: (a) strong & medium line lists with lower energy levels for calculation of selective absorption features that can be used at various temperatures and (b) compressed "super-line" libraries of very weak transitions contributing to the quasi-continuum modelling. Comparisons with raw previously unassigned experimental spectra showed a very good accuracy for integrated absorbance in the entire range of the reported spectra predictions. The data obtained in this work will be made available through the TheoReTS information system (http://theorets.univ-reims.fr, http://theorets.tsu.ru) that contains ab initio born line lists and provides a user-friendly graphical interface for a fast simulation of the CF4 absorption cross-sections and radiance under various temperature conditions from 80 K to 400 K.
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.
| | | | | | | |
Collapse
|
21
|
Chubb KL, Yachmenev A, Tennyson J, Yurchenko SN. Treating linear molecule HCCH in calculations of rotation-vibration spectra. J Chem Phys 2018; 149:014101. [DOI: 10.1063/1.5031844] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Katy L. Chubb
- Department of Physics and Astronomy, University College London, London WC1E 6BT, United Kingdom
| | - Andrey Yachmenev
- Center for Free-Electron Laser Science (CFEL), Deutsches Elektronen-Synchrotron DESY, Notkestrasse 85, 22607 Hamburg,
Germany
| | - Jonathan Tennyson
- Department of Physics and Astronomy, University College London, London WC1E 6BT, United Kingdom
| | - Sergei N. Yurchenko
- Department of Physics and Astronomy, University College London, London WC1E 6BT, United Kingdom
| |
Collapse
|
22
|
Zhao Z, Chen J, Zhang Z, Zhang DH, Wang XG, Carrington T, Gatti F. Computing energy levels of CH4, CHD3, CH3D, and CH3F with a direct product basis and coordinates based on the methyl subsystem. J Chem Phys 2018; 148:074113. [DOI: 10.1063/1.5019323] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Zhiqiang Zhao
- 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
| | - Jun Chen
- 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
- Collaborative Innovation Center of Chemistry for Energy Materials (iChEM), College of Chemistry and Chemical Engineering, Xiamen University, Siming South Road 422, 361005 Xiamen, China
| | - Zhaojun Zhang
- 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
| | - Dong H. Zhang
- 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
| | - 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
| | - Fabien Gatti
- ISMO, Institut des Sciences Moléculaires d’Orsay - UMR 8214 CNRS/Université Paris-Sud, F-91405 Orsay, France
| |
Collapse
|
23
|
Accurate Theoretical Methane Line Lists in the Infrared up to 3000 K and Quasi-continuum Absorption/Emission Modeling for Astrophysical Applications. ACTA ACUST UNITED AC 2017. [DOI: 10.3847/1538-4357/aa8909] [Citation(s) in RCA: 55] [Impact Index Per Article: 7.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
|
24
|
Protasevich AE, Nikitin AV. Matrix elements of vibration kinetic energy operator of tetrahedral molecules in non-orthogonal-dependent coordinates. Mol Phys 2017. [DOI: 10.1080/00268976.2017.1366568] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
Affiliation(s)
- Alexander E. Protasevich
- Laboratory of Atmospheric Absorption Spectroscopy, V.E. Zuev Institute of Atmospheric Optics, SB RAS, 1, Academician Zuev square, 634021, Tomsk, Russian Federation
| | - Andrei V. Nikitin
- Laboratory of Theoretical Spectroscopy, V.E. Zuev Institute of Atmospheric Optics, SB RAS, 1, Academician Zuev square, 634021, Tomsk, Russian Federation
| |
Collapse
|
25
|
Yurchenko SN, Yachmenev A, Ovsyannikov RI. Symmetry-Adapted Ro-vibrational Basis Functions for Variational Nuclear Motion Calculations: TROVE Approach. J Chem Theory Comput 2017; 13:4368-4381. [DOI: 10.1021/acs.jctc.7b00506] [Citation(s) in RCA: 49] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Sergei N. Yurchenko
- Department
of Physics and Astronomy, University College London, London, WC1E 6BT, United Kingdom
| | - Andrey Yachmenev
- Center for Free-Electron Laser Science (CFEL), DESY, Notkestrasse 85, 22607 Hamburg, Germany
| | - Roman I. Ovsyannikov
- Institute of Applied Physics, Russian Academy of Sciences, Ulyanov Street 46, Nizhny Novgorod, Russia 603950
| |
Collapse
|