1
|
Mao Y, Buren B, Yang Z, Chen M. Electronically Nonadiabatic Effects on the Quantum Dynamics of the H a + BeH b+ → Be + + H aH b; H b + BeH a+ Reactions. J Phys Chem A 2022; 126:5574-5581. [PMID: 35948431 DOI: 10.1021/acs.jpca.2c04319] [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
Nonadiabatic effects are ubiquitous and play an important role in many chemical processes. Here, the adiabatic and nonadiabatic quantum scattering calculations of the H + BeH+ reaction are performed using the time-dependent wave packet method based on an accurate diabatic potential energy matrix that includes the lowest two electronic states and their couplings. The resulting integral cross sections reveal that the nonadiabatic effect significantly inhibits the reactivity of the BeH+-depletion channel but enhances that of the H-exchange channel. The vibrational excitation is suppressed, but the translational excitation is promoted for the H2 product in the BeH+-depletion channel when the nonadiabatic coupling is included. However, the nonadiabatic coupling has a mild effect on the H-exchange product-state distribution. When the nonadiabatic effect is considered, the differential cross sections of the H2 product become less polarized because of the formation of an excited-state complex, whereas the corresponding results of the H-exchange channel only present an increase in the magnitude at the backward region.
Collapse
Affiliation(s)
- Ye Mao
- Key Laboratory of Materials Modification by Laser, Electron, and Ion Beams (Ministry of Education), School of Physics, Dalian University of Technology, Dalian 116024, PR China
| | - Bayaer Buren
- Key Laboratory of Materials Modification by Laser, Electron, and Ion Beams (Ministry of Education), School of Physics, Dalian University of Technology, Dalian 116024, PR China
| | - Zijiang Yang
- Key Laboratory of Materials Modification by Laser, Electron, and Ion Beams (Ministry of Education), School of Physics, Dalian University of Technology, Dalian 116024, PR China
| | - Maodu Chen
- Key Laboratory of Materials Modification by Laser, Electron, and Ion Beams (Ministry of Education), School of Physics, Dalian University of Technology, Dalian 116024, PR China
| |
Collapse
|
2
|
Zhang Y, Xu J, Yang H, Xu J. Nonadiabatic dynamics studies of the H( 2S) + RbH(X 1Σ +) reaction: based on new diabatic potential energy surfaces. RSC Adv 2022; 12:19751-19762. [PMID: 35865202 DOI: 10.1039/d2ra03028d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/13/2022] [Accepted: 06/14/2022] [Indexed: 11/21/2022] Open
Abstract
The global diabatic potential energy surfaces (PESs) that correspond to the ground (12A') and first excited states (22A') of the RbH2 system PES are constructed based on 17 786 ab initio points. The neural network method is used to fit the PESs and the topographic features of the new diabatic PESs are discussed in detail. Based on the newly constructed diabatic PESs, the dynamics calculations of the H(2S) + RbH(X1Σ+) → Rb(52S) + H2(X1Σg +)/Rb(52P) + H2(X1Σg +) reactions are performed using the time-dependent wave packet method. The dynamics properties of these two channels such as the reaction probabilities, integral cross sections, and differential cross sections (DCSs) are calculated at state-to-state level of theory. The nonadiabatic effects are discussed in detail, and the results indicate that the adiabatic results are overestimated from the dynamics values. The DCSs of these two channels are forward biased, which indicates that the abstraction mechanism plays a dominant role in the reaction.
Collapse
Affiliation(s)
- Yong Zhang
- NEST Lab, Department of Chemistry, Department of Physics, College of Science, Shanghai University Shanghai 200444 China .,Department of Physics, Tonghua Normal University Tonghua Jilin 134002 China
| | - Jinghua Xu
- Department of Physics, Tonghua Normal University Tonghua Jilin 134002 China
| | - Haigang Yang
- Department of Physics, Tonghua Normal University Tonghua Jilin 134002 China
| | - Jiaqiang Xu
- NEST Lab, Department of Chemistry, Department of Physics, College of Science, Shanghai University Shanghai 200444 China
| |
Collapse
|
3
|
Buren B, Chen M. Stereodynamics-Controlled Product Branching in the Nonadiabatic H + NaD → Na(3s, 3p) + HD Reaction at Low Temperatures. J Phys Chem A 2022; 126:2453-2462. [PMID: 35434992 DOI: 10.1021/acs.jpca.2c00114] [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/30/2022]
Abstract
Nonadiabatic processes play an important role at energies near or higher than conical intersection of adiabatic potential energy surfaces in chemical reactions. In this work, dynamics of the nonadiabatic H + NaD reaction at low temperatures are studied by using the quantum wave packet method based on an improved L-shaped grid. The nonadiabatic H + NaD reaction has two exothermic reaction channels: Na(3s) + HD and Na(3p) + HD; the latter can only occur via nonadiabatic transition. The dynamics results show that the product branching of the H + NaD reaction at collision energies ranging from 20 to 80 cm-1 is controlled by stereodynamics. The Na(3s) and Na(3p) reaction channels occur through collinear collision and side-on collision, respectively. When the collision energy is lower than 20 cm-1, the resonance-mediated reaction mechanism is dominant in both the Na(3s) and Na(3p) reaction channels.
Collapse
Affiliation(s)
- Bayaer Buren
- Key Laboratory of Materials Modification by Laser, Electron, and Ion Beams (Ministry of Education), School of Physics, Dalian University of Technology, Dalian 116024, PR China
| | - Maodu Chen
- Key Laboratory of Materials Modification by Laser, Electron, and Ion Beams (Ministry of Education), School of Physics, Dalian University of Technology, Dalian 116024, PR China
| |
Collapse
|
4
|
Non-adiabatic Couplings Induced Complex-forming Mechanism in the H+MgH +→Mg ++H 2 Reaction. CHINESE J CHEM PHYS 2022. [DOI: 10.1063/1674-0068/cjcp2111237] [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]
|
5
|
Yang Z, Chen H, Mao Y, Chen M. Neural network potential energy surface and quantum dynamics studies for the Ca +( 2S) + H 2 → CaH + + H reaction. Phys Chem Chem Phys 2022; 24:19209-19217. [DOI: 10.1039/d2cp02711a] [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
Reactive collisions of Ca+ ion with H2 molecule play a crucial role in ultracold chemistry, quantum information and other cutting-edge fields, and have been widely concerned experimentally, but the corresponding...
Collapse
|
6
|
Li Y, Liu J, Li J, Zhai Y, Yang J, Qu Z, Li H. A new permutation-symmetry-adapted machine learning diabatization procedure and its application in MgH 2 system. J Chem Phys 2021; 155:214102. [PMID: 34879675 DOI: 10.1063/5.0072004] [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
This work introduces a new permutation-symmetry-adapted machine learning diabatization procedure, termed the diabatization by equivariant neural network (DENN). In this approach, the permutation symmetric and anti-symmetric elements in diabatic potential energy metrics (DPEMs) were simultaneously simulated by the equivariant neural network. The diabatization by deep neural network scheme was adopted for machine learning diabatization, and non-zero diabatic coupling was included to increase accuracy in the near degenerate region. Based on DENN, the global DPEMs for 11A' and 21A' states of MgH2 have been constructed. To the best of our knowledge, these are the first global DPEMs for the MgH2 system. The root-mean-square-errors (RMSEs) for diagonal elements (H11 and H22) and the off-diagonal element (H12) around the conical intersection region were 5.824, 5.307, and 5.796 meV, respectively. The RMSEs of global adiabatic energies for two adiabatic states were 4.613 and 12.755 meV, respectively. The spectroscopic calculations of the 11A' state in the linear HMgH region are in good agreement with the experiment and previous theoretical results. The differences between calculated frequencies and corresponding experiment values are 1.38 and 1.08 cm-1 for anti-symmetric stretching fundamental vibrational frequency and first overtone. The global DPEMs obtained in this work should be useful for further quantum mechanics dynamic simulations on the MgH2 system.
Collapse
Affiliation(s)
- You Li
- Institute of Theoretical Chemistry, College of Chemistry, Jilin University, 2519 Jiefang Road, Changchun 130023, People's Republic of China
| | - Jingmin Liu
- Institute of Theoretical Chemistry, College of Chemistry, Jilin University, 2519 Jiefang Road, Changchun 130023, People's Republic of China
| | - Jiarui Li
- Institute of Theoretical Chemistry, College of Chemistry, Jilin University, 2519 Jiefang Road, Changchun 130023, People's Republic of China
| | - Yu Zhai
- Institute of Theoretical Chemistry, College of Chemistry, Jilin University, 2519 Jiefang Road, Changchun 130023, People's Republic of China
| | - Jitai Yang
- Institute of Theoretical Chemistry, College of Chemistry, Jilin University, 2519 Jiefang Road, Changchun 130023, People's Republic of China
| | - Zexing Qu
- Institute of Theoretical Chemistry, College of Chemistry, Jilin University, 2519 Jiefang Road, Changchun 130023, People's Republic of China
| | - Hui Li
- Institute of Theoretical Chemistry, College of Chemistry, Jilin University, 2519 Jiefang Road, Changchun 130023, People's Republic of China
| |
Collapse
|
7
|
|
8
|
Westermayr J, Marquetand P. Machine Learning for Electronically Excited States of Molecules. Chem Rev 2021; 121:9873-9926. [PMID: 33211478 PMCID: PMC8391943 DOI: 10.1021/acs.chemrev.0c00749] [Citation(s) in RCA: 162] [Impact Index Per Article: 54.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/17/2020] [Indexed: 12/11/2022]
Abstract
Electronically excited states of molecules are at the heart of photochemistry, photophysics, as well as photobiology and also play a role in material science. Their theoretical description requires highly accurate quantum chemical calculations, which are computationally expensive. In this review, we focus on not only how machine learning is employed to speed up such excited-state simulations but also how this branch of artificial intelligence can be used to advance this exciting research field in all its aspects. Discussed applications of machine learning for excited states include excited-state dynamics simulations, static calculations of absorption spectra, as well as many others. In order to put these studies into context, we discuss the promises and pitfalls of the involved machine learning techniques. Since the latter are mostly based on quantum chemistry calculations, we also provide a short introduction into excited-state electronic structure methods and approaches for nonadiabatic dynamics simulations and describe tricks and problems when using them in machine learning for excited states of molecules.
Collapse
Affiliation(s)
- Julia Westermayr
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
| | - Philipp Marquetand
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Vienna
Research Platform on Accelerating Photoreaction Discovery, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Data
Science @ Uni Vienna, University of Vienna, Währinger Strasse 29, 1090 Vienna, Austria
| |
Collapse
|
9
|
Abstract
Electronically excited states of molecules are at the heart of photochemistry, photophysics, as well as photobiology and also play a role in material science. Their theoretical description requires highly accurate quantum chemical calculations, which are computationally expensive. In this review, we focus on not only how machine learning is employed to speed up such excited-state simulations but also how this branch of artificial intelligence can be used to advance this exciting research field in all its aspects. Discussed applications of machine learning for excited states include excited-state dynamics simulations, static calculations of absorption spectra, as well as many others. In order to put these studies into context, we discuss the promises and pitfalls of the involved machine learning techniques. Since the latter are mostly based on quantum chemistry calculations, we also provide a short introduction into excited-state electronic structure methods and approaches for nonadiabatic dynamics simulations and describe tricks and problems when using them in machine learning for excited states of molecules.
Collapse
Affiliation(s)
- Julia Westermayr
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
| | - Philipp Marquetand
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Vienna Research Platform on Accelerating Photoreaction Discovery, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Data Science @ Uni Vienna, University of Vienna, Währinger Strasse 29, 1090 Vienna, Austria
| |
Collapse
|
10
|
Takatsuka K, Arasaki Y. An orbital picture extracted from correlated electronic wavefunctions and application to forbidden reactions: 70 years of the frontier orbital theory. J Chem Phys 2021; 155:064104. [PMID: 34391356 DOI: 10.1063/5.0059370] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023] Open
Abstract
The one-electron picture in molecular electronic state theory, particularly the molecular orbital (MO) theory with the Hartree-Fock approximation, has set a foundation to develop chemical science. Frontier orbital theory, or the theory of HOMO (highest occupied MO)-LUMO (lowest unoccupied MO) interaction, and the conservation rule of orbital symmetry are among the brightest achievements in a molecular orbital picture. After 70 years from the birth of frontier orbital theory, however, electronic wavefunctions treated in current quantum chemistry are often highly correlated and consist of extensive scales of electronic configurations to be more accurate and to cope with far more complicated reactions than concerted reactions. Under such circumstances, the MO approximation itself readily loses its validity, let alone the utter dominance of the HOMO-LUMO interaction. Recently, we have proposed an invariant method to extract general orbitals from such correlated electronic wavefunctions, which we refer to as Energy Natural Orbitals (ENOs) [K. Takatsuka and Y. Arasaki, J. Chem. Phys. 154, 094103 (2021)]. The energies of ENOs are summed exactly to the total electronic energy. The topological (symmetry) properties of a total wavefunction are represented by the relative phases of ENOs along with the continuity and crossing (avoided and conical intersection) among them. Only a small number of ENOs often dominate and characterize chemical reactions. With these properties of ENO, we explore a couple of simple and typical symmetry forbidden reactions, illustrating the effects of electron correlation and degeneracy in relevant ENOs. We propose the notion of "internal conical intersection" among ENOs, which leads to Jahn-Teller effect, pseudo-Jahn-Teller effect, and so on. We dare to explain the primary origin of elementary conical intersections and multidimensional avoided crossing in chemical reactions with the use of the notion of orbital crossing between those of HOMO-HOMO and LUMO-LUMO interactions and so on.
Collapse
Affiliation(s)
- Kazuo Takatsuka
- Fukui Institute for Fundamental Chemistry, Kyoto University, 606-8103 Kyoto, Japan
| | - Yasuki Arasaki
- Fukui Institute for Fundamental Chemistry, Kyoto University, 606-8103 Kyoto, Japan
| |
Collapse
|
11
|
Abstract
Theoretical simulations of electronic excitations and associated processes in molecules are indispensable for fundamental research and technological innovations. However, such simulations are notoriously challenging to perform with quantum mechanical methods. Advances in machine learning open many new avenues for assisting molecular excited-state simulations. In this Review, we track such progress, assess the current state of the art and highlight the critical issues to solve in the future. We overview a broad range of machine learning applications in excited-state research, which include the prediction of molecular properties, improvements of quantum mechanical methods for the calculations of excited-state properties and the search for new materials. Machine learning approaches can help us understand hidden factors that influence photo-processes, leading to a better control of such processes and new rules for the design of materials for optoelectronic applications.
Collapse
|
12
|
Takatsuka K. Electron Dynamics in Molecular Elementary Processes and Chemical Reactions. BULLETIN OF THE CHEMICAL SOCIETY OF JAPAN 2021. [DOI: 10.1246/bcsj.20200388] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Affiliation(s)
- Kazuo Takatsuka
- Fukui Institute for Fundamental Chemistry, Kyoto University, Kyoto 606-8103, Japan
| |
Collapse
|
13
|
Yin Z, Braams BJ, Fu B, Zhang DH. Neural Network Representation of Three-State Quasidiabatic Hamiltonians Based on the Transformation Properties from a Valence Bond Model: Three Singlet States of H3+. J Chem Theory Comput 2021; 17:1678-1690. [DOI: 10.1021/acs.jctc.0c01336] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Zhengxi Yin
- 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, P. R. China
- University of Chinese Academy of Sciences, Beijing 100049, P. R. China
| | - Bastiaan J. Braams
- Centrum Wiskunde & Informatica (CWI), the Dutch National Center for Mathematics and Computer Science, 1098 XG Amsterdam, Netherlands
| | - 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, Dalian 116023, P. R. 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, Dalian 116023, P. R. China
| |
Collapse
|
14
|
Yang Z, Mao Y, Chen M. Quantum Dynamics Studies of the Significant Intramolecular Isotope Effects on the Nonadiabatic Be +( 2P) + HD → BeH +/BeD + + D/H Reaction. J Phys Chem A 2021; 125:235-242. [PMID: 33369408 DOI: 10.1021/acs.jpca.0c09593] [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/30/2022]
Abstract
Quantum time-dependent wave packet dynamics studies on the nonadiabatic Be+(2P) + HD → BeH+/BeD+ + D/H reaction are performed for the first time employing recently constructed diabatic potential energy surfaces. Strong intramolecular isotope effects and unusual results are presented, which are attributed to the dynamic effects of shallow wells induced by avoided crossing on the diagonal V22d surface. The BeH+ + D and BeD+ + H channels are dominated by high-J and low-J partial waves, respectively. The BeD+/BeH+ branching ratio is larger than 10 at low energy and gradually decreases with increasing collision energy. The BeH+ product is primarily distributed at low vibrational states, whereas there exists an obvious population inversion of vibrational states on the BeD+ product. The results of differential cross sections suggest that the formation of the BeH+ + D channel favors a direct reaction process, while the BeD+ + H channel is mainly generated by the complex-forming mechanism.
Collapse
Affiliation(s)
- Zijiang Yang
- Key Laboratory of Materials Modification by Laser, Electron, and Ion Beams (Ministry of Education), School of Physics, Dalian University of Technology, Dalian 116024, PR China
| | - Ye Mao
- Key Laboratory of Materials Modification by Laser, Electron, and Ion Beams (Ministry of Education), School of Physics, Dalian University of Technology, Dalian 116024, PR China
| | - Maodu Chen
- Key Laboratory of Materials Modification by Laser, Electron, and Ion Beams (Ministry of Education), School of Physics, Dalian University of Technology, Dalian 116024, PR China
| |
Collapse
|
15
|
Yin Z, Braams BJ, Guan Y, Fu B, Zhang DH. A fundamental invariant-neural network representation of quasi-diabatic Hamiltonians for the two lowest states of H3. Phys Chem Chem Phys 2021; 23:1082-1091. [DOI: 10.1039/d0cp05047d] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The FI-NN approach is capable of representing highly accurate diabatic PESs with particular and complicated symmetry problems.
Collapse
Affiliation(s)
- Zhengxi Yin
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry
- Dalian Institute of Chemical Physics
- Chinese Academy of Sciences
- Dalian
- P. R. China
| | - Bastiaan J. Braams
- Centrum Wiskunde & Informatica (CWI)
- The Dutch national Center for Mathematics and Computer Science
- The Netherlands
| | - Yafu Guan
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry
- Dalian Institute of Chemical Physics
- Chinese Academy of Sciences
- Dalian
- P. R. 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
- Dalian
- P. R. 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
- Dalian
- P. R. China
| |
Collapse
|
16
|
Yin R, Gao N, Cao J, Li Y, Wang D, Huang X. Global accurate diabatic potential surfaces for the reaction H + Li 2. RSC Adv 2020; 10:39226-39240. [PMID: 35518445 PMCID: PMC9057374 DOI: 10.1039/d0ra05777k] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/02/2020] [Accepted: 09/21/2020] [Indexed: 12/25/2022] Open
Abstract
The adiabatic potential energies for the lowest three states of a Li2H system are calculated with a high level ab initio method (MCSCF/MRCI) with a large basis set (aV5Z). The accurate three dimensional B-spline fitting method is used to map the global adiabatic potential energy surfaces, using the existing adiabatic potential energies, for the lowest two adiabatic states of the title reaction system. The different vibrational states and corresponding energies are studied for the diatomic molecule of reactant and products. In order to clearly understand the nonadiabatic process, the avoided crossing area and conical intersection are carefully studied. For further study of the nonadiabatic dynamic reaction, the diabatic potential energy surfaces are deduced in the present work.
Collapse
Affiliation(s)
- Ruilin Yin
- Institute of Theoretical Chemistry, Jilin University Changchun People's Republic of China
| | - Nan Gao
- Department of Thoracic Surgery, China-Japan Union Hospital of Jilin University Changchun People's Republic of China
| | - Jing Cao
- Institute of Theoretical Chemistry, Jilin University Changchun People's Republic of China
| | - Yanchun Li
- Institute of Theoretical Chemistry, Jilin University Changchun People's Republic of China
| | - Dequan Wang
- Institute of Theoretical Chemistry, Jilin University Changchun People's Republic of China
| | - Xuri Huang
- Institute of Theoretical Chemistry, Jilin University Changchun People's Republic of China
| |
Collapse
|
17
|
Westermayr J, Marquetand P. Machine learning and excited-state molecular dynamics. MACHINE LEARNING-SCIENCE AND TECHNOLOGY 2020. [DOI: 10.1088/2632-2153/ab9c3e] [Citation(s) in RCA: 27] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
|
18
|
Hong Y, Yin Z, Guan Y, Zhang Z, Fu B, Zhang DH. Exclusive Neural Network Representation of the Quasi-Diabatic Hamiltonians Including Conical Intersections. J Phys Chem Lett 2020; 11:7552-7558. [PMID: 32835486 DOI: 10.1021/acs.jpclett.0c02173] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/26/2023]
Abstract
We propose a numerically simple and straightforward, yet accurate and efficient neural networks-based fitting strategy to construct coupled potential energy surfaces (PESs) in a quasi-diabatic representation. The fundamental invariants are incorporated to account for the complete nuclear permutation inversion symmetry. Instead of derivative couplings or interstate couplings, a so-called modified derivative coupling term is fitted by neural networks, resulting in accurate description of near degeneracy points, such as the conical intersections. The adiabatic energies, energy gradients, and derivative couplings are well reproduced, and the vanishing of derivative couplings as well as the isotropic topography of adiabatic and diabatic energies in asymptotic regions are automatically satisfied. All of these features of the coupled global PESs are requisite for accurate dynamics simulations. Our approach is expected to be very useful in developing highly accurate coupled PESs in a quasi-diabatic representation in an efficient machine learning-based way.
Collapse
Affiliation(s)
- Yingyue Hong
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian, P.R. China 116023
- University of Chinese Academy of Sciences, Beijing 100049, People's Republic of China
| | - Zhengxi Yin
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian, P.R. China 116023
- University of Chinese Academy of Sciences, Beijing 100049, People's Republic of China
| | - Yafu Guan
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian, P.R. China 116023
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - 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, P.R. China 116023
| | - 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, Dalian, P.R. China 116023
| | - 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, P.R. China 116023
| |
Collapse
|
19
|
Li W, Wang X, Zhao H, He D. Non-adiabatic dynamics studies of the K(4p 2P) + H 2(X 1Σ) reaction based on new diabatic potential energy surfaces. Phys Chem Chem Phys 2020; 22:16203-16214. [PMID: 32643736 DOI: 10.1039/d0cp02859b] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Global diabatic potential energy surfaces (PESs) of the KH2 system corresponding to the ground (12A') and first excited (22A') states were constructed for the first time. In ab initio calculations, the MRCI-F12 method with AVTZ and def2-QZVP basis sets was adopted and 17 865 ab initio energy points were calculated. The mixing angle, which is used to obtain the diabatic energies, was calculated by the molecular properties of the transition dipole moment. The diabatic PESs were fitted individually by the permutation invariant polynomial neural network method and the topographical features of the diabatic PESs are discussed in detail. The non-adiabatic dynamics studies of the K(4p2P) + H2(v0 = 0, 1, j0 = 0) reaction were carried out using the APH method based on the new diabatic PESs. The collision reaction processes K(4p2P) + H2(v0 = 0, 1, j0 = 0) → H + KH and the quenching processes K(4p2P) + H2(v0 = 0, 1, j0 = 0) → K(4s2S) + H2 were studied at the state-to-state level of theory. For the reaction process, the dynamics results indicated that the vibrational excitation of H2 was significantly more effective at promoting the reaction than the translational energy. In addition, the differential cross-sections were forward-biased scattering, which indicated that the direct abstraction mechanism plays a dominant role in the reaction. For the quenching process, the vibrational excitation of H2 molecules could improve the quenching efficiency obviously.
Collapse
Affiliation(s)
- Wentao Li
- Center for Theoretical and Computational Chemistry and State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China.
| | - Xuemei Wang
- Weifang University of Science and Technology, Shouguang 262700, China
| | - Hailin Zhao
- Center for Theoretical and Computational Chemistry and State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China.
| | - Di He
- School of Physics and Optoelectronics Engineering, Ludong University, Yantai 264025, China
| |
Collapse
|
20
|
Mao Y, Yuan J, Yang Z, Chen M. Quantum dynamics studies of isotope effects in the Mg +(3p) + HD → MgH +/MgD + + D/H insertion reaction. Sci Rep 2020; 10:3410. [PMID: 32098984 PMCID: PMC7042225 DOI: 10.1038/s41598-020-60033-2] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/04/2019] [Accepted: 02/03/2020] [Indexed: 11/17/2022] Open
Abstract
The time-dependent wave packet quantum dynamics studies for the Mg+(3p) + HD → MgH+/MgD+ + D/H diabatic reaction are carried out for the first time on recently developed diabatic YHWCH potential energy surfaces [Phys. Chem. Chem. Phys., 2018, 20, 6638-6647]. The results of reaction probabilities and total integral cross sections show a dramatic preference to the formation of MgD+ over MgH+ owing to the insertion reaction mechanism in the title reaction. The MgD+/MgH+ branching ratio witnesses a monotonic decrease from 10.58 to 3.88 at collision energy range of 0.01 to 0.20 eV, and at the collision energy of 0.114 eV, it is close to the experimental value of 5. The rovibrational state-resolved ICSs of the two channels show the products MgD+ have higher vibrational excitation and hotter rotational state distributions. The opacity function P(J) suggests that the MgH+ + D channel and MgD+ + H channel are dominated by high-b and low-b collisions, respectively. Both forward and backward scattering peaks are found in the differential cross section curves, whereas the angle distributions of products are not strictly forward-backward symmetric because of the short lifetime of the complex in the reaction.
Collapse
Affiliation(s)
- Ye Mao
- Key Laboratory of Materials Modification by Laser, Electron, and Ion Beams (Ministry of Education), School of Physics, Dalian University of Technology, Dalian, 116024, PR China
| | - Jiuchuang Yuan
- Key Laboratory of Materials Modification by Laser, Electron, and Ion Beams (Ministry of Education), School of Physics, Dalian University of Technology, Dalian, 116024, PR China
| | - Zijiang Yang
- Key Laboratory of Materials Modification by Laser, Electron, and Ion Beams (Ministry of Education), School of Physics, Dalian University of Technology, Dalian, 116024, PR China
| | - Maodu Chen
- Key Laboratory of Materials Modification by Laser, Electron, and Ion Beams (Ministry of Education), School of Physics, Dalian University of Technology, Dalian, 116024, PR China.
| |
Collapse
|
21
|
Buren B, Yang Z, Chen M. Dynamics study on the non-adiabatic Na(3p) + HD → NaH/NaD + D/H reaction: insertion-abstraction mechanism. Phys Chem Chem Phys 2020; 22:3633-3642. [PMID: 31998904 DOI: 10.1039/c9cp06026j] [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
Time-dependent wave packet calculations are carried out for two reaction channels of the non-adiabatic Na(3p) + HD → NaH/NaD + D/H reaction. The potential well on the excited state potential energy surface makes the reaction preferable to proceed through the insertion reaction path. The dominance of the NaD + H reaction channel and product rotational state distributions are found to be in agreement with the characteristics of typical adiabatic insertion reactions. However, significant forward scattering peaks in the differential cross sections (DCS) are found to be inconsistent with the forward-backward symmetric scattering characteristic of typical adiabatic insertion reactions, which indicate that the Na(3p) + HD reaction is dominated by a direct reaction mechanism. The comparison between adiabatic and non-adiabatic calculated DCSs reveals that the non-adiabatic couplings in the reaction could reduce the lifetime of the intermediate complex. Finally, the insertion-abstraction mechanism is put forward for the non-adiabatic Na(3p) + HD reaction.
Collapse
Affiliation(s)
- Bayaer Buren
- Key Laboratory of Materials Modification by Laser, Electron, and Ion Beams (Ministry of Education), School of Physics, Dalian University of Technology, Dalian 116024, P. R. China.
| | - Zijiang Yang
- Key Laboratory of Materials Modification by Laser, Electron, and Ion Beams (Ministry of Education), School of Physics, Dalian University of Technology, Dalian 116024, P. R. China.
| | - Maodu Chen
- Key Laboratory of Materials Modification by Laser, Electron, and Ion Beams (Ministry of Education), School of Physics, Dalian University of Technology, Dalian 116024, P. R. China.
| |
Collapse
|
22
|
Buren B, Yang Z, Chen M. Non-adiabatic state-to-state dynamic studies of Na(3p) + H2(v = 1, 2, 3; j = 0) → NaH + H reactions. Chem Phys Lett 2019. [DOI: 10.1016/j.cplett.2019.03.044] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
|
23
|
Yang Z, Wang S, Yuan J, Chen M. Neural network potential energy surface and dynamical isotope effects for the N+(3P) + H2 → NH+ + H reaction. Phys Chem Chem Phys 2019; 21:22203-22214. [DOI: 10.1039/c9cp02798j] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Dynamical isotope effects are calculated for the N+(3P) + H2 → NH+ + H reaction on a new neural network potential energy surface.
Collapse
Affiliation(s)
- Zijiang Yang
- Key Laboratory of Materials Modification by Laser, Electron, and Ion Beams (Ministry of Education)
- School of Physics
- Dalian University of Technology
- Dalian 116024
- P. R. China
| | - Shufen Wang
- Key Laboratory of Materials Modification by Laser, Electron, and Ion Beams (Ministry of Education)
- School of Physics
- Dalian University of Technology
- Dalian 116024
- P. R. China
| | - Jiuchuang Yuan
- State Key Laboratory of Molecular Reaction Dynamics
- Dalian Institute of Chemical Physics
- Chinese Academy of Science
- Dalian 116023
- P. R. China
| | - Maodu Chen
- Key Laboratory of Materials Modification by Laser, Electron, and Ion Beams (Ministry of Education)
- School of Physics
- Dalian University of Technology
- Dalian 116024
- P. R. China
| |
Collapse
|
24
|
Yuan J, Duan Z, Wang S, Liu J, Han K. Significant effects of vibrational excitation of reactant in K + H 2 → H + KH reaction based on a new neural network potential energy surface. Phys Chem Chem Phys 2018; 20:20641-20649. [PMID: 30059113 DOI: 10.1039/c8cp03310b] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
The study of K + H2 collision has a long experimental history, but there have been few theoretical studies due to lack of a global potential energy surface (PES). In this study, a new global PES for the ground state of KH2 system was constructed based on numerous ab initio points, using the permutation invariant polynomial neural network method. The root mean square error (RMSE) of PES is very small (5.64 meV). On the new PES, time-dependent quantum wave packet (TDWP) and quasiclassical trajectory (QCT) calculations were carried out to study the dynamics of the K(2S) + H2(X1Σ+g) → H(2S) + KH(X1Σ+) reaction. Dynamics results show that (i) the K + H2(v = 0) → H + KH reaction scarcely occurred, (ii) the K + H2(v = 1) → H + KH reaction took place in small quantities, and (iii) the K + H2(v = 2) → H + KH reaction occurred in large quantities. This indicates that vibrational energy of the reactant is significantly more effective at promoting the reaction than the translational energy. This characteristic stems from a major physical model in reactive collisions: the vibrationally excited H2 molecule and K atom collide first in a T-geometric configuration and the vibrational motion of the H2 molecule helps separate the two H atoms a large distance after the collision. At a large H-H distance, a broad well exists on the PES, so the heavy K atom could pull back the light H atom to initiate the reaction. Similar to the reactive channel, vibrational excitation of the reactant also has a significant effect on the collision-induced dissociation channel.
Collapse
Affiliation(s)
- Jiuchuang Yuan
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Science, Dalian 116023, P. R. China.
| | | | | | | | | |
Collapse
|
25
|
Yang Z, Yuan J, Wang S, Chen M. Global diabatic potential energy surfaces for the BeH 2 + system and dynamics studies on the Be +( 2P) + H 2(X 1Σ g +) → BeH +(X 1Σ +) + H( 2S) reaction. RSC Adv 2018; 8:22823-22834. [PMID: 35539737 PMCID: PMC9081383 DOI: 10.1039/c8ra04305a] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/21/2018] [Revised: 07/13/2018] [Accepted: 06/14/2018] [Indexed: 11/21/2022] Open
Abstract
The Be+(2P) + H2(X1Σg +) → BeH+(X1Σ+) + H(2S) reaction has great significance for studying diabatic processes and ultracold chemistry. The first global diabatic potential energy surfaces (PESs) which are correlated with the lowest two adiabatic states 12A' and 22A' of the BeH2 + system are constructed by using the neural network method. Ab initio energy points are calculated using the multi-reference configuration interaction method with the Davidson correction and AVQZ basis set. The diabatic energies are obtained from the transformation of ab initio data based on the dipole moment operators. The topographical characteristics of the diabatic PESs are described in detail, and the positions of crossing between the V d 11 and V d 22 are pinpointed. On new diabatic PESs, the time-dependent quantum wave packet method is carried out to study the mechanism of the title reaction. The results of dynamics calculations indicate the reaction has no threshold and the product BeH+ is excited to high vibrational states easily. In addition, the product BeH+ tends to backward scattering at most collision energies.
Collapse
Affiliation(s)
- Zijiang Yang
- Key Laboratory of Materials Modification by Laser, Electron, and Ion Beams (Ministry of Education), School of Physics, Dalian University of Technology Dalian 116024 P. R. China
| | - Jiuchuang Yuan
- Key Laboratory of Materials Modification by Laser, Electron, and Ion Beams (Ministry of Education), School of Physics, Dalian University of Technology Dalian 116024 P. R. China
| | - Shufen Wang
- Key Laboratory of Materials Modification by Laser, Electron, and Ion Beams (Ministry of Education), School of Physics, Dalian University of Technology Dalian 116024 P. R. China
| | - Maodu Chen
- Key Laboratory of Materials Modification by Laser, Electron, and Ion Beams (Ministry of Education), School of Physics, Dalian University of Technology Dalian 116024 P. R. China
| |
Collapse
|