1
|
Zhao H, Sun Z. Higher-Order Split Operator Schemes for Solving Tetratomic Reactions Using the Time-Dependent Wave Packet Method. J Phys Chem A 2024; 128:4911-4922. [PMID: 38847623 DOI: 10.1021/acs.jpca.4c01802] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/21/2024]
Abstract
In this work, using the time-dependent quantum wave packet method, quite a few typical higher-order split operators (HOSOs) were for the first time applied to calculate the tetratomic reactive scattering processes in the hyperspherical coordinate. It was found that the HOSOs were hardly efficient for a tetratomic reaction calculation, unlike those for a triatomic reactive scattering calculation. We proposed an efficient HOSO with a force gradient (denoted as 2G1 in the main text) for efficiently and accurately calculating a tetratomic reaction using the quantum wave packet method. Several typical tetratomic reactions, such as H2 + OH, HF + OH, and H2 + OH+, are calculated for demonstrating the effectiveness of the proposed 2G1 in terms of (product state-resolved) reaction probability and inelastic probability, by comparing with the performance of the previously reported various HOSOs. We suggest that the 2G1 propagator could be applied to efficiently calculate a general tetratomic reaction.
Collapse
Affiliation(s)
- Hailin Zhao
- Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, P.R. China
| | - Zhigang Sun
- Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, P.R. China
| |
Collapse
|
2
|
Ellerbrock R, Johnson KG, Seritan S, Hoppe H, Zhang JH, Lenzen T, Weike T, Manthe U, Martínez TJ. QuTree: A tree tensor network package. J Chem Phys 2024; 160:112501. [PMID: 38497471 DOI: 10.1063/5.0180233] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/08/2023] [Accepted: 02/20/2024] [Indexed: 03/19/2024] Open
Abstract
We present QuTree, a C++ library for tree tensor network approaches. QuTree provides class structures for tensors, tensor trees, and related linear algebra functions that facilitate the fast development of tree tensor network approaches such as the multilayer multiconfigurational time-dependent Hartree approach or the density matrix renormalization group approach and its various extensions. We investigate the efficiency of relevant tensor and tensor network operations and show that the overhead for managing the network structure is negligible, even in cases with a million leaves and small tensors. QuTree focuses on providing simple, high-level routines while retaining easy access to the backend to facilitate novel developments. We demonstrate the capabilities of the package by computing the eigenstates of coupled harmonic oscillator Hamiltonians and performing random circuit simulations on a virtual quantum computer.
Collapse
Affiliation(s)
- Roman Ellerbrock
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, USA
- SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, USA
- Bielefeld University, Universitätsstraße 25, 33615 Bielefeld, Germany
| | - K Grace Johnson
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, USA
- SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, USA
| | - Stefan Seritan
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, USA
- SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, USA
| | - Hannes Hoppe
- Bielefeld University, Universitätsstraße 25, 33615 Bielefeld, Germany
| | - J H Zhang
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, USA
- SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, USA
| | - Tim Lenzen
- Bielefeld University, Universitätsstraße 25, 33615 Bielefeld, Germany
| | - Thomas Weike
- Bielefeld University, Universitätsstraße 25, 33615 Bielefeld, Germany
| | - Uwe Manthe
- Bielefeld University, Universitätsstraße 25, 33615 Bielefeld, Germany
| | - Todd J Martínez
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, USA
- SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, USA
| |
Collapse
|
3
|
Hoppe H, Manthe U. Eigenstate calculation in the state-averaged (multi-layer) multi-configurational time-dependent Hartree approach. J Chem Phys 2024; 160:034104. [PMID: 38230812 DOI: 10.1063/5.0188748] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2023] [Accepted: 12/15/2023] [Indexed: 01/18/2024] Open
Abstract
A new approach for the calculation of eigenstates with the state-averaged (multi-layer) multi-configurational time-dependent Hartree (MCTDH) approach is presented. The approach is inspired by the recent work of Larsson [J. Chem. Phys. 151, 204102 (2019)]. It employs local optimization of the basis sets at each node of the multi-layer MCTDH tree and successive downward and upward sweeps to obtain a globally converged result. At the top node, the Hamiltonian represented in the basis of the single-particle functions (SPFs) of the first layer is diagonalized. Here p wavefunctions corresponding to the p lowest eigenvalues are computed by a block Lanczos approach. At all other nodes, a non-linear operator consisting of the respective mean-field Hamiltonian matrix and a projector onto the space spanned by the respective SPFs is considered. Here, the eigenstate corresponding to the lowest eigenvalue is computed using a short iterative Lanczos scheme. Two different examples are studied to illustrate the new approach: the calculation of the vibrational states of methyl and acetonitrile. The calculations for methyl employ the single-layer MCTDH approach, a general potential energy surface, and the correlation discrete variable representation. A five-layer MCTDH representation and a sum of product-type Hamiltonian are used in the acetonitrile calculations. Very fast convergence and order of magnitude reductions in the numerical effort compared to the previously used block relaxation scheme are found. Furthermore, a detailed comparison with the results of Avila and Carrington [J. Chem. Phys. 134, 054126 (2011)] for acetonitrile highlights the potential problems of convergence tests for high-dimensional systems.
Collapse
Affiliation(s)
- Hannes Hoppe
- 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
|
4
|
Mondal S, Pan H, Liu K. Stretching-mode specificity in the Cl + CH 3D( v1-I, v1-II, and v4 = 1; | jK〉) reactions: dependency on the initial | jK〉 selectivity. Phys Chem Chem Phys 2022; 24:24050-24061. [PMID: 36168830 DOI: 10.1039/d2cp03614b] [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 title reactions were studied at a collisional energy of 5.4 kcal mol-1 in a crossed-beam product-imaging experiment. The dynamics attributes of the dominant ground-state CH2D(00) and the accompanied C-D bend-excited CH2D(61) products were imaged in reactions with totally 16 ro-vibrationally selected states of the CH3D(vi, |jK〉) reagents. We found that all three vibrational excitations yielded marked |jK〉-dependent rate-enhancements in forming the (00, 0/1)s product pairs. Furthermore, for a given rotational |jK〉-mode, a vibrational-mode propensity of v4 > v1-I > v1-II in rate promotion and a clear manifestation of the Fermi-phase-induced interference effect of the latter two were observed. Compared to the reactivity of the rotationless state |jK〉 = |00〉, a minute rotational-excitation of all three stretch-excited CH3D(vi = 1) reagents could yield significantly higher reaction rates for the product pair (00, 0)s, but not so for (00, 1)s. The signals in forming the (61, 0)s pair were clearly notable but smaller than that of the ground-state reaction product pair, (00, 0)g. An opposite propensity of v1-II ≈ v1-I > v4 with a milder dependency on the initial |jK〉-states was observed. The angular distributions of the (00, 0)s pairs were nearly identical for all ro-vibrationally excited reagents, displaying the typical trait for a direct abstraction of the rebound mechanism. Similar distributions were found for the (61, 0)s pairs; yet, both pairs deviated substantially from the peripheral feature of the ground-state reaction pair of (00, 0)g. Those of the (00, 1)s pairs in reactions with v4-excitation featured a prominent forward-peaking distribution-suggestive of a time-delayed, resonance-mediated pathway, again with little dependency on the initial |jK〉-states. As for the reactions with the two Fermi-dyads, v1-I and v1-II, albeit showing globally similar distributions to that for v4, a substantial variation with the initial rotational-mode excitation could be discerned in the forward-peaking features. To unravel the impact of the Fermi-phase on the |jK〉-dependent attributes, we adopted a comparative approach by contrasting the observations in reactions with the Fermi-dyad reagents (the superposition states) to those with the pure-state reagents. Remarkable distinctions are unveiled and elucidated with the unexplained results explicitly pointed out, which call for future theoretical investigations for deeper understanding.
Collapse
Affiliation(s)
- Sohidul Mondal
- Institute of Atomic and Molecular Sciences (IAMS), Academia Sinica, P. O. Box 23-166, Taipei, 10699, Taiwan.
| | - Huilin Pan
- Institute of Atomic and Molecular Sciences (IAMS), Academia Sinica, P. O. Box 23-166, Taipei, 10699, Taiwan. .,Department of Chemistry, Southern University of Science and Technology, Shenzhen, 518055, P. R. China
| | - Kopin Liu
- Institute of Atomic and Molecular Sciences (IAMS), Academia Sinica, P. O. Box 23-166, Taipei, 10699, Taiwan. .,Aerosol Science Research Center, National Sun Yat-sen University, Kaohsiung, 80424, Taiwan.,State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, CAS, Dalian, 116023, P. R. China
| |
Collapse
|
5
|
Han S, Schröder M, Gatti F, Meyer HD, Lauvergnat D, Yarkony DR, Guo H. Representation of Diabatic Potential Energy Matrices for Multiconfiguration Time-Dependent Hartree Treatments of High-Dimensional Nonadiabatic Photodissociation Dynamics. J Chem Theory Comput 2022; 18:4627-4638. [PMID: 35839299 DOI: 10.1021/acs.jctc.2c00370] [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
Conventional quantum mechanical characterization of photodissociation dynamics is restricted by steep scaling laws with respect to the dimensionality of the system. In this work, we examine the applicability of the multi-configurational time-dependent Hartree (MCTDH) method in treating nonadiabatic photodissociation dynamics in two prototypical systems, taking advantage of its favorable scaling laws. To conform to the sum-of-product form, elements of the ab initio diabatic potential energy matrix (DPEM) are re-expressed using the recently proposed Monte Carlo canonical polyadic decomposition method, with enforcement of proper symmetry. The MCTDH absorption spectra and product branching ratios are shown to compare well with those calculated using conventional grid-based methods, demonstrating its promise for treating high-dimensional nonadiabatic photodissociation problems.
Collapse
Affiliation(s)
- Shanyu Han
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| | - Markus Schröder
- Theoretische Chemie, Physikalisch Chemisches Institut, Ruprecht-Karls Universität Heidelberg, D-69120 Heidelberg, Germany
| | - Fabien Gatti
- ISMO, Institut des Sciences Moléculaires d'Orsay─UMR 8214 CNRS/Université Paris-Saclay, F-91405 Orsay, France
| | - Hans-Dieter Meyer
- Theoretische Chemie, Physikalisch Chemisches Institut, Ruprecht-Karls Universität Heidelberg, D-69120 Heidelberg, Germany
| | - David Lauvergnat
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, Orsay 91405, France
| | - 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
|
6
|
Ellerbrock R, Manthe U. A non-hierarchical correlation discrete variable representation. J Chem Phys 2022; 156:134107. [PMID: 35395891 DOI: 10.1063/5.0088509] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
The correlation discrete variable representation (CDVR) facilitates (multi-layer) multi-configurational time-dependent Hartree (MCTDH) calculations with general potentials. It employs a layered grid representation to efficiently evaluate all potential matrix elements appearing in the MCTDH equations of motion. The original CDVR approach and its multi-layer extension show a hierarchical structure: the size of the grids employed at the different layers increases when moving from an upper layer to a lower one. In this work, a non-hierarchical CDVR approach, which uses identically structured quadratures at all layers of the MCTDH wavefunction representation, is introduced. The non-hierarchical CDVR approach crucially reduces the number of grid points required, compared to the hierarchical CDVR, shows superior scaling properties, and yields identical results for all three representations showing the same topology. Numerical tests studying the photodissociation of NOCl and the vibrational states of CH3 demonstrate the accuracy of the non-hierarchical CDVR approach.
Collapse
Affiliation(s)
- Roman Ellerbrock
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, USA
| | - Uwe Manthe
- Theoretische Chemie, Fakultät für Chemie, Universität Bielefeld, Universitätsstr. 25, D-33615 Bielefeld, Germany
| |
Collapse
|
7
|
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
|
8
|
Marjollet A, Inhester L, Welsch R. Initial state-selected scattering for the reactions H + CH4/CHD3 and F + CHD3 employing ring polymer molecular dynamics. J Chem Phys 2022; 156:044101. [DOI: 10.1063/5.0076216] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- A. Marjollet
- Center for Free-Electron Laser Science CFEL, Deutsches Elektronen-Synchrotron DESY, Notkestrasse 85, 22607 Hamburg, Germany
- Department of Physics, Universität Hamburg, Notkestr. 9-11, 22607 Hamburg, Germany
| | - L. Inhester
- Center for Free-Electron Laser Science CFEL, Deutsches Elektronen-Synchrotron DESY, Notkestrasse 85, 22607 Hamburg, Germany
| | - R. Welsch
- Center for Free-Electron Laser Science CFEL, Deutsches Elektronen-Synchrotron DESY, Notkestrasse 85, 22607 Hamburg, Germany
| |
Collapse
|
9
|
A Ten-Dimensional Quantum Dynamics Study of the H+CH 3D→H 2+CH 2D Reaction. CHINESE J CHEM PHYS 2022. [DOI: 10.1063/1674-0068/cjcp2112284] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022]
|
10
|
Weike T, Manthe U. The multi-configurational time-dependent Hartree approach in optimized second quantization: thermal ensembles and statistical sampling. Chem Phys 2021. [DOI: 10.1016/j.chemphys.2021.111413] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
|
11
|
Zhao B. The symmetric C-D stretching spectator mode in the H + CHD 3 → H 2 + CD 3 reaction and its effect on dynamical modeling. Phys Chem Chem Phys 2021; 23:12105-12114. [PMID: 34027536 DOI: 10.1039/d1cp01614h] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
The symmetric C-D stretching mode is a spectator mode in the H + CHD3 → H2 + CD3 reaction. Effects of multiple vibrational excitations of the CHD3 reactant are studied with the quantum transition-state (QTS) framework and an eight-dimensional (8D) model Hamiltonian developed by Palma and Clary. By including many thermal flux eigenstates, results have been obtained up to high energies, allowing the study of the symmetric C-D stretching spectator mode. A new concept of a state-specific thermal flux operator is proposed to analyze the C-D stretching spectator mode in detail, providing a new and insightful venue for studying transition-state control of chemical reactions. Furthermore, as a spectator mode, whether the C-D stretching motion can be excluded in a seven-dimensional (7D) model has not been fully interrogated, although the 7D model is a reasonable approximation and has provided accurate theoretical predictions. By comparing with available results of full-dimensional calculations, both the 7D and 8D models predict reasonably accurate results. However, the 7D model underestimates the mixing of two vibrational states that are in Fermi resonance. Despite its spectator nature, the C-D stretch is important in the dynamical modeling of chemical reaction systems affected by state mixing.
Collapse
Affiliation(s)
- Bin Zhao
- Theoretische Chemie, Fakultät für Chemie, Universität Bielefeld, Universitätsstraße 25, D-33615 Bielefeld, Germany.
| |
Collapse
|
12
|
Weike T, Manthe U. Symmetries in the multi-configurational time-dependent Hartree wavefunction representation and propagation. J Chem Phys 2021; 154:194108. [PMID: 34240912 DOI: 10.1063/5.0054105] [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/14/2022] Open
Abstract
In multi-configurational time-dependent Hartree (MCTDH) approaches, different multi-layered wavefunction representations can be used to represent the same physical wavefunction. Transformations between different equivalent representations of a physical wavefunction that alter the tree structure used in the multi-layer MCTDH wavefunction representation interchange the role of single-particle functions (SPFs) and single-hole functions (SHFs) in the MCTDH formalism. While the physical wavefunction is invariant under these transformations, this invariance does not hold for the standard multi-layer MCTDH equations of motion. Introducing transformed SPFs, which obey normalization conditions typically associated with SHFs, revised equations of motion are derived. These equations do not show the singularities resulting from the inverse single-particle density matrix and are invariant under tree transformations. Based on the revised equations of motion, a new integration scheme is introduced. The scheme combines the advantages of the constant mean-field approach of Beck and Meyer [Z. Phys. D 42, 113 (1997)] and the singularity-free integrator suggested by Lubich [Appl. Math. Res. Express 2015, 311]. Numerical calculations studying the spin boson model in high dimensionality confirm the favorable properties of the new integration scheme.
Collapse
Affiliation(s)
- Thomas Weike
- 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
|
13
|
Zhao B, Manthe U. Direct product-type grid representations for angular coordinates in extended space and their application in the MCTDH approach. J Chem Phys 2021; 154:104115. [PMID: 33722051 DOI: 10.1063/5.0045054] [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/14/2022] Open
Abstract
Multi-configurational time-dependent Hartree (MCTDH) calculations using time-dependent grid representations can be used to accurately simulate high-dimensional quantum dynamics on general ab initio potential energy surfaces. Employing the correlation discrete variable representation, sets of direct product type grids are employed in the calculation of the required potential energy matrix elements. This direct product structure can be a problem if the coordinate system includes polar and azimuthal angles that result in singularities in the kinetic energy operator. In the present work, a new direct product-type discrete variable representation (DVR) for arbitrary sets of polar and azimuthal angles is introduced. It employs an extended coordinate space where the range of the polar angles is taken to be [-π, π]. The resulting extended space DVR resolves problems caused by the singularities in the kinetic energy operator without generating a very large spectral width. MCTDH calculations studying the F·CH4 complex are used to investigate important properties of the new scheme. The scheme is found to allow for more efficient integration of the equations of motion compared to the previously employed cot-DVR approach [G. Schiffel and U. Manthe, Chem. Phys. 374, 118 (2010)] and decreases the required central processing unit times by about an order of magnitude.
Collapse
Affiliation(s)
- 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
|
14
|
Zhang X, Zhang Z, Gatti F, Zhang DH. Full-dimensional quantum dynamics study of isotope effects for the H 2 + NH 2/ND 2/NHD and H 2/D 2/HD + NH 2 reactions. J Chem Phys 2021; 154:074301. [PMID: 33607900 DOI: 10.1063/5.0040002] [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/2022] Open
Abstract
A full-dimensional quantum dynamical study for the bimolecular reactions of hydrogen molecules with amino radicals for different isotopologues is reported. The nonreactive amino radical is described by two Radau vectors that are very close to the valence bond coordinates. Potential-optimized discrete variable representation basis is used for the vibrational coordinates of the amino radical. Starting from the reaction H2 + NH2, we study the isotope effects for the two reagents separately, i.e., H2 + NH2/ND2/NHD and H2/D2/HD + NH2. The effects of different vibrational mode excitations of the reagents on the reactivities are studied. Physical explanations about the isotope effects are also provided thoroughly including the influence of vibrational energy differences between the different isotopologues and the impact of the tunneling effect.
Collapse
Affiliation(s)
- Xiaoren 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
| | - 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-Saclay, 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
|
15
|
Zhao B, Manthe U. Eight-Dimensional Wave Packet Dynamics Within the Quantum Transition-State Framework: State-to-State Reactive Scattering for H2 + CH3 ⇆ H + CH4. J Phys Chem A 2020; 124:9400-9412. [DOI: 10.1021/acs.jpca.0c08461] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Affiliation(s)
- Bin Zhao
- Theoretische Chemie, Fakultät für Chemie, Universität Bielefeld, Universitätsstraße 25, D-33615 Bielefeld, Germany
| | - Uwe Manthe
- Theoretische Chemie, Fakultät für Chemie, Universität Bielefeld, Universitätsstraße 25, D-33615 Bielefeld, Germany
| |
Collapse
|
16
|
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
|
17
|
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
|
18
|
Weike T, Manthe U. The multi-configurational time-dependent Hartree approach in optimized second quantization: Imaginary time propagation and particle number conservation. J Chem Phys 2020; 152:034101. [DOI: 10.1063/1.5140984] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Thomas Weike
- 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
|
19
|
Ellerbrock R, Manthe U, Palma J. A Quasi-Classical Evaluation of the J-Shifting Approximation for the Reactive Cross Sections of F + CHD 3 and F + CH 4. J Phys Chem A 2019; 123:7237-7245. [PMID: 31361132 DOI: 10.1021/acs.jpca.9b06060] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We evaluated the accuracy of the J-shifting approximation to estimate reactant state-selected cross sections for the F+CH4 → HF+CH3 and F+CHD3 → HF+CD3/DF+CHD2 reactions. In particular, we analyzed how the rotational state of methane influences the quality of the approximation. The systems were considered in full dimensionality. Since full-quantum scattering calculations are still unfeasible for these reactions, we employed quasi-classical trajectories (QCT) to calculate the cross sections. The characteristics of the Born-Oppenheimer potential energy surface of these reactions pose a great challenge to the assumptions of the J-shifting approach. In spite of this, we found that it performs well for both reactions if the methane molecule is in the rotational ground state. However, when methane is rotationally excited, the approach affords good results for the F+CH4 system but clearly fails for F+CHD3. The reasons for this failure will be discussed, and a simple procedure to recover good estimators for the cross sections from J = 0 calculations will be introduced.
Collapse
Affiliation(s)
- Roman Ellerbrock
- 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
| | - Juliana Palma
- Departamento de Ciencia y Tecnología, Universidad Nacional de Quilmes, Sáenz Peña 352, Bernal B1876BXD , Consejo Nacional de Investigaciones Científicas y Técnicas , Buenos Aires , Argentina
| |
Collapse
|
20
|
Lenzen T, Eisfeld W, Manthe U. Vibronically and spin-orbit coupled diabatic potentials for X(2P) + CH4→ HX + CH3reactions: Neural network potentials for X = Cl. J Chem Phys 2019; 150:244115. [DOI: 10.1063/1.5109877] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Tim Lenzen
- Theoretische Chemie, Fakultät für Chemie, Universität Bielefeld, Universitätsstraße 25, D-33615 Bielefeld, Germany
| | - Wolfgang Eisfeld
- Theoretische Chemie, Fakultät für Chemie, Universität Bielefeld, Universitätsstraße 25, D-33615 Bielefeld, Germany
| | - Uwe Manthe
- Theoretische Chemie, Fakultät für Chemie, Universität Bielefeld, Universitätsstraße 25, D-33615 Bielefeld, Germany
| |
Collapse
|
21
|
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
|
22
|
Zhao B, Manthe U. Counter-propagating wave packets in the quantum transition state approach to reactive scattering. J Chem Phys 2019; 150:184103. [DOI: 10.1063/1.5097997] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- 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
|
23
|
Lu X, Wang X, Fu B, Zhang D. Theoretical Investigations of Rate Coefficients of H + H2O2 → OH + H2O on a Full-Dimensional Potential Energy Surface. J Phys Chem A 2019; 123:3969-3976. [DOI: 10.1021/acs.jpca.9b02526] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Xiaoxiao Lu
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Zhongshan Road 457, Dalian 116023, China
- Department of Chemical Physics, University of Science and Technology of China, Jinzhai Road 96, Hefei 230026, China
| | - Xingan Wang
- Department of Chemical Physics, University of Science and Technology of China, Jinzhai Road 96, Hefei 230026, China
| | - Bina Fu
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Zhongshan Road 457, Dalian 116023, China
| | - Donghui Zhang
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Zhongshan Road 457, Dalian 116023, China
| |
Collapse
|
24
|
Jasper AW, Sivaramakrishnan R, Klippenstein SJ. Nonthermal rate constants for CH 4 * + X → CH 3 + HX, X = H, O, OH, and O 2. J Chem Phys 2019; 150:114112. [PMID: 30902010 DOI: 10.1063/1.5090394] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023] Open
Abstract
Quasiclassical trajectories are used to compute nonthermal rate constants, k*, for abstraction reactions involving highly-excited methane CH4 * and the radicals H, O, OH, and O2. Several temperatures and internal energies of methane, Evib, are considered, and significant nonthermal rate enhancements for large Evib are found. Specifically, when CH4 * is internally excited close to its dissociation threshold (Evib ≈ D0 = 104 kcal/mol), its reactivity with H, O, and OH is shown to be collision-rate-limited and to approach that of comparably-sized radicals, such as CH3, with k* > 10-10 cm3 molecule-1 s-1. Rate constants this large are more typically associated with barrierless reactions, and at 1000 K, this represents a nonthermal rate enhancement, k*/k, of more than two orders of magnitude relative to thermal rate constants k. We show that large nonthermal rate constants persist even after significant internal cooling, with k*/k > 10 down to Evib ≈ D0/4. The competition between collisional cooling and nonthermal reactivity is studied using a simple model, and nonthermal reactions are shown to account for up to 35%-50% of the fate of the products of H + CH3 = CH4 * under conditions of practical relevance to combustion. Finally, the accuracy of an effective temperature model for estimating k* from k is quantified.
Collapse
Affiliation(s)
- Ahren W Jasper
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Argonne, Illinois 60439, USA
| | - Raghu Sivaramakrishnan
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Argonne, Illinois 60439, USA
| | - Stephen J Klippenstein
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Argonne, Illinois 60439, USA
| |
Collapse
|
25
|
Aieta C, Gabas F, Ceotto M. Parallel Implementation of Semiclassical Transition State Theory. J Chem Theory Comput 2019; 15:2142-2153. [PMID: 30822385 DOI: 10.1021/acs.jctc.8b01286] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
This paper presents the parsctst code, an efficient parallel implementation of the semiclassical transition state theory (SCTST) for reaction rate constant calculations. Parsctst is developed starting from a previously presented approach for the computation of the vibrational density of states of fully coupled anharmonic molecules ( Nguyen et al. Chem. Phys. Lett. 2010 , 499 , 915 ). The parallel implementation makes it practical to tackle reactions involving more than 100 fully coupled anharmonic vibrational degrees of freedom and also includes multidimensional tunneling effects. After describing the pseudocode and demonstrating its computational efficiency, we apply the new code for estimating the rate constant of the proton transfer isomerization reaction of the 2,4,6-tri- tert-butylphenyl to 3,5-di- tert-butylneophyl. Comparison with both theoretical and experimental results is presented. Parsctst code is user-friendly and provides a significant computational time saving compared to serial calculations. We believe that parsctst can boost the application of SCTST as an alternative to the basic transition state theory for accurate kinetics modeling not only in combustion or atmospheric chemistry, but also in organic synthesis, where bigger reactive systems are encountered.
Collapse
Affiliation(s)
- Chiara Aieta
- Dipartimento di Chimica , Università degli Studi di Milano , via C. Golgi 19 , 20133 Milano , Italy
| | - Fabio Gabas
- Dipartimento di Chimica , Università degli Studi di Milano , via C. Golgi 19 , 20133 Milano , Italy
| | - Michele Ceotto
- Dipartimento di Chimica , Università degli Studi di Milano , via C. Golgi 19 , 20133 Milano , Italy
| |
Collapse
|
26
|
Lenzen T, Manthe U. Vibronically and spin-orbit coupled diabatic potentials for X(P) + CH4→ HX + CH3reactions: General theory and application for X(P) = F(2P). J Chem Phys 2019; 150:064102. [DOI: 10.1063/1.5063907] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Tim Lenzen
- Theoretische Chemie, Fakultät für Chemie, Universität Bielefeld, Universitätsstraße 25, D-33615 Bielefeld, Germany
| | - Uwe Manthe
- Theoretische Chemie, Fakultät für Chemie, Universität Bielefeld, Universitätsstraße 25, D-33615 Bielefeld, Germany
| |
Collapse
|
27
|
Li J, Song K, Behler J. A critical comparison of neural network potentials for molecular reaction dynamics with exact permutation symmetry. Phys Chem Chem Phys 2019; 21:9672-9682. [DOI: 10.1039/c8cp06919k] [Citation(s) in RCA: 32] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Several symmetry strategies have been compared in fitting full dimensional accurate potentials for reactive systems based on a neural network approach.
Collapse
Affiliation(s)
- Jun Li
- School of Chemistry and Chemical Engineering, Chongqing University
- Chongqing 401331
- China
- Universität Göttingen, Institut für Physikalische Chemie, Theoretische Chemie
- 37077 Göttingen
| | - Kaisheng Song
- School of Chemistry and Chemical Engineering, Chongqing University
- Chongqing 401331
- China
| | - Jörg Behler
- Universität Göttingen, Institut für Physikalische Chemie, Theoretische Chemie
- 37077 Göttingen
- Germany
| |
Collapse
|
28
|
Welsch R. Kinetic isotope effects in the water forming reaction H2/D2 + OH from rigorous close-coupling quantum dynamics simulations. Phys Chem Chem Phys 2019; 21:17054-17062. [DOI: 10.1039/c9cp02323b] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Rigorous quantum dynamics simulations of thermal rate constants and kinetic isotope effects for the water-forming H2/D2 + OH reaction are presented, which show increased tunneling below 300 K and can serve as benchmarks for approximate methods.
Collapse
Affiliation(s)
- Ralph Welsch
- Center for Free-Electron Laser Science
- DESY
- 22607 Hamburg
- Germany
- The Hamburg Centre for Ultrafast Imaging
| |
Collapse
|
29
|
Lu X, Meng Q, Wang X, Fu B, Zhang DH. Rate coefficients of the H + H2O2→ H2+ HO2reaction on an accurate fundamental invariant-neural network potential energy surface. J Chem Phys 2018; 149:174303. [PMID: 30409010 DOI: 10.1063/1.5063613] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Xiaoxiao Lu
- Department of Chemical Physics, University of Science and Technology of China, Jinzhai Road 96, Hefei 230026, China
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Zhongshan Road 457, Dalian 116023, China
| | - Qingyong Meng
- Department of Applied Chemistry, Northwestern Polytechnical University, Youyi West Road 127, Xi’an 710072, China
| | - Xingan Wang
- Department of Chemical Physics, University of Science and Technology of China, Jinzhai Road 96, Hefei 230026, China
| | - Bina Fu
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Zhongshan Road 457, Dalian 116023, China
| | - 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, Zhongshan Road 457, Dalian 116023, China
| |
Collapse
|
30
|
Welsch R. Low‐Temperature Thermal Rate Constants for the Water Formation Reaction H
2
+OH from Rigorous Quantum Dynamics Calculations. Angew Chem Int Ed Engl 2018; 57:13150-13153. [DOI: 10.1002/anie.201807666] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/05/2018] [Indexed: 11/06/2022]
Affiliation(s)
- Ralph Welsch
- Center for Free-Electron Laser Science, DESY Notkestraße 85 22607 Hamburg Germany
| |
Collapse
|
31
|
Welsch R. Low‐Temperature Thermal Rate Constants for the Water Formation Reaction H
2
+OH from Rigorous Quantum Dynamics Calculations. Angew Chem Int Ed Engl 2018. [DOI: 10.1002/ange.201807666] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Affiliation(s)
- Ralph Welsch
- Center for Free-Electron Laser Science, DESY Notkestraße 85 22607 Hamburg Germany
| |
Collapse
|