1
|
Sathyamurthy N, Mahapatra S. Time-dependent quantum mechanical wave packet dynamics. Phys Chem Chem Phys 2020; 23:7586-7614. [PMID: 33306771 DOI: 10.1039/d0cp03929b] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Starting from a model study of the collinear (H, H2) exchange reaction in 1959, the time-dependent quantum mechanical wave packet (TDQMWP) method has come a long way in dealing with systems as large as Cl + CH4. The fast Fourier transform method for evaluating the second order spatial derivative of the wave function and split-operator method or Chebyshev polynomial expansion for determining the time evolution of the wave function for the system have made the approach highly accurate from a practical point of view. The TDQMWP methodology has been able to predict state-to-state differential and integral reaction cross sections accurately, in agreement with available experimental results for three dimensional (H, H2) collisions, and identify reactive scattering resonances too. It has become a practical computational tool in predicting the observables for many A + BC exchange reactions in three dimensions and a number of larger systems. It is equally amenable to determining the bound and quasi-bound states for a variety of molecular systems. Just as it is able to deal with dissociative processes (without involving basis set expansion), it is able to deal with multi-mode nonadiabatic dynamics in multiple electronic states with equal ease. We present an overview of the method and its strength and limitations, citing examples largely from our own research groups.
Collapse
|
2
|
Mao Y, Yuan J, Yang Z, Chen M. Quantum dynamics studies of isotope effects in the Mg +(3p) + HD → MgH +/MgD + + D/H insertion reaction. Sci Rep 2020; 10:3410. [PMID: 32098984 PMCID: PMC7042225 DOI: 10.1038/s41598-020-60033-2] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/04/2019] [Accepted: 02/03/2020] [Indexed: 11/17/2022] Open
Abstract
The time-dependent wave packet quantum dynamics studies for the Mg+(3p) + HD → MgH+/MgD+ + D/H diabatic reaction are carried out for the first time on recently developed diabatic YHWCH potential energy surfaces [Phys. Chem. Chem. Phys., 2018, 20, 6638-6647]. The results of reaction probabilities and total integral cross sections show a dramatic preference to the formation of MgD+ over MgH+ owing to the insertion reaction mechanism in the title reaction. The MgD+/MgH+ branching ratio witnesses a monotonic decrease from 10.58 to 3.88 at collision energy range of 0.01 to 0.20 eV, and at the collision energy of 0.114 eV, it is close to the experimental value of 5. The rovibrational state-resolved ICSs of the two channels show the products MgD+ have higher vibrational excitation and hotter rotational state distributions. The opacity function P(J) suggests that the MgH+ + D channel and MgD+ + H channel are dominated by high-b and low-b collisions, respectively. Both forward and backward scattering peaks are found in the differential cross section curves, whereas the angle distributions of products are not strictly forward-backward symmetric because of the short lifetime of the complex in the reaction.
Collapse
Affiliation(s)
- Ye Mao
- Key Laboratory of Materials Modification by Laser, Electron, and Ion Beams (Ministry of Education), School of Physics, Dalian University of Technology, Dalian, 116024, PR China
| | - Jiuchuang Yuan
- Key Laboratory of Materials Modification by Laser, Electron, and Ion Beams (Ministry of Education), School of Physics, Dalian University of Technology, Dalian, 116024, PR China
| | - Zijiang Yang
- Key Laboratory of Materials Modification by Laser, Electron, and Ion Beams (Ministry of Education), School of Physics, Dalian University of Technology, Dalian, 116024, PR China
| | - Maodu Chen
- Key Laboratory of Materials Modification by Laser, Electron, and Ion Beams (Ministry of Education), School of Physics, Dalian University of Technology, Dalian, 116024, PR China.
| |
Collapse
|
3
|
Zhang C, Zheng Y, Cao J, Bian W. Quasiclassical trajectory study of the C(1D) + HD reaction. RSC Adv 2017. [DOI: 10.1039/c7ra03966b] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022] Open
Abstract
Isotopic branching ratios are investigated by detailed quasiclassical trajectory calculations on our recent singlet ground and excited potential energy surfaces.
Collapse
Affiliation(s)
- Chunfang Zhang
- CAS Research/Education Center for Excellence in Molecular Sciences
- Institute of Chemistry
- Chinese Academy of Sciences
- Beijing 100190
- China
| | - Yujun Zheng
- School of Physics
- Shandong University
- Jinan 250100
- China
| | - Jianwei Cao
- CAS Research/Education Center for Excellence in Molecular Sciences
- Institute of Chemistry
- Chinese Academy of Sciences
- Beijing 100190
- China
| | - Wensheng Bian
- CAS Research/Education Center for Excellence in Molecular Sciences
- Institute of Chemistry
- Chinese Academy of Sciences
- Beijing 100190
- China
| |
Collapse
|
4
|
Yuan J, Cheng D, Chen M. Time-dependent wave packet and quasiclassical trajectory studies of the Au + HD reaction: competition between the reactive channels. RSC Adv 2014. [DOI: 10.1039/c4ra06297c] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022] Open
|
5
|
Gamallo P, Akpinar S, Defazio P, Petrongolo C. Quantum Dynamics of the Reaction H(2S) + HeH+(X1Σ+) → H2+(X2Σg+) + He(1S) from Cold to Hyperthermal Energies: Time-Dependent Wavepacket Study and Comparison with Time-Independent Calculations. J Phys Chem A 2014; 118:6451-6. [DOI: 10.1021/jp5023289] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Pablo Gamallo
- Departament
de Quı́mica Fı́sica, Institut
de Quı́mica Teòrica i Computacional, Universitat de Barcelona, C/Martı́ i Franquès 1, 08028 Barcelona, Spain
| | - Sinan Akpinar
- Department
of Physics, Firat University, 23169 Elazig, Turkey
| | - Paolo Defazio
- Dipartimento
di Biotecnologie, Chimica, e Farmacia, Università di Siena, Via A. Moro
2, 53100 Siena, Italy
| | - Carlo Petrongolo
- Istituto per i
Processi Chimico-Fisici, Consiglio Nazionale delle Ricerche, Via G. Moruzzi 1, 56124 Pisa, Italy
| |
Collapse
|
6
|
Yang H, Ge M, Zheng Y. Adiabatic wavepacket dynamics study of the N + NH → N2+ H reaction on the ground-state potential energy surface. Chem Phys Lett 2014. [DOI: 10.1016/j.cplett.2013.12.017] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
7
|
Pradhan GB, Balakrishnan N, Kendrick BK. Ultracold collisions of O(1D) and H2: The effects of H2vibrational excitation on the production of vibrationally and rotationally excited OH. J Chem Phys 2013; 138:164310. [DOI: 10.1063/1.4802476] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
|
8
|
|
9
|
GRAY STEPHENK. CHEMICAL REACTION DYNAMICS WITH REAL WAVE PACKETS. JOURNAL OF THEORETICAL & COMPUTATIONAL CHEMISTRY 2012. [DOI: 10.1142/s0219633602000178] [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
An approach to carrying out accurate quantum dynamics simulations of chemical reactions, termed the real wave packet (RWP) method, is outlined. The method focuses on propagation of just the real part of a complex-valued wave packet, halving computational memory and effort requirements in comparison with comparable high accuracy quantum propagation methods. Applications to 3-atom and 4-atom chemical reactions are reviewed. Potential future directions are indicated.
Collapse
Affiliation(s)
- STEPHEN K. GRAY
- Chemistry Division, Argonne National Laboratory, Argonne, Illinois 60439, USA
| |
Collapse
|
10
|
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]
|
11
|
|
12
|
Sun Z, Lin SY, Zheng Y. Adiabatic and non-adiabatic quantum dynamics calculation of O(1D) + D2 → OD + D reaction. J Chem Phys 2011; 135:234301. [DOI: 10.1063/1.3668084] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
13
|
Balucani N, Casavecchia P, Aoiz FJ, Banares L, Castillo JF, Herrero VJ. Dynamics of the O(1D) D2 reaction: A comparison between crossed molecular beam experiments and quasiclassical trajectory calculations on the lowest three potential energy surfaces. Mol Phys 2011. [DOI: 10.1080/149920500058077] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
Affiliation(s)
- N. Balucani
- a Dipartimento di Chimica , Universita di Perugia , Perugia , 06123 , Italy
| | - P. Casavecchia
- a Dipartimento di Chimica , Universita di Perugia , Perugia , 06123 , Italy
| | - F. J. Aoiz
- b Departamento de Quimica Fisica, F. de Quimica , Universidad Complutense , Madrid , 28040 , Spain
| | - L. Banares
- b Departamento de Quimica Fisica, F. de Quimica , Universidad Complutense , Madrid , 28040 , Spain
| | - J. F. Castillo
- b Departamento de Quimica Fisica, F. de Quimica , Universidad Complutense , Madrid , 28040 , Spain
| | - V. J. Herrero
- c Instituto de Estructura de la Materia. C.S.I.C. , Madrid , 28006 , Spain
| |
Collapse
|
14
|
Li YM, Xie ZM, Dang Q. Theoretical Studies of O(1D) + HD (v = 0, j = 0, 1, 2, 3) → OD(H) + H(D) Reaction. CHINESE J CHEM PHYS 2010. [DOI: 10.1088/1674-0068/23/03/310-312] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
|
15
|
|
16
|
Lin SY, Guo H. Adiabatic and Nonadiabatic State-to-State Quantum Dynamics for O(1D) + H2(X1Σg+, υi = ji = 0) → OH(X2Π, υf, jf) + H(2S) Reaction. J Phys Chem A 2009; 113:4285-93. [DOI: 10.1021/jp810948k] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Shi Ying Lin
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131
| |
Collapse
|
17
|
Lin SY, Guo H. Energy dependence of differential and integral cross sections for O(D1)+H2(υi=0,ji=0)→OH(υf,jf)+H reaction. J Chem Phys 2008; 129:124311. [DOI: 10.1063/1.2981063] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
18
|
Bargueño P, González-Lezana T, Larrégaray P, Bonnet L, Rayez JC, Hankel M, Smith SC, Meijer AJHM. Study of the H+O2 reaction by means of quantum mechanical and statistical approaches: the dynamics on two different potential energy surfaces. J Chem Phys 2008; 128:244308. [PMID: 18601333 DOI: 10.1063/1.2944246] [Citation(s) in RCA: 29] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
The possible existence of a complex-forming pathway for the H+O(2) reaction has been investigated by means of both quantum mechanical and statistical techniques. Reaction probabilities, integral cross sections, and differential cross sections have been obtained with a statistical quantum method and the mean potential phase space theory. The statistical predictions are compared to exact results calculated by means of time dependent wave packet methods and a previously reported time independent exact quantum mechanical approach using the double many-body expansion (DMBE IV) potential energy surface (PES) [Pastrana et al., J. Phys. Chem. 94, 8073 (1990)] and the recently developed surface (denoted XXZLG) by Xu et al. [J. Chem. Phys. 122, 244305 (2005)]. The statistical approaches are found to reproduce only some of the exact total reaction probabilities for low total angular momenta obtained with the DMBE IV PES and some of the cross sections calculated at energy values close to the reaction threshold for the XXZLG surface. Serious discrepancies with the exact integral cross sections at higher energy put into question the possible statistical nature of the title reaction. However, at a collision energy of 1.6 eV, statistical rotationally resolved cross sections managed to reproduce the experimental cross sections for the H+O(2)(v=0,j=1)-->OH(v(')=1,j('))+O process reasonably well.
Collapse
Affiliation(s)
- Pedro Bargueño
- Instituto de Fisica Fundamental (CSIC), Serrano 123, 28006 Madrid, Spain
| | | | | | | | | | | | | | | |
Collapse
|
19
|
Yang H, Han KL, Nanbu S, Nakamura H, Balint-Kurti GG, Zhang H, Smith SC, Hankel M. Quantum Mechanical Calculation of Energy Dependence of OCl/OH Product Branching Ratio and Product Quantum State Distributions for the O(1D) + HCl Reaction on All Three Contributing Electronic State Potential Energy Surfaces. J Phys Chem A 2008; 112:7947-60. [DOI: 10.1021/jp803673y] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Huan Yang
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China, Research Institute for Information Technology, Kyushu University, 6-10-1 Hakozaki, Higashi-ku Fukuoka 812-8581, Japan, Institute for Molecular Science, Myodaiji, Okazaki 444-8585, Japan, School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom, and Centre for Computational Molecular Science, Australian Institute for Bioengineering and
| | - Ke-Li Han
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China, Research Institute for Information Technology, Kyushu University, 6-10-1 Hakozaki, Higashi-ku Fukuoka 812-8581, Japan, Institute for Molecular Science, Myodaiji, Okazaki 444-8585, Japan, School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom, and Centre for Computational Molecular Science, Australian Institute for Bioengineering and
| | - Shinkoh Nanbu
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China, Research Institute for Information Technology, Kyushu University, 6-10-1 Hakozaki, Higashi-ku Fukuoka 812-8581, Japan, Institute for Molecular Science, Myodaiji, Okazaki 444-8585, Japan, School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom, and Centre for Computational Molecular Science, Australian Institute for Bioengineering and
| | - Hiroki Nakamura
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China, Research Institute for Information Technology, Kyushu University, 6-10-1 Hakozaki, Higashi-ku Fukuoka 812-8581, Japan, Institute for Molecular Science, Myodaiji, Okazaki 444-8585, Japan, School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom, and Centre for Computational Molecular Science, Australian Institute for Bioengineering and
| | - Gabriel G. Balint-Kurti
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China, Research Institute for Information Technology, Kyushu University, 6-10-1 Hakozaki, Higashi-ku Fukuoka 812-8581, Japan, Institute for Molecular Science, Myodaiji, Okazaki 444-8585, Japan, School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom, and Centre for Computational Molecular Science, Australian Institute for Bioengineering and
| | - Hong Zhang
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China, Research Institute for Information Technology, Kyushu University, 6-10-1 Hakozaki, Higashi-ku Fukuoka 812-8581, Japan, Institute for Molecular Science, Myodaiji, Okazaki 444-8585, Japan, School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom, and Centre for Computational Molecular Science, Australian Institute for Bioengineering and
| | - Sean C. Smith
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China, Research Institute for Information Technology, Kyushu University, 6-10-1 Hakozaki, Higashi-ku Fukuoka 812-8581, Japan, Institute for Molecular Science, Myodaiji, Okazaki 444-8585, Japan, School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom, and Centre for Computational Molecular Science, Australian Institute for Bioengineering and
| | - Marlies Hankel
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China, Research Institute for Information Technology, Kyushu University, 6-10-1 Hakozaki, Higashi-ku Fukuoka 812-8581, Japan, Institute for Molecular Science, Myodaiji, Okazaki 444-8585, Japan, School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom, and Centre for Computational Molecular Science, Australian Institute for Bioengineering and
| |
Collapse
|
20
|
Affiliation(s)
- Aditya N. Panda
- Department of Chemistry, Indian Institute of Technology, Guwahati 781039, India
| |
Collapse
|
21
|
Yang H, Han KL, Nanbu S, Nakamura H, Balint-Kurti GG, Zhang H, Smith SC, Hankel M. Quantum dynamical study of the O(D1)+HCl reaction employing three electronic state potential energy surfaces. J Chem Phys 2008; 128:014308. [DOI: 10.1063/1.2813414] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
22
|
An Important Well Studied Atmospheric Reaction, <mml:math altimg="si1.gif" overflow="scroll" xmlns:xocs="http://www.elsevier.com/xml/xocs/dtd" xmlns:xs="http://www.w3.org/2001/XMLSchema" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xmlns="http://www.elsevier.com/xml/ja/dtd" xmlns:ja="http://www.elsevier.com/xml/ja/dtd" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:tb="http://www.elsevier.com/xml/common/table/dtd" xmlns:sb="http://www.elsevier.com/xml/common/struct-bib/dtd" xmlns:ce="http://www.elsevier.com/xml/common/dtd" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:cals="http://www.elsevier.com/xml/common/cals/dtd" xmlns:sa="http://www.elsevier.com/xml/common/struct-aff/dtd"><mml:mi mathvariant="normal">O</mml:mi><mml:mspace width="0.25em"/><mml:mo stretchy="false">(</mml:mo><mml:mmultiscripts><mml:mi>D</mml:mi><mml:mprescripts/><mml:none/><mml:mn>1</mml:mn></mml:mmultiscripts><mml:mo stretchy="false">)</mml:mo><mml:mo>+</mml:mo><mml:msub><mml:mi mathvariant="normal">H</mml:mi><mml:mn>2</mml:mn></mml:msub></mml:math>. ADVANCES IN QUANTUM CHEMISTRY 2008. [DOI: 10.1016/s0065-3276(07)00203-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|
23
|
Hankel M, Smith SC, Meijer AJHM. State-to-state reaction probabilities for the H+O2(v,j)→O+OH(v′,j′) reaction on three potential energy surfaces. J Chem Phys 2007; 127:064316. [PMID: 17705605 DOI: 10.1063/1.2762220] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We report state-to-state and total reaction probabilities for J=0 and total reaction probabilities for J=2 and 4 for the title reaction, both for ground-state and initially rovibrationally excited reactants. The results for three different potential energy surfaces are compared and contrasted. The potential energy surfaces employed are the DMBE IV surface by Pastrana et al. [J. Phys. Chem. 94, 8073 (1990)], the surface by Troe and Ushakov (TU) [J. Chem. Phys. 115, 3621 (2001)], and the new XXZLG ab initio surface by Xu et al. [J. Chem. Phys. 122, 244305 (2005)]. Our results show that the total reaction probabilities from both the TU and XXZLG surfaces are much smaller in magnitude for collision energies above 1.2 eV compared to the DMBE IV surface. The three surfaces also show different behavior with regards to the effect of initial state excitation. The reactivity is increased on the XXZLG and the TU surfaces and decreased on the DMBE IV surface. Vibrational and rotational product state distributions for the XXZLG and the DMBE IV surface show different behaviors for both types of distributions. Our results show that for energies above 1.25 eV the dynamics on the DMBE IV surface are not statistical. However, there is also evidence that the dynamics on the XXZLG surface are not purely statistical for energies above the onset of the first excited product vibrational state v'=1. The magnitude of the total reaction probability is decreased for J>0 for the DMBE IV and the XXZLG surfaces for ground-state reactants. However, for initially rovibrationally excited reactants, the total reaction probability does not decrease as expected for both surfaces. As a result the total cross section averaged over all Boltzmann accessible rotational states may well be larger than the cross section reported in the literature for j=1.
Collapse
Affiliation(s)
- Marlies Hankel
- Centre for Computational Molecular Science, Australian Institute for Bioengineering and Nanotechnology, The University of Queensland, QLD 4072, Australia.
| | | | | |
Collapse
|
24
|
Chu TS, Han KL, Hankel M, Balint-Kurti GG. Coriolis coupling effects in the calculation of state-to-state integral and differential cross sections for the H+D2 reaction. J Chem Phys 2007; 126:214303. [PMID: 17567192 DOI: 10.1063/1.2735624] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The quantum wavepacket parallel computational code DIFFREALWAVE is used to calculate state-to-state integral and differential cross sections for the title reaction on the BKMP2 surface in the total energy range of 0.4-1.2 eV with D2 initially in its ground vibrational-rotational state. The role of Coriolis couplings in the state-to-state quantum calculations is examined in detail. Comparison of the results from calculations including the full Coriolis coupling and those using the centrifugal sudden approximation demonstrates that both the energy dependence and the angular dependence of the calculated cross sections are extremely sensitive to the Coriolis coupling, thus emphasizing the importance of including it correctly in an accurate state-to-state calculation.
Collapse
Affiliation(s)
- Tian-Shu Chu
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China.
| | | | | | | |
Collapse
|
25
|
Aoiz FJ, Bañares L, Herrero VJ. Dynamics of insertion reactions of H2 molecules with excited atoms. J Phys Chem A 2007; 110:12546-65. [PMID: 17107104 DOI: 10.1021/jp063815o] [Citation(s) in RCA: 82] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Recent progress in the study of insertion reactions of hydrogen molecules with excited atoms is reviewed in this article. In particular, the dynamics of the reaction of O(1D), N(2D), C(1D), and S(1D) with H2 and its isotopomers, which have received a great deal of attention over the past decade, are examined in detail. All of these systems have in common the existence of several potential energy surfaces (PES) correlating with the reagents' states, and consequently, they can give rise to reaction following different adiabatic and nonadiabatic pathways. The main contribution, however, arises from their ground singlet PESs which feature the existence of deep wells with small or null barriers for insertion. Accordingly, these reactions proceed mainly via formation of relatively long-lived collision complexes and display an overall nearly statistical behavior. In spite of their similarities, the various reactions have peculiar characteristics caused by important differences of their respective PESs. The contribution of excited PES to the global reactivity, which has also become an important issue and a challenge both for theory and experiment, is also examined. The different theoretical approaches are discussed in the text, along with the experimental results obtained by a variety of techniques. The recent exact quantum treatments of these reactive systems together with the development of a rigorous statistical model have contributed to a very accurate description which in many cases matches very well the detailed measurements. The quasi-classical trajectory (QCT) method has also provided a fairly accurate description of the reaction dynamics for these systems. In particular, the analysis in terms of collision times has yielded interesting clues about the reaction mechanisms.
Collapse
Affiliation(s)
- F J Aoiz
- Departamento de Química Física, Facultad de Química, Universidad Complutense, 28040 Madrid, Spain.
| | | | | |
Collapse
|
26
|
Hankel M, Smith SC, Allan RJ, Gray SK, Balint-Kurti GG. State-to-state reactive differential cross sections for the H+H2→H2+H reaction on five different potential energy surfaces employing a new quantum wavepacket computer code: DIFFREALWAVE. J Chem Phys 2006; 125:164303. [PMID: 17092069 DOI: 10.1063/1.2358350] [Citation(s) in RCA: 73] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
State-to-state differential cross sections have been calculated for the hydrogen exchange reaction, H+H2-->H2+H, using five different high quality potential energy surfaces with the objective of examining the sensitivity of these detailed cross sections to the underlying potential energy surfaces. The calculations were performed using a new parallel computer code, DIFFREALWAVE. The code is based on the real wavepacket approach of Gray and Balint-Kurti [J. Chem. Phys. 108, 950 (1998)]. The calculations are parallelized over the helicity quantum number Omega' (i.e., the quantum number for the body-fixed z component of the total angular momentum) and wavepackets for each J,Omega' set are assigned to different processors, similar in spirit to the Coriolis-coupled processors approach of Goldfield and Gray [Comput. Phys. Commun. 84, 1 (1996)]. Calculations for J=0-24 have been performed to obtain converged state-to-state differential cross sections in the energy range from 0.4 to 1.2 eV. The calculations employ five different potential energy surfaces, the BKMP2 surface and a hierarchical family of four new ab initio surfaces [S. L. Mielke, et al., J. Chem. Phys. 116, 4142 (2002)]. This family of four surfaces has been calculated using three different hierarchical sets of basis functions and also an extrapolation to the complete basis set limit, the so called CCI surface. The CCI surface is the most accurate surface for the H3 system reported to date. Our calculations of differential cross sections are the first to be reported for the A2, A3, A4, and CCI surfaces. They show that there are some small differences in the cross sections obtained from the five different surfaces, particularly at higher energies. The calculations also show that the BKMP2 performs well and gives cross sections in very good agreement with the results from the CCI surface, displaying only small divergences at higher energies.
Collapse
Affiliation(s)
- Marlies Hankel
- Centre for Computational Molecular Science, The University of Queensland, Queensland QLD 4072, Australia.
| | | | | | | | | |
Collapse
|
27
|
Martínez R, Sierra JD, Gray SK, González M. Time dependent quantum dynamics study of the O++H2(v=0,j=0)→OH++H ion-molecule reaction and isotopic variants (D2,HD). J Chem Phys 2006; 125:164305. [PMID: 17092071 DOI: 10.1063/1.2359727] [Citation(s) in RCA: 35] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The time dependent real wave packet method using the helicity decoupling approximation was used to calculate the cross section evolution with collision energy (excitation function) of the O++H2(v=0,j=0)-->OH++H reaction and its isotopic variants with D2 and HD, using the best available ab initio analytical potential energy surface. The comparison of the calculated excitation functions with exact quantum results and experimental data showed that the present quantum dynamics approach is a very useful tool for the study of the selected and related systems, in a quite wide collision energy interval (approximately 0.0-1.1 eV), involving a much lower computational cost than the quantum exact methods and without a significant loss of accuracy in the cross sections.
Collapse
Affiliation(s)
- Rodrigo Martínez
- Departmento de Química, Universidad de La Rioja, C/Madre de Dios 51, 26006 Logroño, Spain
| | | | | | | |
Collapse
|
28
|
Abstract
This paper is an overview of the theory of reactive scattering, with emphasis on fully quantum mechanical theories that have been developed to describe simple chemical reactions, especially atom-diatom reactions. We also describe related quasiclassical trajectory applications, and in all of this review the emphasis is on methods and applications concerned with state-resolved reaction dynamics. The review first provides an overview of the development of the theory, including a discussion of computational methods based on coupled channel calculations, variational methods, and wave packet methods. Choices of coordinates, including the use of hyperspherical coordinates are discussed, as are basis set and discrete variational representations. The review also summarizes a number of applications that have been performed, especially the two most comprehensively studied systems, H+H2 and F+H2, along with brief discussions of a large number of other systems, including other hydrogen atom transfer reactions, insertion reactions, electronically nonadiabatic reactions, and reactions involving four or more atoms. For each reaction we describe the method used and important new physical insight extracted from the results.
Collapse
Affiliation(s)
- Wenfang Hu
- Department of Chemistry, Northwestern University, Evanston, Illinois 60208-3113, USA
| | | |
Collapse
|
29
|
Balucani N, Capozza G, Leonori F, Segoloni E, Casavecchia P. Crossed molecular beam reactive scattering: from simple triatomic to multichannel polyatomic reactions. INT REV PHYS CHEM 2006. [DOI: 10.1080/01442350600641305] [Citation(s) in RCA: 68] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
|
30
|
Bañares L, Aoiz FJ, González-Lezana T, Herrero VJ, Tanarro I. Influence of rotation and isotope effects on the dynamics of the N(D2)+H2 reactive system and of its deuterated variants. J Chem Phys 2005; 123:224301. [PMID: 16375470 DOI: 10.1063/1.2131075] [Citation(s) in RCA: 45] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Integral cross sections and thermal rate constants have been calculated for the N((2)D)+H(2) reaction and its isotopic variants N((2)D)+D(2) and the two-channel N((2)D)+HD by means of quasiclassical trajectory and statistical quantum-mechanical model methods on the latest ab initio potential-energy surface [T.-S. Ho et al., J. Chem. Phys. 119, 3063 (2003)]. The effect of rotational excitation of the diatom on the dynamics of these reactions has been investigated and interesting discrepancies between the classical and statistical model calculations have been found. Whereas a net effect of reagent rotation on reactivity is always observed in the classical calculations, only a very slight effect is observed in the case of the asymmetric N((2)D)+HD reaction for the statistical quantum-mechanical method. The thermal rate constants calculated on this Potential-Energy Surface using quasiclassical trajectory and statistical model methods are in good agreement with the experimental determinations, although the latter are somewhat larger. A reevaluation of the collinear barrier of the potential surface used in the present study seems timely. Further theoretical and experimental studies are needed for a full understanding of the dynamics of the title reaction.
Collapse
Affiliation(s)
- L Bañares
- Departamento de Química Física, Facultad de Química, Universidad Complutense, 28040 Madrid, Spain.
| | | | | | | | | |
Collapse
|
31
|
Larrégaray P, Bonnet L, Rayez JC. Validity of Phase Space Theory for Atom−Diatom Insertion Reactions. J Phys Chem A 2005; 110:1552-60. [PMID: 16435816 DOI: 10.1021/jp053822x] [Citation(s) in RCA: 45] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Phase space theory (PST) is applied to the calculation of state-resolved integral and differential cross sections for the complex-forming atom-diatom insertion reactions A + H(2) --> AH(2) --> AH + H with A = C((1)D), S((1)D), O((1)D), and N((2)D). In the asymptotic channels, vibration motion is quantized while rotation and translation motions are treated classically. The approach is compared to exact quantum scattering calculations and quantum statistical models. Given the simplicity of PST, the agreement with the previous much more refined treatments is very satisfying. Although PST is a well-established theory, this work is, to our knowledge, the first such systematic comparison of its predictions with accurate quantum scattering and quantum statistical calculations.
Collapse
Affiliation(s)
- P Larrégaray
- Laboratoire de Physico-Chimie Moléculaire, UMR5803, Université Bordeaux1-CNRS, 33405 Talence Cedex, France.
| | | | | |
Collapse
|
32
|
Ying Lin S, Guo H. Quantum statistical and wave packet studies of insertion reactions of S(D1) with H2, HD, and D2. J Chem Phys 2005; 122:074304. [PMID: 15743229 DOI: 10.1063/1.1851500] [Citation(s) in RCA: 66] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A thorough theoretical investigation of the reactions between S(1D) and various hydrogen isotopomers (H2, D2, and HD) has been carried out using a recent ab initio potential energy surface. State-resolved integral and differential cross sections, thermal rate constants, and their dependence on energy or temperature were obtained from quantum mechanical capture probabilities within a statistical model. For comparison, the J=0 reaction probabilities were also computed using an exact wave packet method. The statistical results are in excellent agreement with available exact differential and integral cross sections. The comparison with experimental results shows that the agreement is reasonably good in general, but some significant differences exist, particularly for the SD/SH branching ratio in the S(1D)+HD reaction.
Collapse
Affiliation(s)
- Shi Ying Lin
- Department of Chemistry, University of New Mexico, Albuquerque, New Mexico 87131, USA
| | | |
Collapse
|
33
|
Khachatrian A, Dagdigian PJ. Determination of the internal state distribution of the SD product from the S(1D)+D2 reaction. J Chem Phys 2005; 122:024303. [PMID: 15638582 DOI: 10.1063/1.1827598] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
The S(1D)+D2-->SD+D reaction has been studied through a photolysis-probe experiment in a cell. S(1D) reagent was prepared by 193 nm photolysis of CS2, and the SD(X 2Pi) product was detected by laser fluorescence excitation. The nascent rotational/fine-structure state distribution of the SD(X 2Pi) product was determined. This reaction, previously studied theoretically and in a crossed molecular beam experiment, is known to proceed through formation and decay of a long-lived collision complex involving the deep well in the H2S ground electronic state. The determined SD rotational state distribution in the v=0 vibrational level was found to be approximately statistical, with a small preference for formation of the F1 (Omega=3/2) fine-structure manifold over F2 (Omega=1/2). The branching into the Lambda doublet levels was also investigated, and essentially equal populations of levels of A' and A" symmetry were found. The present results are compared with previous investigations of this reaction and the analogous O(1D)+D2 reaction.
Collapse
Affiliation(s)
- Ani Khachatrian
- Department of Chemistry, The Johns Hopkins University, Baltimore, Maryland 21218-2685, USA
| | | |
Collapse
|
34
|
Lin SY, Guo H. Case Study of a Prototypical Elementary Insertion Reaction: C(1D) + H2 → CH + H. J Phys Chem A 2004. [DOI: 10.1021/jp046039y] [Citation(s) in RCA: 50] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Shi Ying Lin
- Department of Chemistry, University of New Mexico, Albuquerque, New Mexico 87131
| | - Hua Guo
- Department of Chemistry, University of New Mexico, Albuquerque, New Mexico 87131
| |
Collapse
|
35
|
Lin SY, Guo H. Reactions of C(1D) with H2 and its deuterated isotopomers, a wave packet study. J Chem Phys 2004; 121:1285-92. [PMID: 15260670 DOI: 10.1063/1.1764502] [Citation(s) in RCA: 34] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Using a Chebyshev wave packet method, total and state-resolved reaction probabilities (J=0) were calculated for the reactions of C(1D) with various hydrogen isotopomers (H2, D2, and HD, nu i=0, j i=0) on a recent ab initio potential energy surface. For all the isotopic variants, it was found that the initial state specified reaction probabilities have no energy threshold and are strongly oscillatory, indicative of the involvement of long-lived resonances in this barrierless reaction. The J=0 product vibrational and rotational distributions for all three isotopic reactions, and the CH/CD branching ratio for the C+HD reaction, show strong dependence on the collision energy, further underscoring the important role played by the resonances. The generally decaying vibrational distributions and highly excited rotational distributions, which corroborate an insertion mechanism, and the dominance of the CD+H channel in the C+HD reaction are consistent with existing experimental observations. Initial state specified integral cross sections and rate constants were estimated using a capture model. The estimated rate constants were found to be close and in the order kHD>kH2>kD2. Finally, a method to calculate branching ratio in the C+HD reaction is proposed.
Collapse
Affiliation(s)
- Shi Ying Lin
- Department of Chemistry, University of New Mexico, Albuquerque, New Mexico 87131, USA
| | | |
Collapse
|
36
|
Lin SY, Guo H. A wave packet based statistical approach to complex-forming reactions. J Chem Phys 2004; 120:9907-10. [PMID: 15268008 DOI: 10.1063/1.1756584] [Citation(s) in RCA: 65] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A wave packet based statistical model is suggested for complex-forming reactions. This model assumes statistical formation and decay of the long-lived reaction complex and computes reaction cross sections and their energy dependence from capture probabilities. This model is very efficient and reasonably accurate for reactions dominated by long-lived resonances, as confirmed by its application to the C((1)D)+H(2) reaction.
Collapse
Affiliation(s)
- Shi Ying Lin
- Department of Chemistry, University of New Mexico, Albuquerque, NM 87131, USA
| | | |
Collapse
|
37
|
Lin SY, Guo H. Quantum Wave Packet Studies of the C(1D) + H2 → CH + H Reaction: Integral Cross Section and Rate Constant. J Phys Chem A 2004. [DOI: 10.1021/jp031184h] [Citation(s) in RCA: 90] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Shi Ying Lin
- Department of Chemistry, University of New Mexico, Albuquerque, New Mexico 87131
| | - Hua Guo
- Department of Chemistry, University of New Mexico, Albuquerque, New Mexico 87131
| |
Collapse
|
38
|
Lin SY, Guo H. Quantum wave packet study of reactive and inelastic scattering between C(1D) and H2. J Chem Phys 2003. [DOI: 10.1063/1.1624060] [Citation(s) in RCA: 89] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
39
|
Abstract
This review discusses recent quantum scattering calculations on bimolecular chemical reactions in the gas phase. This theory provides detailed and accurate predictions on the dynamics and kinetics of reactions containing three atoms. In addition, the method can now be applied to reactions involving polyatomic molecules. Results obtained with both time-independent and time-dependent quantum dynamical methods are described. The review emphasises the recent development in time-dependent wave packet theories and the applications of reduced dimensionality approaches for treating polyatomic reactions. Calculations on over 40 different reactions are described.
Collapse
|
40
|
Aoiz FJ, Bañares L, Castillo JF, Herrero VJ, Martı́nez-Haya B, Honvault P, Launay JM, Liu X, Lin JJ, Harich SA, Wang CC, Yang X. The O(1D)+H2 reaction at 56 meV collision energy: A comparison between quantum mechanical, quasiclassical trajectory, and crossed beam results. J Chem Phys 2002. [DOI: 10.1063/1.1478693] [Citation(s) in RCA: 68] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
41
|
van Harrevelt R, van Hemert MC, Schatz GC. The CH+H reaction studied with quantum-mechanical and classical trajectory calculations. J Chem Phys 2002. [DOI: 10.1063/1.1459416] [Citation(s) in RCA: 31] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
|
42
|
Takayanagi T. Nonadiabatic quantum reactive scattering calculations for the O(1D)+H2, D2, and HD reactions on the lowest three potential energy surfaces. J Chem Phys 2002. [DOI: 10.1063/1.1434988] [Citation(s) in RCA: 29] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|