1
|
Yang J, Li J, Li J, Li J. Gaussian Process Regression for State-to-State Integral Cross Sections: The Case of the O + O 2 Collision Dissociation Reactions. J Phys Chem A 2024; 128:4966-4975. [PMID: 38869143 DOI: 10.1021/acs.jpca.4c01445] [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
Research on hypersonic vehicles has become increasingly important worldwide in recent years. However, accurately simulating the dynamics of the nonequilibrium high-temperature reactions that are in the hypersonic flow around the vehicles presents a significant challenge as a large number of states and transitions are accessible even for the smallest atom-diatom reaction systems. It is quite difficult, sometimes even impossible, to exhaustively investigate all relevant combinations or determine high-dimensional analytical representations for the state-to-state reaction probabilities. In this study, we used Gaussian process regression (GPR) to fit a model based on only 807 QCT data for training. The confidence interval of the GPR prediction and the Kullback-Leibler (KL) divergence were used to help minimize the sampling amount of data for fitting the converged GPR model. The model aims to predict the state-to-state integral cross section (ICS) of the O + O2 → 3O dissociation reaction under random initial conditions (Et, v, j). In total, it took almost a month to obtain this converged GPR model, but it took only a few seconds to predict the ICS value for any initial condition. For 330 initial conditions not included in the training set, the mean-square error (MSE) between the QCT-calculated ICSs and the GPR-predicted ones is only 0.08 Å2 and the R2 is 0.9986, indicating that the GPR model can replace the direct expensive QCT calculation with high accuracy. Finally, we calculated the equilibrium dissociation rate coefficients based on the StS ICS values predicted by the GPR model, and the results were in good agreement with available experimental and theoretical results. Thus, this study provides an effective and accurate approach to the extensive direct state-to-state reaction dynamic calculations.
Collapse
Affiliation(s)
- Jiawei Yang
- School of Chemistry and Chemical Engineering & Chongqing Key Laboratory of Chemical Theory and Mechanism, Chongqing University, Chongqing 401331, China
| | - Jia Li
- School of Chemistry and Chemical Engineering & Chongqing Key Laboratory of Chemical Theory and Mechanism, Chongqing University, Chongqing 401331, China
| | - Junhong Li
- School of Chemistry and Chemical Engineering & Chongqing Key Laboratory of Chemical Theory and Mechanism, Chongqing University, Chongqing 401331, China
| | - Jun Li
- School of Chemistry and Chemical Engineering & Chongqing Key Laboratory of Chemical Theory and Mechanism, Chongqing University, Chongqing 401331, China
| |
Collapse
|
2
|
Iyengar SS, Ricard TC, Zhu X. Reformulation of All ONIOM-Type Molecular Fragmentation Approaches and Many-Body Theories Using Graph-Theory-Based Projection Operators: Applications to Dynamics, Molecular Potential Surfaces, Machine Learning, and Quantum Computing. J Phys Chem A 2024; 128:466-478. [PMID: 38180503 DOI: 10.1021/acs.jpca.3c05630] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2024]
Abstract
We present a graph-theory-based reformulation of all ONIOM-based molecular fragmentation methods. We discuss applications to (a) accurate post-Hartree-Fock AIMD that can be conducted at DFT cost for medium-sized systems, (b) hybrid DFT condensed-phase studies at the cost of pure density functionals, (c) reduced cost on-the-fly large basis gas-phase AIMD and condensed-phase studies, (d) post-Hartree-Fock-level potential surfaces at DFT cost to obtain quantum nuclear effects, and (e) novel transfer machine learning protocols derived from these measures. Additionally, in previous work, the unifying strategy discussed here has been used to construct new quantum computing algorithms. Thus, we conclude that this reformulation is robust and accurate.
Collapse
Affiliation(s)
- Srinivasan S Iyengar
- Department of Chemistry, Department of Physics, and the Indiana University Quantum Science and Engineering Center (IU-QSEC), Indiana University, 800 E. Kirkwood Avenue, Bloomington, Indiana 47405, United States
| | - Timothy C Ricard
- Department of Chemistry, Department of Physics, and the Indiana University Quantum Science and Engineering Center (IU-QSEC), Indiana University, 800 E. Kirkwood Avenue, Bloomington, Indiana 47405, United States
| | - Xiao Zhu
- Department of Chemistry, Department of Physics, and the Indiana University Quantum Science and Engineering Center (IU-QSEC), Indiana University, 800 E. Kirkwood Avenue, Bloomington, Indiana 47405, United States
| |
Collapse
|
3
|
Iyengar SS, Zhang JH, Saha D, Ricard TC. Graph-| Q⟩⟨ C|: A Quantum Algorithm with Reduced Quantum Circuit Depth for Electronic Structure. J Phys Chem A 2023; 127:9334-9345. [PMID: 37906738 DOI: 10.1021/acs.jpca.3c04261] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/02/2023]
Abstract
The accurate determination of chemical properties is known to have a critical impact on multiple fundamental chemical problems but is deeply hindered by the steep algebraic scaling of electron correlation calculations and the exponential scaling of quantum nuclear dynamics. With the advent of new quantum computing hardware and associated developments in creating new paradigms for quantum software, this avenue has been recognized as perhaps one way to address exponentially complex challenges in quantum chemistry and molecular dynamics. In this paper, we discuss a new approach to drastically reduce the quantum circuit depth (by several orders of magnitude) and help improve the accuracy in the quantum computation of electron correlation energies for large molecular systems. The method is derived from a graph-theoretic approach to molecular fragmentation and enables us to create a family of projection operators that decompose quantum circuits into separate unitary processes. Some of these processes can be treated on quantum hardware and others on classical hardware in a completely asynchronous and parallel fashion. Numerical benchmarks are provided through the computation of unitary coupled-cluster singles and doubles (UCCSD) energies for medium-sized protonated and neutral water clusters using the new quantum algorithms presented here.
Collapse
Affiliation(s)
- Srinivasan S Iyengar
- Department of Chemistry, Department of Physics, and the Indiana University Quantum Science and Engineering Center (IU-QSEC), Indiana University, 800 E. Kirkwood Avenue, Bloomington, Indiana 47405, United States
| | - Juncheng Harry Zhang
- Department of Chemistry, Department of Physics, and the Indiana University Quantum Science and Engineering Center (IU-QSEC), Indiana University, 800 E. Kirkwood Avenue, Bloomington, Indiana 47405, United States
| | - Debadrita Saha
- Department of Chemistry, Department of Physics, and the Indiana University Quantum Science and Engineering Center (IU-QSEC), Indiana University, 800 E. Kirkwood Avenue, Bloomington, Indiana 47405, United States
| | - Timothy C Ricard
- Department of Chemistry, Department of Physics, and the Indiana University Quantum Science and Engineering Center (IU-QSEC), Indiana University, 800 E. Kirkwood Avenue, Bloomington, Indiana 47405, United States
| |
Collapse
|
4
|
Kumar A, DeGregorio N, Ricard T, Iyengar SS. Graph-Theoretic Molecular Fragmentation for Potential Surfaces Leads Naturally to a Tensor Network Form and Allows Accurate and Efficient Quantum Nuclear Dynamics. J Chem Theory Comput 2022; 18:7243-7259. [PMID: 36332133 DOI: 10.1021/acs.jctc.2c00484] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Abstract
Molecular fragmentation methods have revolutionized quantum chemistry. Here, we use a graph-theoretically generated molecular fragmentation method, to obtain accurate and efficient representations for multidimensional potential energy surfaces and the quantum time-evolution operator, which plays a critical role in quantum chemical dynamics. In doing so, we find that the graph-theoretic fragmentation approach naturally reduces the potential portion of the time-evolution operator into a tensor network that contains a stream of coupled lower-dimensional propagation steps to potentially achieve quantum dynamics with reduced complexity. Furthermore, the fragmentation approach used here has previously been shown to allow accurate and efficient computation of post-Hartree-Fock electronic potential energy surfaces, which in many cases has been shown to be at density functional theory cost. Thus, by combining the advantages of molecular fragmentation with the tensor network formalism, the approach yields an on-the-fly quantum dynamics scheme where both the electronic potential calculation and nuclear propagation portion are enormously simplified through a single stroke. The method is demonstrated by computing approximations to the propagator and to potential surfaces for a set of coupled nuclear dimensions within a protonated water wire problem exhibiting the Grotthuss mechanism of proton transport. In all cases, our approach has been shown to reduce the complexity of representing the quantum propagator, and by extension action of the propagator on an initial wavepacket, by several orders, with minimal loss in accuracy.
Collapse
Affiliation(s)
- Anup Kumar
- Department of Chemistry, and the Indiana University Quantum Science and Engineering Center (IU-QSEC), Indiana University, Bloomington, Indiana 47405, United States
| | - Nicole DeGregorio
- Department of Chemistry, and the Indiana University Quantum Science and Engineering Center (IU-QSEC), Indiana University, Bloomington, Indiana 47405, United States
| | - Timothy Ricard
- Department of Chemistry, and the Indiana University Quantum Science and Engineering Center (IU-QSEC), Indiana University, Bloomington, Indiana 47405, United States
| | - Srinivasan S Iyengar
- Department of Chemistry, and the Indiana University Quantum Science and Engineering Center (IU-QSEC), Indiana University, Bloomington, Indiana 47405, United States
| |
Collapse
|
5
|
Zhu X, Iyengar SS. Graph Theoretic Molecular Fragmentation for Multidimensional Potential Energy Surfaces Yield an Adaptive and General Transfer Machine Learning Protocol. J Chem Theory Comput 2022; 18:5125-5144. [PMID: 35994592 DOI: 10.1021/acs.jctc.1c01241] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
Over a series of publications we have introduced a graph-theoretic description for molecular fragmentation. Here, a system is divided into a set of nodes, or vertices, that are then connected through edges, faces, and higher-order simplexes to represent a collection of spatially overlapping and locally interacting subsystems. Each such subsystem is treated at two levels of electronic structure theory, and the result is used to construct many-body expansions that are then embedded within an ONIOM-scheme. These expansions converge rapidly with many-body order (or graphical rank) of subsystems and have been previously used for ab initio molecular dynamics (AIMD) calculations and for computing multidimensional potential energy surfaces. Specifically, in all these cases we have shown that CCSD and MP2 level AIMD trajectories and potential surfaces may be obtained at density functional theory cost. The approach has been demonstrated for gas-phase studies, for condensed phase electronic structure, and also for basis set extrapolation-based AIMD. Recently, this approach has also been used to derive new quantum-computing algorithms that enormously reduce the quantum circuit depth in a circuit-based computation of correlated electronic structure. In this publication, we introduce (a) a family of neural networks that act in parallel to represent, efficiently, the post-Hartree-Fock electronic structure energy contributions for all simplexes (fragments), and (b) a new k-means-based tessellation strategy to glean training data for high-dimensional molecular spaces and minimize the extent of training needed to construct this family of neural networks. The approach is particularly useful when coupled cluster accuracy is desired and when fragment sizes grow in order to capture nonlocal interactions accurately. The unique multidimensional k-means tessellation/clustering algorithm used to determine our training data for all fragments is shown to be extremely efficient and reduces the needed training to only 10% of data for all fragments to obtain accurate neural networks for each fragment. These fully connected dense neural networks are then used to extrapolate the potential energy surface for all molecular fragments, and these are then combined as per our graph-theoretic procedure to transfer the learning process to a full system energy for the entire AIMD trajectory at less than one-tenth the cost as compared to a regular fragmentation-based AIMD calculation.
Collapse
Affiliation(s)
- Xiao Zhu
- Department of Chemistry and Department of Physics, Indiana University, 800 E. Kirkwood Avenue, Bloomington 47405, Indiana, United States
| | - Srinivasan S Iyengar
- Department of Chemistry and Department of Physics, Indiana University, 800 E. Kirkwood Avenue, Bloomington 47405, Indiana, United States
| |
Collapse
|
6
|
Zhang T, Lu Y, Cheng X. State-Specific Dynamic Study of the Exchange and Dissociation Reaction for O(3P) and O2($${}^{3}\Sigma _{g}^{ - }$$) Collision by Quasi-Classical Trajectory. RUSSIAN JOURNAL OF PHYSICAL CHEMISTRY A 2022. [DOI: 10.1134/s0036024422040331] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
7
|
Kumar A, DeGregorio N, Iyengar SS. Graph-Theory-Based Molecular Fragmentation for Efficient and Accurate Potential Surface Calculations in Multiple Dimensions. J Chem Theory Comput 2021; 17:6671-6690. [PMID: 34623129 DOI: 10.1021/acs.jctc.1c00065] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
We present a multitopology molecular fragmentation approach, based on graph theory, to calculate multidimensional potential energy surfaces in agreement with post-Hartree-Fock levels of theory but at the density functional theory cost. A molecular assembly is coarse-grained into a set of graph-theoretic nodes that are then connected with edges to represent a collection of locally interacting subsystems up to an arbitrary order. Each of the subsystems is treated at two levels of electronic structure theory, the result being used to construct many-body expansions that are embedded within an ONIOM scheme. These expansions converge rapidly with the many-body order (or graphical rank) of subsystems and capture many-body interactions accurately and efficiently. However, multiple graphs, and hence multiple fragmentation topologies, may be defined in molecular configuration space that may arise during conformational sampling or from reactive, bond breaking and bond formation, events. Obtaining the resultant potential surfaces is an exponential scaling proposition, given the number of electronic structure computations needed. We utilize a family of graph-theoretic representations within a variational scheme to obtain multidimensional potential surfaces at a reduced cost. The fast convergence of the graph-theoretic expansion with increasing order of many-body interactions alleviates the exponential scaling cost for computing potential surfaces, with the need to only use molecular fragments that contain a fewer number of quantum nuclear degrees of freedom compared to the full system. This is because the dimensionality of the conformational space sampled by the fragment subsystems is much smaller than the full molecular configurational space. Additionally, we also introduce a multidimensional clustering algorithm, based on physically defined criteria, to reduce the number of energy calculations by orders of magnitude. The molecular systems benchmarked include coupled proton motion in protonated water wires. The potential energy surfaces and multidimensional nuclear eigenstates obtained are shown to be in very good agreement with those from explicit post-Hartree-Fock calculations that become prohibitive as the number of quantum nuclear dimensions grows. The developments here provide a rigorous and efficient alternative to this important chemical physics problem.
Collapse
Affiliation(s)
- Anup Kumar
- Department of Chemistry and Department of Physics, Indiana University, 800 E. Kirkwood Avenue, Bloomington, Indiana 47405, United States
| | - Nicole DeGregorio
- Department of Chemistry and Department of Physics, Indiana University, 800 E. Kirkwood Avenue, Bloomington, Indiana 47405, United States
| | - Srinivasan S Iyengar
- Department of Chemistry and Department of Physics, Indiana University, 800 E. Kirkwood Avenue, Bloomington, Indiana 47405, United States
| |
Collapse
|
8
|
Zhang JH, Ricard TC, Haycraft C, Iyengar SS. Weighted-Graph-Theoretic Methods for Many-Body Corrections within ONIOM: Smooth AIMD and the Role of High-Order Many-Body Terms. J Chem Theory Comput 2021; 17:2672-2690. [PMID: 33891416 DOI: 10.1021/acs.jctc.0c01287] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
We present a weighted-graph-theoretic approach to adaptively compute contributions from many-body approximations for smooth and accurate post-Hartree-Fock (pHF) ab initio molecular dynamics (AIMD) of highly fluxional chemical systems. This approach is ONIOM-like, where the full system is treated at a computationally feasible quality of treatment (density functional theory (DFT) for the size of systems considered in this publication), which is then improved through a perturbative correction that captures local many-body interactions up to a certain order within a higher level of theory (post-Hartree-Fock in this publication) described through graph-theoretic techniques. Due to the fluxional and dynamical nature of the systems studied here, these graphical representations evolve during dynamics. As a result, energetic "hops" appear as the graphical representation deforms with the evolution of the chemical and physical properties of the system. In this paper, we introduce dynamically weighted, linear combinations of graphs, where the transition between graphical representations is smoothly achieved by considering a range of neighboring graphical representations at a given instant during dynamics. We compare these trajectories with those obtained from a set of trajectories where the range of local many-body interactions considered is increased, sometimes to the maximum available limit, which yields conservative trajectories as the order of interactions is increased. The weighted-graph approach presents improved dynamics trajectories while only using lower-order many-body interaction terms. The methods are compared by computing dynamical properties through time-correlation functions and structural distribution functions. In all cases, the weighted-graph approach provides accurate results at a lower cost.
Collapse
Affiliation(s)
- Juncheng Harry Zhang
- Department of Chemistry and Department of Physics, Indiana University, 800 East Kirkwood Avenue, Bloomington, Indiana 47405, United States
| | - Timothy C Ricard
- Department of Chemistry and Department of Physics, Indiana University, 800 East Kirkwood Avenue, Bloomington, Indiana 47405, United States
| | - Cody Haycraft
- Department of Chemistry and Department of Physics, Indiana University, 800 East Kirkwood Avenue, Bloomington, Indiana 47405, United States
| | - Srinivasan S Iyengar
- Department of Chemistry and Department of Physics, Indiana University, 800 East Kirkwood Avenue, Bloomington, Indiana 47405, United States
| |
Collapse
|
9
|
Pan TJ, Stephani KA. Rovibrationally state-specific collision model for the O 2(Σg-3) + O(P3) system in DSMC. J Chem Phys 2021; 154:104306. [PMID: 33722033 DOI: 10.1063/5.0027411] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A rovibrationally state-specific collision model for the O2(Σg-3)+O(P3) system is presented for direct simulation Monte Carlo, including rotation-vibration-translation energy transfer, exchange, dissociation, and recombination processes. The two-step binary collision approach is employed to model recombination reactions. Two available cross section databases by Andrienko/Boyd and Esposito/Capitelli are employed for the rovibrationally resolved model (rv-STS) and vibrationally resolved model (v-STS), respectively. The difference between rv-STS and v-STS comes from two contributions: the multisurface factor of dissociation (fMS) and the rotational averaging process. The dissociation cross section with the constant fMS is typically larger than with the variable fMS, especially for the low vibrational energy states. On the other hand, the cross sections resulting from the rotationally averaged database are found to underpredict the dissociation rate coefficient at low temperatures. In the rovibrational heating case, the rv-STS predicts faster relaxation than the v-STS, which also shows a lower quasi-steady-state temperature than v-STS. In the rovibrational cooling case, the rv-STS shows a faster relaxation than v-STS, which also presents a thermal non-equilibrium between rovibrational and translational mode during the cooling process.
Collapse
Affiliation(s)
- Tzu-Jung Pan
- Center for Hypersonics and Entry Systems Studies, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| | - Kelly A Stephani
- Center for Hypersonics and Entry Systems Studies, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| |
Collapse
|
10
|
Geistfeld E, Schwartzentruber TE. QCT calculations of O 2 + O collisions: Comparison to molecular beam experiments. J Chem Phys 2020; 153:184302. [DOI: 10.1063/5.0024870] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- E. Geistfeld
- Department of Aerospace Engineering and Mechanics, University of Minnesota, Minneapolis, Minnesota 55455, USA
| | - T. E. Schwartzentruber
- Department of Aerospace Engineering and Mechanics, University of Minnesota, Minneapolis, Minnesota 55455, USA
| |
Collapse
|
11
|
Venturi S, Sharma MP, Lopez B, Panesi M. Data-Inspired and Physics-Driven Model Reduction for Dissociation: Application to the O 2 + O System. J Phys Chem A 2020; 124:8359-8372. [PMID: 32886505 DOI: 10.1021/acs.jpca.0c04516] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
This work presents an in-depth discussion on the nonequilibrium dissociation of O2 molecules colliding with O atoms, combining quasi-classical trajectory calculations, master equation, and dimensionality reduction. A rovibrationally resolved database for all of the elementary collisional processes is constructed by including all nine adiabatic electronic states of O3 in the QCT calculations. A detailed analysis of the ab initio data set reveals that for a rovibrational level, the probability of dissociating is mostly dictated by its deficit in internal energy compared to the centrifugal barrier. Because of the assumption of rotational equilibrium, the conventional vibrational-specific calculations fail to characterize such a dependence. Based on this observation, a new physics-based grouping strategy for application to coarse-grained models is proposed. By relying on a hybrid technique made of rovibrationally resolved excitation coupled to coarse-grained dissociation, the new approach is compared to the vibrational-specific model and the direct solution of the rovibrational state-to-state master equation. Simulations are performed in a zero-dimensional isothermal and isochoric chemical reactor for a wide range of temperatures (1500-20,000 K). The study shows that the main contribution to the model inadequacy of vibrational-specific approaches originates from the incapability of characterizing dissociation, rather than the energy transfers. Even when constructed with only twenty groups, the new reduced-order model outperforms the vibrational-specific one in predicting all of the QoIs related to dissociation kinetics. At the highest temperature, the accuracy in the mole fraction is improved by 2000%.
Collapse
Affiliation(s)
- S Venturi
- University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
| | - M P Sharma
- University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
| | - B Lopez
- University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
| | - M Panesi
- University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
| |
Collapse
|
12
|
Ricard TC, Iyengar SS. Efficient and Accurate Approach To Estimate Hybrid Functional and Large Basis-Set Contributions to Condensed-Phase Systems and Molecule–Surface Interactions. J Chem Theory Comput 2020; 16:4790-4812. [DOI: 10.1021/acs.jctc.9b01089] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Affiliation(s)
- Timothy C. Ricard
- Department of Chemistry and Department of Physics, Indiana University, 800 E. Kirkwood Avenue, Bloomington, Indiana 47405, United States
| | - Srinivasan S. Iyengar
- Department of Chemistry and Department of Physics, Indiana University, 800 E. Kirkwood Avenue, Bloomington, Indiana 47405, United States
| |
Collapse
|
13
|
Venturi S, Jaffe RL, Panesi M. Bayesian Machine Learning Approach to the Quantification of Uncertainties on Ab Initio Potential Energy Surfaces. J Phys Chem A 2020; 124:5129-5146. [DOI: 10.1021/acs.jpca.0c02395] [Citation(s) in RCA: 28] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Affiliation(s)
- S. Venturi
- University of Illinois at Urbana−Champaign, Urbana, Illinois 61801, United States
| | - R. L. Jaffe
- NASA Ames Research Center, Moffett Field, California 94035-1000, United States
| | - M. Panesi
- University of Illinois at Urbana−Champaign, Urbana, Illinois 61801, United States
| |
Collapse
|
14
|
Andrienko DA. The importance of O 3 excited potential energy surfaces in O 2-O high-temperature kinetics. J Chem Phys 2020; 152:044305. [PMID: 32007085 DOI: 10.1063/1.5142191] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022] Open
Abstract
The mechanism of vibrational relaxation and dissociation in the O2-O system at elevated temperatures is investigated by means of molecular dynamics. The most recent O3 potential energy surfaces (PESs), obtained from the first principles quantum mechanical calculations [Varga et al., J. Chem. Phys. 147, 154312 (2017)], are used to derive a complete set of state-specific rate coefficients of vibrational energy transfer and dissociation. Unlike most of the previous efforts that utilize only the lowest and supposedly most reactive 11A' O3 PES [A. Varandas and A. Pais, Mol. Phys. 65, 843 (1988)], this paper demonstrates the necessity to account for a complete ensemble of all excited O3 PESs that correlate with O2(X) and O(3P) when high-temperature kinetics is of interest. At the same time, it is found that the Varandas 11A' O3 PES adequately describes vibrational energy transfer and dissociating dynamics when compared to the most recent 11A' O3 PES by Varga et al. [J. Chem Phys. 147, 154312 (2017)]. The differences between this new dataset and previous rate coefficients are quantified by the master equation model.
Collapse
Affiliation(s)
- Daniil A Andrienko
- Department of Aerospace Engineering, Texas A&M University, H. R. Bright Building, 701 Ross Street, College Station, Texas 77840, USA
| |
Collapse
|
15
|
Pan TJ, Wilson TJ, Stephani KA. Vibrational state-specific model for dissociation and recombination of the O2( 3Σg−)+O( 3P) system in DSMC. J Chem Phys 2019; 150:074305. [DOI: 10.1063/1.5035283] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022] Open
Affiliation(s)
- Tzu-Jung Pan
- Department of Mechanical Science and Engineering, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801-2946, USA
| | - Taiyo J. Wilson
- Department of Mechanical Science and Engineering, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801-2946, USA
| | - Kelly A. Stephani
- Department of Mechanical Science and Engineering, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801-2946, USA
| |
Collapse
|
16
|
Quantum scattering theory for collisional energy transfer. ACTA ACUST UNITED AC 2019. [DOI: 10.1016/b978-0-444-64207-3.00002-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register]
|
17
|
Oblapenko GP. Calculation of Vibrational Relaxation Times Using a Kinetic Theory Approach. J Phys Chem A 2018; 122:9615-9625. [PMID: 30482021 DOI: 10.1021/acs.jpca.8b09897] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
In the present work, a method for computation of vibrational relaxation times based on a kinetic theory definition is utilized to calculate vibrational relaxation times of molecules present in air (N2, O2, and NO) in collisions with air species particles. An overview of available experiment VT relaxation time measurements, as well as quasi-classical trajectory (QCT) calculation results, and various empirical models, is given. Different inelastic cross-section models are used for the computation of the relaxation times, and their parameters are adjusted to fit the available experimental data and QCT results. It is shown that the proposed method of calculation can give a quantitative and qualitative agreement with the available data in a wide range of temperatures; the obtained interaction parameters may be used not only for vibrational relaxation time calculation within a multitemperature framework but also for development of state-specific models for use in CFD and DSMC codes.
Collapse
Affiliation(s)
- G P Oblapenko
- Saint Petersburg State University , 7/9 Universitetskaya nab. , 199034 St. Petersburg , Russia
| |
Collapse
|
18
|
Ricard TC, Iyengar SS. Efficiently Capturing Weak Interactions in ab Initio Molecular Dynamics with on-the-Fly Basis Set Extrapolation. J Chem Theory Comput 2018; 14:5535-5552. [DOI: 10.1021/acs.jctc.8b00803] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Affiliation(s)
- Timothy C. Ricard
- Department of Chemistry and Department of Physics, Indiana University, 800 E. Kirkwood Avenue, Bloomington, Indiana 47405, United States
| | - Srinivasan S. Iyengar
- Department of Chemistry and Department of Physics, Indiana University, 800 E. Kirkwood Avenue, Bloomington, Indiana 47405, United States
| |
Collapse
|
19
|
Pogosbekyan MY, Sergievskaya AL. Simulation of the Oxygen Dissociation Reaction under Thermally Nonequilibrium Conditions: Models, Trajectory Calculations, and the Experiment. RUSSIAN JOURNAL OF PHYSICAL CHEMISTRY B 2018. [DOI: 10.1134/s1990793118020239] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
20
|
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
|
21
|
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
|
22
|
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
|
23
|
Andrienko DA, Boyd ID. Thermal relaxation of molecular oxygen in collisions with nitrogen atoms. J Chem Phys 2016; 145:014309. [PMID: 27394110 DOI: 10.1063/1.4955199] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Daniil A. Andrienko
- Department of Aerospace Engineering, University of Michigan, 1320 Beal Ave., Ann Arbor, Michigan 48108, USA
| | - Iain D. Boyd
- Department of Aerospace Engineering, University of Michigan, 1320 Beal Ave., Ann Arbor, Michigan 48108, USA
| |
Collapse
|
24
|
Kulakhmetov M, Gallis M, Alexeenko A. Ab initio-informed maximum entropy modeling of rovibrational relaxation and state-specific dissociation with application to the O2 + O system. J Chem Phys 2016; 144:174302. [PMID: 27155635 DOI: 10.1063/1.4947590] [Citation(s) in RCA: 64] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022] Open
Affiliation(s)
- Marat Kulakhmetov
- Aeronautics and Astronautics, Purdue University, West Lafayette, Indiana 47907, USA
| | - Michael Gallis
- Engineering Sciences Center, Sandia National Laboratories, Albuquerque, New Mexico 87185, USA
| | - Alina Alexeenko
- Aeronautics and Astronautics, Purdue University, West Lafayette, Indiana 47907, USA
| |
Collapse
|
25
|
Andrienko DA, Boyd ID. Rovibrational energy transfer and dissociation in O2-O collisions. J Chem Phys 2016; 144:104301. [PMID: 26979687 DOI: 10.1063/1.4943114] [Citation(s) in RCA: 90] [Impact Index Per Article: 11.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A set of state-specific transition rates for each rovibrational level is generated for the O2(X(3)Σ(g)(-))-O(3)P system using the quasi-classical trajectory method at temperatures observed in hypersonic flows. A system of master equations describes the relaxation of the rovibrational ensemble to thermal equilibrium under ideal heat bath conditions at a constant translational temperature. Vibrational and rotational relaxation times, obtained from the average internal energies, exhibit a pattern inherent in a chemically reactive collisional pair. An intrinsic feature of the O3 molecular system with a large attractive potential is a weak temperature dependence of the rovibrational transition rates. For this reason, the quasi-steady vibrational and rotational temperatures experience a maximum at increasing translational temperature. The energy rate coefficients, that characterize the average loss of internal energy due to dissociation, quickly diminish at high temperatures, compared to other molecular systems.
Collapse
Affiliation(s)
- Daniil A Andrienko
- Department of Aerospace Engineering, University of Michigan, 1320 Beal Ave., Ann Arbor, Michigan 48108, USA
| | - Iain D Boyd
- Department of Aerospace Engineering, University of Michigan, 1320 Beal Ave., Ann Arbor, Michigan 48108, USA
| |
Collapse
|
26
|
Reactivity and Relaxation of Vibrationally/Rotationally Excited Molecules with Open Shell Atoms. ACTA ACUST UNITED AC 2016. [DOI: 10.1007/978-1-4419-8185-1_2] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/22/2023]
|
27
|
Andrienko D, Boyd ID. Investigation of oxygen vibrational relaxation by quasi-classical trajectory method. Chem Phys 2015. [DOI: 10.1016/j.chemphys.2015.07.023] [Citation(s) in RCA: 39] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
28
|
Teixeira OBM, Mota VC, Garcia de la Vega JM, Varandas AJC. Single-Sheeted Double Many-Body Expansion Potential Energy Surface for Ground-State ClO2. J Phys Chem A 2014; 118:4851-62. [DOI: 10.1021/jp503744x] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- O. B. M. Teixeira
- Departamento
de Química Física Aplicada, Universidad Autónoma de Madrid, 28049 Madrid, Spain
| | - V. C. Mota
- Departamento
de Física, Universidade Federal do Espirito Santo, 29075-910 Vitória, Brazil
| | - J. M. Garcia de la Vega
- Departamento
de Química Física Aplicada, Universidad Autónoma de Madrid, 28049 Madrid, Spain
| | - A. J. C. Varandas
- Departamento
de Física, Universidade Federal do Espirito Santo, 29075-910 Vitória, Brazil
- Centro de Química,
and Departamento de Química, Universidade de Coimbra, 3004-535 Coimbra, Portugal
| |
Collapse
|
29
|
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
|
30
|
Li J, Varandas AJC. Accurate ab-Initio-Based Single-Sheeted DMBE Potential-Energy Surface for Ground-State N2O. J Phys Chem A 2012; 116:4646-56. [DOI: 10.1021/jp302173h] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Jing Li
- Departamento de Química, Universidade de Coimbra, 3004-535 Coimbra, Portugal
| | | |
Collapse
|
31
|
|
32
|
Lino da Silva M, Loureiro J, Guerra V. A multiquantum dataset for vibrational excitation and dissociation in high-temperature O2–O2 collisions. Chem Phys Lett 2012. [DOI: 10.1016/j.cplett.2012.01.074] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/14/2022]
|
33
|
Mota VC, Caridade PJSB, Varandas AJC. Ab Initio-Based Global Double Many-Body Expansion Potential Energy Surface for the First 2A″ Electronic State of NO2. J Phys Chem A 2012; 116:3023-34. [DOI: 10.1021/jp300031q] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- V. C. Mota
- Departamento de Química, Universidade de Coimbra, 3004-535 Coimbra, Portugal
| | | | - A. J. C. Varandas
- Departamento de Química, Universidade de Coimbra, 3004-535 Coimbra, Portugal
| |
Collapse
|
34
|
Rampino S, Faginas Lago N, Laganà A, Huarte-Larrañaga F. An extension of the grid empowered molecular simulator to quantum reactive scattering. J Comput Chem 2011; 33:708-14. [DOI: 10.1002/jcc.22878] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/26/2011] [Accepted: 10/28/2011] [Indexed: 11/10/2022]
|
35
|
Song YZ, Varandas AJC. Accurate Double Many-Body Expansion Potential Energy Surface for Ground-State HS2 Based on ab Initio Data Extrapolated to the Complete Basis Set Limit. J Phys Chem A 2011; 115:5274-83. [DOI: 10.1021/jp201980m] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Y. Z. Song
- Departamento de Química, Universidade de Coimbra, 3004-535 Coimbra, Portugal
| | - A. J. C. Varandas
- Departamento de Química, Universidade de Coimbra, 3004-535 Coimbra, Portugal
| |
Collapse
|
36
|
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
|
37
|
BOGGIO-PASQUA M, VORONIN AI, HALVICK PH, RAYEZ JC, VARANDAS AJC. Coupled ab initio potential energy surfaces for the two lowest 2A′ electronic states of the C2H molecule. Mol Phys 2009. [DOI: 10.1080/00268970009483396] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
Affiliation(s)
- M. BOGGIO-PASQUA
- a Laboratoire de Physico-Chimie Moléculaire—UMR 5803/CNRS , Université Bordeaux 1 , 33405 , Talence Cedex, France
| | - A. I. VORONIN
- a Laboratoire de Physico-Chimie Moléculaire—UMR 5803/CNRS , Université Bordeaux 1 , 33405 , Talence Cedex, France
| | - PH. HALVICK
- a Laboratoire de Physico-Chimie Moléculaire—UMR 5803/CNRS , Université Bordeaux 1 , 33405 , Talence Cedex, France
| | - J.-C. RAYEZ
- a Laboratoire de Physico-Chimie Moléculaire—UMR 5803/CNRS , Université Bordeaux 1 , 33405 , Talence Cedex, France
| | - A. J. C. VARANDAS
- b Departamento de Química , Universidade de Coimbra , 3049 , Coimbra Codex , Portugal
| |
Collapse
|
38
|
Rampino S, Skouteris D, Laganà A. The O + O2 reaction: quantum detailed probabilities and thermal rate coefficients. Theor Chem Acc 2009. [DOI: 10.1007/s00214-009-0524-1] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
|
39
|
Esposito F, Armenise I, Capitta G, Capitelli M. O–O2 state-to-state vibrational relaxation and dissociation rates based on quasiclassical calculations. Chem Phys 2008. [DOI: 10.1016/j.chemphys.2008.04.004] [Citation(s) in RCA: 101] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
|
40
|
Esposito F, Capitelli M. The relaxation of vibrationally excited O2 molecules by atomic oxygen. Chem Phys Lett 2007. [DOI: 10.1016/j.cplett.2007.06.099] [Citation(s) in RCA: 46] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
41
|
Charlo D, Clary DC. Quantum-mechanical calculations on pressure and temperature dependence of three-body recombination reactions: application to ozone formation rates. J Chem Phys 2006; 120:2700-7. [PMID: 15268414 DOI: 10.1063/1.1635361] [Citation(s) in RCA: 47] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
Abstract
A quantum-mechanical model is designed for the calculation of termolecular association reaction rate coefficients in the low-pressure fall-off regime. The dynamics is set up within the energy transfer mechanism and the kinetic scheme is the steady-state approximation. We applied this model to the formation of ozone O + O2 + M --> O3 + M for M = Ar, making use of semiquantitative potential energy surfaces. The stabilization process is treated by means of the vibrational close-coupling infinite order sudden scattering theory. Major approximations include the neglect of the O3 vibrational bending mode and rovibrational couplings. We calculated individual isotope-specific rate constants and rate constant ratios over the temperature range 10-1000 K and the pressure fall-off region 10(-7)-10(2) bar. The present results show a qualitative and semiquantitative agreement with available experiments, particularly in the temperature region of atmospheric interest.
Collapse
Affiliation(s)
- David Charlo
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
| | | |
Collapse
|
42
|
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
|
43
|
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
|
44
|
Varandas AJC, Rodrigues SPJ. New Double Many-Body Expansion Potential Energy Surface for Ground-State HCN from a Multiproperty Fit to Accurate ab Initio Energies and Rovibrational Calculations. J Phys Chem A 2005; 110:485-93. [PMID: 16405320 DOI: 10.1021/jp051434p] [Citation(s) in RCA: 41] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
An accurate single-sheeted double many-body expansion potential energy surface has been obtained for the ground electronic state of the hydrogen cyanide molecule via a multiproperty fit to ab initio energies and rovibrational data. This includes 106 rovibrational levels and 2313 discrete points, which are fit with a rmsd of 4 cm(-1) and 2.42 kcal mol(-1), respectively, and seven zero first-derivatives that are reproduced at three stationary points. Since the potential also describes accurately the appropriate asymptotic limits at the various dissociation channels, it is commended both for the simulation of rovibrational spectra and reaction dynamics.
Collapse
Affiliation(s)
- A J C Varandas
- Departamento de Química, Universidade de Coimbra 3004-535 Coimbra, Portugal
| | | |
Collapse
|
45
|
Tashiro M, Schinke R. The effect of spin–orbit coupling in complex forming O(3P) +O2 collisions. J Chem Phys 2003. [DOI: 10.1063/1.1616919] [Citation(s) in RCA: 41] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022] Open
|
46
|
Fleurat-Lessard P, Grebenshchikov SY, Schinke R, Janssen C, Krankowsky D. Isotope dependence of the O+O2 exchange reaction: Experiment and theory. J Chem Phys 2003. [DOI: 10.1063/1.1595091] [Citation(s) in RCA: 63] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022] Open
|
47
|
Yeh KL, Xie D, Zhang DH, Lee SY, Schinke R. Time-Dependent Wave Packet Study of the O + O2 (v = 0, j = 0) Exchange Reaction. J Phys Chem A 2003. [DOI: 10.1021/jp034471u] [Citation(s) in RCA: 27] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Affiliation(s)
- Ka-Lo Yeh
- Department of Computational Science The National University of Singapore, Singapore 119260
| | - Daiqian Xie
- Department of Computational Science The National University of Singapore, Singapore 119260
| | - Dong H. Zhang
- Department of Computational Science The National University of Singapore, Singapore 119260
| | - Soo-Y Lee
- Department of Chemistry The National University of Singapore, Singapore 119260
| | - Reinhard Schinke
- Max-Planck-Institut für Strömungsforschung, D-37073 Göttingen, Germany
| |
Collapse
|
48
|
|
49
|
Fleurat-Lessard P, Grebenshchikov SY, Siebert R, Schinke R, Halberstadt N. Theoretical investigation of the temperature dependence of the O+O2 exchange reaction. J Chem Phys 2003. [DOI: 10.1063/1.1525255] [Citation(s) in RCA: 66] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
50
|
Baker TA, Gellene GI. Classical and quasi-classical trajectory calculations of isotope exchange and ozone formation proceeding through O+O2 collision complexes. J Chem Phys 2002. [DOI: 10.1063/1.1508373] [Citation(s) in RCA: 28] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
|