1
|
Ge F, Wang R, Qu C, Zheng P, Nandi A, Conte R, Houston PL, Bowman JM, Dral PO. Tell Machine Learning Potentials What They Are Needed For: Simulation-Oriented Training Exemplified for Glycine. J Phys Chem Lett 2024; 15:4451-4460. [PMID: 38626460 DOI: 10.1021/acs.jpclett.4c00746] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/18/2024]
Abstract
Machine learning potentials (MLPs) are widely applied as an efficient alternative way to represent potential energy surfaces (PESs) in many chemical simulations. The MLPs are often evaluated with the root-mean-square errors on the test set drawn from the same distribution as the training data. Here, we systematically investigate the relationship between such test errors and the simulation accuracy with MLPs on an example of a full-dimensional, global PES for the glycine amino acid. Our results show that the errors in the test set do not unambiguously reflect the MLP performance in different simulation tasks, such as relative conformer energies, barriers, vibrational levels, and zero-point vibrational energies. We also offer an easily accessible solution for improving the MLP quality in a simulation-oriented manner, yielding the most precise relative conformer energies and barriers. This solution also passed the stringent test by diffusion Monte Carlo simulations.
Collapse
Affiliation(s)
- Fuchun Ge
- State Key Laboratory of Physical Chemistry of Solid Surfaces, College of Chemistry and Chemical Engineering, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and Innovation Laboratory for Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen, Fujian 361005, China
| | - Ran Wang
- State Key Laboratory of Physical Chemistry of Solid Surfaces, College of Chemistry and Chemical Engineering, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and Innovation Laboratory for Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen, Fujian 361005, China
| | - Chen Qu
- Independent Researcher, Toronto, Ontario M9B0E3, Canada
| | - Peikun Zheng
- State Key Laboratory of Physical Chemistry of Solid Surfaces, College of Chemistry and Chemical Engineering, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and Innovation Laboratory for Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen, Fujian 361005, China
| | - Apurba Nandi
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
- Department of Physics and Materials Science, University of Luxembourg, Luxembourg City L-1511, Luxembourg
| | - Riccardo Conte
- Dipartimento di Chimica, Università degli Studi di Milano, via Golgi 19, 20133 Milano, Italy
| | - Paul L Houston
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States
| | - Joel M Bowman
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| | - Pavlo O Dral
- State Key Laboratory of Physical Chemistry of Solid Surfaces, College of Chemistry and Chemical Engineering, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and Innovation Laboratory for Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen, Fujian 361005, China
| |
Collapse
|
2
|
Barbe A, Mikhailenko S, Starikova E, Tyuterev V. High Resolution Infrared Spectroscopy in Support of Ozone Atmospheric Monitoring and Validation of the Potential Energy Function. Molecules 2022; 27:911. [PMID: 35164172 PMCID: PMC8838290 DOI: 10.3390/molecules27030911] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/30/2021] [Revised: 01/21/2022] [Accepted: 01/22/2022] [Indexed: 11/16/2022] Open
Abstract
The first part of this review is a brief reminder of general information concerning atmospheric ozone, particularly related to its formation, destruction, observations of its decrease in the stratosphere, and its increase in the troposphere as a result of anthropogenic actions and solutions. A few words are said about the abandonment of the Airbus project Alliance, which was expected to be the substitute of the supersonic Concorde. This project is over due to the theoretical evaluation of the impact of a fleet in the stratosphere and has been replaced by the A380, which is now operating. The largest part is devoted to calculations and observations of the transitions in the infrared range and their applications for the atmosphere based both on effective models (Hamiltonian, symmetry rules, and dipole moments) and ab initio calculations. The complementarities of the two approaches are clearly demonstrated, particularly for the creation of an exhaustive line list consisting of more than 300,000 lines reaching experimental accuracies (from 0.00004 to 0.001 cm-1) for positions and a sub percent for the intensities in the 10 microns region. This contributes to definitively resolving the issue of the observed discrepancies between line intensity data in different spectral regions: between the infrared and ultraviolet ranges, on the one hand, and between 10 and 5 microns on the other hand. The following section is devoted to the application of recent work to improve the knowledge about the behavior of potential function at high energies. A controversial issue related to the shape of the potential function in the transition state range near the dissociation is discussed.
Collapse
Affiliation(s)
- Alain Barbe
- Groupe de Spectrométrie Moléculaire et Atmosphérique, UMR CNRS 7331, Université de Reims, UFR Sciences Exactes et Naturelles, CEDEX02, BP 1039-51687 Reims, France;
| | - Semen Mikhailenko
- Laboratory of Theoretical Spectroscopy, V.E. Zuev Institute of Atmospheric Optics SB RAS, 634055 Tomsk, Russia; (S.M.); (E.S.)
- Climate and Environmental Physics Laboratory, Ural Federal University, 19, Mira av., 620002 Yekaterinburg, Russia
| | - Evgeniya Starikova
- Laboratory of Theoretical Spectroscopy, V.E. Zuev Institute of Atmospheric Optics SB RAS, 634055 Tomsk, Russia; (S.M.); (E.S.)
| | - Vladimir Tyuterev
- Groupe de Spectrométrie Moléculaire et Atmosphérique, UMR CNRS 7331, Université de Reims, UFR Sciences Exactes et Naturelles, CEDEX02, BP 1039-51687 Reims, France;
- Laboratory of Theoretical Spectroscopy, V.E. Zuev Institute of Atmospheric Optics SB RAS, 634055 Tomsk, Russia; (S.M.); (E.S.)
- Laboratory of Quantum Mechanics of Molecules and Radiative Processes, Tomsk State University, 634050 Tomsk, Russia
| |
Collapse
|
3
|
Hui WCH, Lemke KH. The ozone–water complex: CCSD(T)/CBS structures and anharmonic vibrational spectroscopy of O 3(H 2O) n, ( n = 1 − 2). J Chem Phys 2020; 153:084302. [DOI: 10.1063/5.0015597] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Wallace C. H. Hui
- Department of Earth Sciences, University of Hong Kong, Pokfulam Road, Pok Fu Lam, Hong Kong, SAR
| | - Kono H. Lemke
- Department of Earth Sciences, University of Hong Kong, Pokfulam Road, Pok Fu Lam, Hong Kong, SAR
| |
Collapse
|
4
|
Tyuterev VG, Barbe A, Jacquemart D, Janssen C, Mikhailenko SN, Starikova EN. Ab initio predictions and laboratory validation for consistent ozone intensities in the MW, 10 and 5 μm ranges. J Chem Phys 2019; 150:184303. [DOI: 10.1063/1.5089134] [Citation(s) in RCA: 28] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Vl. G. Tyuterev
- Tomsk State Research University, TSU, Tomsk 634050, Russia
- GSMA UMR CNRS 7331, UFR Sciences, Université de Reims, BP 1039, 51687 Reims, France
| | - A. Barbe
- GSMA UMR CNRS 7331, UFR Sciences, Université de Reims, BP 1039, 51687 Reims, France
| | - D. Jacquemart
- MONARIS, Sorbonne Université, CNRS, 75252 Paris, France
| | - C. Janssen
- LERMA-IPSL, Sorbonne Université, Observatoire de Paris, Université PSL, CNRS, 75252 Paris, France
| | - S. N. Mikhailenko
- V.E. Zuev Institute of Atmospheric Optics, SB RAS, Tomsk 634055, Russia
| | - E. N. Starikova
- V.E. Zuev Institute of Atmospheric Optics, SB RAS, Tomsk 634055, Russia
| |
Collapse
|
5
|
Mankodi TK, Bhandarkar UV, Puranik BP. Dissociation cross section for high energy O 2-O 2 collisions. J Chem Phys 2018; 148:144305. [PMID: 29655354 DOI: 10.1063/1.5020125] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Collision-induced dissociation cross section database for high energy O2-O2 collisions (up to 30 eV) is generated and published using the quasiclassical trajectory method on the singlet, triplet, and quintet spin ground state O4 potential energy surfaces. At equilibrium conditions, these cross sections predict reaction rate coefficients that match those obtained experimentally. The main advantage of the cross section database based on ab initio computations is in the study of complex flows with high degree of non-equilibrium. Direct simulation Monte Carlo simulations using the reactive cross section databases are carried out for high enthalpy hypersonic oxygen flow over a cylinder at rarefied ambient conditions. A comparative study with the phenomenological total collision energy chemical model is also undertaken to point out the difference and advantage of the reported ab initio reaction model.
Collapse
Affiliation(s)
- T K Mankodi
- Department of Mechanical Engineering, Indian Institute of Technology Bombay, Mumbai, India
| | - U V Bhandarkar
- Department of Mechanical Engineering, Indian Institute of Technology Bombay, Mumbai, India
| | - B P Puranik
- Department of Mechanical Engineering, Indian Institute of Technology Bombay, Mumbai, India
| |
Collapse
|
6
|
Varga Z, Paukku Y, Truhlar DG. Potential energy surfaces for O + O2 collisions. J Chem Phys 2017; 147:154312. [PMID: 29055336 DOI: 10.1063/1.4997169] [Citation(s) in RCA: 55] [Impact Index Per Article: 7.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Zoltan Varga
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, USA
| | - Yuliya Paukku
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, USA
| | - Donald G. Truhlar
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, USA
| |
Collapse
|
7
|
Mankodi TK, Bhandarkar UV, Puranik BP. Dissociation cross sections for N 2 + N → 3N and O 2 + O → 3O using the QCT method. J Chem Phys 2017; 146:204307. [PMID: 28571362 DOI: 10.1063/1.4983813] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
Abstract
Cross sections for the homo-nuclear atom-diatom collision induced dissociations (CIDs): N2 + N and O2 + O are calculated using Quasi-Classical Trajectory (QCT) method on ab initio Potential Energy Surfaces (PESs). A number of studies for these reactions carried out in the past focused on the CID cross section values generated using London-Eyring-Polanyi-Sato PES and seldom listed the CID cross section data. A highly accurate CASSCF-CASPT2 N3 and a new O3 global PES are used for the present QCT analysis and the CID cross section data up to 30 eV relative energy are also published. In addition, an interpolating scheme based on spectroscopic data is introduced that fits the CID cross section for the entire ro-vibrational spectrum using QCT data generated at chosen ro-vibrational levels. The rate coefficients calculated using the generated CID cross section compare satisfactorily with the existing experimental and theoretical results. The CID cross section data generated will find an application in the development of a more precise chemical reaction model for Direct Simulation Monte Carlo code simulating hypersonic re-entry flows.
Collapse
Affiliation(s)
- Tapan K Mankodi
- Department of Mechanical Engineering, Indian Institute of Technology Bombay, Mumbai, India
| | - Upendra V Bhandarkar
- Department of Mechanical Engineering, Indian Institute of Technology Bombay, Mumbai, India
| | - Bhalchandra P Puranik
- Department of Mechanical Engineering, Indian Institute of Technology Bombay, Mumbai, India
| |
Collapse
|
8
|
Tyuterev VG, Kochanov RV, Tashkun SA. Accurateab initiodipole moment surfaces of ozone: First principle intensity predictions for rotationally resolved spectra in a large range of overtone and combination bands. J Chem Phys 2017; 146:064304. [DOI: 10.1063/1.4973977] [Citation(s) in RCA: 41] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
9
|
Xu L, Wang F, Huang J, Yang C, Yu L, Fan Y. l -Proline and thiourea co-catalyzed condensation of acetone. Tetrahedron 2016. [DOI: 10.1016/j.tet.2016.05.039] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
|
10
|
Dawes R, Lolur P, Li A, Jiang B, Guo H. Communication: An accurate global potential energy surface for the ground electronic state of ozone. J Chem Phys 2013; 139:201103. [DOI: 10.1063/1.4837175] [Citation(s) in RCA: 93] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
11
|
Towards insight into properties and stabilities of complexes of ozone with CO2, CS2 and SCO species. Struct Chem 2013. [DOI: 10.1007/s11224-013-0365-3] [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]
|
12
|
Tyuterev VG, Kochanov RV, Tashkun SA, Holka F, Szalay PG. New analytical model for the ozone electronic ground state potential surface and accurate ab initio vibrational predictions at high energy range. J Chem Phys 2013; 139:134307. [DOI: 10.1063/1.4821638] [Citation(s) in RCA: 91] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
13
|
Bozkaya U. Orbital-optimized third-order Møller-Plesset perturbation theory and its spin-component and spin-opposite scaled variants: Application to symmetry breaking problems. J Chem Phys 2011; 135:224103. [DOI: 10.1063/1.3665134] [Citation(s) in RCA: 51] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
14
|
Le HM, Dinh TS, Le HV. Molecular Dynamics Investigations of Ozone on an Ab Initio Potential Energy Surface with the Utilization of Pattern-Recognition Neural Network for Accurate Determination of Product Formation. J Phys Chem A 2011; 115:10862-70. [DOI: 10.1021/jp206531s] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Hung M. Le
- Faculty of Materials Science, College of Science, Vietnam National University, Ho Chi Minh City, Vietnam 750000
| | - Thach S. Dinh
- Faculty of Materials Science, College of Science, Vietnam National University, Ho Chi Minh City, Vietnam 750000
| | - Hieu V. Le
- Faculty of Materials Science, College of Science, Vietnam National University, Ho Chi Minh City, Vietnam 750000
| |
Collapse
|
15
|
Dawes R, Lolur P, Ma J, Guo H. Communication: Highly accurate ozone formation potential and implications for kinetics. J Chem Phys 2011; 135:081102. [DOI: 10.1063/1.3632055] [Citation(s) in RCA: 78] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
16
|
Han H, Suo B, Xie D, Lei Y, Wang Y, Wen Z. Electronic structure calculations of low-lying electronic states of O3. Phys Chem Chem Phys 2011; 13:2723-31. [DOI: 10.1039/c0cp01300e] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
17
|
Holka F, Szalay PG, Müller T, Tyuterev VG. Toward an Improved Ground State Potential Energy Surface of Ozone. J Phys Chem A 2010; 114:9927-35. [DOI: 10.1021/jp104182q] [Citation(s) in RCA: 58] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Filip Holka
- Faculty of Materials Science and Technology, Slovak University of Technology, Paulínska 16, 917 24 Trnava, Slovak Republic, Laboratory of Theoretical Chemistry, Institute of Chemistry,Eötvös University, P.O. Box 32, H-1518, Budapest 112, Hungary, Institute of Advanced Simulation, Jülich Supercomputer Centre, Research Centre Jülich GmbH D-52425 Jülich, Germany, and Groupe de Spectromé trie Moléculaire et Atmosphérique, Reims University, F-51687, Reims Cedex 2, France
| | - Péter G. Szalay
- Faculty of Materials Science and Technology, Slovak University of Technology, Paulínska 16, 917 24 Trnava, Slovak Republic, Laboratory of Theoretical Chemistry, Institute of Chemistry,Eötvös University, P.O. Box 32, H-1518, Budapest 112, Hungary, Institute of Advanced Simulation, Jülich Supercomputer Centre, Research Centre Jülich GmbH D-52425 Jülich, Germany, and Groupe de Spectromé trie Moléculaire et Atmosphérique, Reims University, F-51687, Reims Cedex 2, France
| | - Thomas Müller
- Faculty of Materials Science and Technology, Slovak University of Technology, Paulínska 16, 917 24 Trnava, Slovak Republic, Laboratory of Theoretical Chemistry, Institute of Chemistry,Eötvös University, P.O. Box 32, H-1518, Budapest 112, Hungary, Institute of Advanced Simulation, Jülich Supercomputer Centre, Research Centre Jülich GmbH D-52425 Jülich, Germany, and Groupe de Spectromé trie Moléculaire et Atmosphérique, Reims University, F-51687, Reims Cedex 2, France
| | - Vladimir G. Tyuterev
- Faculty of Materials Science and Technology, Slovak University of Technology, Paulínska 16, 917 24 Trnava, Slovak Republic, Laboratory of Theoretical Chemistry, Institute of Chemistry,Eötvös University, P.O. Box 32, H-1518, Budapest 112, Hungary, Institute of Advanced Simulation, Jülich Supercomputer Centre, Research Centre Jülich GmbH D-52425 Jülich, Germany, and Groupe de Spectromé trie Moléculaire et Atmosphérique, Reims University, F-51687, Reims Cedex 2, France
| |
Collapse
|
18
|
Abstract
The ground and low-lying states of ozone (O(3)) have been studied by multireference variational methods and large basis sets. We have constructed potential energy curves along the bending coordinate for (1,2) (1)A('), (1,2) (1)A("), (1,2) (3)A('), and (1,2) (3)A(") symmetries, optimizing at the same time the symmetric stretching coordinate. Thirteen minima have been located whose geometrical and energetic characteristics are in very good agreement with existing experimental data. Special emphasis has been given to the interpretation of the chemical bond through valence-bond-Lewis diagrams; their appropriate use captures admirably the bonding nature of the O(3) molecule. The biradical character of its ground state, adopted long ago by the scientific community, does not follow from a careful analysis of its wave function.
Collapse
Affiliation(s)
- Apostolos Kalemos
- Laboratory of Physical Chemistry, Department of Chemistry, National and Kapodistrian University of Athens, P.O. Box 64 004, Zografou, Athens 157 10, Greece.
| | | |
Collapse
|
19
|
Rosmus P, Linguerri R, Komiha N. First-principle computations of rotational-vibrational transition probabilities. Mol Phys 2008. [DOI: 10.1080/00268970802054040] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
|
20
|
Campargue A, Barbe A, De Backer-Barilly MR, Tyuterev VG, Kassi S. The near infrared spectrum of ozone by CW-cavity ring down spectroscopy between 5850 and 7000 cm−1: new observations and exhaustive review. Phys Chem Chem Phys 2008; 10:2925-46. [DOI: 10.1039/b719773j] [Citation(s) in RCA: 64] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
21
|
Zúñiga J, Picón JAG, Bastida A, Requena A. Optimal internal coordinates, vibrational spectrum, and effective Hamiltonian for ozone. J Chem Phys 2007; 126:244305. [PMID: 17614547 DOI: 10.1063/1.2743441] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
In this paper the authors use the optimal internal vibrational coordinates previously determined for the electronic ground state of the ozone molecule to study the vibrational spectrum of the molecule employing the second empirical potential energy surface calculated by Tyuterev et al. [Chem. Phys. Lett. 316, 271 (2000)]. First, the authors compute variationally all the bound vibrational energy levels of the molecule up to the dissociation limit and state the usefulness of the optimal coordinates in this respect, which allows us to converge all the bound levels using relatively small anharmonic basis sets. By analyzing the expansion coefficients of the wave functions, they show then that a large portion of the vibrational spectrum of O3 can be structured in nearly separable polyadic groups characterized by the polyad quantum number N=n1+n2+n(theta) corresponding to the optimal internal coordinates. Accordingly, they determine an internal effective vibrational Hamiltonian for O3 by fitting the effective Hamiltonian parameters to the experimental vibrational frequencies, using as input parameters in the fit those extracted from an analytical second-order Van Vleck perturbation theory calculation. It is finally shown that the internal effective Hamiltonian thus obtained accurately describes the vibrational spectrum of ozone in the low and medium energy regimes.
Collapse
Affiliation(s)
- José Zúñiga
- Departamento de Química Física, Universidad de Murcia, 30100 Murcia, Spain.
| | | | | | | |
Collapse
|
22
|
Evangelista FA, Allen WD, Schaefer HF. Coupling term derivation and general implementation of state-specific multireference coupled cluster theories. J Chem Phys 2007; 127:024102. [PMID: 17640114 DOI: 10.1063/1.2743014] [Citation(s) in RCA: 242] [Impact Index Per Article: 14.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Simple closed-form expressions are derived for the "same vacuum" renormalization terms that arise in state-specific multireference coupled cluster (MRCC) theories. Explicit equations are provided for these coupling terms through the triple excitation level of MRCC theory, and a general expression is included for arbitrary-order excitations. The first production-level code (PSIMRCC) for state-specific and rigorously size-extensive Mukherjee multireference coupled cluster singles and doubles (MkCCSD) computations has been written. This code is also capable of evaluating analogous Brillouin-Wigner multireference energies (BWCCSD), including a posteriori size-extensivity corrections. Using correlation-consistent basis sets (cc-pVXZ, X=D,T,Q), MkCCSD and BWCCSD were tested and compared on two classic multireference problems: (1) the dissociation potential curve of molecular fluorine (F(2)) and (2) the structure and vibrational frequencies of ozone. Comparison with experimental data shows that the Mukherjee method is generally superior to the Brillouin-Wigner theory in predicting energies, structures, and vibrational frequencies. Particularly accurate results for F(2) are obtained by applying the MkCCSD method with localized molecular orbitals. Although the MkCCSD theory greatly improves upon single-reference CCSD for the geometric parameters and a(1) vibrational frequencies of ozone, the antisymmetric stretching frequency omega(3)(b(2)) remains pathological and cannot be properly treated without the inclusion of connected triple excitations. Finally, preliminary multireference MkCCSD results are reported for the singlet-triplet splittings in ortho-, meta-, and para-benzyne, coming within 1.5 kcal mol(-1) of experiment in all cases.
Collapse
|
23
|
|
24
|
Grebenshchikov SY, Schinke R, Qu ZW, Zhu H. Absorption spectrum and assignment of the Chappuis band of ozone. J Chem Phys 2006; 124:204313. [PMID: 16774338 DOI: 10.1063/1.2196881] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
New global diabatic potential energy surfaces of the electronic states 1B1 and 1A2 of ozone and the non-adiabatic coupling surface between them are constructed from electronic structure calculations. These surfaces are used to study the visible photodissociation in the Chappuis band by means of quantum mechanical calculations. The calculated absorption spectrum and its absolute intensity are in good agreement with the experimental results. A vibrational assignment of the diffuse structures in the Chappuis band system is proposed on the basis of the nodal structures of the underlying resonance states.
Collapse
Affiliation(s)
- S Yu Grebenshchikov
- Max-Planck-Institut für Dynamik und Selbstorganisation, D-37073 Göttingen, Germany.
| | | | | | | |
Collapse
|
25
|
Schinke R, Grebenshchikov SY, Ivanov MV, Fleurat-Lessard P. DYNAMICAL STUDIES OF THE OZONE ISOTOPE EFFECT: A Status Report. Annu Rev Phys Chem 2006; 57:625-61. [PMID: 16599823 DOI: 10.1146/annurev.physchem.57.032905.104542] [Citation(s) in RCA: 85] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
▪ Abstract Dynamical studies of the recombination of O and O2 to form ozone are reviewed. The focus is the intriguing isotope dependence of the recombination rate coefficient as observed by Mauersberger and coworkers in the last decade. The key quantity for understanding of this dependence appears to be the difference of zero-point energies of the two fragmentation channels to which excited ozone can dissociate, i.e., X + YZ ← XYZ* → XY + Z, where X, Y, and Z stand for the three isotopes of oxygen. Besides the isotope dependence, the variation of the recombination rate coefficient with pressure and temperature is also addressed. Despite the numerous approaches of recent years, the recombination of ozone is far from being satisfactorily explained; there are still several essential questions to be solved by detailed theoretical analysis. We mainly discuss—and critically assess—the results of our own investigations of the ozone kinetics. The work of other research groups is also evaluated.
Collapse
Affiliation(s)
- R Schinke
- Max-Planck-Institut für Dynamik und Selbstorganisation, D-37073 Göttingen, Germany.
| | | | | | | |
Collapse
|
26
|
Hino O, Kinoshita T, Chan GKL, Bartlett RJ. Tailored coupled cluster singles and doubles method applied to calculations on molecular structure and harmonic vibrational frequencies of ozone. J Chem Phys 2006; 124:114311. [PMID: 16555892 DOI: 10.1063/1.2180775] [Citation(s) in RCA: 70] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
To assess the separation of dynamic and nondynamic correlations and orbital choice, we calculate the molecular structure and harmonic vibrational frequencies of ozone with the recently developed tailored coupled cluster singles and doubles method (TCCSD). We employ the Hartree-Fock and complete active space (CAS) self-consistent field (SCF) orbitals to perform TCCSD calculations. When using the Hartree-Fock orbitals, it is difficult to reproduce the experimental vibrational frequency of the asymmetric stretching mode. On the other hand, the TCCSD based on the CASSCF orbitals in a correlation consistent polarized valence triple zeta basis yields excellent results with the two symmetric vibrations differing from the experimental harmonic values by 2 cm(-1) and the asymmetric vibration differing by 9 cm(-1).
Collapse
Affiliation(s)
- Osamu Hino
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853-1301, USA.
| | | | | | | |
Collapse
|
27
|
Lin SY, Guo H. Quantum Statistical Study of O + O2 Isotopic Exchange Reactions: Cross Sections and Rate Constants. J Phys Chem A 2005; 110:5305-11. [PMID: 16623456 DOI: 10.1021/jp0556299] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Using a wave packet based statistical model, we compute cross sections and thermal rate constants for various isotopic variants of the O + O2 exchange reaction on a recently modified ab initio potential energy surface. The calculation predicts a highly excited rotational distribution and relatively cold vibrational distribution for the diatomic product. A small but important threshold effect was identified for the (16)O + 18O2 reaction, which is suggested to contribute to the experimentally observed negative temperature dependence of the rate ratio, k(18O + 16O2)/k(16O + 18O2). Despite reasonable agreement with quasiclassical trajectory results, however, the calculated thermal rate constants are smaller than experimental measurements by a factor from 2 to 5. The experimentally observed negative temperature dependence of the rate constants is not reproduced. Possible reasons for the theory-experiment discrepancies are discussed.
Collapse
Affiliation(s)
- Shi Ying Lin
- Department of Chemistry, University of New Mexico, Albuquerque, New Mexico 87131, USA
| | | |
Collapse
|
28
|
Elliott R, Compton R, Levis R, Matsika S. Excited Electronic States of the Cyclic Isomers of O3 and SO2. J Phys Chem A 2005; 109:11304-11. [PMID: 16331915 DOI: 10.1021/jp054455s] [Citation(s) in RCA: 18] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
The low-lying electronic states of O3 and SO2 in their bent and cyclic isomers up to about 10 eV are calculated using the multireference configuration interaction (MRCI) method with a standard Gaussian correlation consistent polarized triple-zeta (cc-pVTZ) basis set. The vertical excitation energies, electron configurations, and oscillator strengths of these states are reported. The molecular orbital structures and excited states of the cyclic isomers are discussed in relation to the bent ones. Coherent anti-Stokes Raman spectroscopy (CARS) schemes for detecting the synthesis of the cyclic isomers are suggested.
Collapse
Affiliation(s)
- Ruth Elliott
- Department of Chemistry and Center for Advanced Photonics Research, Temple University, Philadelphia, Pennsylvania 19122, USA
| | | | | | | |
Collapse
|
29
|
Qu ZW, Zhu H, Grebenshchikov SY, Schinke R. The photodissociation of ozone in the Hartley band: A theoretical analysis. J Chem Phys 2005; 123:074305. [PMID: 16229568 DOI: 10.1063/1.2001650] [Citation(s) in RCA: 40] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Three-dimensional diabatic potential energy surfaces for the lowest four electronic states of ozone with 1A' symmetry-termed X, A, B, and R-are constructed from electronic structure calculations. The diabatization is performed by reassigning corresponding energy points. Although approximate, these diabatic potential energy surfaces allow one to study the uv photodissociation of ozone on a level of theory not possible before. In the present work photoexcitation in the Hartley band and subsequent dissociation into the singlet channel, O3X+hnu-->O(1D)+O2(a 1Deltag), are investigated by means of quantum mechanical and classical trajectory calculations using the diabatic potential energy surface of the B state. The calculated low-resolution absorption spectrum as well as the vibrational and rotational state distributions of O2(a 1Deltag) are in good agreement with available experimental results.
Collapse
Affiliation(s)
- Z-W Qu
- Max-Planck-Institut für Dynamik und Selbstorganisation, D-37073 Göttingen, Germany
| | | | | | | |
Collapse
|
30
|
Baloïtcha E, Balint-Kurti GG. Theory of the photodissociation of ozone in the Hartley continuum: Potential energy surfaces, conical intersections, and photodissociation dynamics. J Chem Phys 2005; 123:014306. [PMID: 16035834 DOI: 10.1063/1.1903947] [Citation(s) in RCA: 42] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Ab initio potential energy and transition dipole moment surfaces are presented for the five lowest singlet even symmetry electronic states of ozone. The surfaces are calculated using the complete active space self consistent field method followed by contracted multireference configuration interaction (MRCI) calculations. A slightly reduced augmented correlation consistent valence triple-zeta orbital basis set is used. The ground and excited state energies of the molecule have been computed at 9282 separate nuclear geometries. Cuts through the potential energy surfaces, which pass through the geometry of the minimum of the ground electronic state, show several closely avoided crossings. Close examination, and higher level calculations, very strongly suggests that some of these seemingly avoided crossings are in fact associated with non-symmetry related conical intersections. Diabatic potential energy and transition dipole moment surfaces are created from the computed ab initio adiabatic MRCI energies and transition dipole moments. The transition dipole moment connecting the ground electronic state to the diabatic B state surface is by far the strongest. Vibrational-rotational wavefunctions and energies are computed using the ground electronic state. The energy level separations compare well with experimentally determined values. The ground vibrational state wavefunction is then used, together with the diabatic B<--X transition dipole moment surface, to form an initial wavepacket. The analysis of the time-dependent quantum dynamics of this wavepacket provides the total and partial photodissociation cross sections for the system. Both the total absorption cross section and the predicted product quantum state distributions compare well with experimental observations. A discussion is also given as to how the observed alternation in product diatom rotational state populations might be explained.
Collapse
Affiliation(s)
- Ezinvi Baloïtcha
- Center for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom.
| | | |
Collapse
|
31
|
Zúñiga J, Picón JAG, Bastida A, Requena A. On the use of optimal internal vibrational coordinates for symmetrical bent triatomic molecules. J Chem Phys 2005; 122:224319. [PMID: 15974680 DOI: 10.1063/1.1929738] [Citation(s) in RCA: 13] [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 use of generalized internal coordinates for the variational calculation of excited vibrational states of symmetrical bent triatomic molecules is considered with applications to the SO2, O3, NO2, and H2O molecules. These coordinates depend on two external parameters which can be properly optimized. We propose a simple analytical method to determine the optimal internal coordinates for this kind of molecules based on the minimization with respect to the external parameters of the zero-point energy, assuming only quadratic terms in the Hamiltonian and no quadratic coupling between the optimal coordinates. The optimal values of the parameters thus obtained are shown to agree quite well with those that minimize the sum of a number of unconverged energies of the lowest vibrational states, computed variationally using a small basis function set. The unconverged variational calculation uses a basis set consisting of the eigenfunctions of the uncoupled anharmonic internal coordinate Hamiltonian. Variational calculations of the excited vibrational states for the four molecules considered carried out with an increasing number of basis functions, also evidence the excellent convergence properties of the optimal internal coordinates versus those provided by other normal and local coordinate systems.
Collapse
Affiliation(s)
- José Zúñiga
- Departmento de Química Física, Universidad de Murcia, 30100 Murcia, Spain.
| | | | | | | |
Collapse
|
32
|
|
33
|
Baloïtcha E, Balint-Kurti GG. Theory of the photodissociation of ozone in the Hartley continuum; effect of vibrational excitation and O(1D) atom velocity distribution. Phys Chem Chem Phys 2005; 7:3829-33. [PMID: 16358032 DOI: 10.1039/b511640f] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The effect of vibrational excitation on the photodissociation cross section of ozone in the Hartley continuum is examined. The calculations make use of newly computed potential energy and transition dipole moment surfaces. The initial vibrational states of the ozone are computed using grid based techniques and the first few ab initio computed vibrational energy level spacings agree to within 10 cm(-1) with experimental values. The computed total absorption cross sections arising from different initial vibrational states of ozone are discussed in the light of the nature of the transition dipole moment surface. The computed cross section for excitation from the ground vibrational-rotational state is in good agreement with the experimentally measured cross section. Excitation of the asymmetric stretching vibration of ozone has a marked effect on both the form and magnitude of the photodissociation cross section. The velocity distributions of highly reactive O(1D) atoms arising from the photodissociation process in different wavelength ranges is also presented. The results show that the O(1D) atoms travel with a most probable translational velocity of 2.030 km s(-1) corresponding to a translational energy of 0.342 eV or 33.0 kJ mol(-1).
Collapse
Affiliation(s)
- Ezinvi Baloïtcha
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol, UK BS8 1TS.
| | | |
Collapse
|
34
|
Diehr M, Rosmus * P, Carter S, Knowles PJ. Theoretical rovibrational line intensities in the electronic ground state of ozone. Mol Phys 2004. [DOI: 10.1080/00268970410001722993] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
35
|
|
36
|
Schinke R, Fleurat-Lessard P. The transition-state region of the O(3P)+O2(3Σg−) potential energy surface. J Chem Phys 2004; 121:5789-93. [PMID: 15367004 DOI: 10.1063/1.1784776] [Citation(s) in RCA: 37] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
New electronic structure calculations for the transition-state region of the lowest ozone potential energy surface are reported. A two-dimensional potential energy surface in the asymptotic channel is calculated with the O(2) bond distance being fixed. The calculations are performed at the multireference average quadratic coupled cluster level of theory using full-valence complete active space self-consistent field wave functions and the augmented correlation consistent polarized V6Z atomic basis set. The general shape of the potential energy surface as predicted in earlier studies, that is, a narrow transition state below the O+O(2) asymptote, is confirmed by the present calculations. The transition state is 181 cm(-1) below the asymptote and 72 cm(-1) above the van der Waals-like minimum. The changes in the O+O(2)-->O(3) (*) capture cross section and rate constant when the new potential energy surface is employed are investigated by means of classical trajectory calculations.
Collapse
Affiliation(s)
- Reinhard Schinke
- Max-Planck-Institut für Strömungsforschung, Bunsenstrasse 10, D-37073 Göttingen, Federal Republic of Germany.
| | | |
Collapse
|
37
|
Lu SI. Diffusion quantum Monte Carlo for equilibrium structures and harmonic frequencies of ethane and ozone molecules. J Chem Phys 2004; 120:10423-5. [PMID: 15268070 DOI: 10.1063/1.1738114] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Application of the Ornstein-Uhlenbeck diffusion quantum Monte Carlo method in combination with a trial wave function constructed from the floating spherical Gaussian orbitals and spherical Gaussian geminals to studies on the equilibrium structures and harmonic frequencies of ethane and ozone is presented. These Monte Carlo computed results are compared with those of experiments as well as the coupled cluster methods with the correlation consistent basis sets for the two molecules. For ozone, we also compare the Monte Carlo results with the results from multireference calculations.
Collapse
Affiliation(s)
- Shih-I Lu
- Department of Applied Chemistry, Fooyin University, 151 Chinhsueh Road, Ta-Liao Hsiang, Kaohsiung Hsien, 831 Taiwan.
| |
Collapse
|
38
|
|
39
|
Venayagamoorthy M, Ford T. Ab initio molecular orbital studies of the vibrational spectra of some van der Waals complexes. Part 2. Complexes of ozone with carbon dioxide, carbonyl sulphide and carbon disulphide. J Mol Struct 2003. [DOI: 10.1016/s0022-2860(02)00670-1] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|
40
|
Rosmus P, Palmieri P, Schinke R. The asymptotic region of the potential energy surfaces relevant for the O(3P)+O2(X 3Σg−)⇌O3 reaction. J Chem Phys 2002. [DOI: 10.1063/1.1491396] [Citation(s) in RCA: 47] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
41
|
Charlo D, Clary DC. Quantum-mechanical calculations on termolecular association reactions XY+Z+M→XYZ+M: Application to ozone formation. J Chem Phys 2002. [DOI: 10.1063/1.1485069] [Citation(s) in RCA: 58] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
42
|
Siebert R, Fleurat-Lessard P, Schinke R, Bittererová M, Farantos SC. The vibrational energies of ozone up to the dissociation threshold: Dynamics calculations on an accurate potential energy surface. J Chem Phys 2002. [DOI: 10.1063/1.1473664] [Citation(s) in RCA: 159] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
43
|
|
44
|
Xie D, Guo H, Peterson KA. Ab initio characterization of low-lying triplet state potential-energy surfaces and vibrational frequencies in the Wulf band of ozone. J Chem Phys 2001. [DOI: 10.1063/1.1417502] [Citation(s) in RCA: 17] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|