1
|
Czakó G, Gruber B, Papp D, Tajti V, Tasi DA, Yin C. First-principles mode-specific reaction dynamics. Phys Chem Chem Phys 2024; 26:15818-15830. [PMID: 38639072 DOI: 10.1039/d4cp00417e] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/20/2024]
Abstract
Controlling the outcome of chemical reactions by exciting specific vibrational and/or rotational modes of the reactants is one of the major goals of modern reaction dynamics studies. In the present Perspective, we focus on first-principles vibrational and rotational mode-specific dynamics computations on reactions of neutral and anionic systems beyond six atoms such as X + C2H6 [X = F, Cl, OH], HX + C2H5 [X = Br, I], OH- + CH3I, and F- + CH3CH2Cl. The dynamics simulations utilize high-level ab initio analytical potential energy surfaces and the quasi-classical trajectory method. Besides initial state specificity and the validity of the Polanyi rules, mode-specific vibrational-state assignment for polyatomic product species using normal-mode analysis and Gaussian binning is also discussed and compared with experiment.
Collapse
Affiliation(s)
- Gábor Czakó
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group, Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged, Rerrich Béla tér 1, Szeged H-6720, Hungary.
| | - Balázs Gruber
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group, Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged, Rerrich Béla tér 1, Szeged H-6720, Hungary.
| | - Dóra Papp
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group, Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged, Rerrich Béla tér 1, Szeged H-6720, Hungary.
| | - Viktor Tajti
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group, Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged, Rerrich Béla tér 1, Szeged H-6720, Hungary.
| | - Domonkos A Tasi
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group, Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged, Rerrich Béla tér 1, Szeged H-6720, Hungary.
| | - Cangtao Yin
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group, Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged, Rerrich Béla tér 1, Szeged H-6720, Hungary.
| |
Collapse
|
2
|
Espinosa-Garcia J, Rangel C, Corchado JC. Role of the Vibrational and Translational Energies in the CN(v)+C 2H 6(ν 1, ν 2, ν 5 and ν 9) Reactions. A Theoretical QCT Study. Chemphyschem 2024; 25:e202300997. [PMID: 38421195 DOI: 10.1002/cphc.202300997] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/27/2023] [Revised: 02/27/2024] [Accepted: 02/28/2024] [Indexed: 03/02/2024]
Abstract
Quasi-classical trajectory (QCT) calculations were conducted on the newly developed full-dimensional potential energy surface, PES-2023, to analyse two critical aspects: the influence of vibrational versus translational energy in promoting reactivity, and the impact of vibrational excitation within similar vibrational modes. The former relates to Polanyi's rules, while the latter concerns mode selectivity. Initially, the investigation revealed that independent vibrational excitation by a single quantum of ethane's symmetric and asymmetric stretching modes (differing by only 15 cm-1) yielded comparable dynamics, reaction cross-sections, HCN(v) vibrational product distributions, and scattering distributions. This observation dismisses any significant mode selectivity. Moreover, an equivalent amount of energy provided as translational energy (at total energies of 9.6 and 20.0 kcal mol-1) gave rise to slightly lower reactivity compared to the same amount of energy provided as vibrational energy. This effect is more evident at low energies, presenting a counterintuitive scenario in an 'early transition state' reaction. These findings challenge the straightforward application of Polanyi's rules in polyatomic systems. Regarding CN(v) vibrational excitation, our calculations reveal that the reaction cross-section remains practically unaffected by this vibrational excitation, suggesting that the CN stretching mode is a spectator mode. The results were rationalized by considering several factors: the strong coupling between different vibrational modes, and between vibrational modes and the reaction coordinate; and a significant vibrational energy redistribution within the ethane reactant before collision. This redistribution creates an unphysical energy flow, resulting in loss of adiabaticity and vibrational memory before the reactants' collision. These theoretical findings require future confirmation through experimental or theoretical quantum mechanical studies, which are currently unavailable.
Collapse
Affiliation(s)
- Joaquin Espinosa-Garcia
- Área de Quimica Fisica and Instituto de Computación Científica Avanzada de Extremadura, Universidad de Extremadura, Badajoz, Spain
| | - Cipriano Rangel
- Área de Quimica Fisica and Instituto de Computación Científica Avanzada de Extremadura, Universidad de Extremadura, Badajoz, Spain
| | - Jose C Corchado
- Área de Quimica Fisica and Instituto de Computación Científica Avanzada de Extremadura, Universidad de Extremadura, Badajoz, Spain
| |
Collapse
|
3
|
Wu Y, Hu Y, Li Z, Ma J. Molecular Dynamics Simulation of Thermal Nonequilibrium and Chemical Reaction Processes in Hydrogen Combustion. J Phys Chem A 2024; 128:2643-2655. [PMID: 38530707 DOI: 10.1021/acs.jpca.3c08131] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/28/2024]
Abstract
Using reactive force field (ReaxFF) and molecular dynamics simulation, we investigate the combustion process of hydrogen-oxygen systems in initial thermal nonequilibrium states with different translational and rovibrational temperatures for oxygen. The system studied in this work contains 300 oxygen molecules and 700 hydrogen molecules with a density of 7 times the air density. For this system, the characteristic relaxation times of oxygen and hydrogen vibrational energies are 0.173 and 0.249 ns, respectively. 0.6% of hydrogen undergoes a chemical reaction with oxygen during the thermal nonequilibrium relaxation stage. For the distribution of translational energy and vibrational energy of oxygen in the thermal nonequilibrium state, the maximum mean error of the statistical distribution in the simulation and the Boltzmann distribution at temperature calculated from the average kinetic energy of molecules is about 2.25 × 10-5. At the same time, it was observed in the simulation that many-body interactions play a certain role in the combustion process. Furthermore, we compare the ignition time and temperature rise behavior of different combustion mechanisms and molecular dynamics simulations starting from the thermal equilibrium state. These results will provide meaningful references for the construction of thermal nonequilibrium combustion chemical reaction mechanisms.
Collapse
Affiliation(s)
- Yimiao Wu
- Institute of Atomic and Molecular Physics, Sichuan University, Chengdu 610065, P. R. China
| | - Yongxin Hu
- Institute of Atomic and Molecular Physics, Sichuan University, Chengdu 610065, P. R. China
| | - Zhiwei Li
- Institute of Atomic and Molecular Physics, Sichuan University, Chengdu 610065, P. R. China
| | - Jianyi Ma
- Institute of Atomic and Molecular Physics, Sichuan University, Chengdu 610065, P. R. China
| |
Collapse
|
4
|
Huber ME, Lewis TWR, Meta M, Ard SG, Liu Y, Sweeny BC, Guo H, Ončák M, Shuman NS, Meyer J. Ta + and Nb + + CO 2: intersystem crossing in ion-molecule reactions. Phys Chem Chem Phys 2024; 26:8670-8680. [PMID: 38437035 DOI: 10.1039/d3cp05549c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/06/2024]
Abstract
The reactions of Ta+ and Nb+ with CO2 proceed only by a highly efficient oxygen atom transfer reaction to the respective oxide at room temperature in the gas phase. Although the product spin states are not determined, thermochemistry dictates that they must be different from ground state quintet Ta+ and Nb+, implying that intersystem crossing (ISC) has occurred. Recent reactive scattering experiments found dominant indirect dynamics for the reaction with Ta+ hinting at a bottleneck along the reaction path. The question on the nature of the bottleneck, whether it involves a crossing point or a transition state, could not be finally answered because theory located both close to each other. Here, we aim at shedding further light onto the impact of intersystem crossing on the reaction dynamics and ultimately the reactivity of transition metal ion reactions in the gas phase. We employ a combination of thermal kinetics for Ta+ and Nb+ with CO2 using a selected-ion flow tube (SIFT) apparatus and differential scattering cross sections for Nb+ + CO2 from crossed-beam velocity map imaging. The reaction with niobium again shows dominant indirect dynamics and in general very similar dynamics compared to Ta+ + CO2. At thermal energies, both reactions show sub-collisional rate constants with small negative temperature dependencies. Experiments are complemented by high level quantum chemical calculations of the minimum energy pathway. Statistical modelling well-reproduces the experimental thermal rate constants, and suggests that the Nb+ reaction is rate-limited by the intersystem crossing at thermal energies.
Collapse
Affiliation(s)
- Maximilian E Huber
- RPTU Kaiserslautern-Landau, Fachbereich Chemie und Forschungszentrum OPTIMAS, Erwin-Schrödinger Str. 52, 67663 Kaiserslautern, Germany.
| | - Tucker W R Lewis
- Air Force Research Laboratory, Space Vehicles Directorate, Kirtland Air Force Base, New Mexico 87117, USA.
| | - Marcel Meta
- RPTU Kaiserslautern-Landau, Fachbereich Chemie und Forschungszentrum OPTIMAS, Erwin-Schrödinger Str. 52, 67663 Kaiserslautern, Germany.
| | - Shaun G Ard
- Air Force Research Laboratory, Space Vehicles Directorate, Kirtland Air Force Base, New Mexico 87117, USA.
| | - Yang Liu
- Department of Chemistry and Chemical Biology, Center for Computational Chemistry, University of New Mexico, Albuquerque, New Mexico 87131, USA
| | - Brendan C Sweeny
- Air Force Research Laboratory, Space Vehicles Directorate, Kirtland Air Force Base, New Mexico 87117, USA.
| | - Hua Guo
- Department of Chemistry and Chemical Biology, Center for Computational Chemistry, University of New Mexico, Albuquerque, New Mexico 87131, USA
| | - Milan Ončák
- Universität Innsbruck, Institut für Ionenenphysik und Angewandte Physik, Technikerstra. 25, 6020 Innsbruck, Austria
| | - Nicholas S Shuman
- Air Force Research Laboratory, Space Vehicles Directorate, Kirtland Air Force Base, New Mexico 87117, USA.
| | - Jennifer Meyer
- RPTU Kaiserslautern-Landau, Fachbereich Chemie und Forschungszentrum OPTIMAS, Erwin-Schrödinger Str. 52, 67663 Kaiserslautern, Germany.
| |
Collapse
|
5
|
Rawat AMS, Sahoo J, Mahapatra S. Combined Quantum Mechanical and Quasi-Classical State-to-State Dynamical Study on the Isotopic Effect in H/D + LiH +/LiD + → H 2/HD/D 2 + Li + Reactions. J Phys Chem A 2023; 127:10733-10746. [PMID: 38096485 DOI: 10.1021/acs.jpca.3c04131] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2023]
Abstract
Coriolis-coupled quantum mechanical (QM-CC) and quasi-classical trajectory (QCT) calculations are carried out to investigate the dynamics of the H(D) + LiH+(v = 0, j = 0) → H2(HD) (v', j') + Li+ reactions on the ground electronic state potential energy surface reported by Martinazzo et al. (Martinazzo et al., J. Chem. Phys. 2003, 119, 11241). The QM-CC and QCT results at the initial state-selected and state-to-state levels are used to investigate the validity and accuracy of the QCT method for these exoergic barrierless reactions. Furthermore, the QCT method is used to understand the isotopic effects on reaction observables like total and state-to-state integral cross section, differential cross section, product energy disposal, and rate constants of H(D) + LiH+(v = 0, j = 0) → H2(HD) (v', j') + Li+ and H(D) + LiD+(v = 0, j = 0) → HD(D2) (v', j') + Li+ reactions. Attempts are also made to understand the impact of the isotopic substitution on the reaction mechanism. It is observed that QM-CC and QCT results closely follow each other at the initial state-selected and state-to-state levels. Noticeable kinematic effects of reagents on the reactivity and mechanism of the reactions are also observed.
Collapse
Affiliation(s)
| | - Jayakrushna Sahoo
- School of Chemistry, University of Hyderabad, Hyderabad 500046, India
| | - S Mahapatra
- School of Chemistry, University of Hyderabad, Hyderabad 500046, India
| |
Collapse
|
6
|
Pan H, Zhao B, Guo H, Liu K. State-to-State Dynamics in Mode-Selective Polyatomic Reactions. J Phys Chem Lett 2023; 14:10412-10419. [PMID: 37955874 DOI: 10.1021/acs.jpclett.3c02853] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2023]
Abstract
Chemical reactions are intrinsically quantum mechanical transformations of reactants to products. Recent experimental and theoretical advances have enabled the exploration of reaction dynamics with a quantum state resolution for both reactants and products. To this end, reactions involving more than three atoms are of particular interest, because they exhibit rich dynamics concerning the role of different reactant modes in controlling reactivity and product energy disposal. A clear understanding of the state-to-state dynamics requires new paradigms. In this Perspective, we examine some new concepts that have emerged from recent state-to-state studies of polyatomic reactions and illustrate the key role played by the transition state.
Collapse
Affiliation(s)
- Huilin Pan
- Department of Chemistry, Southern University of Science and Technology, Shenzhen 518055, P. R. China
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, CAS, Dalian 116023, P. R. China
| | - Bin Zhao
- Department of Chemistry, Southern University of Science and Technology, Shenzhen 518055, P. R. China
| | - Hua Guo
- Department of Chemistry and Chemical Biology, Center for Computational Chemistry, University of New Mexico, Albuquerque, New Mexico 87131, United States
| | - Kopin Liu
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, CAS, Dalian 116023, P. R. China
- Aerosol Science Research Center, National Sun Yat-sen University, Kaohsiung 80424, Taiwan
- Institute of Atomic and Molecular Sciences (IAMS), Academia Sinica, Taipei 10617, Taiwan
| |
Collapse
|
7
|
Yin C, Czakó G. Full-dimensional automated potential energy surface development and detailed dynamics for the CH 2OO + NH 3 reaction. Phys Chem Chem Phys 2023; 25:26917-26922. [PMID: 37787004 DOI: 10.1039/d3cp03469k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/04/2023]
Abstract
With the help of the ROBOSURFER program package, a global full-dimensional potential energy surface (PES) for the reaction of the Criegee intermediate, CH2OO, with the NH3 molecule is developed iteratively using different ab initio methods and the monomial symmetrization fitting approach. The final permutationally-invariant analytical PES is constructed based on 23447 geometries and the corresponding ManyHF-based CCSD(T)-F12b/cc-pVTZ-F12 energies. The accuracy of the PES is confirmed by the excellent agreement of its stationary-point properties and one-dimensional potential energy curves compared with the corresponding ab initio data. The reaction probabilities and integral cross sections are calculated for the ground-state and several vibrationally excited-state reactions by quasi-classical trajectory simulations. Remarkable is that the maximum impact parameter b where reactivity vanishes is almost independent of collision energy ranging from 1 to 40 kcal mol-1, and the reaction probability increases with increasing collision energy for this negative-barrier reaction. At the same time, a slight mode-specificity effect is observed. In addition, the deuterium effect is investigated and the sudden vector projection is discussed.
Collapse
Affiliation(s)
- Cangtao Yin
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group, Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged, Rerrich Béla tér 1, Szeged H-6720, Hungary.
| | - Gábor Czakó
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group, Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged, Rerrich Béla tér 1, Szeged H-6720, Hungary.
| |
Collapse
|
8
|
Song H, Guo H. Theoretical Insights into the Dynamics of Gas-Phase Bimolecular Reactions with Submerged Barriers. ACS PHYSICAL CHEMISTRY AU 2023; 3:406-418. [PMID: 37780541 PMCID: PMC10540288 DOI: 10.1021/acsphyschemau.3c00009] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 03/08/2023] [Revised: 06/02/2023] [Accepted: 06/06/2023] [Indexed: 10/03/2023]
Abstract
Much attention has been paid to the dynamics of both activated gas-phase bimolecular reactions, which feature monotonically increasing integral cross sections and Arrhenius kinetics, and their barrierless capture counterparts, which manifest monotonically decreasing integral cross sections and negative temperature dependence of the rate coefficients. In this Perspective, we focus on the dynamics of gas-phase bimolecular reactions with submerged barriers, which often involve radicals or ions and are prevalent in combustion, atmospheric chemistry, astrochemistry, and plasma chemistry. The temperature dependence of the rate coefficients for such reactions is often non-Arrhenius and complex, and the corresponding dynamics may also be quite different from those with significant barriers or those completely dominated by capture. Recent experimental and theoretical studies of such reactions, particularly at relatively low temperatures or collision energies, have revealed interesting dynamical behaviors, which are discussed here. The new knowledge enriches our understanding of the dynamics of these unusual reactions.
Collapse
Affiliation(s)
- 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
| | - Hua Guo
- Department
of Chemistry and Chemical Biology, University
of New Mexico, Albuquerque, New Mexico 87131, United States
| |
Collapse
|
9
|
Gruber B, Tajti V, Czakó G. Vibrational Mode-Specific Dynamics of the OH + C 2H 6 Reaction. J Phys Chem A 2023; 127:7364-7372. [PMID: 37620310 PMCID: PMC10493966 DOI: 10.1021/acs.jpca.3c04328] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/27/2023] [Revised: 08/04/2023] [Indexed: 08/26/2023]
Abstract
We investigate the effects of the initial vibrational excitations on the dynamics of the OH + C2H6 → H2O + C2H5 reaction using the quasi-classical trajectory method and a full-dimensional analytical ab initio potential energy surface. Excitation of the initial CH, CC, and OH stretching modes enhances, slightly inhibits, and does not affect the reactivity, respectively. Translational energy activates the early-barrier title reaction more efficiently than OH and CC stretching excitations, in accord with the Polanyi rules whereas CH stretching modes have similar or higher efficacy than translation, showing that these rules are not always valid in polyatomic processes. Scattering angle, initial attack angle, and product translational energy distributions show the dominance of direct stripping with increasing collision energy, side-on OH and isotropic C2H6 attack preferences, and substantial reactant-product translational energy transfer without any significant mode specificity. The reactant vibrational excitation energy of OH and C2H6 flows into the H2O and C2H5 product vibrations, respectively, whereas product rotations are not affected. The computed mode-specific H2O vibrational distributions show that initial OH excitation appears in the asymmetric stretching vibration of the H2O product and allow comparison with experiments.
Collapse
Affiliation(s)
- Balázs Gruber
- MTA-SZTE Lendület Computational
Reaction Dynamics Research Group, Interdisciplinary Excellence Centre
and Department of Physical Chemistry and Materials Science, Institute
of Chemistry, University of Szeged, Rerrich Béla tér 1, Szeged H-6720, Hungary
| | - Viktor Tajti
- MTA-SZTE Lendület Computational
Reaction Dynamics Research Group, Interdisciplinary Excellence Centre
and Department of Physical Chemistry and Materials Science, Institute
of Chemistry, University of Szeged, Rerrich Béla tér 1, Szeged H-6720, Hungary
| | - Gábor Czakó
- MTA-SZTE Lendület Computational
Reaction Dynamics Research Group, Interdisciplinary Excellence Centre
and Department of Physical Chemistry and Materials Science, Institute
of Chemistry, University of Szeged, Rerrich Béla tér 1, Szeged H-6720, Hungary
| |
Collapse
|
10
|
Oda K, Tsutsumi T, Keshavamurthy S, Furuya K, Armentrout PB, Taketsugu T. Dynamically Hidden Reaction Paths in the Reaction of CF 3+ + CO. ACS PHYSICAL CHEMISTRY AU 2022; 2:388-398. [PMID: 36193292 PMCID: PMC9524575 DOI: 10.1021/acsphyschemau.2c00012] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 12/03/2022]
Abstract
![]()
Reaction paths on
a potential energy surface are widely used in
quantum chemical studies of chemical reactions. The recently developed
global reaction route mapping (GRRM) strategy automatically constructs
a reaction route map, which provides a complete picture of the reaction.
Here, we thoroughly investigate the correspondence between the reaction
route map and the actual chemical reaction dynamics for the CF3+ + CO reaction studied by guided ion beam tandem
mass spectrometry (GIBMS). In our experiments, FCO+, CF2+, and CF+ product ions were observed,
whereas if the collision partner is N2, only CF2+ is observed. Interestingly, for reaction with CO, GRRM-predicted
reaction paths leading to the CF+ + F2CO product
channel are found at a barrier height of about 2.5 eV, whereas the
experimentally obtained threshold for CF+ formation was
7.48 ± 0.15 eV. In other words, the ion was not obviously observed
in the GIBMS experiment, unless a much higher collision energy than
the requisite energy threshold was provided. On-the-fly molecular
dynamics simulations revealed a mechanism that hides these reaction
paths, in which a non-statistical energy distribution at the first
collisionally reached transition state prevents the reaction from
proceeding along some reaction paths. Our results highlight the existence
of dynamically hidden reaction paths that may be inaccessible in experiments
at specific energies and hence the importance of reaction dynamics
in controlling the destinations of chemical reactions.
Collapse
Affiliation(s)
- Kohei Oda
- Graduate School of Chemical Sciences and Engineering, Hokkaido University, Sapporo 060-0810, Japan
| | - Takuro Tsutsumi
- Department of Chemistry, Faculty of Science, Hokkaido University, Sapporo 060-0810, Japan
| | - Srihari Keshavamurthy
- Department of Chemistry, Faculty of Science, Hokkaido University, Sapporo 060-0810, Japan
- Department of Chemistry, Indian Institute of Technology, Kanpur 208 016, India
| | - Kenji Furuya
- Faculty of Arts and Science, Kyushu University, Motooka, Fukuoka 819-0395, Japan
- Department of Molecular and Material Sciences, Kyushu University, Kasuga, Fukuoka 816-8580, Japan
- Department of Chemistry, University of Utah, Salt Lake City 84112, United States
| | - P. B. Armentrout
- Department of Chemistry, University of Utah, Salt Lake City 84112, United States
| | - Tetsuya Taketsugu
- Department of Chemistry, Faculty of Science, Hokkaido University, Sapporo 060-0810, Japan
- Institute for Chemical Reaction Design and Discovery (WPI-ICReDD), Hokkaido University, Sapporo 001-0021, Japan
| |
Collapse
|
11
|
Goswami S, San Vicente Veliz JC, Upadhyay M, Bemish RJ, Meuwly M. Quantum and quasi-classical dynamics of the C( 3P) + O 2( 3Σ-g) → CO( 1Σ +) + O( 1D) reaction on its electronic ground state. Phys Chem Chem Phys 2022; 24:23309-23322. [PMID: 36165004 PMCID: PMC9533374 DOI: 10.1039/d2cp02840a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The dynamics of the C(3P) + O2(3Σ−g) → CO(1Σ+) + O(1D) reaction on its electronic ground state is investigated by using time-dependent wave packet propagation (TDWP) and quasi-classical trajectory (QCT) simulations. For the moderate collision energies considered (Ec = 0.001 to 0.4 eV, corresponding to a range from 10 K to 4600 K) the total reaction probabilities from the two different treatments of the nuclear dynamics agree very favourably. The undulations present in P(E) from the quantum mechanical treatment can be related to stabilization of the intermediate CO2 complex with lifetimes on the 0.05 ps time scale. This is also confirmed from direct analysis of the TDWP simulations and QCT trajectories. Product diatom vibrational and rotational level resolved state-to-state reaction probabilities from TDWP and QCT simulations agree well except for the highest product vibrational states (v′ ≥ 15) and for the lowest product rotational states (j′ ≤ 10). Opening of the product vibrational level CO(v′ = 17) requires ∼0.2 eV from QCT and TDWP simulations with O2(j = 0) and decreases to 0.04 eV if all initial rotational states are included in the QCT analysis, compared with Ec > 0.04 eV obtained from experiments. It is thus concluded that QCT simulations are suitable for investigating and realistically describe the C(3P) + O2(3Σ−g) → CO(1Σ+) + O(1D) reaction down to low collision energies when compared with results from a quantum mechanical treatment using TDWPs. The dynamics of the C(3P) + O2(3Σ−g) → CO(1Σ+) + O(1D) reaction on its electronic ground state is investigated by using time-dependent wave packet propagation (TDWP) and quasi-classical trajectory (QCT) simulations.![]()
Collapse
Affiliation(s)
- Sugata Goswami
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland.
| | | | - Meenu Upadhyay
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland.
| | - Raymond J Bemish
- Air Force Research Laboratory, Space Vehicles Directorate, Kirtland AFB, New Mexico 87117, USA
| | - Markus Meuwly
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland. .,Department of Chemistry, Brown University, RI, USA
| |
Collapse
|
12
|
Han S, Zhao B, Conte R, Malbon CL, Bowman JM, Yarkony DR, Guo H. Nonadiabatic Reactive Quenching of OH( A2Σ +) by H 2: Origin of High Vibrational Excitation in the H 2O Product. J Phys Chem A 2022; 126:6944-6952. [PMID: 36137233 DOI: 10.1021/acs.jpca.2c05704] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The nonadiabatic dynamics of the reactive quenching channel of the OH(A2Σ+) + H2/D2 collisions is investigated with a semiclassical surface hopping method, using a recently developed four-state diabatic potential energy matrix (DPEM). In agreement with experimental observations, the H2O/HOD products are found to have significant vibrational excitation. Using a Gaussian binning method, the H2O vibrational state distribution is determined. The preferential energy disposal into the product vibrational modes is rationalized by an extended Sudden Vector Projection model, in which the h and g vectors associated with the conical intersection are found to have large projections with the product normal modes. However, our calculations did not find significant insertion trajectories, suggesting the need for further improvement of the DPEM.
Collapse
Affiliation(s)
- Shanyu Han
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| | - Bin Zhao
- Department of Chemistry, Southern University of Science and Technology, Shenzhen, 518055, China
| | - Riccardo Conte
- Dipartimento di Chimica, Università degli Studi di Milano, via Golgi 19, 20133 Milano, Italy
| | - Christopher L Malbon
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Joel M Bowman
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| | - David R Yarkony
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| |
Collapse
|
13
|
Ma D, Ma J. Full-dimensional quantum mechanical calculations for the tunneling behavior of HOCO dissociation to H + CO 2. Phys Chem Chem Phys 2022; 24:15321-15329. [PMID: 35703166 DOI: 10.1039/d1cp04269f] [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 tunneling behavior during HOCO dissociation to H + CO2 was investigated by full-dimensional quantum mechanical calculations based on an accurate global potential energy surface. The tunneling lifetimes for the low-lying 1500 vibrational states were calculated using the low-storage filter diagonalization method after a 1 million-step Chebyshev propagation. In the calculated energy range, the lifetimes of different vibrational states with similar energy are found to differ by 3-4 orders of magnitude, and the lower limit for these tunneling lifetimes is consistent with the experimental results reported by Continetti et al. For the given vibrational progressions, the lifetime of the vibrational states decreases with the increasing energy level, which is consistent with the results of 1D simulation calculations reported by Bowman, but the declining curve obviously fluctuates, and the declining slope is significantly different from that obtained by 1D simulation. Due to a difference in the effective barrier width, the mode-specific behavior of the tunneling effect is manifested in that the O-C-O' and H-O-C bend modes lead to the largest enhancement and an inhibitory effect on the tunneling process, respectively.
Collapse
Affiliation(s)
- Dandan Ma
- Institute of Atomic and Molecular Physics, Sichuan University, Chengdu, Sichuan, 610065, China.
| | - Jianyi Ma
- Institute of Atomic and Molecular Physics, Sichuan University, Chengdu, Sichuan, 610065, China.
| |
Collapse
|
14
|
Yin R, Jiang B, Guo H. Mechanism and Dynamics of CO 2 Formation in Formic Acid Decomposition on Pt Surfaces. ACS Catal 2022. [DOI: 10.1021/acscatal.2c01598] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Rongrong Yin
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| | - Bin Jiang
- Department of Chemical Physics, School of Chemistry and Materials Science, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| |
Collapse
|
15
|
Gao D, Xin X, Wang D, Szabó P, Lendvay G. Theoretical dynamics studies of the CH 3 + HBr → CH 4 + Br reaction: integral cross sections, rate constants and microscopic mechanism. Phys Chem Chem Phys 2022; 24:10548-10560. [PMID: 35445671 DOI: 10.1039/d2cp00066k] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023]
Abstract
Quantum and quasi-classical dynamics calculations have been performed for the reaction of HBr with CH3. The accurate ab initio-based potential energy surface function developed earlier for this reaction displays a potential well corresponding to a reactant complex and a submerged potential barrier. The integral cross sections were calculated on this potential energy surface using both a six-degree-of-freedom reduced dimensional quantum dynamics and the quasi-classical trajectory method and very good agreement was found between the two approaches. The cross sections were found to diverge when the collision energy decreases, indicating that the reactant attraction is responsible for the dynamics at low collision energy. The quantum mechanical and the quasi-classical rate constants also agree very well and almost exactly reproduce the experimental results at low temperatures up to 540 K. The negative activation energy observed experimentally is confirmed by the calculations and is a consequence of the long-range attraction between the reactants. From the classical trajectories mechanistic details have been extracted. It is found that at very low collision energy, the reacting system crosses the potential barrier because the forces within the complex guide them, although some 30% is reflected from the product side of the barrier. When the collision energy increases, the system does not follow the most favorable path and the reactants are, with increasing probability, reflected from the repulsive walls of the nonreactive parts of the reactants, providing a picture beyond the decreasing excitation function.
Collapse
Affiliation(s)
- Delu Gao
- College of Physics and Electronics, Shandong Normal University, Jinan 250014, Shandong, China.
| | - Xin Xin
- College of Physics and Electronics, Shandong Normal University, Jinan 250014, Shandong, China.
| | - Dunyou Wang
- College of Physics and Electronics, Shandong Normal University, Jinan 250014, Shandong, China.
| | - Péter Szabó
- Faculté des Sciences, des Technologies et de Médecine, Département Physique et sciences des matériaux, Campus Limpertsberg, Université du Luxembourg 162 A, avenue de la, Faïencerie L-1511, Luxembourg
| | - György Lendvay
- Institute of Materials and Environmental Chemistry, Research Centre for Natural Sciences, Magyar Tudósok krt. 2, H-1117 Budapest, Hungary. .,Department of General and Inorganic Chemistry, University of Pannonia, Egyetem utca 10, Veszprém, H-8200, Hungary
| |
Collapse
|
16
|
Ellerbrock R, Zhao B, Manthe U. Vibrational control of the reaction pathway in the H + CHD 3 → H 2 + CD 3 reaction. SCIENCE ADVANCES 2022; 8:eabm9820. [PMID: 35353570 PMCID: PMC8967217 DOI: 10.1126/sciadv.abm9820] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 10/26/2021] [Accepted: 02/08/2022] [Indexed: 06/14/2023]
Abstract
An accurate full-dimensional quantum state-to-state simulation of the six-atom title reaction based on first-principles theory is reported. Counterintuitive effects are found: Increasing the energy in the reactant's CD3 umbrella vibration reduces the energy in the corresponding product vibration. An in-depth analysis reveals the crucial role of the effective dynamical transition state: Its geometry is controlled by the vibrational states of the reactants and subsequently controls the quantum state distribution of the products. This finding enables generalizing the concept of transition state control of chemical reactions to the quantum state-specific level.
Collapse
Affiliation(s)
- Roman Ellerbrock
- Theoretische Chemie, Fakultät für Chemie, Universität Bielefeld, Universitätsstr. 25, D-33615 Bielefeld, Germany
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, CA 94305, USA
- SLAC National Accelerator Laboratory, Menlo Park, CA 94025, USA
| | - Bin Zhao
- Theoretische Chemie, Fakultät für Chemie, Universität Bielefeld, Universitätsstr. 25, D-33615 Bielefeld, Germany
| | - Uwe Manthe
- Theoretische Chemie, Fakultät für Chemie, Universität Bielefeld, Universitätsstr. 25, D-33615 Bielefeld, Germany
| |
Collapse
|
17
|
Chen Q, Zhang S, Hu X, Xie D, Guo H. Reaction Pathway Control via Reactant Vibrational Excitation and Impact on Product Vibrational Distributions: The O + HO 2 → OH + O 2 Atmospheric Reaction. J Phys Chem Lett 2022; 13:1872-1878. [PMID: 35175051 DOI: 10.1021/acs.jpclett.2c00053] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
Chemical reactions often have multiple pathways, the control of which is of fundamental and practical importance. In this Letter, we examine the dynamics of the O + HO2 → OH + O2 reaction, which plays an important role in atmospheric chemistry, using quasi-classical trajectories on a recently developed full-dimensional potential energy surface (PES). This reaction has two pathways leading to the same products: the H abstraction pathway (Oa + HObOc → OaH + ObOc) and the O abstraction pathway (Oa + HObOc → ObH + OaOc). Under thermal conditions, the reaction is dominated by the latter channel, which is barrierless, leading to vibrational excitation of the O2 product. However, we demonstrate that excitation of the HO2 reactant in its O-H (v1) vibrational mode results in dramatic switching of the reaction pathway to the activated H abstraction channel, which leads to a highly excited OH product vibrational state distribution. The implications of such dynamical effects in the atmospheric chemistry are discussed.
Collapse
Affiliation(s)
- Qixin Chen
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Shuwen Zhang
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Xixi Hu
- Kuang Yaming Honors School, Institute for Brain Sciences, Nanjing University, Nanjing 210023, China
| | - Daiqian Xie
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| |
Collapse
|
18
|
Song H, Xie W, Zhang C, Yang M. Toward a Comprehensive Understanding of Mode-Specific Dynamics of Polyatomic Reactions: A Full-Dimensional Quantum Dynamics Study of the H + NH 3 Reaction. J Phys Chem A 2022; 126:663-669. [PMID: 35080397 DOI: 10.1021/acs.jpca.1c08399] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Mode specificity not only sheds light on reaction dynamics but also opens the door for chemical reaction control. This work reports a state-of-the-art full-dimensional quantum dynamics study on the prototypical hydrogen abstraction reaction of hydrogen with ammonia, which serves as a benchmark for advancing our fundamental understanding of polyatomic reaction dynamics. By taking advantage of the (3 + 1) Radau-Jacobi coordinates, the bond-specific probabilities are resolved with the reactant NH3 initiated from either a non-degenerate or degenerate stretching vibrational state. The observed different atom-specific abstraction probabilities from individual states of the degenerate pair are rationalized in the local mode representation according to the different vibrational energy deposited in each N-H bond. It is verified that the three H atoms in NH3 have the same abstraction probability as that from the degenerate pair and the linear combination of the degenerate pair gives the same reaction probability. In addition, the symmetric and asymmetric stretching modes of the reactant NH3 have comparable efficacies on driving the reaction.
Collapse
Affiliation(s)
- Hongwei Song
- 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
| | - Weiyu Xie
- Institute of Chemical Materials, China Academy of Engineering Physics, Mianyang 621900, China
| | - Chaoyang Zhang
- Institute of Chemical Materials, China Academy of Engineering Physics, Mianyang 621900, China
| | - Minghui Yang
- 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.,Wuhan National Laboratory for Optoelectronics, Huazhong University of Science and Technology, Wuhan 430071, China
| |
Collapse
|
19
|
Dissection of the Multichannel Reaction O( 3P) + C 2H 2: Differential Cross-Sections and Product Energy Distributions. MOLECULES (BASEL, SWITZERLAND) 2022; 27:molecules27030754. [PMID: 35164017 PMCID: PMC8838145 DOI: 10.3390/molecules27030754] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/25/2021] [Revised: 01/17/2022] [Accepted: 01/19/2022] [Indexed: 11/17/2022]
Abstract
The O(3P) + C2H2 reaction plays an important role in hydrocarbon combustion. It has two primary competing channels: H + HCCO (ketenyl) and CO + CH2 (triplet methylene). To further understand the microscopic dynamic mechanism of this reaction, we report here a detailed quasi-classical trajectory study of the O(3P) + C2H2 reaction on the recently developed full-dimensional potential energy surface (PES). The entrance barrier TS1 is the rate-limiting barrier in the reaction. The translation of reactants can greatly promote reactivity, due to strong coupling with the reaction coordinate at TS1. The O(3P) + C2H2 reaction progress through a complex-forming mechanism, in which the intermediate HCCHO lives at least through the duration of a rotational period. The energy redistribution takes place during the creation of the long-lived high vibrationally (and rotationally) excited HCCHO in the reaction. The product energy partitioning of the two channels and CO vibrational distributions agree with experimental data, and the vibrational state distributions of all modes of products present a Boltzmann-like distribution.
Collapse
|
20
|
Mallick S, Kumar P. Effect of microsolvation on the mode specificity of the OH˙(H 2O) + HCl reaction. Phys Chem Chem Phys 2021; 23:25246-25255. [PMID: 34734608 DOI: 10.1039/d1cp01300a] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The present study investigates the mode specificity in the microsolvated OH˙(H2O) + HCl reaction using on-the-fly direct dynamics simulation. To the best of our knowledge, this is the first study which aims to gain insights into the effect of microsolvation on the mode selectivity. Our investigation reveals that, similar to the gas phase OH˙ + HCl reaction, the microsolvated reaction is also predominantly affected by the vibrational excitation of the HCl mode, whereas the OH vibrational mode behaves as a spectator. Interestingly, in contrast to the behavior of the bare reaction, the integral cross section at the ground state of the microsolvated reaction decreases with an increase in translational energy. However, for the vibrational excited states, the reactivity of the microsolvated reaction is found to be higher than that of the bare reaction within the selected range of translational energies.
Collapse
Affiliation(s)
- Subhasish Mallick
- Department of Chemistry, Malaviya National Institute of Technology Jaipur, Jaipur, 302017, India.
| | - Pradeep Kumar
- Department of Chemistry, Malaviya National Institute of Technology Jaipur, Jaipur, 302017, India.
| |
Collapse
|
21
|
Papp D, Czakó G. Vibrational mode-specific dynamics of the F( 2P 3/2) + C 2H 6 → HF + C 2H 5 reaction. J Chem Phys 2021; 155:154302. [PMID: 34686045 DOI: 10.1063/5.0069658] [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/14/2022] Open
Abstract
We investigate the competing effect of vibrational and translational excitation and the validity of the Polanyi rules in the early- and negative-barrier F(2P3/2) + C2H6 → HF + C2H5 reaction by performing quasi-classical dynamics simulations on a recently developed full-dimensional multi-reference analytical potential energy surface. The effect of five normal-mode excitations of ethane on the reactivity, the mechanism, and the post-reaction energy flow is followed through a wide range of collision energies. Promoting effects of vibrational excitations and interaction time, related to the slightly submerged barrier, are found to be suppressed by the early-barrier-induced translational enhancement, in contrast to the slightly late-barrier Cl + C2H6 reaction. The excess vibrational energy mostly converts into ethyl internal excitation while collision energy is transformed into product separation. The substantial reaction energy excites the HF vibration, which tends to show mode-specificity and translational energy dependence as well. With increasing collision energy, direct stripping becomes dominant over the direct rebound and indirect mechanisms, being basically independent of reactant excitation.
Collapse
Affiliation(s)
- Dóra Papp
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group, Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged, Rerrich Béla tér 1, Szeged H-6720, Hungary
| | - Gábor Czakó
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group, Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged, Rerrich Béla tér 1, Szeged H-6720, Hungary
| |
Collapse
|
22
|
Csorba B, Szabó P, Góger S, Lendvay G. The Role of Zero-Point Vibration and Reactant Attraction in Exothermic Bimolecular Reactions with Submerged Potential Barriers: Theoretical Studies of the R + HBr → RH + Br (R = CH 3, HO) Systems. J Phys Chem A 2021; 125:8386-8396. [PMID: 34543008 PMCID: PMC8488937 DOI: 10.1021/acs.jpca.1c05839] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
The dynamics of the reactions CH3 + HBr → CH4 + Br and HO + HBr → H2O + Br have been studied using the quasiclassical trajectory method to explore the interplay of the vibrational excitation of the breaking bond and the potential energy surface characterized by a prereaction van der Waals well and a submerged barrier to reaction. The attraction between the reactants is favorable for the reaction, because it brings together the reactants without any energy investment. The reaction can be thought to be controlled by capture. The trajectory calculations indeed provide excitation functions typical to capture: the reaction cross sections diverge when the collision energy is reduced toward zero. Excitation of reactant vibration accelerates both reactions. The barrier on the potential surface is so early that the coupling between the degrees of freedom at the saddle point geometry is negligible. However, the trajectory calculations show that when the breaking bond is stretched at the time of the encounter, an attractive force arises, as if the radical approached a HBr molecule whose bond is partially broken. As a result, the dynamics of the reaction are controlled more by the temporary "dynamical", vibrationally induced than by the "static" van der Waals attraction even when the reactants are in vibrational ground state. The cross sections are shown to drop to very small values when the amplitude of the breaking bond's vibration is artificially reduced, which provides an estimate of the reactivity due to the "static" attraction. Without zero-point vibration these reactions would be very slow, which is a manifestation of a unique quantum effect. Reactions where the reactivity is determined by dynamical factors such as the vibrationally enhanced attraction are found to be beyond the range of applicability of Polanyi's rules.
Collapse
Affiliation(s)
- Benjámin Csorba
- Institute of Materials and Environmental Chemistry, Research Centre for Natural Sciences, Magyar tudósok krt. 2., H-1117 Budapest, Hungary
| | - Péter Szabó
- Institute of Materials and Environmental Chemistry, Research Centre for Natural Sciences, Magyar tudósok krt. 2., H-1117 Budapest, Hungary
| | - Szabolcs Góger
- Institute of Materials and Environmental Chemistry, Research Centre for Natural Sciences, Magyar tudósok krt. 2., H-1117 Budapest, Hungary
| | - György Lendvay
- Institute of Materials and Environmental Chemistry, Research Centre for Natural Sciences, Magyar tudósok krt. 2., H-1117 Budapest, Hungary.,Center for Natural Sciences, Faculty of Engineering, University of Pannonia, Egyetem u. 10. Veszprém, 8200 Hungary
| |
Collapse
|
23
|
Shi W, Wang K, Zhang P, Yu L, Li A. Mode-specific dynamics in multichannel reaction NH + + H 2. Phys Chem Chem Phys 2021; 23:20352-20358. [PMID: 34490857 DOI: 10.1039/d1cp02603h] [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 vibrational- and rotational-mode specificity of the multichannel NH+ + H2 reaction is studied on a recently constructed ab initio-based global potential energy surface using an initial state selected quasi-classical trajectory method, and the trajectories are analyzed using an isometric feature mapping and k-means approach. All excitation modes promote two reactions (R1: NH'+ + H2 → NH+ + HH' and R4: NH'+ + H2 → NH2+ + H') where both NH and HH bonds are broken, but reduce the reactivity of the proton-transfer reaction R2 (NH'+ + H2 → N + H'H2+) at low collision energies. For the hydrogen-transfer reaction R3 (NH'+ + H2 → HNH'+ + H), the rotational excitation of NH+ enhances the reactivity remarkably, while its vibrational excitation has an inhibiting effect on the reaction. The trajectory analyses show that the vibrational and rotational excitations of NH+ make R3 tend to go over a submerged saddle point instead of extracting hydrogen atoms directly. On the other hand, the motions of the H2 reactant facilitate the enhancement of the reactivity but they do not affect the mechanism of R3. In addition, the results suggest that the coupling of the isometric feature mapping and the k-means approach in the trajectory analysis is an appropriate tool for reaction-dynamics studies.
Collapse
Affiliation(s)
- Weiliang Shi
- Key Laboratory of Synthetic and Natural Functional Molecule of the Ministry of Education, College of Chemistry and Materials Science, Northwest University, Xi'an, 710127, P. R. China.
| | - Kun Wang
- Key Laboratory of Synthetic and Natural Functional Molecule of the Ministry of Education, College of Chemistry and Materials Science, Northwest University, Xi'an, 710127, P. R. China.
| | - Ping Zhang
- Key Laboratory of Synthetic and Natural Functional Molecule of the Ministry of Education, College of Chemistry and Materials Science, Northwest University, Xi'an, 710127, P. R. China.
| | - Le Yu
- Key Laboratory of Synthetic and Natural Functional Molecule of the Ministry of Education, College of Chemistry and Materials Science, Northwest University, Xi'an, 710127, P. R. China.
| | - Anyang Li
- Key Laboratory of Synthetic and Natural Functional Molecule of the Ministry of Education, College of Chemistry and Materials Science, Northwest University, Xi'an, 710127, P. R. China.
| |
Collapse
|
24
|
Papp D, Li J, Guo H, Czakó G. Vibrational mode-specificity in the dynamics of the Cl + C 2H 6 → HCl + C 2H 5 reaction. J Chem Phys 2021; 155:114303. [PMID: 34551541 DOI: 10.1063/5.0062677] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Abstract
We report a detailed dynamics study on the mode-specificity of the Cl + C2H6 → HCl + C2H5 H-abstraction reaction. We perform quasi-classical trajectory simulations using a recently developed high-level ab initio full-dimensional potential energy surface by exciting five different vibrational modes of ethane at four collision energies. We find that all the studied vibrational excitations, except that of the CC-stretching mode, clearly promote the title reaction, and the vibrational enhancements are consistent with the predictions of the Sudden Vector Projection (SVP) model, with the largest effect caused by the CH-stretching excitations. Intramolecular vibrational redistribution is also monitored for the differently excited ethane molecule. Our results indicate that the mechanism of the reaction changes with increasing collision energy, with no mode-specificity at high energies. The initial translational energy mostly converts into product recoil, while a significant part of the excess vibrational energy remains in the ethyl radical. An interesting competition between translational and vibrational energies is observed for the HCl vibrational distribution: the effect of exciting the low-frequency ethane modes, having small SVP values, is suppressed by translational excitation, whereas a part of the excess vibrational energy pumped into the CH-stretching modes (larger SVP values) efficiently flows into the HCl vibration.
Collapse
Affiliation(s)
- Dóra Papp
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group, Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged, Rerrich Béla tér 1, Szeged H-6720, Hungary
| | - Jun Li
- School of Chemistry and Chemical Engineering and Chongqing Key Laboratory of Theoretical and Computational Chemistry, Chongqing University, Chongqing 401331, China
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, USA
| | - Gábor Czakó
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group, Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged, Rerrich Béla tér 1, Szeged H-6720, Hungary
| |
Collapse
|
25
|
Xin R, Xiang H, Tian L, Li Y, Song H. Kinetic and Dynamic Studies of the F( 2P) + ND 3 → DF + ND 2 Reaction. J Phys Chem A 2021; 125:8025-8032. [PMID: 34478289 DOI: 10.1021/acs.jpca.1c06515] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The fast F reaction with NH3 poses a big challenge to experimental studies because of secondary chemical and collisional reactions. The quasi-classical trajectory method is utilized to investigate the mode specificity, product energy disposal, and temperature dependence of the thermal rate coefficient of F + ND3 → DF + ND2 on a recently developed potential energy surface. The effect of isotopic substitution is explored by comparing the F + ND3 reaction with the F + NH3 reaction. The computed results permit a better understanding of the F + ammonia reaction. The DF vibrational state has a Λ-type distribution, in accordance with the experimental measurement by the fast flow reactor technique. The product ND2 is dominantly populated in the ground state, and a considerable amount of ND2 is produced in the fundamental states of the bending mode. The similar vibrational state distributions of HF and NH2 in the F + NH3 reaction indicate a weak isotopic substitution effect on the product energy disposal. Exciting the umbrella mode of ND3 suppresses the reaction at low energies below 5 kcal mol-1, in sharp contrast to the observation in the F + NH3 reaction. These dynamical behaviors can be partially explained by the sudden vector projection model. In addition, the thermal rate coefficient of F + ND3 shows no temperature dependence in the range between 150 and 2000 K. There exists an inverse kinetic isotope effect at temperatures from 150 to 1500 K.
Collapse
Affiliation(s)
- Rong Xin
- 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.,College of Physical Science and Technology, Huazhong Normal University, Wuhan 430079, China
| | - Haipan Xiang
- 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.,College of Physical Science and Technology, Huazhong Normal University, Wuhan 430079, China
| | - Li Tian
- 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.,College of Physical Science and Technology, Huazhong Normal University, Wuhan 430079, 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
|
26
|
Wester R. Fifty years of nucleophilic substitution in the gas phase. MASS SPECTROMETRY REVIEWS 2021; 41:627-644. [PMID: 34060119 PMCID: PMC9291629 DOI: 10.1002/mas.21705] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 03/31/2021] [Revised: 05/05/2021] [Accepted: 05/05/2021] [Indexed: 06/12/2023]
Abstract
Bimolecular nucleophilic substitution ( S N 2 ) reactions have become a model system for the investigation of structure-reactivity relationships, stereochemistry, solvent influences, and detailed atomistic dynamics. In this review, the progress during five decades of experimental and theoretical research on gas phase S N 2 reactions is discussed. Many advancements of the employed methods have led to a tremendous increase in our understanding of the properties and the dynamics of these reactions. For reactions involving six atoms a quantitative agreement of the differential reactive scattering cross sections has already been achieved, in the future it is expected that even larger polyatomic reactions systems become tractable. Furthermore, studies with higher precision, improved reactant control, and a more accurate theoretical treatment of quantum effects are envisioned.
Collapse
Affiliation(s)
- Roland Wester
- Institut für Ionenphysik und Angewandte PhysikUniversität InnsbruckTechnikerstraße 256020 InnsbruckAustria
| |
Collapse
|
27
|
Sahoo J, Rawat AMS, Mahapatra S. Theoretical Study of the Energy Disposal Mechanism and the State-Resolved Quantum Dynamics of the H + LiH + → H 2 + Li + Reaction. J Phys Chem A 2021; 125:3387-3397. [PMID: 33876630 DOI: 10.1021/acs.jpca.1c01811] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Despite several studies in the literature, the detailed quantum state-to-state level mechanism of the astrophysically important exoergic barrierless H + LiH+ → H2 + Li+ reaction is yet to be understood. In this work, we have investigated the energy disposal mechanism of the reaction in terms of integral reaction cross section, product internal state distributions, differential cross section, and rate constant. Fully converged and Coriolis coupled quantum mechanical calculations based on a time-dependent wave packet method have been performed at the state-to-state level on the ab initio electronic ground state potential energy surface (PES) constructed by Martinazzo et al. (J. Chem. Phys. 2003, 119, 11241-11248). The agreement between the present quantum mechanical and previous quasi-classical results is found even at very low relative translational energies of reagents. A non-statistical inverse Boltzmann vibrational distribution for the product is found. This is attributed to the "attractive" nature of the underlying PES, which facilitates the excess energy release mostly as product vibration (60-80%). The energy disposal in products is found to be unaffected by the rovibrational excitation of the reagent diatom due to the weak coupling between the vibrational modes of the reagent and the product. The mild effect of collision energy on the product energy disposal is ascribed to the effective coupling between the translational modes of the reagent and the product. It is found that the collisions lead to the formation of the product H2 in its rovibrationally excited levels.
Collapse
Affiliation(s)
- Jayakrushna Sahoo
- School of Chemistry, University of Hyderabad, Hyderabad 500046, India
| | | | - S Mahapatra
- School of Chemistry, University of Hyderabad, Hyderabad 500046, India
| |
Collapse
|
28
|
Chen Q, Hu X, Guo H, Xie D. Insights into the Formation of Hydroxyl Radicals with Nonthermal Vibrational Excitation in the Meinel Airglow. J Phys Chem Lett 2021; 12:1822-1828. [PMID: 33577325 DOI: 10.1021/acs.jpclett.1c00159] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
To understand night time airglow in the Meinel bands and heat conversion from the highly excited OH radicals in the upper atmosphere via the important atmospheric reaction H + O3 → OH + O2, we report here a quasi-classical trajectory study of the reaction dynamics on a recently developed full-dimensional potential energy surface (PES). Our results indicate that the reaction energy of this highly exoergic reaction is almost exclusively channeled into the vibration of the OH product, underscoring an extreme departure from the statistical limit. The calculated OH vibrational distribution is highly inverted and peaks near the highest accessible vibrational state, in excellent agreement with experimental observations, validating the accuracy of the PES. More importantly, the dynamical origin of the nonthermal excitation of the OH vibrational mode is identified by its large projection onto the reaction coordinate at a small potential barrier in the entrance channel, which controls the energy flow into various degrees of freedom in the products.
Collapse
Affiliation(s)
- Qixin Chen
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Xixi Hu
- Kuang Yaming Honors School, Institute for Brain Sciences, Jiangsu Key Laboratory of Vehicle Emissions Control, Center of Modern Analysis, Nanjing University, Nanjing 210023, China
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| | - Daiqian Xie
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| |
Collapse
|
29
|
Quasi-Classical Trajectory Study of the CN + NH 3 Reaction Based on a Global Potential Energy Surface. Molecules 2021; 26:molecules26040994. [PMID: 33668582 PMCID: PMC7918900 DOI: 10.3390/molecules26040994] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/26/2021] [Revised: 02/10/2021] [Accepted: 02/11/2021] [Indexed: 11/17/2022] Open
Abstract
Based on a combination of valence-bond and molecular mechanics functions which were fitted to high-level ab initio calculations, we constructed an analytical full-dimensional potential energy surface, named PES-2020, for the hydrogen abstraction title reaction for the first time. This surface is symmetrical with respect to the permutation of the three hydrogens in ammonia, it presents numerical gradients and it improves the description presented by previous theoretical studies. In order to analyze its quality and accuracy, stringent tests were performed, exhaustive kinetics and dynamics studies were carried out using quasi-classical trajectory calculations, and the results were compared with the available experimental evidence. Firstly, the properties (geometry, vibrational frequency and energy) of all stationary points were found to reasonably reproduce the ab initio information used as input; due to the complicated topology with deep wells in the entrance and exit channels and a “submerged” transition state, the description of the intermediate complexes was poorer, although it was adequate to reasonably simulate the kinetics and dynamics of the title reaction. Secondly, in the kinetics study, the rate constants simulated the experimental data in the wide temperature range of 25–700 K, improving the description presented by previous theoretical studies. In addition, while previous studies failed in the description of the kinetic isotope effects, our results reproduced the experimental information. Finally, in the dynamics study, we analyzed the role of the vibrational and rotational excitation of the CN(v,j) reactant and product angular scattering distribution. We found that vibrational excitation by one quantum slightly increased reactivity, thus reproducing the only experimental measurement, while rotational excitation strongly decreased reactivity. The scattering distribution presented a forward-backward shape, associated with the presence of deep wells along the reaction path. These last two findings await experimental confirmation.
Collapse
|
30
|
Tian L, Song H, Yang M. Effects of bending excitation on the reaction dynamics of fluorine atoms with ammonia. Phys Chem Chem Phys 2021; 23:2715-2722. [PMID: 33491710 DOI: 10.1039/d0cp05790h] [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
Vibrational excitation has been established as an efficient way to control the chemical reaction outcome. Stretching vibration of polyatomic molecules is believed to be efficient to promote abstraction reactions since energy is placed directly into the breaking bond. In this work, we report on a counterexample showing that exciting the low-frequency umbrella bending mode of ammonia enhances its reaction with fluorine atoms much more than exciting the high-frequency symmetric or asymmetric stretching mode over a wide range of collision energy, validated using both quasiclassical trajectory simulations and full-dimensional quantum dynamics calculations under the centrifugal-sudden approximation. This interesting mode-specific reaction dynamic originates from the increased chance of capturing the fluorine atom by ammonia due to the enlarged attractive interaction between them and the enhancement of the direct stripping reaction mediated by two submerged barriers.
Collapse
Affiliation(s)
- Li Tian
- 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. and 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, Wuhan Institute of Physics and Mathematics, Innovation Academy for Precision Measurement Science and Technology, Chinese Academy of Sciences, Wuhan 430071, China.
| | - Minghui Yang
- 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
|
31
|
Xin R, Pan M, Song H, Yang M. Mode- and Bond-Selected Reaction of H with Local Mode Molecule HDS. J Phys Chem A 2020; 124:10162-10170. [PMID: 33252233 DOI: 10.1021/acs.jpca.0c09415] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Understanding mode- and bond-selected dynamics of elementary chemical reactions is of central importance in molecular reaction dynamics. The initial state-selected time-dependent wave packet method is employed to study the mode and bond selectivity, isotopic branching ratio, and temperature dependence of rate constants of the two-channel reaction of H with local mode molecule HDS. For the abstraction channel, fundamental excitation of the HS (DS) bond of the reactant HDS significantly enhances the H-abstraction (D-abstraction) reaction, whose efficacy is higher than the same amount of translational energy except at low energies just above the energy threshold. This is in sharp contrast to the prediction of Polanyi rules: translational energy is more efficient than vibrational energy in enhancing a reaction with an early barrier. The recent sudden vector projection model is then applied to rationalize the observed mode specificity, which, however, shows that the translational mode vector has a larger coupling with the reaction coordinate than the stretching vector of the active bond, implying a reversed relative efficacy on promoting the reaction as well. In contrast, the mode and bond specificity for the exchange channel is not as strong as for the abstraction channel due to the regulation of the shallow well along the reaction path.
Collapse
Affiliation(s)
- Rong Xin
- 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.,College of Physical Science and Technology, Huazhong Normal University, Wuhan 430079, China
| | - 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.,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, Wuhan Institute of Physics and Mathematics, Innovation Academy for Precision Measurement Science and Technology, Chinese Academy of Sciences, Wuhan 430071, China
| | - Minghui Yang
- 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
|
32
|
Goswami S, Sahoo J, Paul SK, Rao TR, Mahapatra S. Effect of Reagent Vibration and Rotation on the State-to-State Dynamics of the Hydrogen Exchange Reaction, H + H 2 → H 2 + H. J Phys Chem A 2020; 124:9343-9359. [PMID: 33124827 DOI: 10.1021/acs.jpca.0c06707] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
State-to-state dynamics of the benchmark hydrogen exchange reaction H + H2 (v = 0-4, j = 0-3) → H2 (v', j') + H is investigated with the aid of the real wave packet approach of Gray and Balint-Kurti (J. Chem. Phys. 1998, 108, 950-962) and electronic ground BKMP2 potential energy surface of Boothroyd et al. (J. Chem. Phys. 1996, 104, 7139-7152). Initial state-selected and product state-resolved reaction probabilities, integral cross section, and product diatom vibrational and rotational level populations at a few collision energies are reported to elucidate the energy disposal mechanism. State-specific thermal rate constants are also calculated and compared with the available literature results. Coriolis coupling terms of the nuclear Hamiltonian are included, and calculations are parallelized over the helicity quantum number, Ω'. Attempts are made, in particular, to study the effect of reagent vibrational and rotational excitations on the dynamical attributes. It is found that the calculations become computationally expensive with reagent vibrational and rotational excitation. Reagent vibrational excitation is found to enhance the reactivity and has significant impact on the energy disposal to the vibrational and rotational degrees of freedom of the product. The interplay of reagent translational and vibrational energy on the product vibrational distribution unfolds an important aspect of the energy disposal mechanism. The effect of reagent rotation on the state-to-state dynamics is found not to be very significant, and the weak effect turns out to be specific to v'.
Collapse
Affiliation(s)
- Sugata Goswami
- School of Chemistry, University of Hyderabad, Hyderabad 500046, India
| | - Jayakrushna Sahoo
- School of Chemistry, University of Hyderabad, Hyderabad 500046, India
| | - Suranjan K Paul
- School of Chemistry, University of Hyderabad, Hyderabad 500046, India
| | - T Rajagopala Rao
- School of Chemistry, University of Hyderabad, Hyderabad 500046, India
| | - S Mahapatra
- School of Chemistry, University of Hyderabad, Hyderabad 500046, India
| |
Collapse
|
33
|
Li J, Zhao B, Xie D, Guo H. Advances and New Challenges to Bimolecular Reaction Dynamics Theory. J Phys Chem Lett 2020; 11:8844-8860. [PMID: 32970441 DOI: 10.1021/acs.jpclett.0c02501] [Citation(s) in RCA: 32] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
Dynamics of bimolecular reactions in the gas phase are of foundational importance in combustion, atmospheric chemistry, interstellar chemistry, and plasma chemistry. These collision-induced chemical transformations are a sensitive probe of the underlying potential energy surface(s). Despite tremendous progress in past decades, our understanding is still not complete. In this Perspective, we survey the recent advances in theoretical characterization of bimolecular reaction dynamics, stimulated by new experimental observations, and identify key new challenges.
Collapse
Affiliation(s)
- Jun Li
- School of Chemistry and Chemical Engineering & Chongqing Key Laboratory of Theoretical and Computational Chemistry, Chongqing University, Chongqing 401331, China
| | - Bin Zhao
- Theoretische Chemie, Fakultät für Chemie, Universität Bielefeld, Universitätsstraße 25, D-33615 Bielefeld, Germany
| | - Daiqian Xie
- Institute of Theoretical and Computational Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| |
Collapse
|
34
|
|
35
|
Roy S, K. J. N, Tiwari N, Tiwari AK. Energetics and dynamics of CH4 and H2O dissociation on metal surfaces. INT REV PHYS CHEM 2020. [DOI: 10.1080/0144235x.2020.1765598] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
Affiliation(s)
- Sudipta Roy
- Department of Chemical Sciences, Indian Institute of Science Education and Research Kolkata, Mohanpur, India
| | - Nayanthara K. J.
- Department of Chemical Sciences, Indian Institute of Science Education and Research Kolkata, Mohanpur, India
| | - Nidhi Tiwari
- Department of Chemical Sciences, Indian Institute of Science Education and Research Kolkata, Mohanpur, India
| | - Ashwani K. Tiwari
- Department of Chemical Sciences, Indian Institute of Science Education and Research Kolkata, Mohanpur, India
| |
Collapse
|
36
|
Michaelsen T, Bastian B, Strübin P, Meyer J, Wester R. Proton transfer dynamics modified by CH-stretching excitation. Phys Chem Chem Phys 2020; 22:12382-12388. [PMID: 32319988 DOI: 10.1039/d0cp00727g] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Gaining insight how specific rovibrational states influence reaction kinetics and dynamics is a fundamental goal of physical chemistry. Purely statistical approaches often fail to predict the influence of a specific state on the reaction outcome, evident in a great number of both experimental and theoretical studies. Most detailed insight in atomistic reaction mechanisms is achieved using accurate collision experiments and high level dynamics calculations. For ion-molecule reactions such experiments are scarce. Here we show the influence of symmetric CH-stretching vibration on the rate and dynamics of proton transfer in the reaction of F- + CH3I. We find a pronounced shift in the reaction dynamics for excited reactions from indirect to preferred direct dynamics at higher collision energy. Moreover, excited reactions occur at larger impact parameters. Finally, we compare vibrational excitation with collision energy and find that vibration is overall more efficient in promoting reactivity, which agrees with recent theoretical calculations.
Collapse
Affiliation(s)
- Tim Michaelsen
- Institut für Ionenphysik und Angewandte Physik, Universität Innsbruck, Technikerstraße 25/3, 6020 Innsbruck, Austria.
| | - Björn Bastian
- Institut für Ionenphysik und Angewandte Physik, Universität Innsbruck, Technikerstraße 25/3, 6020 Innsbruck, Austria.
| | - Patrick Strübin
- Institut für Ionenphysik und Angewandte Physik, Universität Innsbruck, Technikerstraße 25/3, 6020 Innsbruck, Austria.
| | - Jennifer Meyer
- Institut für Ionenphysik und Angewandte Physik, Universität Innsbruck, Technikerstraße 25/3, 6020 Innsbruck, Austria.
| | - Roland Wester
- Institut für Ionenphysik und Angewandte Physik, Universität Innsbruck, Technikerstraße 25/3, 6020 Innsbruck, Austria.
| |
Collapse
|
37
|
Michaelsen T, Bastian B, Ayasli A, Strübin P, Meyer J, Wester R. Influence of Vibrational Excitation on the Reaction of F - with CH 3I: Spectator Mode Behavior, Enhancement, and Suppression. J Phys Chem Lett 2020; 11:4331-4336. [PMID: 32383877 PMCID: PMC7277560 DOI: 10.1021/acs.jpclett.0c01095] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/08/2020] [Accepted: 05/04/2020] [Indexed: 06/10/2023]
Abstract
Detailed insight into chemical reaction dynamics can be obtained by probing the effect of mode-specific vibrational excitation. Suppression or enhancement of reactivity is possible as is already known from the Polanyi rules. In the reaction F- + CH3I, we found vibrational enhancement, suppression, and spectator mode dynamics in the four different reaction channels. For this system we have probed the influence of symmetric CH-stretching vibration over a collision energy range of 0.7-2.3 eV. Proton transfer is significantly enhanced, while for the nucleophilic substitution channel the spectator mode dynamics at lower collision energies unexpectedly move toward enhancement at higher collision energies. In contrast, for two halide abstraction channels, forming FI- and FHI-, we found an overall suppression, which stems mainly from a suppression of the FHI- product. We compare these results to quasiclassical trajectory calculations and with the sudden vector projection model.
Collapse
Affiliation(s)
- Tim Michaelsen
- Institut für Ionenphysik und Angewandte Physik, Universität Innsbruck, Technikerstraße 25/3, 6020 Innsbruck, Austria
| | - Björn Bastian
- Institut für Ionenphysik und Angewandte Physik, Universität Innsbruck, Technikerstraße 25/3, 6020 Innsbruck, Austria
| | - Atilay Ayasli
- Institut für Ionenphysik und Angewandte Physik, Universität Innsbruck, Technikerstraße 25/3, 6020 Innsbruck, Austria
| | - Patrick Strübin
- Institut für Ionenphysik und Angewandte Physik, Universität Innsbruck, Technikerstraße 25/3, 6020 Innsbruck, Austria
| | - Jennifer Meyer
- Institut für Ionenphysik und Angewandte Physik, Universität Innsbruck, Technikerstraße 25/3, 6020 Innsbruck, Austria
| | - Roland Wester
- Institut für Ionenphysik und Angewandte Physik, Universität Innsbruck, Technikerstraße 25/3, 6020 Innsbruck, Austria
| |
Collapse
|
38
|
Marjollet A, Welsch R. Nuclear quantum effects in state-selective scattering from ring polymer molecular dynamics. J Chem Phys 2020; 152:194113. [DOI: 10.1063/5.0004179] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Adrien Marjollet
- Center for Free-Electron Laser Science, DESY, Notkestrasse 85, 22607 Hamburg, Germany
- Department of Physics, Universität Hamburg, Jungiusstrasse 9, 20355 Hamburg, Germany
| | - Ralph Welsch
- Center for Free-Electron Laser Science, DESY, Notkestrasse 85, 22607 Hamburg, Germany
- The Hamburg Centre for Ultrafast Imaging, Luruper Chaussee 149, 22761 Hamburg, Germany
| |
Collapse
|
39
|
Abstract
Nonstatistical dynamics is important for many chemical reactions. The Rice-Ramsperger-Kassel-Marcus (RRKM) theory of unimolecular kinetics assumes a reactant molecule maintains a statistical microcanonical ensemble of vibrational states during its dissociation so that its unimolecular dynamics are time independent. Such dynamics results when the reactant's atomic motion is chaotic or irregular. Intrinsic non-RRKM dynamics occurs when part of the reactant's phase space consists of quasiperiodic/regular motion and a bottleneck exists, so that the unimolecular rate constant is time dependent. Nonrandom excitation of a molecule may result in short-time apparent non-RRKM dynamics. For rotational activation, the 2J + 1 K levels for a particular J may be highly mixed, making K an active degree of freedom, or K may be a good quantum number and an adiabatic degree of freedom. Nonstatistical dynamics is often important for bimolecular reactions and their intermediates and for product-energy partitioning of bimolecular and unimolecular reactions. Post–transition state dynamics is often highly complex and nonstatistical.
Collapse
Affiliation(s)
- Bhumika Jayee
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409, USA
| | - William L. Hase
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409, USA
| |
Collapse
|
40
|
Liu Y, Song H, Xie D, Li J, Guo H. Mode Specificity in the OH + HO2 → H2O + O2 Reaction: Enhancement of Reactivity by Exciting a Spectator Mode. J Am Chem Soc 2020; 142:3331-3335. [DOI: 10.1021/jacs.9b12467] [Citation(s) in RCA: 27] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Yang Liu
- School of Chemistry and Chemical Engineering, Chongqing University, Chongqing 401331, 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
| | - Daiqian Xie
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, 210093, China
| | - Jun Li
- School of Chemistry and Chemical Engineering, Chongqing University, Chongqing 401331, China
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| |
Collapse
|
41
|
Bal KM, Bogaerts A, Neyts EC. Ensemble-Based Molecular Simulation of Chemical Reactions under Vibrational Nonequilibrium. J Phys Chem Lett 2020; 11:401-406. [PMID: 31865709 DOI: 10.1021/acs.jpclett.9b03356] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
We present an approach to incorporate the effect of vibrational nonequilibrium in molecular dynamics (MD) simulations. A perturbed canonical ensemble, in which selected modes are excited to higher temperature while all others remain equilibrated at low temperature, is simulated by applying a specifically tailored bias potential. Our method can be readily applied to any (classical or quantum mechanical) MD setup at virtually no additional computational cost and allows the study of reactions of vibrationally excited molecules in nonequilibrium environments such as plasmas. In combination with enhanced sampling methods, the vibrational efficacy and mode selectivity of vibrationally stimulated reactions can then be quantified in terms of chemically relevant observables, such as reaction rates and apparent free energy barriers. We first validate our method for the prototypical hydrogen exchange reaction and then show how it can capture the effect of vibrational excitation on a symmetric SN2 reaction and radical addition on CO2.
Collapse
Affiliation(s)
- Kristof M Bal
- Research Group PLASMANT, Department of Chemistry , University of Antwerp , Universiteitsplein 1 , 2610 Antwerp , Belgium
| | - Annemie Bogaerts
- Research Group PLASMANT, Department of Chemistry , University of Antwerp , Universiteitsplein 1 , 2610 Antwerp , Belgium
| | - Erik C Neyts
- Research Group PLASMANT, Department of Chemistry , University of Antwerp , Universiteitsplein 1 , 2610 Antwerp , Belgium
| |
Collapse
|
42
|
Karmakar S, Keshavamurthy S. Intramolecular vibrational energy redistribution and the quantum ergodicity transition: a phase space perspective. Phys Chem Chem Phys 2020; 22:11139-11173. [DOI: 10.1039/d0cp01413c] [Citation(s) in RCA: 26] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
The onset of facile intramolecular vibrational energy flow can be related to features in the connected network of anharmonic resonances in the classical phase space.
Collapse
Affiliation(s)
- Sourav Karmakar
- Department of Chemistry
- Indian Institute of Technology
- Kanpur
- India
| | | |
Collapse
|
43
|
Jambrina PG, Zanchet A, Menéndez M, Herrero VJ, Aoiz FJ. Unexpected dynamical effects change the lambda-doublet propensity in the tunneling region for the O( 3P) + H 2 reaction. Phys Chem Chem Phys 2019; 21:25389-25396. [PMID: 31709441 DOI: 10.1039/c9cp04690a] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
One of the most relevant features of the O(3P) + H2 reaction is that it occurs on two different potential energy surfaces (PESs) of symmetries A' and A'' that correlate reactants and products. The respective saddle points, which correspond to a collinear arrangement, are the same for both PESs, whilst the barrier height rises more abruptly on the 3A' PES than on the 3A'' PES. Accordingly, the reactivity on the 3A'' PES should be always higher than on the 3A' PES. In this work, we present accurate quantum-scattering calculations showing that this is not always the case for rotationless reactants, where dynamical factors near the reaction threshold cause the 3A' PES to dominate at energies around the barrier. Further calculation of cross sections and Λ-doublet populations has allowed us to establish how the reaction mechanism changes from the deep tunneling regime to hyperthermal energies.
Collapse
Affiliation(s)
- P G Jambrina
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad de Salamanca, 37003, Salamanca, Spain
| | - A Zanchet
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad de Salamanca, 37003, Salamanca, Spain and Departamento de Química Física, Facultad de Química, Universidad Complutense de Madrid (Unidad Asociada CSIC), 28040 Madrid, Spain.
| | - M Menéndez
- Departamento de Química Física, Facultad de Química, Universidad Complutense de Madrid (Unidad Asociada CSIC), 28040 Madrid, Spain.
| | - V J Herrero
- Instituto de Estructura de la Materia, IEM-CSIC c/Serrano 123, 28006 Madrid, Spain
| | - F J Aoiz
- Departamento de Química Física, Facultad de Química, Universidad Complutense de Madrid (Unidad Asociada CSIC), 28040 Madrid, Spain.
| |
Collapse
|
44
|
QCT study of the vibrational and translational role in the H + C2H6(ν1, ν2, ν5, ν7, ν9 and ν10) reactions. Theor Chem Acc 2019. [DOI: 10.1007/s00214-019-2504-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|
45
|
Zhang Z, Gatti F, Zhang DH. Full dimensional quantum mechanical calculations of the reaction probability of the H + NH3 collision based on a mixed Jacobi and Radau description. J Chem Phys 2019; 150:204301. [DOI: 10.1063/1.5096047] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Zhaojun Zhang
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, People’s Republic of China
| | - Fabien Gatti
- ISMO, Institut des Sciences Moléculaires d’Orsay—UMR 8214 CNRS/Université Paris-Sud, F-91405 Orsay, France
| | - Dong H. Zhang
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, People’s Republic of China
| |
Collapse
|
46
|
Jiang B, Guo H. Dynamics in reactions on metal surfaces: A theoretical perspective. J Chem Phys 2019; 150:180901. [PMID: 31091904 DOI: 10.1063/1.5096869] [Citation(s) in RCA: 31] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Recent advances in theoretical characterization of reaction dynamics on metal surfaces are reviewed. It is shown that the widely available density functional theory of metals and their interactions with molecules have enabled first principles theoretical models for treating surface reaction dynamics. The new theoretical tools include methods to construct high-dimensional adiabatic potential energy surfaces, to characterize nonadiabatic processes within the electronic friction models, and to describe dynamics both quantum mechanically and classically. Three prototypical surface reactions, namely, dissociative chemisorption, Eley-Rideal reactions, and recombinative desorption, are surveyed with a focus on some representative examples. While principles governing gas phase reaction dynamics may still be applicable, the presence of the surface introduces a higher level of complexity due to strong interaction between the molecular species and metal substrate. Furthermore, most of these reactive processes are impacted by energy exchange with surface phonons and/or electron-hole pair excitations. These theoretical studies help to interpret and rationalize experimental observations and, in some cases, guide experimental explorations. Knowledge acquired in these fundamental studies is expected to impact many practical problems in a wide range of interfacial processes.
Collapse
Affiliation(s)
- 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 230026, China
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, USA
| |
Collapse
|
47
|
Ghosh S, Ray D, Tiwari AK. Effects of alloying on mode-selectivity in H2O dissociation on Cu/Ni bimetallic surfaces. J Chem Phys 2019; 150:114702. [DOI: 10.1063/1.5085696] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Smita Ghosh
- Department of Chemical Sciences, Indian Institute of Science Education and Research Kolkata, Mohanpur 741246, India
| | - Dhiman Ray
- Department of Chemistry, University of California Irvine, Irvine, California 92617, USA
| | - Ashwani K. Tiwari
- Department of Chemical Sciences, Indian Institute of Science Education and Research Kolkata, Mohanpur 741246, India
| |
Collapse
|
48
|
Unke OT, Brickel S, Meuwly M. Sampling reactive regions in phase space by following the minimum dynamic path. J Chem Phys 2019; 150:074107. [DOI: 10.1063/1.5082885] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023] Open
Affiliation(s)
- Oliver T. Unke
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| | - Sebastian Brickel
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| | - Markus Meuwly
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| |
Collapse
|
49
|
Gudem M, Hazra A. Mechanism of the Chemiluminescent Reaction between Nitric Oxide and Ozone. J Phys Chem A 2019; 123:715-722. [PMID: 30380861 DOI: 10.1021/acs.jpca.8b08812] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The gas phase reaction of nitric oxide with ozone to give chemiluminescence is used extensively for detection of nitrogen oxides. The molecular mechanism of chemiluminescence in this reaction is not known. So far, the only chemiluminescent systems studied in depth are certain cycloperoxides, which emit light following decomposition. Given our understanding of the mechanism of chemiluminescence in those molecules, one would expect by extension that in the NO + O3 reaction the chemiluminescent species (NO2 in this case) is formed in the excited state through a reaction pathway that diverges from the ground state pathway near the transition state. A systematic search for such a pathway leads us to conclude that such a mechanism is unlikely. Instead, our study suggests that chemiluminescence in the NO + O3 reaction is due to emission from the NO2 vibronic states associated with the ground (X̃ 2A1) and first excited (à 2B2) electronic states, which are populated in the nascent NO2 produced in the reaction. The vibronic coupling between the X̃ 2A1 and à 2B2 states of NO2 is due to a conical intersection (CI), which is geometrically and energetically close to the à 2B2 minimum energy geometry and only 1.3 eV higher than ground state NO2. Further, the CI is 1.2 eV lower than the energy of the NO + O3 reactants and therefore thermodynamically accessible following the reaction. An analysis of the product energy distribution indicates that the major fraction of the reaction energy is channeled into the vibrational modes of NO2, sufficient to populate the vibronic states of NO2 around the X̃/à CI. These vibronic states show dipole-allowed emission in a frequency range that is consistent with the observed broad chemiluminescence spectrum.
Collapse
Affiliation(s)
- Mahesh Gudem
- Department of Chemistry , Indian Institute of Science Education and Research Pune , Dr. Homi Bhabha Road , Pune 411008 , Maharashtra , India
| | - Anirban Hazra
- Department of Chemistry , Indian Institute of Science Education and Research Pune , Dr. Homi Bhabha Road , Pune 411008 , Maharashtra , India
| |
Collapse
|
50
|
State-to-state dynamics of the Cl(2P) + C2H6(ν5, ν1 = 0, 1) → HCl(v′, j′) + C2H5 hydrogen abstraction reactions. Theor Chem Acc 2019. [DOI: 10.1007/s00214-019-2416-3] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|