1
|
Chang H, Li W, Sun Z. New Diabatic Potential Energy Surfaces for the Li + H 2 Reaction and Time-Dependent Quantum Wave Packet Studies. J Phys Chem A 2024; 128:4412-4424. [PMID: 38787593 DOI: 10.1021/acs.jpca.4c00539] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/25/2024]
Abstract
New global diabatic potential energy surfaces (DPESs) for the ground (12A') and first excited (22A') states for the Li + H2 system were developed, with more than 30,000 energy points at the IC-MRCI+Q level of theory, utilizing the aug-cc-pV5Z basis set for the H atoms and the cc-pCV5Z basis set for the Li atom, fitted by a single neural network (NN) with symmetry. Product state-resolved quantum dynamics calculations of the nonadiabatic reaction Li (2P) + H2 (X 1 ∑g+, v0 = 0, j0 = 0) → LiH (X 1∑+) + H(2S) were carried out using these new DPESs and also the previous HYLC-DPESs. The numerical results suggested that our newly constructed DPESs provided an accurate description of the LiH2 system.
Collapse
Affiliation(s)
- Hanwen Chang
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Wentao Li
- Weifang University of Science and Technology, Shouguang 262700, China
| | - Zhigang Sun
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China
| |
Collapse
|
2
|
Chen H, Buren B, Yang Z, Chen M. An effective approximation of Coriolis coupling in reactive scattering: application to the time-dependent wave packet calculations. Phys Chem Chem Phys 2023; 25:22927-22940. [PMID: 37591811 DOI: 10.1039/d3cp00530e] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/19/2023]
Abstract
Coriolis coupling plays a crucial role in reactive scattering, but dynamics calculations including the complete Coriolis coupling significantly increase the difficulty of numerical evolution due to the corresponding expensive matrix processing. The coupled state approximation that completely ignores the off-diagonal Coriolis coupling saves computational cost significantly but its error is usually unacceptable. In this paper, an improved coupled state approximation inspired by recently published results [D. Yang, X. Hu, D. H. Zhang and D. Xie, J. Chem. Phys., 2018, 148, 084101.] of the inelastic scattering problem is extended to deal with the reactive scattering. The calculations using the time-dependent wave packet method reveal that the new method can accurately reproduce the rigorous results of the H + HD (j0 < 3) → D + H2 reaction and immensely improve the computational efficiency. Additionally, we extend the new method to the non-adiabatic Li(2p) + H2 (v0 = 0, j0 = 0, 1) → H + LiH reaction, showcasing its advantages of low computational cost and high accuracy.
Collapse
Affiliation(s)
- Hanghang 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.
| | - Bayaer Buren
- School of Science, Shenyang University of Technology, Shenyang 110870, 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
|
3
|
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
|
4
|
Accurate Adiabatic and Diabatic Potential Energy Surfaces for the Reaction of He + H 2. BIOMED RESEARCH INTERNATIONAL 2022; 2022:7552881. [PMID: 35757471 PMCID: PMC9225863 DOI: 10.1155/2022/7552881] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 04/10/2022] [Accepted: 05/12/2022] [Indexed: 11/17/2022]
Abstract
The accurate adiabatic and diabatic potential energy surfaces, which are for the two lowest states of He + H2, are presented in this study. The Molpro 2012 software package is used, and the large basis sets (aug-cc-pV5Z) are selected. The high-level MCSCF/MRCI method is employed to calculate the adiabatic potential energy points of the title reaction system. The triatomic reaction system is described by Jacobi coordinates, and the adiabatic potential energy surfaces are fitted accurately using the B-spline method. The equilibrium structures and electronic energies for the H2 are provided, and the corresponding different levels of vibrational energies of the ground state are deduced. To better express the diabatic process of the whole reaction, avoid crossing points being calculated and conical intersection also being optimized. Meanwhile, the diabatic potential energy surfaces of the reaction process are constructed. This study will be helpful for the analysis of histopathology and for the study in biological and medical mechanisms.
Collapse
|
5
|
He D, Li W, Wang M. A study on the non-adiabatic dynamics of the Li(2p) + H2 → Li(2 s) + H2 quenching reaction calculated by time-dependent wavepacket method. Chem Phys Lett 2021. [DOI: 10.1016/j.cplett.2021.138910] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
|
6
|
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
|
7
|
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
|
8
|
Bai Y, Buren B, Yang Z, Chen M. Time-dependent wave packet dynamics study of non-adiabatic Li(2p) + HD → LiH/LiD + D/H reaction in a diabatic representation. Chem Phys Lett 2021. [DOI: 10.1016/j.cplett.2020.138279] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
9
|
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
|
10
|
Sathyamurthy N, Mahapatra S. Time-dependent quantum mechanical wave packet dynamics. Phys Chem Chem Phys 2020; 23:7586-7614. [PMID: 33306771 DOI: 10.1039/d0cp03929b] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Starting from a model study of the collinear (H, H2) exchange reaction in 1959, the time-dependent quantum mechanical wave packet (TDQMWP) method has come a long way in dealing with systems as large as Cl + CH4. The fast Fourier transform method for evaluating the second order spatial derivative of the wave function and split-operator method or Chebyshev polynomial expansion for determining the time evolution of the wave function for the system have made the approach highly accurate from a practical point of view. The TDQMWP methodology has been able to predict state-to-state differential and integral reaction cross sections accurately, in agreement with available experimental results for three dimensional (H, H2) collisions, and identify reactive scattering resonances too. It has become a practical computational tool in predicting the observables for many A + BC exchange reactions in three dimensions and a number of larger systems. It is equally amenable to determining the bound and quasi-bound states for a variety of molecular systems. Just as it is able to deal with dissociative processes (without involving basis set expansion), it is able to deal with multi-mode nonadiabatic dynamics in multiple electronic states with equal ease. We present an overview of the method and its strength and limitations, citing examples largely from our own research groups.
Collapse
|
11
|
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
|
12
|
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]
|
13
|
Li W, Sun J, He D. Non-adiabatic dynamics studies of the H( 2S) + LiH(X 1Σ +) reaction by time-dependent wave packet method. Phys Chem Chem Phys 2020; 22:17587-17596. [PMID: 32716453 DOI: 10.1039/d0cp01803a] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Non-adiabatic dynamics studies of the H + LiH (v0 = 0, j0 = 0) reaction are carried out based on the diabatic potential energy surfaces (PESs) reported by He and co-workers (Sci. Rep., 2016, 6, 25083) in the collision energy range from 0.001 to 1.0 eV. The H(2S) + LiH(X1Σ+) → Li(22S)/Li(22P) + H2(X1Σg+) depletion reactions and the H'(2S) + LiH(X1Σ+) → H(2S) + LiH'(X1Σ+) exchange reaction are studied at the state-to-state level of theory. The dynamics properties, such as reaction probability, integral cross section, differential cross section, the ro-vibrational state distributions of the product and specific-state rate constant are calculated. In addition, the dynamics values of the H(2S) + LiH(X1Σ+) → Li(22S) + H2(X1Σg+) depletion reaction are compared with available theoretical results, which are based on the adiabatic PESs. Large discrepancies are found between the adiabatic and diabatic values, especially at low J values, which indicate that the non-adiabatic effect is very great and cannot be simply ignored. Furthermore, the deviations between adiabatic and diabatic values gradually decrease as the collision energy increases.
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.
| | | | | |
Collapse
|
14
|
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
|
15
|
Song J, Zhu Z. Dynamics studies of the Li(2S) + H2(X1Σg+) → LiH (X1Σ+) + H(2S) reaction by time-dependent wave packet and quasi-classical trajectory methods. COMPUT THEOR CHEM 2020. [DOI: 10.1016/j.comptc.2020.112703] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
16
|
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
|
17
|
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]
|
18
|
Wang S, Yang Z, Yuan J, Chen M. New diabatic potential energy surfaces of the NaH 2 system and dynamics studies for the Na(3p) + H 2 → NaH + H reaction. Sci Rep 2018; 8:17960. [PMID: 30568250 PMCID: PMC6299287 DOI: 10.1038/s41598-018-35987-z] [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: 06/13/2018] [Accepted: 11/14/2018] [Indexed: 12/04/2022] Open
Abstract
The Na(3p) + H2(X1Σg+) → NaH(X1Σ+) + H(2S) reaction plays an important role in the field of diabatic reaction dynamics. A set of new diabatic potential energy surfaces (PESs) of the NaH2 system are structured, which include the diabatic coupling between the lowest two adiabatic states. The electronic structure calculations are performed on the multi-reference configuration interaction level with the cc-pwCVQZ and aug-cc-PVQZ basis sets for Na and H atoms. 32402 geometries are chosen to construct the diabatic data by a unitary transformation based on the molecular property method. The diabatic matrix elements of [Formula: see text], [Formula: see text] and [Formula: see text] ([Formula: see text]) are fitted by the artificial neural network model. The spectroscopic constants of diatoms obtained from the present PESs are consistent with the experimental data. The topographical and intersection characteristics of the [Formula: see text] and [Formula: see text] surfaces are discussed. Based on the new PESs, the time-dependent quantum wave packet calculations are carried out to study the reaction mechanism of the title reaction in detail.
Collapse
Affiliation(s)
- 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
| | - 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
- 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
|
19
|
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
|
20
|
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
|
21
|
Fu L, Wang D, Huang X. Accurate potential energy surfaces for the first two lowest electronic states of the Li (2p) + H2 reaction. RSC Adv 2018; 8:15595-15602. [PMID: 35539505 PMCID: PMC9080088 DOI: 10.1039/c8ra02504e] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/23/2018] [Accepted: 04/15/2018] [Indexed: 01/17/2023] Open
Abstract
The accuracy of three-dimensional adiabatic and diabatic potential energy surfaces is calculated using ab initio methods and is numerically fitted for the two lowest electronic states 1 and 22A′ of the LiH2 system, which are very important for the Li (2p) + H2 reaction. The finite difference method is performed to generate the mixing angles, which are used to educe the diabatic potential from the adiabatic potential. The accurate conical intersection (CI) is studied in this work with three different basis sets. The energy of the conical intersection is slightly lower (nearly 0.12 eV) than that of the perpendicular intermediate on the first excited state. By analyzing the potential energy surfaces in this work we can suggest that the most possible reaction pathway for the title reaction is Li (2p) + H2 → LiH2 (22A′) (C2v) → CI → LiH2 (12A′) (C2v) → LiH⋯H → LiH (X1∑g+) + H. The conical intersection and (22A′) intermediate may play a vital role in the title reaction. Accurate diabatic potential energy surfaces for the Li (2p) + H2 → LiH (X) + H reaction are produced.![]()
Collapse
Affiliation(s)
- Liwei Fu
- Laboratory of Theoretical and Computational Chemistry
- Institute of Theoretical Chemistry
- Jilin University
- Changchun
- People’s Republic of China
| | - Dequan Wang
- Laboratory of Theoretical and Computational Chemistry
- Institute of Theoretical Chemistry
- Jilin University
- Changchun
- People’s Republic of China
| | - Xuri Huang
- Laboratory of Theoretical and Computational Chemistry
- Institute of Theoretical Chemistry
- Jilin University
- Changchun
- People’s Republic of China
| |
Collapse
|
22
|
Yuan J, He D, Wang S, Chen M, Han K. Diabatic potential energy surfaces of MgH2+ and dynamic studies for the Mg+(3p) + H2 → MgH+ + H reaction. Phys Chem Chem Phys 2018; 20:6638-6647. [DOI: 10.1039/c7cp08679b] [Citation(s) in RCA: 27] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The global diabatic potential energy surfaces for the Mg+(3p) + H2 → MgH+ + H reaction are structured for the first time.
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
| | - Di He
- School of Physics and Optoelectronic Engineering
- Ludong University
- Yantai 264025
- 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
| | - 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
| | - Keli Han
- State Key Laboratory of Molecular Reaction Dynamics
- Dalian Institute of Chemical Physics
- Chinese Academy of Science
- Dalian 116023
- P. R. China
| |
Collapse
|
23
|
Wang S, Yuan J, Li H, Chen M. A neural network potential energy surface for the NaH 2 system and dynamics studies on the H( 2S) + NaH(X 1Σ +) → Na( 2S) + H 2(X 1Σ g+) reaction. Phys Chem Chem Phys 2017; 19:19873-19880. [PMID: 28597884 DOI: 10.1039/c7cp02153d] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
In order to study the dynamics of the reaction H(2S) + NaH(X1Σ+) → Na(2S) + H2(X1Σg+), a new potential energy surface (PES) for the ground state of the NaH2 system is constructed based on 35 730 ab initio energy points. Using basis sets of quadruple zeta quality, multireference configuration interaction calculations with Davidson correction were carried out to obtain the ab initio energy points. The neural network method is used to fit the PES, and the root mean square error is very small (0.00639 eV). The bond lengths, dissociation energies, zero-point energies and spectroscopic constants of H2(X1Σg+) and NaH(X1Σ+) obtained on the new NaH2 PES are in good agreement with the experiment data. On the new PES, the reactant coordinate-based time-dependent wave packet method is applied to study the reaction dynamics of H(2S) + NaH(X1Σ+) → Na(2S) + H2(X1Σg+), and the reaction probabilities, integral cross-sections (ICSs) and differential cross-sections (DCSs) are obtained. There is no threshold in the reaction due to the absence of an energy barrier on the minimum energy path. When the collision energy increases, the ICSs decrease from a high value at low collision energy. The DCS results show that the angular distribution of the product molecules tends to the forward direction. Compared with the LiH2 system, the NaH2 system has a larger mass and the PES has a larger well at the H-NaH configuration, which leads to a higher ICS value in the H(2S) + NaH(X1Σ+) → Na(2S) + H2(X1Σg+) reaction. Because the H(2S) + NaH(X1Σ+) → Na(2S) + H2(X1Σg+) reaction releases more energy, the product molecules can be excited to a higher vibrational state.
Collapse
Affiliation(s)
- Shufen Wang
- Key Laboratory of Materials Modification by Laser, Electron, and Ion Beams (Ministry of Education), School of Physics and Optoelectronic Technology, Dalian University of Technology, Dalian, 116024, P. R. China.
| | | | | | | |
Collapse
|
24
|
He D, Yuan J, Chen M. Influence of rovibrational excitation on the non-diabatic state-to-state dynamics for the Li(2p) + H 2 → LiH + H reaction. Sci Rep 2017; 7:3084. [PMID: 28596568 PMCID: PMC5465068 DOI: 10.1038/s41598-017-03274-y] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/05/2016] [Accepted: 04/26/2017] [Indexed: 12/03/2022] Open
Abstract
The non-adiabatic state-to-state dynamics of the Li(2p) + H2 → LiH + H reaction has been studied using the time-dependent wave packet method, based on a set of diabatic potential energy surfaces recently developed by our group. Integral cross sections (ICSs) can be increase more than an order of magnitude by the vibrational excitation of H2, whereas the ICSs are barely affected by the rotational excitation of H2. Moreover, ICSs of the title reaction with vibrationally excited H2 decrease rapidly with increasing collision energy, which is a typical feature of non-threshold reaction. This phenomenon implies that the title reaction can transformed from an endothermic to an exothermic reaction by vibrational excitation of H2. With the increase of the collision energy, the sideways and backward scattered tendencies of LiH for the Li(2p) + H2(v = 0, j = 0, 1) → LiH + H reactions are enhanced slightly, while the backward scattering tendency of LiH for the Li(2p) + H2(v = 1, j = 0) → LiH + H reaction becomes remarkably weakened. For the reaction with vibrationally excited H2 molecule, both direct and indirect reaction mechanism exist simultaneously.
Collapse
Affiliation(s)
- Di He
- Key Laboratory of Materials Modification by Laser, Electron, and Ion Beams (Ministry of Education), School of Physics and Optoelectronic Technology, 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 and Optoelectronic Technology, 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 and Optoelectronic Technology, Dalian University of Technology, Dalian, 116024, P. R. China.
| |
Collapse
|
25
|
Wang S, He D, Li W, Chen M. A global potential energy surface and dynamics study of the Au + + H 2 → H + Au +H reaction. RSC Adv 2017. [DOI: 10.1039/c7ra05223e] [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] Open
Abstract
A global potential energy surface (PES) of the ground state of the Au+H2 system was constructed using a neural network method with permutation invariant polynomials.
Collapse
Affiliation(s)
- Shufen Wang
- Key Laboratory of Materials Modification by Laser, Electron, and Ion Beams (Ministry of Education)
- School of Physics and Optoelectronic Technology
- Dalian University of Technology
- Dalian 116024
- PR China
| | - Di He
- Key Laboratory of Materials Modification by Laser, Electron, and Ion Beams (Ministry of Education)
- School of Physics and Optoelectronic Technology
- Dalian University of Technology
- Dalian 116024
- PR China
| | - Wentao Li
- Key Laboratory of Materials Modification by Laser, Electron, and Ion Beams (Ministry of Education)
- School of Physics and Optoelectronic Technology
- 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 and Optoelectronic Technology
- Dalian University of Technology
- Dalian 116024
- PR China
| |
Collapse
|
26
|
He D, Yuan J, Li H, Chen M. A new potential energy surface of LiHCl system and dynamic studies for the Li( 2S) + HCl(X 1Σ +) → LiCl(X 1Σ +) + H( 2S) reaction. J Chem Phys 2016; 145:234312. [PMID: 28010078 DOI: 10.1063/1.4972229] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A new global potential energy surface (PES) is constructed for the ground state of LiHCl system based on high-quality ab initio energy points calculated using multi-reference configuration interaction calculations with the Davidson correction. The AVQZ and WCVQZ basis sets are employed for H and Li atoms, respectively. To compensate the relativistic effects of heavy element, the AWCVQZ-DK basis set is employed for Cl atom. The neural network method is used for fitting the PES, and the root mean square error is small (1.36 × 10-2 eV). The spectroscopic constants of the diatoms obtained from the new PES agree well with experimental data. The geometric characteristics of the transition state and the complex are examined and compared with the previous theoretical values. To study the reaction dynamics of the Li(2S) + HCl(X1Σ+) → LiCl(X1Σ+) + H(2S) reaction, quantum reactive scattering dynamics calculations using collection reactant-coordinate-based wave packet method are conducted based on the new PES. The results of the reaction probabilities indicate that a small barrier exists along the reaction path as observed from the PES. The integral cross section curves reveal that the product molecule LiCl is easily excited. In addition, the reaction is dominated by forward scattering, and similar pattern is observed from Becker's experiment.
Collapse
Affiliation(s)
- Di He
- Key Laboratory of Materials Modification by Laser, Electron, and Ion Beams (Ministry of Education), School of Physics and Optoelectronic Technology, Dalian University of Technology, Dalian 116024, People's Republic of China
| | - Jiuchuang Yuan
- Key Laboratory of Materials Modification by Laser, Electron, and Ion Beams (Ministry of Education), School of Physics and Optoelectronic Technology, Dalian University of Technology, Dalian 116024, People's Republic of China
| | - Huixing Li
- Key Laboratory of Materials Modification by Laser, Electron, and Ion Beams (Ministry of Education), School of Physics and Optoelectronic Technology, Dalian University of Technology, Dalian 116024, People's Republic of China
| | - Maodu Chen
- Key Laboratory of Materials Modification by Laser, Electron, and Ion Beams (Ministry of Education), School of Physics and Optoelectronic Technology, Dalian University of Technology, Dalian 116024, People's Republic of China
| |
Collapse
|