1
|
Mazo-Sevillano PD, Aguado A, Goicoechea JR, Roncero O. Quantum study of the CH3+ photodissociation in full-dimensional neural network potential energy surfaces. J Chem Phys 2024; 160:184307. [PMID: 38738612 DOI: 10.1063/5.0206895] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/05/2024] [Accepted: 04/22/2024] [Indexed: 05/14/2024] Open
Abstract
C H 3 + , a cornerstone intermediate in interstellar chemistry, has recently been detected for the first time by using the James Webb Space Telescope. The photodissociation of this ion is studied here. Accurate explicitly correlated multi-reference configuration interaction ab initio calculations are done, and full-dimensional potential energy surfaces are developed for the three lower electronic states, with a fundamental invariant neural network method. The photodissociation cross section is calculated using a full-dimensional quantum wave packet method in heliocentric Radau coordinates. The wave packet is represented in angular and radial grids, allowing us to reduce the number of points physically accessible, requiring to push up the spurious states appearing when evaluating the angular kinetic terms, through projection technique. The photodissociation spectra, when employed in astrochemical models to simulate the conditions of the Orion bar, result in a lesser destruction of CH3+ compared to that obtained when utilizing the recommended values in the kinetic database for astrochemistry.
Collapse
Affiliation(s)
- Pablo Del Mazo-Sevillano
- Unidad Asociada UAM-IFF-CSIC, Departamento de Química Física Aplicada, Facultad de Ciencias M-14, Universidad Autónoma de Madrid, 28049 Madrid, Spain
| | - Alfredo Aguado
- Unidad Asociada UAM-IFF-CSIC, Departamento de Química Física Aplicada, Facultad de Ciencias M-14, Universidad Autónoma de Madrid, 28049 Madrid, Spain
| | - Javier R Goicoechea
- Instituto de Física Fundamental (IFF-CSIC), C.S.I.C., Serrano 123, 28006 Madrid, Spain
| | - Octavio Roncero
- Instituto de Física Fundamental (IFF-CSIC), C.S.I.C., Serrano 123, 28006 Madrid, Spain
| |
Collapse
|
2
|
Murakami T, Ibuki S, Hashimoto Y, Kikuma Y, Takayanagi T. Dynamics study of the post-transition-state-bifurcation process of the (HCOOH)H + → CO + H 3O +/HCO + + H 2O dissociation: application of machine-learning techniques. Phys Chem Chem Phys 2023; 25:14016-14027. [PMID: 37161528 DOI: 10.1039/d3cp00252g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/11/2023]
Abstract
The process of protonated formic acid dissociating from the transition state was studied using ring-polymer molecular dynamics (RPMD), classical MD, and quasi-classical trajectory (QCT) simulations. Temperature had a strong influence on the branching fractions for the HCO+ + H2O and CO + H3O+ dissociation channels. The RPMD and classical MD simulations showed similar behavior, but the QCT dynamics were significantly different owing to the excess energies in the quasi-classical trajectories. Machine-learning analysis identified several key features in the phase information of the vibrational motions at the transition state. We found that the initial configuration and momentum of a hydrogen atom connected to a carbon atom and the shrinking coordinate of the CO bond at the transition state play a role in the dynamics of HCO+ + H2O production.
Collapse
Affiliation(s)
- Tatsuhiro Murakami
- Department of Chemistry, Saitama University, Shimo-Okubo 255, Sakura-ku, Saitama City, Saitama, 338-8570, Japan.
- Department of Materials & Life Sciences, Faculty of Science & Technology, Sophia University, 7-1 Kioicho, Chiyoda-ku, Tokyo, 102-8554, Japan
| | - Shunichi Ibuki
- Department of Chemistry, Saitama University, Shimo-Okubo 255, Sakura-ku, Saitama City, Saitama, 338-8570, Japan.
| | - Yu Hashimoto
- Department of Chemistry, Saitama University, Shimo-Okubo 255, Sakura-ku, Saitama City, Saitama, 338-8570, Japan.
| | - Yuya Kikuma
- Department of Chemistry, Saitama University, Shimo-Okubo 255, Sakura-ku, Saitama City, Saitama, 338-8570, Japan.
| | - Toshiyuki Takayanagi
- Department of Chemistry, Saitama University, Shimo-Okubo 255, Sakura-ku, Saitama City, Saitama, 338-8570, Japan.
| |
Collapse
|
3
|
del Mazo-Sevillano P, Félix-González D, Aguado A, Sanz-Sanz C, Kwon DH, Roncero O. Vibrational, non-adiabatic and isotopic effects in the dynamics of the H 2 + H 2+ → H 3+ + H reaction: application to plasma modelling. Mol Phys 2023. [DOI: 10.1080/00268976.2023.2183071] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/04/2023]
Affiliation(s)
- P. del Mazo-Sevillano
- Department of Mathematics and Computer Science, FU Berlin, Berlin, Germany
- Unidad Asociada UAM-IFF-CSIC, Departamento de Química Física Aplicada, Facultad de Ciencias M-14, Universidad Autónoma de Madrid, Madrid, Spain
| | - D. Félix-González
- Unidad Asociada UAM-IFF-CSIC, Departamento de Química Física Aplicada, Facultad de Ciencias M-14, Universidad Autónoma de Madrid, Madrid, Spain
| | - A. Aguado
- Unidad Asociada UAM-IFF-CSIC, Departamento de Química Física Aplicada, Facultad de Ciencias M-14, Universidad Autónoma de Madrid, Madrid, Spain
| | - C. Sanz-Sanz
- Unidad Asociada UAM-IFF-CSIC, Departamento de Química Física Aplicada, Facultad de Ciencias M-14, Universidad Autónoma de Madrid, Madrid, Spain
| | - D.-H. Kwon
- Nuclear Physics Application Research Division, Korea Atomic Energy Research Institute, Daejeon, Republic of Korea
| | - O. Roncero
- Instituto de Física Fundamental, IFF-CSIC, Madrid, Spain
| |
Collapse
|
4
|
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
|
5
|
Murakami T, Iida R, Hashimoto Y, Takahashi Y, Takahashi S, Takayanagi T. Ring-Polymer Molecular Dynamics and Kinetics for the H – + C 2H 2 → H 2 + C 2H – Reaction Using the Full-Dimensional Potential Energy Surface. J Phys Chem A 2022; 126:9244-9258. [DOI: 10.1021/acs.jpca.2c05851] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Affiliation(s)
- Tatsuhiro Murakami
- Department of Chemistry, Saitama University, Shimo-Okubo 255, Sakura-ku, Saitama City, Saitama338-8570, Japan
- Department of Materials & Life Sciences, Faculty of Science & Technology, Sophia University, 7-1 Kioicho, Chiyoda-ku, Tokyo102-8554, Japan
| | - Ryusei Iida
- Department of Chemistry, Saitama University, Shimo-Okubo 255, Sakura-ku, Saitama City, Saitama338-8570, Japan
| | - Yu Hashimoto
- Department of Chemistry, Saitama University, Shimo-Okubo 255, Sakura-ku, Saitama City, Saitama338-8570, Japan
| | - Yukinobu Takahashi
- Department of Chemistry, Saitama University, Shimo-Okubo 255, Sakura-ku, Saitama City, Saitama338-8570, Japan
| | - Soma Takahashi
- Department of Chemistry, Saitama University, Shimo-Okubo 255, Sakura-ku, Saitama City, Saitama338-8570, Japan
| | - Toshiyuki Takayanagi
- Department of Chemistry, Saitama University, Shimo-Okubo 255, Sakura-ku, Saitama City, Saitama338-8570, Japan
| |
Collapse
|
6
|
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
|
7
|
Xiang H, Tian L, Li Y, Song H. Energy- and Local-Gradient-Based Neural Network Method for Accurately Describing Long-Range Interaction: Application to the H 2 + CO + Reaction. J Phys Chem A 2022; 126:352-363. [PMID: 34989591 DOI: 10.1021/acs.jpca.1c09719] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The long-range interaction plays an important role in theoretically describing ion-molecule reaction. However, most energy-based neural network fitting methods usually introduce spurious long-range interactions. In this work, we propose an energy- and local-gradient-based neural network (ELGNN) method to fit potential energy surfaces (PESs). K-means clustering is employed to divide the whole configuration space into three regions: reactant asymptotic region, interaction region, and product asymptotic region. In the interaction region, only the energies of sampled points are computed, while in the asymptotic regions, the gradients of partially sampled configurations are calculated as well, and both the energies and energy gradients (if necessary) are used to fit long-range interactions. These regions are joined together by switching functions. The ELGNN method is first applied to fit the PES of the H2 + CO+ reaction, which has significant long-range interactions. It is found that the ELGNN method works better than the energy-based NN method in describing long-range interactions. The dynamics and kinetics of the reaction are then investigated on the new PES.
Collapse
Affiliation(s)
- Haipan Xiang
- College of Physical Science and Technology, Huazhong Normal University, Wuhan 430079, China.,State Key Laboratory of Magnetic Resonance and Atomic and Molecular Physics, Innovation Academy for Precision Measurement Science and Technology, Chinese Academy of Sciences, Wuhan 430071, China
| | - Li Tian
- College of Physical Science and Technology, Huazhong Normal University, Wuhan 430079, China.,State Key Laboratory of Magnetic Resonance and Atomic and Molecular Physics, Innovation Academy for Precision Measurement Science and Technology, Chinese Academy of Sciences, Wuhan 430071, China
| | - Yong Li
- College of Physical Science and Technology, Huazhong Normal University, Wuhan 430079, China
| | - Hongwei Song
- State Key Laboratory of Magnetic Resonance and Atomic and Molecular Physics, Innovation Academy for Precision Measurement Science and Technology, Chinese Academy of Sciences, Wuhan 430071, China
| |
Collapse
|
8
|
Zhang RM, Chen W, Truhlar DG, Xu X. Master Equation Study of Hydrogen Abstraction from HCHO by OH Via a Chemically Activated Intermediate. Faraday Discuss 2022; 238:431-460. [DOI: 10.1039/d2fd00024e] [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
The abstraction reaction of hydrogen from formaldehyde by OH radical plays an important role in formaldehyde oxidation. The reaction involves a bimolecular association to form a chemically activated hydrogen-bonded reaction...
Collapse
|
9
|
Xu X, Chen J, Lu X, Fang W, Liu S, Zhang DH. Strong non-Arrhenius behavior at low temperatures in the OH + HCl → H 2O + Cl reaction due to resonance induced quantum tunneling. Chem Sci 2022; 13:7955-7961. [PMID: 35865883 PMCID: PMC9258319 DOI: 10.1039/d2sc01958b] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/05/2022] [Accepted: 06/11/2022] [Indexed: 11/21/2022] Open
Abstract
The OH + HCl reaction possesses many Feshbach resonances trapped in the hydrogen bond well in the entrance channel, which substantially enhance the reaction rates at low temperatures.
Collapse
Affiliation(s)
- Xin Xu
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, Liaoning, China
| | - Jun Chen
- State Key Laboratory of Structure Chemistry, Fujian Institute of Research on the Structure of Matter, Chinese Academy of Sciences, Fuzhou 350002, Fujian, China
| | - Xiaoxiao Lu
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, Liaoning, China
| | - Wei Fang
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, Liaoning, China
| | - Shu Liu
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, Liaoning, China
| | - Dong H. Zhang
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, Liaoning, China
| |
Collapse
|
10
|
Saito K, Hashimoto Y, Takayanagi T. Ring-Polymer Molecular Dynamics Calculations of Thermal Rate Coefficients and Branching Ratios for the Interstellar H 3+ + CO → H 2 + HCO +/HOC + Reaction and Its Deuterated Analogue. J Phys Chem A 2021; 125:10750-10756. [PMID: 34918514 DOI: 10.1021/acs.jpca.1c09160] [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
The reaction between H3+ and CO is important in understanding the H3+ destruction mechanism in the interstellar medium. In this work, thermal rate coefficients for the H3+ + CO and D3+ + CO reactions are calculated using ring-polymer molecular dynamics (RPMD) on a high-level machine-learning potential energy surface. The RPMD results agree well with the classical molecular dynamics results, where nuclear quantum effects are completely ignored, whereas the agreement between the RPMD results and the previous quasi-classical trajectory is good only at low temperatures. The calculated [HCO+]/[HOC+] product branching ratios decrease as the temperature increases, and the product branching is exclusively determined by the initial collisional orientation, which governs the formation of an ion-dipole complex, H3+···CO or H3+···OC, that dissociates into H2 + HCO+ or H2 + HOC+, respectively, via a direct mechanism. However, the contribution of the indirect mechanism via the rearrangement between H3+···CO and H3+···OC increases as the temperature increases, although its absolute fraction is small.
Collapse
Affiliation(s)
- Kohei Saito
- Department of Chemistry, Saitama University, Shimo-Okubo 255, Sakura-ku, Saitama 338-8570, Japan
| | - Yu Hashimoto
- Department of Chemistry, Saitama University, Shimo-Okubo 255, Sakura-ku, Saitama 338-8570, Japan
| | - Toshiyuki Takayanagi
- Department of Chemistry, Saitama University, Shimo-Okubo 255, Sakura-ku, Saitama 338-8570, Japan
| |
Collapse
|
11
|
Pan M, Xiang H, Li Y, Song H. Study on the kinetics and dynamics of the H 2 + NH 2- reaction on a high-level ab initio potential energy surface. Phys Chem Chem Phys 2021; 23:17848-17855. [PMID: 34612274 DOI: 10.1039/d1cp02423j] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Gas-phase ion-molecule reactions play major roles in many fields of chemistry and physics. The reaction of an amino radical anion with a hydrogen molecule is one of the simplest proton transfer reactions involving anions. A globally accurate full-dimensional potential energy surface (PES) for the NH2- + H2 reaction is developed by the fundamental invariant-neural network method, resulting in a root mean square error of 0.116 kcal mol-1. Quasi-classical trajectory calculations are then carried out on the newly developed PES to give integral cross sections, differential cross sections and thermal rate coefficients. This reaction has two reaction channels, proton transfer and hydrogen exchange. The reactivity of the proton transfer channel is about one or two orders of magnitude stronger than that of the hydrogen exchange channel in the energy range studied. Vibrational excitation of H2 promotes the proton transfer reaction, while fundamental excitation of each vibrational mode of NH2- has a negligible effect. In addition, the theoretical rate coefficients of the proton transfer reaction on the PES show inverse temperature dependence from 150 to 750 K, in accordance with the available experimental results.
Collapse
Affiliation(s)
- Mengyi Pan
- State Key Laboratory of Magnetic Resonance and Atomic and Molecular Physics, Wuhan Institute of Physics and Mathematics, Innovation Academy for Precision Measurement Science and Technology, Chinese Academy of Sciences, Wuhan 430071, China.
| | | | | | | |
Collapse
|
12
|
Guo H, Worth G, Domcke W. Quantum dynamics with ab initio potentials. J Chem Phys 2021; 155:080401. [PMID: 34470339 DOI: 10.1063/5.0066234] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, USA
| | - Graham Worth
- Department of Chemistry, University College London, London WC1H 0AJ, United Kingdom
| | - Wolfgang Domcke
- Department of Chemistry, Technical University of Munich, D-85747 Garching, Germany
| |
Collapse
|
13
|
Li C, Liu Q, Zhang L, Li Y, Jiang B. Ring polymer molecular dynamics in gas-surface reactions: tests on initial sampling and potential energy landscape. Mol Phys 2021. [DOI: 10.1080/00268976.2021.1941367] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
Affiliation(s)
- Chen Li
- Hefei National Laboratory for Physical Science at the Microscale, Department of Chemical Physics, Key Laboratory of Surface and Interface Chemistry and Energy Catalysis of Anhui Higher Education Institutes, University of Science and Technology of China, Hefei, People’s Republic of China
| | - Qinghua Liu
- Hefei National Laboratory for Physical Science at the Microscale, Department of Chemical Physics, Key Laboratory of Surface and Interface Chemistry and Energy Catalysis of Anhui Higher Education Institutes, University of Science and Technology of China, Hefei, People’s Republic of China
| | - Liang Zhang
- Hefei National Laboratory for Physical Science at the Microscale, Department of Chemical Physics, Key Laboratory of Surface and Interface Chemistry and Energy Catalysis of Anhui Higher Education Institutes, University of Science and Technology of China, Hefei, People’s Republic of China
| | - Yongle Li
- Department of Physics, International Center of Quantum and Molecular Structures and Shanghai Key Laboratory of High Temperature Superconductors, Shanghai University, Shanghai, People’s Republic of China
| | - Bin Jiang
- Hefei National Laboratory for Physical Science at the Microscale, Department of Chemical Physics, Key Laboratory of Surface and Interface Chemistry and Energy Catalysis of Anhui Higher Education Institutes, University of Science and Technology of China, Hefei, People’s Republic of China
| |
Collapse
|