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
|
Li J, Tu Z, Xiang H, Li Y, Song H. Theoretical studies on the kinetics and dynamics of the BeH + + H 2O reaction: comparison with the experiment. Phys Chem Chem Phys 2023; 25:20997-21005. [PMID: 37503894 DOI: 10.1039/d3cp02322b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/29/2023]
Abstract
The reaction of BeH+ with background gaseous H2O may play a role in qubit loss for quantum information processing with Be+ as trapped ions, and yet its reaction mechanism has not been well understood until now. In this work, a globally accurate, full-dimensional ground-state potential energy surface (PES) for the BeH+ + H2O reaction was constructed by fitting a total of 170 438 ab initio energy points at the level of RCCSD(T)-F12/aug-cc-pVTZ using the fundamental invariant-neural network method. The total root-mean-square error of the final PES was 0.178 kcal mol-1. For comparison, quasi-classical trajectory calculations were carried out on the PES at an experimental temperature of 150 K. The obtained thermal rate constant and product branching ratio of the BeD+ + H2O reaction agreed quite well with experimental results. In addition, the vibrational state distributions and energy disposals of the products were calculated and rationalized using the sudden vector projection model.
Collapse
Affiliation(s)
- Jiaqi Li
- College of Physical Science and Technology, Huazhong Normal University, Wuhan 430079, China.
- State Key Laboratory of Magnetic Resonance and Atomic and Molecular Physics, Innovation Academy for Precision Measurement Science and Technology, Chinese Academy of Sciences, Wuhan 430071, China.
| | - Zhao Tu
- 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.
- School of Chemical and Environmental Engineering, Hubei Minzu University, Enshi 445000, 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.
- School of Physics and Electronics, Hunan University, Changsha 410082, 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
|
3
|
Pavković N, Milovanović B, Stanojević A, Etinski M, Petković M. Proton leap: shuttling of protons onto benzonitrile. Phys Chem Chem Phys 2022; 24:3958-3969. [PMID: 35099492 DOI: 10.1039/d1cp04338b] [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 detailed description of chemical transformations in the interstellar medium allows deciphering the origin of a number of small and medium - sized organic molecules. We present density functional theory analysis of proton transfer from the trihydrogen cation and the ethenium cation to benzonitrile, a recently discovered species in the Taurus Molecular Cloud 1. Detailed energy transformations along the reaction paths were analyzed using the interacting quantum atoms methodology, which elucidated how the proton carrier influences the lightness to deliver the proton to benzonitrile's nitrogen atom. The proton carriers' deformation energy represents the largest destabilizing effect, whereas a proton's promotion energy, the benzonitrile-proton Coulomb attraction, as well as non-classical benzonitrile-proton and carrier-proton interaction are the dominant stabilizing energy components. As two ion-molecule reactions proceed without energy barriers, rate constants were estimated using the classical capture theory and were found to be an order of magnitude larger for the reaction with the trihydrogen cation compared to that with the ethenium cation (∼10-8 and 10-9 cm3 s-1, respectively). These results were obtained both with quantum chemical and ab initio molecular dynamics simulations (the latter performed at 10 K and 100 K), confirming that up to 100 K both systems choose energetically undemanding routes by tracking the corresponding minimum energy paths. A concept of a turning point is introduced, which is an equivalent to the transition state in barrierless reactions.
Collapse
Affiliation(s)
- Nemanja Pavković
- University of Belgrade - Faculty of Physical Chemistry, Studentski trg 12-16, 11158 Belgrade, Serbia.
| | - Branislav Milovanović
- University of Belgrade - Faculty of Physical Chemistry, Studentski trg 12-16, 11158 Belgrade, Serbia.
| | - Ana Stanojević
- University of Belgrade - Faculty of Physical Chemistry, Studentski trg 12-16, 11158 Belgrade, Serbia.
| | - Mihajlo Etinski
- University of Belgrade - Faculty of Physical Chemistry, Studentski trg 12-16, 11158 Belgrade, Serbia.
| | - Milena Petković
- University of Belgrade - Faculty of Physical Chemistry, Studentski trg 12-16, 11158 Belgrade, Serbia.
| |
Collapse
|
4
|
Product Vibrational State Distributions of the F CH 3OH Reaction on a Full-Dimensional Accurate Potential Energy Surface. CHINESE J CHEM PHYS 2022. [DOI: 10.1063/1674-0068/cjcp2111252] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022]
|
5
|
Zhu Y, Li R, Song H. Kinetic and dynamic studies of the NH 2+ + H 2 reaction on a high-level ab initio potential energy surface. Phys Chem Chem Phys 2022; 24:25663-25672. [DOI: 10.1039/d2cp03859e] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Abstract
The dynamics and kinetics of the NH2+ + H2 reaction are investigated on a newly developed ab initio potential energy surface using the quasi-classical trajectory method.
Collapse
Affiliation(s)
- Yongfa Zhu
- School of Chemistry and Chemical Engineering, Hubei Polytechnic University, Huangshi 435003, China
| | - Rui Li
- School of Chemistry and Chemical Engineering, Hubei Polytechnic University, Huangshi 435003, 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
| |
Collapse
|
6
|
Saito K, Hashimoto Y, Takayanagi T. Ring-Polymer Molecular Dynamics Calculations of Thermal Rate Coefficients and Branching Ratios for the Interstellar H 3+ + CO → H 2 + HCO +/HOC + Reaction and Its Deuterated Analogue. J Phys Chem A 2021; 125:10750-10756. [PMID: 34918514 DOI: 10.1021/acs.jpca.1c09160] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The reaction between H3+ and CO is important in understanding the H3+ destruction mechanism in the interstellar medium. In this work, thermal rate coefficients for the H3+ + CO and D3+ + CO reactions are calculated using ring-polymer molecular dynamics (RPMD) on a high-level machine-learning potential energy surface. The RPMD results agree well with the classical molecular dynamics results, where nuclear quantum effects are completely ignored, whereas the agreement between the RPMD results and the previous quasi-classical trajectory is good only at low temperatures. The calculated [HCO+]/[HOC+] product branching ratios decrease as the temperature increases, and the product branching is exclusively determined by the initial collisional orientation, which governs the formation of an ion-dipole complex, H3+···CO or H3+···OC, that dissociates into H2 + HCO+ or H2 + HOC+, respectively, via a direct mechanism. However, the contribution of the indirect mechanism via the rearrangement between H3+···CO and H3+···OC increases as the temperature increases, although its absolute fraction is small.
Collapse
Affiliation(s)
- Kohei Saito
- Department of Chemistry, Saitama University, Shimo-Okubo 255, Sakura-ku, Saitama 338-8570, Japan
| | - Yu Hashimoto
- Department of Chemistry, Saitama University, Shimo-Okubo 255, Sakura-ku, Saitama 338-8570, Japan
| | - Toshiyuki Takayanagi
- Department of Chemistry, Saitama University, Shimo-Okubo 255, Sakura-ku, Saitama 338-8570, Japan
| |
Collapse
|