1
|
Wang J, Xie C, Hu X, Guo H, Xie D. Impact of Geometric Phase on Dynamics of Complex-Forming Reactions: H + O 2 → OH + O. J Phys Chem Lett 2024; 15:4237-4243. [PMID: 38602563 DOI: 10.1021/acs.jpclett.4c00789] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/12/2024]
Abstract
Reaction dynamics on the ground electronic state might be significantly influenced by conical intersections (CIs) via the geometric phase (GP), as demonstrated for activated reactions (i.e., the H + H2 exchange reaction). However, there have been few investigations of GP effects in complex-forming reactions. Here, we report a full quantum dynamical study of an important reaction in combustion (H + O2 → OH + O), which serves as a proving ground for studying GP effects therein. The results reveal significant differences in reaction probabilities and differential cross sections (DCSs) obtained with and without GP, underscoring its strong impact. However, the GP effects are less pronounced for the reaction integral cross sections, apparently due to the integral of the DCS over the scattering angle. Further analysis indicated that the cross section has roughly the same contributions from the two topologically distinct paths around the CI, namely, the direct and looping paths.
Collapse
Affiliation(s)
- Junyan Wang
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, Jiangsu 210023, China
| | - Changjian Xie
- Institute of Modern Physics, Shaanxi Key Laboratory for Theoretical Physics Frontiers, Northwest University, Xi'an, Shaanxi 710127, 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, Jiangsu 210023, China
- Hefei National Laboratory, Hefei, Anhui 230088, China
| | - Hua Guo
- Department of Chemistry and Chemical Biology, Center for Computational Chemistry, 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, Jiangsu 210023, China
- Hefei National Laboratory, Hefei, Anhui 230088, China
| |
Collapse
|
2
|
Wang J, An F, Chen J, Hu X, Guo H, Xie D. Accurate Full-Dimensional Global Diabatic Potential Energy Matrix for the Two Lowest-Lying Electronic States of the H + O 2 ↔ HO + O Reaction. J Chem Theory Comput 2023; 19:2929-2938. [PMID: 37161259 DOI: 10.1021/acs.jctc.3c00291] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/11/2023]
Abstract
A new and more accurate diabatic potential energy matrix (DPEM) is developed for the two lowest-lying electronic states of HO2, covering both the strong interaction region and reaction asymptotes. The ab initio calculations were performed at the Davidson corrected multireference configuration interaction level with the augmented correlation-consistent polarized valence quintuple-zeta basis set (MRCI+Q/AV5Z). The accuracy of the electronic structure calculations is validated by excellent agreement with the experimental HO2 equilibrium geometry, fundamental vibrational frequencies, and H + O2 ↔ OH + O reaction energy. Through the combination of an electronic angular momentum-method and a configuration interaction vector-based method, the mixing angle between the first two 2A″ states of HO2 was successfully determined. Elements of the 2×2 DPEM were fit to neural networks with a proper account of the complete nuclear permutation inversion symmetry of HO2. The DPEM correctly predicted the properties of conical intersection seams at linear and T-shape geometries, thus providing a reliable platform for studying both the spectroscopy of HO2 and the nonadiabatic dynamics for the H + O2 ↔ OH + O reaction.
Collapse
Affiliation(s)
- Junyan Wang
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Feng An
- Research Center for Graph Computing, Zhejiang Lab, Hangzhou 311121, China
| | - Junjie 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, Jiangsu Key Laboratory of Vehicle Emissions Control, Center of Modern Analysis, Nanjing University, Nanjing 210023, China
- Hefei National Laboratory, Hefei 230088, 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
- Hefei National Laboratory, Hefei 230088, China
| |
Collapse
|
3
|
Kendrick BK. Quantum reactive scattering calculations for the cold and ultracold Li + LiNa → Li 2 + Na reaction. J Chem Phys 2021; 154:124303. [PMID: 33810695 DOI: 10.1063/5.0045712] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
A first-principles based quantum dynamics study of the Li + LiNa(v = 0, j = 0) → Li2(v', j') + Na reaction is reported for collision energies spanning the ultracold (1 nK) to cold (1 K) regimes. A full-dimensional ab initio potential energy surface for the ground electronic state of Li2Na is utilized that includes an accurate treatment of the long-range interactions. The Li + LiNa reaction is barrierless and exoergic and exhibits a deep attractive potential well that supports complex formation. Thus, significant reactivity occurs even for collision temperatures approaching absolute zero. The reactive scattering calculations are based on a numerically exact time-independent quantum dynamics methodology in hyperspherical coordinates. Total and rotationally resolved rate coefficients are reported at 56 collision energies and include all contributing partial waves. Several shape resonances are observed in many of the rotationally resolved rate coefficients and a small resonance feature is also reported in the total rate coefficient near 50 mK. Of particular interest, the angular distributions or differential cross sections are reported as a function of both the collision energy and scattering angle. Unique quantum fingerprints (bumps, channels, and ripples) are observed in the angular distributions for each product rotational state due to quantum interference and shape resonance contributions. The Li + LiNa reaction is under active experimental investigation so that these intriguing features could be verified experimentally when sufficient product state resolution becomes feasible for collision energies below 1 K.
Collapse
Affiliation(s)
- Brian K Kendrick
- Theoretical Division (T-1, MS B221), Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| |
Collapse
|
4
|
Kendrick BK, Li H, Li M, Kotochigova S, Croft JFE, Balakrishnan N. Non-adiabatic quantum interference in the ultracold Li + LiNa → Li 2 + Na reaction. Phys Chem Chem Phys 2021; 23:5096-5112. [PMID: 33576359 DOI: 10.1039/d0cp05499b] [Citation(s) in RCA: 17] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
Electronically non-adiabatic effects play an important role in many chemical reactions. However, how these effects manifest in cold and ultracold chemistry remains largely unexplored. Here for the first time we present from first principles the non-adiabatic quantum dynamics of the reactive scattering between ultracold alkali-metal LiNa molecules and Li atoms. We show that non-adiabatic dynamics induces quantum interference effects that dramatically alter the ultracold rotationally resolved reaction rate coefficients. The interference effect arises from the conical intersection between the ground and an excited electronic state that is energetically accessible even for ultracold collisions. These unique interference effects might be exploited for quantum control applications such as a quantum molecular switch. The non-adiabatic dynamics are based on full-dimensional ab initio potential energy surfaces for the two electronic states that includes the non-adiabatic couplings and an accurate treatment of the long-range interactions. A statistical analysis of rotational populations of the Li2 product reveals a Poisson distribution implying the underlying classical dynamics are chaotic. The Poisson distribution is robust and amenable to experimental verification and appears to be a universal property of ultracold reactions involving alkali metal dimers.
Collapse
Affiliation(s)
- Brian K Kendrick
- Theoretical Division (T-1, MS B221), Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA.
| | - Hui Li
- Department of Physics, Temple University, Philadelphia, PA 19122, USA
| | - Ming Li
- Department of Physics, Temple University, Philadelphia, PA 19122, USA
| | | | - James F E Croft
- Dodd-Walls Centre for Photonic and Quantum Technologies, Dunedin 9054, New Zealand and Department of Physics, University of Otago, Dunedin 9054, New Zealand
| | | |
Collapse
|
5
|
Huang J, Kendrick BK, Zhang DH. Mechanistic Insights into Ultracold Chemical Reactions under the Control of the Geometric Phase. J Phys Chem Lett 2021; 12:2160-2165. [PMID: 33626281 DOI: 10.1021/acs.jpclett.1c00133] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
Ultracold chemical reactions involve collision temperatures approaching absolute zero, and for molecular systems that exhibit a barrierless and exoergic reaction path significant reactivity can occur. In addition, many molecules contain a conical intersection, and the associated geometric phase has been shown to significantly alter the outcome of ultracold reactions. Here we report a quantum dynamics study for the ultracold O + OH → H + O2 reaction. An analysis of the scattering wave functions reveals explicitly the nature of the quantum interference between the direct and looping reaction pathways around the conical intersection and thus illustrates how the reaction proceeds under the control of the geometric phase for the first time. The wave function analysis should generalize to other ultracold reactions that contain a conical intersection. Our findings indicate that quantum control techniques such as an optical lattice trap or the initial state orientation may be effective in controlling the reactivity.
Collapse
Affiliation(s)
- Jiayu Huang
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China
| | - Brian K Kendrick
- Theoretical Division (T-1, MS B221), Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Dong H Zhang
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China
| |
Collapse
|
6
|
Abstract
The prospect of cooling matter down to temperatures that are close to absolute zero raises intriguing questions about how chemical reactivity changes under these extreme conditions. Although some types of chemical reaction still occur at 1 μK, they can no longer adhere to the conventional picture of reactants passing over an activation energy barrier to become products. Indeed, at ultracold temperatures, the system enters a fully quantum regime, and quantum mechanics replaces the classical picture of colliding particles. In this Review, we discuss recent experimental and theoretical developments that allow us to explore chemical reactions at temperatures that range from 100 K to 500 nK. Although the field is still in its infancy, exceptional control has already been demonstrated over reactivity at low temperatures.
Collapse
|
7
|
Kendrick BK. Nonadiabatic Ultracold Quantum Reactive Scattering of Hydrogen with Vibrationally Excited HD( v = 5-9). J Phys Chem A 2019; 123:9919-9933. [PMID: 31647679 DOI: 10.1021/acs.jpca.9b07318] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The results from electronically non-adiabatic and adiabatic quantum reactive scattering calculations are presented for the H + HD(v = 5-9) → H + HD(v', j') reaction at ultracold collision energies from 10 nK to 60 K. Several experimentally verifiable signatures of the geometric phase are reported in the total and vibrationally and rotationally resolved rate coefficients. Most notable is the predicted 2 orders of magnitude enhancement of the rotationally resolved ultracold rates of odd symmetry relative to those of even symmetry. Prominent shape resonances appear at higher collision energies (100 mK to 20 K), which could be measured experimentally. Significant geometric phase effects are also reported on the resonance energies and lifetimes. In particular, an enhancement (suppression) of the l = 1 (l = 2) shape resonances for HD(v = 5, 6) is predicted for even symmetry relative to those of odd symmetry.
Collapse
Affiliation(s)
- Brian K Kendrick
- Theoretical Division , Los Alamos National Laboratory , Group T-1, Mail Stop B221, Los Alamos , New Mexico 87544 , United States
| |
Collapse
|
8
|
Abstract
The role of the geometric phase effect in chemical reaction dynamics has long been a topic of active experimental and theoretical investigations. The topic has received renewed interest in recent years in cold and ultracold chemistry where it was shown to play a decisive role in state-to-state chemical dynamics. We provide a brief review of these developments focusing on recent studies of O + OH and hydrogen exchange in the H + H 2 and D + HD reactions at cold and ultracold temperatures. Non-adiabatic effects in ultracold chemical dynamics arising from the conical intersection between two electronic potential energy surfaces are also briefly discussed. By taking the hydrogen exchange reaction as an illustrative example it is shown that the inclusion of the geometric phase effect captures the essential features of non-adiabatic dynamics at collision energies below the conical intersection.
Collapse
|
9
|
Kendrick BK. Non-adiabatic quantum reactive scattering calculations for the ultracold hydrogen exchange reaction: H + H2(v=4-8,j=0) → H + H2(v′,j′). Chem Phys 2018. [DOI: 10.1016/j.chemphys.2018.07.001] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
|
10
|
Zuo JX, Hu XX, Xie DQ. Quantum Dynamics of Oxyhydrogen Complex-Forming Reactions for the HO2 and HO3 Systems. CHINESE J CHEM PHYS 2018. [DOI: 10.1063/1674-0068/31/cjcp1804060] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022]
Affiliation(s)
- Jun-xiang Zuo
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210093, China
| | - Xi-xi Hu
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210093, China
| | - Dai-qian Xie
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210093, China
| |
Collapse
|
11
|
Kendrick BK. Non-adiabatic quantum reactive scattering in hyperspherical coordinates. J Chem Phys 2018; 148:044116. [DOI: 10.1063/1.5014989] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Brian K. Kendrick
- Theoretical Division (T-1, MS B221), Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| |
Collapse
|
12
|
Croft JFE, Makrides C, Li M, Petrov A, Kendrick BK, Balakrishnan N, Kotochigova S. Universality and chaoticity in ultracold K+KRb chemical reactions. Nat Commun 2017; 8:15897. [PMID: 28722014 PMCID: PMC5524979 DOI: 10.1038/ncomms15897] [Citation(s) in RCA: 34] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/01/2017] [Accepted: 05/10/2017] [Indexed: 11/17/2022] Open
Abstract
A fundamental question in the study of chemical reactions is how reactions proceed at a collision energy close to absolute zero. This question is no longer hypothetical: quantum degenerate gases of atoms and molecules can now be created at temperatures lower than a few tens of nanokelvin. Here we consider the benchmark ultracold reaction between, the most-celebrated ultracold molecule, KRb and K. We map out an accurate ab initio ground-state potential energy surface of the K2Rb complex in full dimensionality and report numerically-exact quantum-mechanical reaction dynamics. The distribution of rotationally resolved rates is shown to be Poissonian. An analysis of the hyperspherical adiabatic potential curves explains this statistical character revealing a chaotic distribution for the short-range collision complex that plays a key role in governing the reaction outcome. Studying chemical reactions near zero temperature in detail is challenging both in theory and practice. Here the authors report an explicit quantum mechanical study of the benchmark ultracold reaction between a K atom and a KRb molecule, important for future controlled chemistry experiments.
Collapse
Affiliation(s)
- J F E Croft
- Department of Chemistry, University of Nevada, Las Vegas, Nevada 89154, USA
| | - C Makrides
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, USA
| | - M Li
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, USA
| | - A Petrov
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, USA.,NRC 'Kurchatov Institute' PNPI, Gatchina, Leningrad District 188300 Russia.,Division of Quantum Mechanics, St Petersburg State University, 7/9 Universitetskaya nab., St Petersburg 199034, Russia
| | - B K Kendrick
- Theoretical Division (T-1, MS B221), Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| | - N Balakrishnan
- Department of Chemistry, University of Nevada, Las Vegas, Nevada 89154, USA
| | - S Kotochigova
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, USA
| |
Collapse
|
13
|
Wu Y, Zhang C, Ma H. Ab initio conical intersections for the Si( 1D) + H 2 reaction system: a lowest five singlet states study. RSC Adv 2017. [DOI: 10.1039/c7ra01021d] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022] Open
Abstract
Conical intersections and geometric phase effects of the Si(1D) + H2 system were clarified intuitively, and important features of them are revealed.
Collapse
Affiliation(s)
- Yanan Wu
- Beijing National Laboratory for Molecular Sciences
- Institute of Chemistry
- Chinese Academy of Sciences
- Beijing 100190
- China
| | - Chunfang Zhang
- Beijing National Laboratory for Molecular Sciences
- Institute of Chemistry
- Chinese Academy of Sciences
- Beijing 100190
- China
| | - Haitao Ma
- Beijing National Laboratory for Molecular Sciences
- Institute of Chemistry
- Chinese Academy of Sciences
- Beijing 100190
- China
| |
Collapse
|
14
|
Kendrick BK, Hazra J, Balakrishnan N. Geometric phase effects in the ultracold H + H2reaction. J Chem Phys 2016; 145:164303. [DOI: 10.1063/1.4966037] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- B. K. Kendrick
- Theoretical Division (T-1, MS B221), Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| | - Jisha Hazra
- Department of Chemistry, University of Nevada, Las Vegas, Nevada 89154, USA
| | - N. Balakrishnan
- Department of Chemistry, University of Nevada, Las Vegas, Nevada 89154, USA
| |
Collapse
|
15
|
Balakrishnan N. Perspective: Ultracold molecules and the dawn of cold controlled chemistry. J Chem Phys 2016; 145:150901. [DOI: 10.1063/1.4964096] [Citation(s) in RCA: 157] [Impact Index Per Article: 19.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- N. Balakrishnan
- Department of Chemistry, University of Nevada, Las Vegas, Nevada 89154, USA
| |
Collapse
|
16
|
Kendrick BK, Hazra J, Balakrishnan N. Geometric Phase Appears in the Ultracold Hydrogen Exchange Reaction. PHYSICAL REVIEW LETTERS 2015; 115:153201. [PMID: 26550721 DOI: 10.1103/physrevlett.115.153201] [Citation(s) in RCA: 41] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/23/2015] [Indexed: 06/05/2023]
Abstract
Quantum reactive scattering calculations for the hydrogen exchange reaction H+H_{2}(v=4,j=0)→H+H_{2}(v^{'}, j^{'}) and its isotopic analogues are reported for ultracold collision energies. Because of the unique properties associated with ultracold collisions, it is shown that the geometric phase effectively controls the reactivity. The rotationally resolved rate coefficients computed with and without the geometric phase are shown to differ by up to 4 orders of magnitude. The effect is also significant in the vibrationally resolved and total rate coefficients. The dynamical origin of the effect is discussed and the large geometric phase effect reported here might be exploited to control the reactivity through the application of external fields or by the selection of a particular nuclear spin state.
Collapse
Affiliation(s)
- B K Kendrick
- Theoretical Division (T-1, MS B221), Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| | - Jisha Hazra
- Department of Chemistry, University of Nevada, Las Vegas, Nevada 89154, USA
| | - N Balakrishnan
- Department of Chemistry, University of Nevada, Las Vegas, Nevada 89154, USA
| |
Collapse
|