1
|
Braunstein M, Bonnet L. An efficient algorithm for capturing quantum effects in classical reactive scattering: application to D + H+3 → H 2D + + H. Phys Chem Chem Phys 2023; 25:1602-1605. [PMID: 36541279 DOI: 10.1039/d2cp05108g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Motivated by a recent semiclassical analysis of chemical reaction thresholds [Bonnet et al., J. Chem. Phys., 2022, 157, 094114], we present an efficient algorithm for including zero-point energy (ZPE) effects in classical reactive scattering. The algorithm is an extension of the quasi-classical trajectory (QCT) Gaussian binning method. We apply it to the astrophysically important D + H+3 reaction, where there are significant quantum effects and where application of other methods is problematic [Braunstein et al., Phys. Chem. Chem. Phys., 2022, 24, 5489]. The rate constants computed with the new, general algorithm closely match recent Ring Polymer Molecular Dynamics (RPMD) [Bulut et al., J. Phys. Chem. A, 2019, 123, 8766] and experimentally derived [Bowen et al., J. Chem. Phys., 2021, 154, 084307] ones spanning ∼4 orders of magnitude from 70 to 1500 K.
Collapse
Affiliation(s)
- Matthew Braunstein
- Spectral Sciences Incorporated, 4 Fourth Avenue, Burlington, MA 01824, USA.
| | - Laurent Bonnet
- CNRS, Université de Bordeaux, ISM, UMR 5255, F-33400 Talence, France.
| |
Collapse
|
2
|
Braunstein M, Bonnet L, Roncero O. Capturing quantum effects with quasi-classical trajectories in the D + H+3 → H 2D + + H reaction. Phys Chem Chem Phys 2022; 24:5489-5505. [PMID: 35171152 DOI: 10.1039/d1cp04244k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We present quasi-classical trajectory (QCT) cross sections, rate constants, and product state distributions for the D + H+3 → H2D+ + H reaction. Using the same H+4 potential surface, the rate constants obtained from several QCT-based methods correcting for zero-point effects by Gaussian binning the product H2D+ are compared to ring polymer molecular dynamics (RPMD) rate constants [Bulut et al., J. Phys. Chem. A, 2019, 123, 8766] which include quantum effects and to recent experimentally derived rate constants [Bowen et al., J. Chem. Phys., 2021, 154, 084307]. QCT with standard binning predicts rate constants that increase slowly as the temperature decreases from 1500 to 100 K. In contrast, the RPMD rate constants decrease rapidly with decreasing temperature. By 100 K, the QCT standard binning rate constant is more than 3 orders of magnitude larger than the RPMD rate constant. We show that QCT with Gaussian binning and proper normalization captures the zero-point effects and reproduces the RPMD rate constants over a large temperature range. Furthermore, the simple technique of counting only reactive trajectories with vibrational energy above the product zero-point energy matches the RPMD results well down to ∼300 K. The present Gaussian binned rate constants are in fair agreement with new experimentally derived rate constants from 100 to 1500 K. However, because the Gaussian binned rate constants do not include tunneling, important at lower temperatures, and the RPMD and experimentally derived rate constants have significant differences, the roles of the competing effects of zero-point energy, internal excitation of the H+3, and quantum tunneling are not simple and require further study for a consistent picture of the dynamics. Since rate constants for complex forming reactions, such as the title reaction, are difficult to converge with RPMD, alternative QCT-based methods, which include quantum effects and in addition provide product state distributions as described here, are highly desirable.
Collapse
Affiliation(s)
- Matthew Braunstein
- Spectral Sciences Incorporated, 4 Fourth Avenue, Burlington, MA 01824, USA.
| | - Laurent Bonnet
- CNRS, Université de Bordeaux, ISM, UMR 5255, F-33400 Talence, France
| | - Octavio Roncero
- Instituto de Fisica Fundamental (IFF-CSIC), C.S.I.C., Serrano 123, 28006 Madrid, Spain
| |
Collapse
|
3
|
Liu A, Li J. Dynamical investigations of the O( 3P) + H 2O reaction at high collision energies on an accurate full-dimensional potential energy surface. Mol Phys 2021. [DOI: 10.1080/00268976.2021.1944686] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
Affiliation(s)
- Aike Liu
- Chemical Synthesis and Pollution Control Key Laboratory of Sichuan Province, Institute of Applied Chemistry, College of Chemistry and Chemical Engineering, China West Normal University, Nanchong, People’s Republic of China
| | - Jun Li
- School of Chemistry and Chemical Engineering & Chongqing Key Laboratory of Theoretical and Computational Chemistry, Chongqing University, Chongqing, People’s Republic of China
| |
Collapse
|
4
|
Li J, Li J. A Full-Dimensional Potential Energy Surface and Dynamics of the Multichannel Reaction between H and HO 2. J Phys Chem A 2021; 125:1540-1552. [PMID: 33591185 DOI: 10.1021/acs.jpca.0c11213] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
In addition to its vital significance in combustion and atmospheric chemistry, the reaction between H' and HO2 on the ground triplet state represents a prototype with multiple product channels, including H2 + O2, OH + OH, O + H2O, and H + H'O2. In this work, a full-dimensional accurate potential energy surface (PES) for the title reaction was developed to provide reliable descriptions for all dynamically relevant regions. Using this PES, we adopted the quasi-classical trajectory approach to study the corresponding reaction dynamics, including the reactivity of each product channel and the associated product branching ratio, the product energy distributions, product angular distributions, and associated microscopic mechanisms. For representing distributions of the product energies, such as product translational energy as well as product rotational and vibrational energies, both the traditional histogram and the kernel density estimation (KDE) methods were used and compared. It seems that the features of the resulting distributions in this work are very similar to each other among different methods. The KDE method is suggested for statistics, particularly for those populations with small oscillations in the histogram plot.
Collapse
Affiliation(s)
- Jia Li
- School of Chemistry and Chemical Engineering & Chongqing Key Laboratory of Theoretical and Computational Chemistry, Chongqing University, Chongqing 401331, China
| | - Jun Li
- School of Chemistry and Chemical Engineering & Chongqing Key Laboratory of Theoretical and Computational Chemistry, Chongqing University, Chongqing 401331, China
| |
Collapse
|
5
|
Masoumpour MS, Daryanavard M. Ab Initio-Based Global Potential Energy Surface and Reaction Dynamics for H 2S + C. J Phys Chem A 2020; 124:7901-7910. [DOI: 10.1021/acs.jpca.0c05856] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
| | - Marzieh Daryanavard
- Department of chemistry, Estahban Higher Education Center, Estahban 74519-44655, Iran
| |
Collapse
|
6
|
da Silva RS, Ballester MY. A theoretical study of energy transfer in Ar( 1S) + SO 2( X ̃ 1 A ') collisions: Cross sections and rate coefficients for vibrational transitions. J Chem Phys 2018; 149:144309. [PMID: 30316261 DOI: 10.1063/1.5051349] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Vibrational transitions, induced by collisions between rare-gas atoms and molecules, play a key role in many problems of interest in physics and chemistry. A theoretical investigation of the translation-to-vibration (T-V) energy transfer process in argon atom and sulfur dioxide molecule collisions is presented here. For such a purpose, the framework of the quasi-classical trajectory (QCT) methodology was followed over the range of translational energies 2 ≤ Etr/kcal mol-1 ≤ 100. A new realistic potential energy surface (PES) for the ArSO2 system was developed using pairwise addition for the four-body energy term within the double many-body expansion. The topological features of the obtained function are compared with a previous one reported by Hippler et al. [J. Phys. Chem. 90, 6158 (1986)]. To test the accuracy of the PES, additional coupled cluster singles and doubles method with a perturbative contribution of connected triples calculations were carried out for the global minimum configuration. From dynamical calculations, the cross sections for the T-V excitation process indicate a barrier-type mechanism due to strong repulsive interactions between SO2 molecules and the Ar atom. Corrections to zero-point energy leakage in QCT were carried out using vibrational energy quantum mechanical threshold of the complex and variations. Rate coefficients and cross sections are calculated for some vibrational transitions using pseudo-quantization approaches of the vibrational energy of products. Main attributes of the title molecular collision are discussed and compared with available information in the literature.
Collapse
Affiliation(s)
- Ramon S da Silva
- Departamento de Física, Universidade Federal de Juiz de Fora, Juiz de Fora, MG 36036-330, Brazil
| | - Maikel Y Ballester
- Departamento de Física, Universidade Federal de Juiz de Fora, Juiz de Fora, MG 36036-330, Brazil
| |
Collapse
|
7
|
Zheng R, Zhu Y, Song H. Mode specific dynamics in bond selective reaction O′(3P) + HOD → O′H + OD/O′D + OH. J Chem Phys 2018; 149:054304. [DOI: 10.1063/1.5037492] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Rui Zheng
- School of Mathematics and Statistics, North China University of Water Resources and Electric Power, Zhengzhou 450011, China
- State Key Laboratory of Magnetic Resonance and Atomic and Molecular Physics, Wuhan Institute of Physics and Mathematics, Chinese Academy of Sciences, Wuhan 430071, China
| | - Yongfa Zhu
- State Key Laboratory of Magnetic Resonance and Atomic and Molecular Physics, Wuhan Institute of Physics and Mathematics, Chinese Academy of Sciences, Wuhan 430071, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Hongwei Song
- State Key Laboratory of Magnetic Resonance and Atomic and Molecular Physics, Wuhan Institute of Physics and Mathematics, Chinese Academy of Sciences, Wuhan 430071, China
| |
Collapse
|
8
|
Masoumpour MS, Mousavipour SH. A Theoretical Study on the Dynamics of O(3P) + H2S(1A1) Reaction on an Interpolated Potential Energy Surface. BULLETIN OF THE CHEMICAL SOCIETY OF JAPAN 2017. [DOI: 10.1246/bcsj.20170107] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
|
9
|
Mazarei E, Mousavipour SH. A Theoretical Study on the Dynamics of the Reaction of CH Radicals with Water. J Phys Chem A 2017; 121:8033-8047. [DOI: 10.1021/acs.jpca.7b05504] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Elham Mazarei
- Department of Chemistry,
College of Science, Shiraz University, Shiraz, Iran
| | | |
Collapse
|
10
|
Li J, Dawes R, Guo H. An accurate multi-channel multi-reference full-dimensional global potential energy surface for the lowest triplet state of H2O2. Phys Chem Chem Phys 2016; 18:29825-29835. [DOI: 10.1039/c6cp06232f] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
A full dimensional potential energy surface for the lowest triplet state of H2O2 was developed at the MRCI-F12 level.
Collapse
Affiliation(s)
- Jun Li
- School of Chemistry and Chemical Engineering
- Chongqing University
- Chongqing 401331
- China
- Department of Chemistry and Chemical Biology
| | - Richard Dawes
- Department of Chemistry
- Missouri University of Science and Technology
- Rolla
- USA
| | - Hua Guo
- Department of Chemistry and Chemical Biology
- University of New Mexico
- Albuquerque
- USA
| |
Collapse
|
11
|
Li J, Jiang B, Song H, Ma J, Zhao B, Dawes R, Guo H. From ab Initio Potential Energy Surfaces to State-Resolved Reactivities: X + H2O ↔ HX + OH [X = F, Cl, and O(3P)] Reactions. J Phys Chem A 2015; 119:4667-87. [DOI: 10.1021/acs.jpca.5b02510] [Citation(s) in RCA: 85] [Impact Index Per Article: 9.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Jun Li
- Department
of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
- School of Chemistry
and Chemical Engineering, Chongqing University, Chongqing 400044, China
| | - Bin Jiang
- Department
of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| | - Hongwei Song
- Department
of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| | - Jianyi Ma
- Institute of Atomic
and Molecular Physics, Sichuan University, Chengdu, Sichuan 610065, China
| | - Bin Zhao
- Department
of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| | - Richard Dawes
- Department
of Chemistry, Missouri University of Science and Technology, Rolla, Missouri 65409, United States
| | - Hua Guo
- Department
of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| |
Collapse
|
12
|
Braunstein M, Conforti PF. Classical Dynamics of H2O Vibrational Self-Relaxation. J Phys Chem A 2015; 119:3311-22. [DOI: 10.1021/jp5115527] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Matthew Braunstein
- Spectral Sciences Incorporated, 4 Fourth Avenue, Burlington, Massachusetts 01803, United States
| | - Patrick F. Conforti
- Spectral Sciences Incorporated, 4 Fourth Avenue, Burlington, Massachusetts 01803, United States
| |
Collapse
|
13
|
Hama T, Watanabe N. Surface Processes on Interstellar Amorphous Solid Water: Adsorption, Diffusion, Tunneling Reactions, and Nuclear-Spin Conversion. Chem Rev 2013; 113:8783-839. [DOI: 10.1021/cr4000978] [Citation(s) in RCA: 211] [Impact Index Per Article: 19.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Tetsuya Hama
- Institute of Low Temperature
Science, Hokkaido University, N19W8 Kita-ku, Sapporo, Hokkaido 060-0819, Japan
| | - Naoki Watanabe
- Institute of Low Temperature
Science, Hokkaido University, N19W8 Kita-ku, Sapporo, Hokkaido 060-0819, Japan
| |
Collapse
|
14
|
Jiang B, Guo H. Control of Mode/Bond Selectivity and Product Energy Disposal by the Transition State: X + H2O (X = H, F, O(3P), and Cl) Reactions. J Am Chem Soc 2013; 135:15251-6. [DOI: 10.1021/ja408422y] [Citation(s) in RCA: 115] [Impact Index Per Article: 10.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Bin Jiang
- Department of Chemistry and
Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| | - Hua Guo
- Department of Chemistry and
Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| |
Collapse
|
15
|
Li J, Guo H. A new ab initio based global HOOH(13A″) potential energy surface for the O(3P) + H2O(X1A1) ↔ OH(X2Π) + OH(X2Π) reaction. J Chem Phys 2013; 138:194304. [DOI: 10.1063/1.4804418] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
|