1
|
Peng Y, Zhang H. Mechanism and Kinetics of Methane Combustion. Part II: Potential Energy Surface for Hydrogen-Abstraction Reaction of CH 4 + O( 3P). J Phys Chem A 2022; 126:1946-1959. [PMID: 35298157 DOI: 10.1021/acs.jpca.1c10860] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Methane combustion plays an important role in various fields such as combustion chemistry and atmospheric chemistry of the stratosphere. Highly accurate study of its initial reaction remains a key challenge. Here, through extensive studies with a state-of-the-art ab initio and neural network method, we present a potential energy surface of the O(3P) + CH4 → OH + CH3 reaction on the ground state 13A and the first excited state 23A. In this work, the energies of 10 167 points covering all important regions are obtained with state-averaged complete active space self-consistent field calculations and then fitted using the Levenberg-Marquardt algorithm with a root-mean-square error of 0.391 and 0.442 kcal/mol for the 13A and 23A states, respectively. This study explores the characteristics of the radical van der Waals (VdW) complex and reveals a detailed mechanism of the methane combustion initial reaction. Within the scope of this mechanism, this surface gives a fairly accurate description of the regions around the saddle point, conical intersection, and vdW wells in the entrance for efficient computational simulations. As a theoretical study on a prototypical polyatomic reaction, it is hopeful that this work will modify our understanding of the primary process in hydrocarbon combustion.
Collapse
Affiliation(s)
- Ya Peng
- Department of Engineering Physics, Tsinghua University, Beijing 100084, P.R. China
| | - Hui Zhang
- Department of Engineering Physics, Tsinghua University, Beijing 100084, P.R. China
| |
Collapse
|
2
|
Roncero O, Zanchet A, Aguado A. Low temperature reaction dynamics for CH 3OH + OH collisions on a new full dimensional potential energy surface. Phys Chem Chem Phys 2018; 20:25951-25958. [PMID: 30294740 PMCID: PMC6290987 DOI: 10.1039/c8cp04970j] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
Abstract
Is the rise of the rate constant measured in laval expansion experiments of OH with organic molecules at low temperatures due to the reaction between the reactants or due to the formation of complexes with the buffer gas? This question has importance for understanding the evolution of prebiotic molecules observed in different astrophysical objects. Among these molecules methanol is one of the most widely observed, and its reaction with OH has been studied by several groups showing a fast increase in the rate constant under 100 K. Transition state theory doesn't reproduce this behavior and here dynamical calculations are performed on a new full dimensional potential energy surface developed for this purpose. The calculated classical reactive cross sections show an increase at low collision energies due to a complex forming mechanism. However, the calculated rate constant at temperatures below 100 K remains lower than the observed one. Quantum effects are likely responsible for the measured behavior at low temperatures.
Collapse
Affiliation(s)
- Octavio Roncero
- Instituto de Física Fundamental (IFF-CSIC), C.S.I.C., Serrano 123, Madrid 28006, Spain.
| | | | | |
Collapse
|
3
|
Xu T, Wu H, Zhang LL, Wang XL, Zhao J, Meng QT. State-to-state quantum dynamics of F(2P)+HO(2Π)→O(3P)+HF(1Σ+) reaction on 13A′ and 23A′′ surfaces. Mol Phys 2018. [DOI: 10.1080/00268976.2018.1512719] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
Affiliation(s)
- Ting Xu
- School of Physics and Electronics, Shandong Normal University, Jinan, People’s Republic of China
| | - Hui Wu
- College of Science, Jiangnan University, Wuxi, People’s Republic of China
| | - Lu-Lu Zhang
- School of Physics and Electronics, Shandong Normal University, Jinan, People’s Republic of China
| | - Xian-Long Wang
- School of Physics and Electronics, Shandong Normal University, Jinan, People’s Republic of China
| | - Juan Zhao
- College of Science, Shandong Jiaotong University, Jinan, People’s Republic of China
| | - Qing-Tian Meng
- School of Physics and Electronics, Shandong Normal University, Jinan, People’s Republic of China
| |
Collapse
|
4
|
Zanchet A, Del Mazo P, Aguado A, Roncero O, Jiménez E, Canosa A, Agúndez M, Cernicharo J. Full dimensional potential energy surface and low temperature dynamics of the H 2CO + OH → HCO + H 2O reaction. Phys Chem Chem Phys 2018; 20:5415-5426. [PMID: 28959812 PMCID: PMC6031300 DOI: 10.1039/c7cp05307j] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
A new method is proposed to analytically represent the potential energy surface of reactions involving polyatomic molecules capable of accurately describing long-range interactions and saddle points, needed to describe low-temperature collisions. It is based on two terms, a reactive force field term and a many-body term. The reactive force field term accurately describes the fragments, long-range interactions among them and the saddle points for reactions. The many-body term increases the desired accuracy everywhere else. This method has been applied to the OH + H2CO → H2O + HCO reaction, giving a barrier of 27.4 meV. The simulated classical rate constants with this potential are in good agreement with recent experimental results [Ocaña et al., Astrophys. J., 2017, submitted], showing an important increase at temperatures below 100 K. The reaction mechanism is analyzed in detail here, and explains the observed behavior at low energy by the formation of long-lived collision complexes, with roaming trajectories, with a capture observed for very long impact parameters, >100 a.u., determined by the long-range dipole-dipole interaction.
Collapse
Affiliation(s)
- Alexandre Zanchet
- Instituto de Física Fundamental, Consejo Superior de Investigaciones Científicas, c/Serrano 123, 28006 Madrid, Spain.
| | | | | | | | | | | | | | | |
Collapse
|
5
|
Ocaña AJ, Jiménez E, Ballesteros B, Canosa A, Antiñolo M, Albaladejo J, Agúndez M, Cernicharo J, Zanchet A, del Mazo P, Roncero O, Aguado A. Is the gas-phase OH+H 2CO reaction a source of HCO in interstellar cold dark clouds? A kinetic, dynamic and modelling study. THE ASTROPHYSICAL JOURNAL 2017; 850:28. [PMID: 29880977 PMCID: PMC5988043 DOI: 10.3847/1538-4357/aa93d9] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
Chemical kinetics of neutral-neutral gas-phase reactions at ultralow temperatures is a fascinating research subject with important implications on the chemistry of complex organic molecules in the interstellar medium (T∼10-100K). Scarce kinetic information is currently available for this kind of reactions at T<200 K. In this work we use the CRESU (Cinétique de Réaction en Ecoulement Supersonique Uniforme, which means Reaction Kinetics in a Uniform Supersonic Flow) technique to measure for the first time the rate coefficients (k) of the gas-phase OH+H2CO reaction between 22 and 107 K. k values greatly increase from 2.1×10-11 cm3 s-1 at 107 K to 1.2×10-10 cm3 s-1 at 22 K. This is also confirmed by quasi-classical trajectories (QCT) at collision energies down to 0.1 meV performed using a new full dimension and ab initio potential energy surface, recently developed which generates highly accurate potential and includes long range dipole-dipole interactions. QCT calculations indicate that at low temperatures HCO is the exclusive product for the OH+H2CO reaction. In order to revisit the chemistry of HCO in cold dense clouds, k is reasonably extrapolated from the experimental results at 10K (2.6×10-10 cm3 s-1). The modeled abundances of HCO are in agreement with the observations in cold dark clouds for an evolving time of 105-106 yrs. The different sources of production of HCO are presented and the uncertainties in the chemical networks discussed. This reaction can be expected to be a competitive process in the chemistry of prestellar cores. The present reaction is shown to account for a few percent of the total HCO production rate. Extensions to photodissociation regions and diffuse clouds environments are also commented.
Collapse
Affiliation(s)
- A. J. Ocaña
- Departamento de Química Física. Facultad de Ciencias y Tecnologías Químicas. Universidad de Castilla-La Mancha. Avda. Camilo José Cela 1B. 13071, Ciudad Real, Spain
| | - E. Jiménez
- Departamento de Química Física. Facultad de Ciencias y Tecnologías Químicas. Universidad de Castilla-La Mancha. Avda. Camilo José Cela 1B. 13071, Ciudad Real, Spain
- Instituto de Investigación en Combustión y Contaminación Atmosférica. Universidad de Castilla-La Mancha. Camino de Moledores s/n. 13071, Ciudad Real, Spain
| | - B. Ballesteros
- Departamento de Química Física. Facultad de Ciencias y Tecnologías Químicas. Universidad de Castilla-La Mancha. Avda. Camilo José Cela 1B. 13071, Ciudad Real, Spain
- Instituto de Investigación en Combustión y Contaminación Atmosférica. Universidad de Castilla-La Mancha. Camino de Moledores s/n. 13071, Ciudad Real, Spain
| | - A. Canosa
- Institut de Physique de Rennes, UMR 6251 CNRS-Université de Rennes 1. Campus de Beaulieu, Bât 11C, 263 Av. Général Leclerc, 35042, Rennes, France
| | - M. Antiñolo
- Instituto de Investigación en Combustión y Contaminación Atmosférica. Universidad de Castilla-La Mancha. Camino de Moledores s/n. 13071, Ciudad Real, Spain
| | - J. Albaladejo
- Departamento de Química Física. Facultad de Ciencias y Tecnologías Químicas. Universidad de Castilla-La Mancha. Avda. Camilo José Cela 1B. 13071, Ciudad Real, Spain
- Instituto de Investigación en Combustión y Contaminación Atmosférica. Universidad de Castilla-La Mancha. Camino de Moledores s/n. 13071, Ciudad Real, Spain
| | - M. Agúndez
- Instituto de Ciencia de Materiales de Madrid. Consejo Superior de Investigaciones Científicas. C/ Sor Juana Inés de la Cruz, 3. 28049, Cantoblanco, Madrid, Spain
| | - J. Cernicharo
- Instituto de Ciencia de Materiales de Madrid. Consejo Superior de Investigaciones Científicas. C/ Sor Juana Inés de la Cruz, 3. 28049, Cantoblanco, Madrid, Spain
| | - A. Zanchet
- Instituto de Física Fundamental, CSIC, C/ Serrano, 123, 28006 Madrid, Spain
| | - P. del Mazo
- Instituto de Física Fundamental, CSIC, C/ Serrano, 123, 28006 Madrid, Spain
| | - O. Roncero
- Instituto de Física Fundamental, CSIC, C/ Serrano, 123, 28006 Madrid, Spain
| | - A. Aguado
- Departamento de Química Física Aplicada (UAM), Unidad Asociada IFF-CSIC, Facultad de Ciencias C-XIV, Universidad Autónoma de Madrid, 28049, Madrid, Spain
| |
Collapse
|
6
|
Product rotational angular momentum polarization of H+FCl (v=0-5; j=0, 3, 6, 9) → HF+Cl and HCl+F at Erel=0.5-20 kcal mol(-1). J Mol Model 2015; 21:309. [PMID: 26564146 DOI: 10.1007/s00894-015-2827-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/06/2015] [Accepted: 09/22/2015] [Indexed: 10/22/2022]
Abstract
The rotational angular momentum polarizations of product molecules of the title reactions on the ground potential energy surface 1 (2)A' of DHTSN [Deskevic et al. J Chem Phys 2006, 124, 224303] have been studied using the quasi-classical trajectory method. Reaction dynamic results of the HF product channel comparing with another channel of HCl with 100,000 trajectories can be accurately resolved. We show the value of the polar p(ϑr) in the range of 0° ≤ ϑr ≤ 180(°), azimuthal p(φr) in the range of 0° ≤ φr ≤ 360(°), and dihedral p(ϑr, φr) in the ranges of 0(°) ≤ ϑr ≤ 180(°) and 0(°) ≤ φr ≤ 360(°); the angular distributions of the product molecules HF and HCl at relative Erel = 0.5, 1, 2, 5, 10, 15, and 20 kcal mol(-1); and four polarization-dependent differential cross sections (PDDCSs) of HF and HCl at Erel = 0.5, 1, 2, 5, 10, and 15 kcal mol(-1). p(φr) distributions at v = 0-5, and j = 0, 3, 6, 9 at every Erel are plotted cylindrically together. The stereo dynamic transformation reaction dependent upon the rovibrational states of the reactant molecule FCl and its relative translational energies around 0.5-5 kcal mol(-1) can be significantly differentiated. Translational and rovibrational enhancements of the title reactions on both early barrier potential energy surfaces have been shown in great detail and clarified. Reaction mechanisms of forward and backward scattering of the product molecules HF and HCl, respectively, have been obtained. Graphical Abstract H + FCl → either HF + Cl (left) or HCl + F (right) is moving along a trajectory on the respective PES.
Collapse
|
7
|
González-Lavado E, Corchado JC, Espinosa-Garcia J. The hydrogen abstraction reaction O(3P) + CH4: a new analytical potential energy surface based on fit to ab initio calculations. J Chem Phys 2014; 140:064310. [PMID: 24527918 DOI: 10.1063/1.4864358] [Citation(s) in RCA: 36] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Based exclusively on high-level ab initio calculations, a new full-dimensional analytical potential energy surface (PES-2014) for the gas-phase reaction of hydrogen abstraction from methane by an oxygen atom is developed. The ab initio information employed in the fit includes properties (equilibrium geometries, relative energies, and vibrational frequencies) of the reactants, products, saddle point, points on the reaction path, and points on the reaction swath, taking especial caution respecting the location and characterization of the intermediate complexes in the entrance and exit channels. By comparing with the reference results we show that the resulting PES-2014 reproduces reasonably well the whole set of ab initio data used in the fitting, obtained at the CCSD(T) = FULL/aug-cc-pVQZ//CCSD(T) = FC/cc-pVTZ single point level, which represents a severe test of the new surface. As a first application, on this analytical surface we perform an extensive dynamics study using quasi-classical trajectory calculations, comparing the results with recent experimental and theoretical data. The excitation function increases with energy (concave-up) reproducing experimental and theoretical information, although our values are somewhat larger. The OH rotovibrational distribution is cold in agreement with experiment. Finally, our results reproduce experimental backward scattering distribution, associated to a rebound mechanism. These results lend confidence to the accuracy of the new surface, which substantially improves the results obtained with our previous surface (PES-2000) for the same system.
Collapse
Affiliation(s)
| | - Jose C Corchado
- Departamento de Química Física, Universidad de Extremadura, 06071 Badajoz, Spain
| | | |
Collapse
|
8
|
Zhao D, He X, Guo W. Stereodynamics investigation of F + HO → HF + O(1D) on the ground singlet potential energy surface by means of the quasi-classical trajectory method. CAN J CHEM 2014. [DOI: 10.1139/cjc-2013-0401] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
Abstract
The stereodynamics calculation of F + HO → HF + O(1D) was carried out using the quasi-classical trajectory method on the 11A′ potential energy surface provided by Gomez-Carrasco et al. (Chem. Phys. Lett. 2007, 435, 188). The effect of the collision energy, isotopic substitution, and different initial ro-vibrational states on the reaction is discussed. It is found that for the initial ground state of HO (v = 0, j = 0), the degree of the forward scattering and the product polarizations remarkably change as the collision energy varies. Isotopic effect leads to the increase of alignment and decrease of orientation of product rotational angular momentum. Moreover, the P(θr) distribution and P(φr) distribution change noticeably by varying the initial vibrational number. The initial vibrational excitation plays a more important role in the enhancement of alignment and orientation distribution of j′ for the title reaction. Although the influence of the initial rotational excitation effect on the aligned and oriented distribution of product is not stronger than that of the initial vibrational excitation effect, the initial rotational excitation makes the alignment of the product rotational angular momentum decrease to some extent. The probabilities show that the reactivity of the title reaction strongly depends on the initial vibrational state.
Collapse
Affiliation(s)
- Dan Zhao
- School of Chemical Engineering, Dalian University of Technology, Dalian 116023, China
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China
| | - Xiaohu He
- School of Chemical Engineering, Dalian University of Technology, Dalian 116023, China
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China
| | - Wei Guo
- School of Electric Engineering, University of South China, Hengyang 421001, China
| |
Collapse
|
9
|
Koner D, Panda AN. Quantum dynamical study of the He + NeH+ reaction on a new analytical potential energy surface. J Phys Chem A 2013; 117:13070-8. [PMID: 24256154 DOI: 10.1021/jp408550c] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
An analytical potential energy surface (PES) for the ground state of the [HeHNe](+) system has been constructed from a set of 19,605 ab initio data points, obtained from coupled cluster singles and doubles with perturbative triples correction calculations and the aug-cc-pVQZ basis set. The PES is based on the many-body expansion form proposed by Aguado and Paniagua (J. Chem. Phys. 1992, 96, 1265), and it has a root-mean-square error of 0.03 kcal/mol. The minimum energy pathways (MEPs) for different Ne-H-He angles are calculated, and it is found that the MEP for 180° (linear) goes through the deepest potential energy well. Preliminary quantum dynamical studies are performed for the He + NeH(+) (v = 0-2, j = 0-3) → HeH(+) + Ne reaction in the 0.0-0.5 eV collision energy range. Quantum calculations are carried out using a time-dependent wave packet method within the centrifugal sudden approximation. Reaction probabilities exhibit strong oscillatory behavior arising because of the metastable [HeHNe](+). Vibrational excitation has been found to enhance the reaction cross sections.
Collapse
Affiliation(s)
- Debasish Koner
- Department of Chemistry, Indian Institute of Technology Guwahati , Guwahati 781039, India
| | | |
Collapse
|
10
|
Gómez-Carrasco S, Bulut N, Bañares L, Roncero O. Wave packet calculations on nonadiabatic effects for the O(3P)+HF(1Σ+) reaction under hyperthermal conditions. J Chem Phys 2012; 137:114309. [PMID: 22998264 DOI: 10.1063/1.4753811] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
We present wave packet calculations of total and state-to-state reaction probabilities and integral cross sections for the nonadiabatic dynamics of the O((3)P)+HF → F((2)P)+OH((2)Π) reaction at hyperthermal collision energies ranging from 1.2 to 2.4 eV. The validity of the centrifugal sudden approximation is discussed for the title reaction and a comprehensive investigation of the influence of nonadiabatic effects on the dynamics of this reactive system at high (hyperthermal) collision energies is presented. In general, nonadiabatic effects are negligible for averaged observables, such as total reaction probabilities and integral cross sections, but they are clearly observed in detailed observables such as rotationally state-resolved reaction probabilities. A critical discussion of nonadiabatic effects on the dynamics of the title reaction is carried out by comparing with the reverse reaction and the characteristics of the adiabatic and diabatic potential energy surfaces involved.
Collapse
Affiliation(s)
- S Gómez-Carrasco
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad de Salamanca, Plaza de los Caídos s/n, 37008 Salamanca, Spain
| | | | | | | |
Collapse
|
11
|
ZHAO JUAN. EFFECT OF THE REACTANT ALIGNMENT ON STEREODYNAMICS FOR THE REACTION F + OH. JOURNAL OF THEORETICAL & COMPUTATIONAL CHEMISTRY 2012. [DOI: 10.1142/s0219633611006621] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
Quasi-classical trajectory (QCT) calculations are performed for the reaction F + OH → HF + O based on the adiabatic potential-energy surface (PES) of the 13A″ triplet state. The reaction probability as a function of incident angle has been presented. The differential cross sections (DCSs), the distribution of angle between k and j′, P(θr) and the distribution of dihedral angle denoting k – k′ – j′ correlation, P(ϕr) have also been calculated at the different incident angles in the center-of-mass (CM) frame, respectively. The effects of reactant alignment on stereodynamics of the reaction are firstly revealed, which provides the theoretical foundation for the related experiment and enriches the theories of the stereodynamics.
Collapse
Affiliation(s)
- JUAN ZHAO
- School of Pharmaceutical Science and Technology, Dalian University of Technology, Dalian 116024, P. R. China
| |
Collapse
|
12
|
Zanchet A, González-Lezana T, Roncero O, Jorfi M, Honvault P, Hankel M. An accurate study of the dynamics of the C+OH reaction on the second excited 14A″ potential energy surface. J Chem Phys 2012; 136:164309. [DOI: 10.1063/1.4705426] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
13
|
YUE XIANFANG, GENG JIANPEI. STATE-TO-STATE REACTION DYNAMICS OF THE D(H) + FO → OD(OH) + F REACTIONS ON THE LOWEST TWO TRIPLET STATES. JOURNAL OF THEORETICAL & COMPUTATIONAL CHEMISTRY 2012. [DOI: 10.1142/s0219633612500034] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
Quasiclassical trajectory (QCT) calculations have been performed to investigate the reaction dynamics of the title reactions at a state-to-state level. The recently developed adiabatic potential energy surfaces (PESs) of the two triplet electronic states of 13A″ and 13A′ (Gómez-Carrasco S. et al., J Chem Phys121:4605, 2004; J Chem Phys123:114310, 2005) are employed in the present QCT calculations. Product vibrational and rotational state distributions have been calculated at three collision energies of 0.5, 0.7 and 0.9 eV. The product vibrational state distributions are found to be Gaussian distributions. Both the vibrational and rotational state distributions depend clearly on the isotope substitution, electronic PESs and collision energies. The observed phenomena are probably attributed to the competition of direct and indirect reactions.
Collapse
Affiliation(s)
- XIAN-FANG YUE
- Department of Physics and Information Engineering, Jining University, Jining, 273155, P. R., China
| | - JIAN-PEI GENG
- School of Physics, Shandong University, Jinan 250100, P. R. China
| |
Collapse
|
14
|
Zhao D, Chu T, Hao C. Vector correlations in the F + HO → HF + O reaction and its isotopic variant. J Mol Model 2012; 18:3283-9. [PMID: 22252836 DOI: 10.1007/s00894-012-1351-1] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/11/2011] [Accepted: 01/01/2012] [Indexed: 11/29/2022]
Abstract
The F + H(D)O → HF(DF) + O reactions have been studied using quasi-classical trajectory (QCT) calculation method, based on the three different potential energy surfaces (PESs) of Gomez-Carrasco et al. (J Chem Phys 2004, 121:4605; J Chem Phys 2005, 123:114310; Chem Phys Lett 2007, 435:188). Facilitated with the analysis of the QCT results, the pictures for product scattering and product polarizations have been presented to investigate the vector correlations in the two reactions, with effects of isotope substitution and electronic state as well as collision energy being revealed at a chemical stereodynamical level.
Collapse
Affiliation(s)
- Dan Zhao
- School of Chemical Engineering, Dalian University of Technology, Dalian 116023, China
| | | | | |
Collapse
|
15
|
Sun Z, Yang W, Zhang DH. Higher-order split operator schemes for solving the Schrödinger equation in the time-dependent wave packet method: applications to triatomic reactive scattering calculations. Phys Chem Chem Phys 2012; 14:1827-45. [DOI: 10.1039/c1cp22790d] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
16
|
Wang Y, Tian B, Qu L, Chen J, Li H. Quasiclassical Trajectory Calculations for the Reaction Ne + H 2+→ NeH ++ H. B KOREAN CHEM SOC 2011. [DOI: 10.5012/bkcs.2011.32.12.4210] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
|
17
|
Bulut N, Roncero O, Jorfi M, Honvault P. Accurate time dependent wave packet calculations for the N + OH reaction. J Chem Phys 2011; 135:104307. [DOI: 10.1063/1.3633240] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
18
|
A quasiclassical trajectory analysis of stereodynamics of the H + FCl (v = 0 – 3, j = 0 – 3) → HCl + F reaction. J CHEM SCI 2011. [DOI: 10.1007/s12039-011-0063-x] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/17/2022]
|
19
|
The rate constants of the H+FO reaction and its isotope variant on two electronic states. Chem Phys 2011. [DOI: 10.1016/j.chemphys.2011.01.010] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
20
|
Zhao J. Quasi-classical trajectory study of the O(1D) + HF reaction dynamics on 11A′ potential energy surface. CAN J CHEM 2011. [DOI: 10.1139/v11-059] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
Abstract
The quasi-classical trajectory (QCT) calculations for the title reaction were carried out using the recently developed, accurate potential energy surface (PES) of the [Formula: see text] singlet state of the OHF system The integral cross section and the product rotational alignment factor [Formula: see text] were calculated as a function of collision energy. In addition, I discovered the effect of isotopic substitution on stereodynamics for the title reaction, and the influence of the rotation excitation of the reagent on stereodynamics is also presented. Both the scalar and vector properties of the reaction O(1D) + HF → OH + F(2P) are studied in this paper. It was found that the reaction is mainly controlled by an indirect reaction mechanism, and that the deep noncollinear insertion HOF well has a great impact on the dynamics of the reaction. The conclusions drawn in this paper will draw from references to similar reactions, and provide a theoretical foundation for related experiments.
Collapse
Affiliation(s)
- Juan Zhao
- School of Pharmaceutical Science and Technology, Dalian University of Technology, Dalian, 116024, China
| |
Collapse
|
21
|
Kurban M, Gogtas F, Karabulut E, Tutuk R. Quantum dynamics of reaction on the state. Mol Phys 2011. [DOI: 10.1080/00268976.2011.554327] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
22
|
Gogtas F, Tutuk R, Kurban M. Time-dependent quantum study of H(2S) + FO(2Π) → OH(2Π) + F(2P) reaction on the 13A′ and 13A″ states. J Comput Chem 2010; 31:2607-11. [DOI: 10.1002/jcc.21555] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
|
23
|
Zanchet A, González-Lezana T, Aguado A, Gómez-Carrasco S, Roncero O. Nonadiabatic state-to-state reactive collisions among open shell reactants with conical intersections: the OH((2)Pi) + F((2)P) example. J Phys Chem A 2010; 114:9733-42. [PMID: 20465247 DOI: 10.1021/jp101914a] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Accurate wave packet calculations on the OH((2)Pi) + F((2)P) → O((3)P) + HF((1)Sigma(+)) reactive collisions are performed using a recently proposed coupled diabatic states. Adiabatic and nonadiabatic dynamics are compared in detail, analyzing the final state distribution of products. It is found that with the new surfaces a significant increase of the rate constant is obtained, with noticeable nonadiabatic effects. The inclusion of the spin-orbit splittings for the calculation of the electronic partition function produces an important increase of the reaction rate constants, yielding a rather good agreement with the experimental results. It is also concluded that spin-orbit couplings are also necessary in the entrance channel to describe this reaction.
Collapse
Affiliation(s)
- Alexandre Zanchet
- Unidad Asociada UAM-CSIC, Instituto de Física Fundamental, CSIC, Serrano 123, 28006 Madrid, Spain
| | | | | | | | | |
Collapse
|
24
|
Zhao D, Zhang TY, Chu TS. A quasi-classical trajectory (QCT) study of the H + OF reaction stereodynamics. CAN J CHEM 2010. [DOI: 10.1139/v10-063] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
Abstract
Based on the global three-dimensional adiabatic potential surface of the 13A′ excited electronic state (J. Chem. Phys. 2005, 123, 114310) of the OHF system, we investigated the H + OF → OH + F/HF + O reaction stereodynamics by using the quasi-classical trajectory (QCT) method. The four polarization-dependent differential cross sections (PDDCSs) and the three angular distributions P(θr), P([Formula: see text]), P(θr, [Formula: see text]) were calculated at a low collision energy of 0.48 eV for both product channels. The results indicated that the products are backward-scattering on the triplet state, and the product rotational angular momenta are aligned or oriented. Moreover, there is a remarkable difference between the polarization behaviors of the two product channels. Product orientation exhibited by the OH + F channel is found to be absent in the HF + O channel at this collision energy albeit the latter shows stronger alignment along the direction perpendicular to the reagent relative velocity k than OH + F.
Collapse
Affiliation(s)
- Dan Zhao
- School of Chemical Engineering, Dalian University of Technology, Dalian 116023, China
- Institute for Computational Sciences and Engineering, Qingdao University, Qingdao 266071, China
| | - Tian Yu Zhang
- School of Chemical Engineering, Dalian University of Technology, Dalian 116023, China
- Institute for Computational Sciences and Engineering, Qingdao University, Qingdao 266071, China
| | - Tian Shu Chu
- School of Chemical Engineering, Dalian University of Technology, Dalian 116023, China
- Institute for Computational Sciences and Engineering, Qingdao University, Qingdao 266071, China
| |
Collapse
|
25
|
Han B, Zong F, Wang C, Ma W, Zhou J. Product polarization on the 3A″ electronic state in the H+FO reaction and its isotope variant. Chem Phys 2010. [DOI: 10.1016/j.chemphys.2010.06.028] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
26
|
Quasi-classical trajectory approach to the stereo-dynamics of the reaction F + HO → HF + O. Sci China Chem 2010. [DOI: 10.1007/s11426-010-0043-x] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
27
|
Chu T. Quantum mechanics and quasiclassical study of the H/D+FO --> OH/OD+F, HF/DF+O reactions: chemical stereodynamics. J Comput Chem 2009; 31:1385-96. [PMID: 19859917 DOI: 10.1002/jcc.21423] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Abstract
The time-dependent quantum wave packet and the quasi-classical trajectory (QCT) calculations for the title reactions are carried out using three recent-developed accurate potential energy surfaces of the 1(1)A', 1(3)A', and 1(3)A'' states. The two commonly used polarization-dependent differential cross sections, dsigma(00)/domega(t), dsigma(20)/domega(t), with omega(t) being the polar coordinates of the product velocity omega', and the three angular distributions, P(theta(r)), P(Phi(r)), and P(theta(r),Phi(r)), with theta(r), Phi(r) being the polar angles of the product angular momentum, are generated in the center-of-mass frame using the QCT method to gain insight into the alignment and the orientation of the product molecules. Influences of the potential energy surface, the collision energy, and the isotope mass on the stereodynamics are shown and discussed. Validity of the QCT calculation has been examined and proved in the comparison with the quantum wave packet calculation.
Collapse
Affiliation(s)
- Tianshu Chu
- Institute for Computational Sciences and Engineering, Qingdao University, Qingdao, 266071, China.
| |
Collapse
|
28
|
|
29
|
Chu TS, Zhang H, Yuan SP, Fu AP, Si HZ, Tian FH, Duan YB. A Comparative Study of the H + FO (v = 0, j = 0) → (OH + F)/(HF + O) Reaction from QM and QCT Methods. J Phys Chem A 2009; 113:3470-5. [DOI: 10.1021/jp9003505] [Citation(s) in RCA: 34] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- T. S. Chu
- Institute for Computational Sciences and Engineering, Qingdao University, Shandong 266071, People’s Republic of China, and Centre for Computational Molecular Science, Chemistry Building, The University of Queensland, Qld 4072, Brisbane, Australia
| | - H. Zhang
- Institute for Computational Sciences and Engineering, Qingdao University, Shandong 266071, People’s Republic of China, and Centre for Computational Molecular Science, Chemistry Building, The University of Queensland, Qld 4072, Brisbane, Australia
| | - S. P. Yuan
- Institute for Computational Sciences and Engineering, Qingdao University, Shandong 266071, People’s Republic of China, and Centre for Computational Molecular Science, Chemistry Building, The University of Queensland, Qld 4072, Brisbane, Australia
| | - A. P. Fu
- Institute for Computational Sciences and Engineering, Qingdao University, Shandong 266071, People’s Republic of China, and Centre for Computational Molecular Science, Chemistry Building, The University of Queensland, Qld 4072, Brisbane, Australia
| | - H. Z. Si
- Institute for Computational Sciences and Engineering, Qingdao University, Shandong 266071, People’s Republic of China, and Centre for Computational Molecular Science, Chemistry Building, The University of Queensland, Qld 4072, Brisbane, Australia
| | - F. H. Tian
- Institute for Computational Sciences and Engineering, Qingdao University, Shandong 266071, People’s Republic of China, and Centre for Computational Molecular Science, Chemistry Building, The University of Queensland, Qld 4072, Brisbane, Australia
| | - Y. B. Duan
- Institute for Computational Sciences and Engineering, Qingdao University, Shandong 266071, People’s Republic of China, and Centre for Computational Molecular Science, Chemistry Building, The University of Queensland, Qld 4072, Brisbane, Australia
| |
Collapse
|
30
|
|
31
|
|
32
|
Gogtas F. Time-dependent quantum study of the kinetics of the H(2S) + FO(2II) → OH(2II) + F(2P) reaction. J Comput Chem 2008; 29:1889-94. [DOI: 10.1002/jcc.20948] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
|
33
|
Gómez-Carrasco S, Aguado A, Paniagua M, Roncero O. Coupled diabatic potential energy surfaces for studying the nonadiabatic dynamics at conical intersections in angular resolved photodetachment simulations of OHF--->OHF+e-. J Chem Phys 2007; 125:164321. [PMID: 17092087 DOI: 10.1063/1.2363988] [Citation(s) in RCA: 31] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
An energy-based method is proposed for the diabatization of the OH(2Pi)+F(2P)-->O(3P)+HF(1Sigma+) reaction. It is demonstrated that the diabatic representation obtained is regularized, i.e., the residual derivative couplings do not present singularities at the conical intersections appearing along the reaction path. This method only requires the knowledge of the 1,2 3A" and 1 3A' eigenvalues and does not require any adjustable parameter. Thus, many convergence problems arising in other derivative-based diabatization methods are avoided, and the description of the configuration space along the reaction path is enormously simplified. Three-dimensional coupled diabatic energy surfaces are obtained by an interpolation procedure using approximately 4000 accurate ab initio points. The angular resolved photodetachment cross sections are obtained in the diabatic and adiabatic representations using a wave packet method. An excellent agreement is obtained with recent experimental data [D. M. Neumark, Phys. Chem. Chem. Phys. 7, 433 (2005)] for high electron kinetic energies where only the triplet electronic states contribute.
Collapse
Affiliation(s)
- Susana Gómez-Carrasco
- Unidad Asociada UAM-CSIC, Instituto de Matemáticas y Física Fundamental, CSIC, Serrano 123, 28006 Madrid, Spain
| | | | | | | |
Collapse
|
34
|
Transition state spectroscopy of open shell systems: Angle-resolved photodetachment spectra for the adiabatic singlet states of OHF. J Photochem Photobiol A Chem 2007. [DOI: 10.1016/j.jphotochem.2007.01.027] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
|
35
|
Gómez-Carrasco S, Luz Hernández M, Alvariño JM. Quantum and quasiclassical state-selected O(1D) + HF reaction dynamics and kinetics on a new MRCI ground singlet potential energy surface. Chem Phys Lett 2007. [DOI: 10.1016/j.cplett.2006.12.073] [Citation(s) in RCA: 29] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
|