1
|
Shu Y, Zhao H, Sun Z. The Interaction-Asymptotic Region Decomposition Method in Jacobi Coordinates: Triatomic Reactive Scatterings. J Phys Chem A 2024; 128:7982-7993. [PMID: 39256192 DOI: 10.1021/acs.jpca.4c04101] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/12/2024]
Abstract
The interaction-asymptotic region decomposition (IARD) technique has been proven to be a good solution to the long-standing coordinate problem in reactive scattering calculations. In this work, the IARD technique was further developed using Jacobi coordinates for the interaction region, instead of the previously used hyperspherical coordinates. Although the Jacobi coordinate may not be as optimal as the hyperspherical coordinates for describing the interaction region in reactive scattering processes, it has simpler kinetic operators and provides a more physically intuitive picture. By developing an intermediate interpolation method, which could efficiently transform the overlapped wave functions from the asymptotic regions to the interaction region, the new implementation of the IARD technique for triatomic reactive scatterings is similarly efficient and accurate. The differential cross sections of the H+H2, and product state-resolved reaction probabilities of the F+HD and16O+36O2 reactions, which involve products of extremely low translational energy and are challenging for a single coordinate-based method, were calculated as numerical examples to show the ability of the new method.
Collapse
Affiliation(s)
- Yajian Shu
- 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
| | - Hailin Zhao
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China
| | - Zhigang Sun
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China
| |
Collapse
|
2
|
Zhao H, Sun Z. An improved method for reactive scatterings in ultra-cold conditions using the time-dependent approach. Phys Chem Chem Phys 2024; 26:22790-22797. [PMID: 39163005 DOI: 10.1039/d4cp02340d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/21/2024]
Abstract
A new efficient method for considering the long-range effect of reactive scattering processes in ultra-cold conditions has been developed using the time-dependent quantum wave packet theory, where the initial wave packet could be placed at a position near the interaction region. This is in contrast to previous methods, where the initial wave packet has to be placed far from the interaction region. The new method reduces the numerical effort significantly. Typical reactions, such as S(1D) + H2, D+ + H2, and 7Li + 7Li2 (v0 = 1, j0 = 0), under cold or ultra-cold conditions, are used to demonstrate the numerical efficiency of the new method.
Collapse
Affiliation(s)
- Hailin Zhao
- Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, P. R. China.
| | - Zhigang Sun
- Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, P. R. China.
| |
Collapse
|
3
|
Li W, Dong B, Niu X, Wang M, Zhang Y. Non-adiabatic dynamics studies of the C+(2P1/2, 3/2) + H2 reaction: Based on global diabatic potential energy surfaces of CH2. J Chem Phys 2024; 161:074302. [PMID: 39145555 DOI: 10.1063/5.0223199] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/12/2024] [Accepted: 07/31/2024] [Indexed: 08/16/2024] Open
Abstract
Global diabatic potential energy surfaces (PESs) of CH2+ are constructed using the neural network method with a specific function based on 18 213 ab initio points. The multi-reference configuration interaction method with the aug-cc-pVQZ basis set is adopted to perform the ab initio calculations. The topographical properties of the diabatic PESs are examined in detail. In general, the diabatic PESs provide an accurate quasi-diabatic representation. To validate the diabatic PESs, the dynamics studies of the C+(2P1/2, 3/2) + H2 (v0 = 0, j0 = 0) → H + CH+(X1Σ+) reaction are performed using the time-dependent wave packet method. The reaction probabilities, integral cross sections, differential cross sections, and rate constants are calculated and compared with the experimental and theoretical results. Non-adiabatic dynamics results are in good agreement with experimental data. In addition, the non-adiabatic effect in the C+(2P1/2, 3/2) + H2 reaction is significant due to the non-adiabatic results being obviously larger than adiabatic values. The reasonable non-adiabatic dynamics results indicate that present diabatic PESs can be recommended for any type of dynamics study.
Collapse
Affiliation(s)
- Wentao Li
- Weifang University of Science and Technology, Shouguang 262700, China
| | - Bin Dong
- Weifang University of Science and Technology, Shouguang 262700, China
| | - Xianghong Niu
- School of Science, Nanjing University of Posts and Telecommunications, Nanjing 210003, China
| | - Meishan Wang
- College of Integrated Ciruits, Ludong University, Yantai 264025, People's Republic of China
| | - Yong Zhang
- Department of Physics, Tonghua Normal University, Tonghua, Jilin 134002, China
| |
Collapse
|
4
|
Li W, Zhang Z, Niu X, He D, Xing W, Zhang Y. Diabatic Potential Energy Surfaces of SrH 2+ and Dynamics Studies of the Sr +(5s 2S) + H 2 Reaction. J Phys Chem A 2024; 128:6677-6684. [PMID: 39093206 DOI: 10.1021/acs.jpca.4c03648] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/04/2024]
Abstract
Based on the ab initio energy points of both ground and excited states, a neural network fitting method combined with a specific function was successfully used to construct the diabatic potential energy surfaces (PESs) of the SrH2+ system. The topographical features of the diabatic PESs were examined in detail. The results indicate that the nonadiabatic transition characteristics between ground and excited states are accurately described by the newly constructed diabatic PESs. To verify the validity and applicability of the diabatic PESs, as well as the nonadiabatic effects during the reaction process, the quantum dynamics studies of the Sr+(5s2S) + H2 reaction were performed based on both adiabatic and diabatic PESs. The dynamics results indicate that adiabatic dynamics results are dozens of times larger than those of nonadiabatic. This illustrates the significant effect of nonadiabaticity, indicating that adiabatic dynamics results often overestimate the actual values. The integral cross sections (ICSs) were calculated and compared with the experimental data. The diabatic ICSs are in good agreement with the experimental results. The reasonable dynamics results indicate that the newly constructed diabatic PESs are suitable for the relevant dynamics studies.
Collapse
Affiliation(s)
- Wentao Li
- Weifang University of Science and Technology, Shouguang 262700, China
| | - Zhijun Zhang
- Weifang University of Science and Technology, Shouguang 262700, China
| | - Xianghong Niu
- School of Science, Nanjing University of Posts and Telecommunications, Nanjing 210003, China
| | - Di He
- Weifang University of Science and Technology, Shouguang 262700, China
| | - Wei Xing
- College of Physics and Electronic Engineering, Xinyang Normal University, Xinyang 464000, China
| | - Yong Zhang
- Department of Physics, Tonghua Normal University, Tonghua ,Jilin134002, China
| |
Collapse
|
5
|
Li W, Liang Y, Niu X, He D, Xing W, Zhang Y. Construction of diabatic potential energy surfaces for the SiH2+ system and dynamics studies of the Si+(2P1/2, 3/2) + H2 reaction. J Chem Phys 2024; 161:044310. [PMID: 39051835 DOI: 10.1063/5.0219621] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/18/2024] [Accepted: 07/11/2024] [Indexed: 07/27/2024] Open
Abstract
The construction of diabatic potential energy surfaces (PESs) for the SiH2+ system, related to the ground (12A') and excited states (22A'), has been successfully achieved. This was accomplished by utilizing high-level ab initio energy points, employing a neural network fitting method in conjunction with a specifically designed function. The newly constructed diabatic PESs are carefully examined for dynamics calculations of the Si+(2P1/2, 3/2) + H2 reaction. Through time-dependent quantum wave packet calculations, the reaction probabilities, integral cross sections (ICSs), and differential cross sections (DCSs) of the Si+(2P1/2, 3/2) + H2 reaction were reported. The dynamics results indicate that the total ICS is in excellent agreement with experimental data within the collision energy range studied. The results also indicate that the SiH+ ion is hardly formed via the Si+(2P3/2) + H2 reaction. The results from the DCSs suggest that the "complex-forming" reaction mechanism predominates in the low collision energy region. Conversely, the forward abstraction reaction mechanism is dominant in the high collision energy region.
Collapse
Affiliation(s)
- Wentao Li
- Weifang University of Science and Technology, Shouguang 262700, China
| | - Yongping Liang
- Weifang University of Science and Technology, Shouguang 262700, China
| | - Xianghong Niu
- School of Science, Nanjing University of Posts and Telecommunications, Nanjing 210003, China
| | - Di He
- Weifang University of Science and Technology, Shouguang 262700, China
| | - Wei Xing
- College of Physics and Electronic Engineering, Xinyang Normal University, Xinyang 464000, People's Republic of China
| | - Yong Zhang
- Department of Physics, Tonghua Normal University, Tonghua, Jilin 134002, China
| |
Collapse
|
6
|
Mao Y, Chen H, Yang Z, Buren B, Chen M. Stereodynamic control of nonadiabatic processes in low-energy Be +( 2P) + H 2 ( v = 0, j = 2) collisions. Phys Chem Chem Phys 2024; 26:19812-19821. [PMID: 38988212 DOI: 10.1039/d4cp01996b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/12/2024]
Abstract
Controlling the relative arrangement of colliding molecules is crucial for determining the dynamical outcomes of chemical processes and has emerged as a hot spot of experimental research. Here, the quantum scattering calculations are conducted to investigate the stereodynamic control in collisions between Be+(2P) and H2 (v = 0, j = 2), which undergo nonadiabatic transitions to the electronic ground state. Stereodynamic preparation is achieved by controlling the initial alignment of the H2 bond axis relative to the scattering frame. For product BeH+ in the reactive process, the differential cross sections (DCSs) are significantly enhanced in the forward and sideways hemispheres when the alignment angle β is 60°. For the product H2 in the quenching channel, the β = 0° preparation can result in a more than one-fold increase in the DCS at a polar scattering angle of 0°. Furthermore, varying the alignment angle β also has noteworthy effects on the rotational-state distributions of BeH+ products. Specifically, β = 0° preparation can induce the disappearance of the bimodal distribution of rotational states at a collision energy of 0.05 eV.
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, P. R. China.
| | - 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, 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.
| | - Bayaer Buren
- School of Science, Shenyang University of Technology, Shenyang 110870, 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
|
7
|
Yang Z, Cao F, Cheng H, Liu S, Sun J. A Globally Accurate Neural Network Potential Energy Surface and Quantum Dynamics Studies on Be +( 2S) + H 2/D 2 → BeH +/BeD + + H/D Reactions. Molecules 2024; 29:3436. [PMID: 39065017 DOI: 10.3390/molecules29143436] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/07/2024] [Revised: 07/18/2024] [Accepted: 07/20/2024] [Indexed: 07/28/2024] Open
Abstract
Chemical reactions between Be+ ions and H2 molecules have significance in the fields of ultracold chemistry and astrophysics, but the corresponding dynamics studies on the ground-state reaction have not been reported because of the lack of a global potential energy surface (PES). Herein, a globally accurate ground-state BeH2+ PES is constructed using the neural network model based on 18,657 ab initio points calculated by the multi-reference configuration interaction method with the aug-cc-PVQZ basis set. On the newly constructed PES, the state-to-state quantum dynamics calculations of the Be+(2S) + H2(v0 = 0; j0 = 0) and Be+(2S) + D2(v0 = 0; j0 = 0) reactions are performed using the time-dependent wave packet method. The calculated results suggest that the two reactions are dominated by the complex-forming mechanism and the direct abstraction process at relatively low and high collision energies, respectively, and the isotope substitution has little effect on the reaction dynamics characteristics. The new PES can be used to further study the reaction dynamics of the BeH2+ system, such as the effects of rovibrational excitations and alignment of reactant molecules, and the present dynamics data could provide an important reference for further experimental studies at a finer level.
Collapse
Affiliation(s)
- Zijiang Yang
- School of Physics and Electronic Technology, Liaoning Normal University, Dalian 116029, China
| | - Furong Cao
- School of Physics and Electronic Technology, Liaoning Normal University, Dalian 116029, China
| | - Huiying Cheng
- School of Physics and Electronic Technology, Liaoning Normal University, Dalian 116029, China
| | - Siwen Liu
- School of Physics and Electronic Technology, Liaoning Normal University, Dalian 116029, China
| | - Jingchang Sun
- School of Physics and Electronic Technology, Liaoning Normal University, Dalian 116029, China
| |
Collapse
|
8
|
Chen H, Mao Y, Yang Z, Chen M. A Neural Network Potential Energy Surface and Quantum Dynamics Study of Ca( 1S) + H 2( v 0 = 0, j 0 = 0) → CaH + H Reaction. ACS OMEGA 2024; 9:30804-30812. [PMID: 39035896 PMCID: PMC11256353 DOI: 10.1021/acsomega.4c03465] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 04/10/2024] [Revised: 06/12/2024] [Accepted: 06/18/2024] [Indexed: 07/23/2024]
Abstract
The reactive collision between Ca and H2 molecules has attracted great interest experimentally due to the key role of the product CaH molecule in the field of astrophysics and cold molecules. However, quantum dynamics calculations for this system have not been reported due to the lack of a global potential energy surface (PES). Herein, a globally accurate PES of the ground-state CaH2 is developed by combining 11365 high-level ab initio points and permutation invariant polynomial neural network method. Based on the newly constructed PES, the state-to-state quantum dynamics calculations for the Ca(1S) + H2 (v 0 = 0, j 0 = 0) → CaH + H reaction are carried out using the time-dependent wave packet method. The dynamic results reveal that the reaction follows the complex-forming mechanism near the reactive threshold, whereas both the indirect insertion mechanism and direct abstraction mechanism have effects at higher collision energies. The newly constructed PES can be used to further study the influence of isotope substitution, rovibrational excitation, and spatial orientation of reactant molecules on reaction dynamics.
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, 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, China
| | - Zijiang Yang
- School
of Physics and Electronic Technology, Liaoning
Normal University, Dalian 116029, 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, China
| |
Collapse
|
9
|
Buren B, Zhang J, Li Y. Quantum Dynamics Studies of the Li + Na 2 ( V = 0, j = 0) → Na + NaLi Reaction on a New Neural Network Potential Energy Surface. J Phys Chem A 2024; 128:5115-5127. [PMID: 38889710 DOI: 10.1021/acs.jpca.4c01891] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/20/2024]
Abstract
The ultracold reaction offers a unique opportunity to elucidate the intricate microscopic mechanism of chemical reactions, and the Na2Li system serves as a pivotal reaction system in the investigation of ultracold reactions. In this work, a high-precision potential energy surface (PES) of the Na2Li system is constructed based on high-level ab initio energy points and the neural network (NN) method, and a proper asymptotic functional form is adopted for the long-range interaction, which is suitable for the study of cold or ultracold collisions. Based on the new NN PES, the dynamics of the Li + Na2 (v = 0, j = 0) → Na + NaLi reaction are studied in the collision energy range of 10-7 to 80 cm-1. In the high collision energy range of 8 to 80 cm-1, the dynamics of the reaction is studied using the time-dependent wave packet method and the statistical quantum mechanical (SQM) method. Comparing the results of the two methods, it is found that the SQM method provides a rough description of the product ro-vibrational state distribution but overestimates the integral cross-section values. With the decrease of collision energy, the reaction differential cross section gradually changes from forward-backward symmetric scattering to predominant forward scattering. In the low collision energy range from 10-7 to 8 cm-1, the SQM method is used to study the reaction dynamics, and the rate constant in the Wigner threshold region is estimated to be 2.87 × 10-10 cm3/s.
Collapse
Affiliation(s)
- Bayaer Buren
- School of Science, Shenyang University of Technology, Shenyang 110870, China
| | - Jiapeng Zhang
- Department of Physics, Liaoning University, Shenyang 110036, China
| | - Yongqing Li
- Department of Physics, Liaoning University, Shenyang 110036, China
| |
Collapse
|
10
|
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
|
11
|
Chang H, Li W, Sun Z. New Potential Energy Surface for the H + Cl 2 Reaction and Quantum Dynamics Studies. J Phys Chem A 2024; 128:4425-4438. [PMID: 38805307 DOI: 10.1021/acs.jpca.4c00542] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/30/2024]
Abstract
The reaction of H + Cl2 → HCl + Cl plays a crucial role in various fields. However, no previous study has investigated this reaction using accurate quantum mechanical methods. In this paper, we construct a global potential energy surface (PES) using the neural network method with more than 20,000 ab initio energies obtained by the MRCI-F12+Q method with the aug-cc-pV5Z basis and extrapolated to the complete basis set limit. The spin-orbit coupling of the Cl atom has been considered in the PES. With this new PES, product state-resolved quantum dynamics calculations for the H + Cl2 (v0 = 0, j0 = 0-2) → HCl + Cl reaction was carried out. Numerical results show that the initial rotational excitation of the Cl2 has negligible effects on the reactivity. Product state-resolved integral cross sections (ICS) and rate constants reveal that the HCl is most favorably produced in its v' = 2 vibrational state. The calculated product vibrational state-resolved and total reaction rate constants suggest that the new global PES is accurate enough, as compared with the available experimental measurements.
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
|
12
|
Wang R, Sun Z, Alexander MH. Development of the Time-Independent Methods for the Cl + H 2/F + HD Reaction Using Hyper-Spherical Coordinates Including (Full) Spin-Orbit Characteristics. J Chem Theory Comput 2024; 20:3449-3461. [PMID: 38691764 DOI: 10.1021/acs.jctc.4c00161] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/03/2024]
Abstract
Recently, a combined study of high-resolution molecular crossed beam experiment and accurate full-dimensional time-dependent theory, including full spin-orbit characteristics on the effect of electronic spin and orbital angular momenta in the F + HD reaction, was reported by some of us, focusing on the partial wave resonance phenomenon (Science 2021, 371, 936-940). It revealed that the time-dependent theory could explain all of the details observed in the high-resolution experiment. Here, we develop two time-independent close-coupling methods using hyperspherical coordinates, including the two-state model, where only a part of the spin-orbit characteristics is considered, and the six-state model, where the full spin-orbit characteristics is considered. With these two newly developed theoretical models and the adiabatic theoretical model, the detailed reaction dynamics of the F + HD (v = 0, j = 0) reaction and the Cl + H2 (v = 0, j = 0) reaction are investigated and compared. Some of the results are compared with the time-dependent quantum wave packet theory and the experimental observations, and good agreements have been obtained, which suggests the validity of the pure-procession approximation in the six-state model using different theoretical methods. This work demonstrates the ability of the reactive scattering theory including full spin-orbit characteristics for describing the reactions of a halogen atom plus hydrogen molecule and its isotopologues.
Collapse
Affiliation(s)
- Ransheng Wang
- 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, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Zhigang Sun
- 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, China
| | - Millard H Alexander
- Department of Chemistry and Biochemistry and Institute for Physical Science and Technology, University of Maryland, College Park, Maryland 20742, United States
| |
Collapse
|
13
|
Wang R, Zhao H, Sun Z. Reactant-Product Decoupling Technique Using the Intermediate Coordinate Method. J Phys Chem A 2024; 128:3726-3741. [PMID: 38666315 DOI: 10.1021/acs.jpca.4c01148] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/12/2024]
Abstract
Although the reactant-product decoupling (RPD) technique was proposed over two decades ago, it remains an efficient approach for calculating product state-resolved information on some simple direct reactions using the quantum wave packet method. In the past, usually the RPD technique employed the collocation method to transform the wave function between reactant and product arrangements, which requires quite large computational efforts. In this work, the intermediate coordinate (IC) method is employed to realize the RPD technique. Numerical examples demonstrate that this new IC RPD (IRPD) technique has superior computational efficiency compared with the original method employing the collocation method. Especially, the new IRPD technique significantly saves disk space and computer memory. To illustrate the features of our new method, the total reaction probabilities of the H + H2, H + Br2, and F + H2 reactions with J = 0 and the differential cross sections of the H + H2 and F + H2 reactions at a series of collision energy are calculated and presented. With this efficient and effective new RPD technique, the Li + HF reaction, which involves sharp resonances with long-range wave functions in the van der Waals wells in both the reactant and product arrangements, is also calculated with several J at the product state-resolved level to reveal the ability of the RPD technique for describing resonance wave functions. With these numerical examples, it is found that, for the reaction with resonances, the RPD approach should be applied carefully. Otherwise, it is very possible that the resonances could disappear with the application of the RPD technique.
Collapse
Affiliation(s)
- Ransheng Wang
- 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, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Hailin Zhao
- 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, China
| | - Zhigang Sun
- 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, China
| |
Collapse
|
14
|
Wang R, Shi H, Xie T, Sun Z. Significant Isotope Effects from the Nonadiabatic Couplings in the Cl( 2P) + HD( v = 0, j = 0) Reaction. J Phys Chem A 2024; 128:3301-3310. [PMID: 38648526 DOI: 10.1021/acs.jpca.4c00830] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/25/2024]
Abstract
The impact of non-Born-Oppenheimer couplings on the isotopic effects in the reaction of the Cl(2P) atom with the HD (v = 0, j = 0) molecule is investigated with our recently developed nonadiabatic time-independent quantum scattering methods, where the full open-shell characteristics are included in the six-state model, and also with the recently developed two-state model solving by time-independent methods, where part of the open-shell characteristic is included. The same reaction is also calculated with the simple adiabatic model using the lowest adiabatic potential energy surface. Compared with the results from different models, it is found that the reactivity of the Cl + HD → HCl + D channel is significantly overestimated in the adiabatic model. In contrast, the reactivity of the other channel agrees well with the nonadiabatic models. This is due to the van der Waals well in the reactant channel being changed a lot by including the nonadiabatic couplings. These quantum dynamics calculations suggest that sometimes the adiabatic model should be used with caution; otherwise, it may result in significant deviations for some reactions.
Collapse
Affiliation(s)
- Ransheng Wang
- 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, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Haimei Shi
- 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, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Ting Xie
- 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, China
| | - Zhigang Sun
- 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, China
| |
Collapse
|
15
|
Chang H, Li W, Sun Z. Product State-Resolved Reactive Scattering Studies of the H + Cl 2 ( v0 = 1-3, j0 = 0) → HCl + Cl Reaction by the Time-Dependent Wave Packet Method. J Phys Chem A 2024; 128:2997-3006. [PMID: 38593417 DOI: 10.1021/acs.jpca.4c01062] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/11/2024]
Abstract
The typical hydrogen atom plus halogen molecule reaction H + Cl2 → HCl + Cl has implications across many fields. In this paper, product state-resolved quantum dynamics calculations for the vibrationally excited reaction H + Cl2 (v0 = 1-3, j0 = 0) → HCl + Cl were conducted using the time-dependent wave packet method on a newly developed accurate potential energy surface. Numerical results indicate that the initial vibrational excitation of Cl2 does enhance the reactivity for this early barrier reaction, although less than the enhancement of the translational energy. The calculated product vibrational state-resolved integral cross sections and rate constants reveal that the product vibrational state distribution and the initial vibrational state of Cl2 are highly correlated. The thermal rate constant in the temperature range from 100 to 1000 K was given and is found to be in reasonable agreement with the experimental measurements.
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
|
16
|
Zhao H, Sun Z. Theoretical Development of the Interaction-Asymptotic Region Decomposition Method for Tetratomic Reactive Scattering. J Chem Theory Comput 2024; 20:1802-1810. [PMID: 38262035 DOI: 10.1021/acs.jctc.3c01368] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2024]
Abstract
An accurate and efficient time-dependent wave packet method is proposed for solving the product state-resolved reaction probabilities of the tetratomic reactive system. In this method, the entire scattering process is divided into the interaction region and multiple asymptotic regions, sharing the same spirit as the interaction-asymptotic region decomposition (IARD) approach in a triatomic reactive scattering process. The hyperspherical coordinate is adopted in the interaction region, while the corresponding Jacobi coordinate is employed in each asymptotic region. Therefore, in this IARD method, the "coordinate problem", the difficulty of expressing the wave function in the entire region using a single coordinate system, can be effectively avoided, and only a very small number of the grid points (or the basis functions) are required. For the numerical illustration, the typical tetratomic reaction H2 + OH with zero total angular momentum is calculated, and compared with other quantum wave packet methods. Our proposed IARD method for the tetratomic reactive system is much more efficient and accurate.
Collapse
Affiliation(s)
- Hailin Zhao
- Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian, P.R. China 116023
| | - Zhigang Sun
- Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian, P.R. China 116023
| |
Collapse
|
17
|
He D, Li W, Li Q, Chen S, Wang L, Liu Y, Wang M. The impact of non-adiabatic effects on reaction dynamics: a study based on the adiabatic and non-adiabatic potential energy surfaces of CaH 2. Phys Chem Chem Phys 2023; 25:22744-22754. [PMID: 37605513 DOI: 10.1039/d3cp02416d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/23/2023]
Abstract
The two-state non-adiabatic potential energy matrices of the CaH2+ system are calculated via a diabatization approach by using a neural network model. Subsequently, the adiabatic and non-adiabatic potential energy surfaces (PESs) are constructed based on these non-adiabatic potential energy matrices. Furthermore, based on the adiabatic and non-adiabatic PESs, the Ca+(4s2S) + H2(X1Σ+g) → H(2S) + CaH+(X1Σ+) reaction is studied using the time-dependent wave packet method. Comparative analysis of the experimental and theoretical integral reaction cross-sections (ICSs) indicates that the maximum deviation between the results obtained from the adiabatic PES and the corresponding experimental value is 12.7 bohr2; in contrast, the maximum discrepancy between the theoretical result derived from the non-adiabatic PES and the experimental value is merely 0.42 bohr2. The potential well along the reaction path acts as a 'filter', selectively guiding intermediates with longer lifetimes in the potential well back to the reactant channel. This phenomenon indicates that the non-adiabatic effects significantly influence the reaction dynamics of the CaH2+ system.
Collapse
Affiliation(s)
- Di He
- School of Physics and Optoelectronics Engineering, Ludong University, Yantai 264025, China.
| | - Wentao Li
- Weifang University of Science and Technology, Shouguang 262700, China
| | - Quanjiang Li
- School of Physics and Optoelectronics Engineering, Ludong University, Yantai 264025, China.
| | - Shenghui Chen
- School of Physics and Optoelectronics Engineering, Ludong University, Yantai 264025, China.
| | - Li Wang
- School of Physics and Optoelectronics Engineering, Ludong University, Yantai 264025, China.
| | - Yanli Liu
- School of Physics and Optoelectronics Engineering, Ludong University, Yantai 264025, China.
| | - Meishan Wang
- School of Physics and Optoelectronics Engineering, Ludong University, Yantai 264025, China.
| |
Collapse
|
18
|
Yang Z, Chen H, Buren B, Chen M. Globally Accurate Gaussian Process Potential Energy Surface and Quantum Dynamics Studies on the Li( 2S) + Na 2 → LiNa + Na Reaction at Low Collision Energies. Molecules 2023; 28:2938. [PMID: 37049701 PMCID: PMC10096016 DOI: 10.3390/molecules28072938] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/03/2023] [Revised: 03/22/2023] [Accepted: 03/23/2023] [Indexed: 03/29/2023] Open
Abstract
The LiNa2 reactive system has recently received great attention in the experimental study of ultracold chemical reactions, but the corresponding theoretical calculations have not been carried out. Here, we report the first globally accurate ground-state LiNa2 potential energy surface (PES) using a Gaussian process model based on only 1776 actively selected high-level ab initio training points. The constructed PES had high precision and strong generalization capability. On the new PES, the quantum dynamics calculations on the Li(2S) + Na2(v = 0, j = 0) → LiNa + Na reaction were carried out in the 0.001-0.01 eV collision energy range using an improved time-dependent wave packet method. The calculated results indicate that this reaction is dominated by a complex-forming mechanism at low collision energies. The presented dynamics data provide guidance for experimental research, and the newly constructed PES could be further used for ultracold reaction dynamics calculations on this reactive system.
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, China
| | - 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, China
| | - Bayaer Buren
- School of Science, Shenyang University of Technology, Shenyang 110870, 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, China
| |
Collapse
|
19
|
Zhu Z, Feng Y, Li W. The state-to-state dynamics of the N + NH( 3Σ −) → N 2(X 1 Σ g+) + H reaction: based on a new global potential energy surface. Mol Phys 2023. [DOI: 10.1080/00268976.2022.2162455] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023]
Affiliation(s)
- Ziliang Zhu
- Shandong Peninsula Engineering Research Center of Comprehensive Brine Utilization, Weifang University of Science and Technology, Shouguang, People’s Republic of China
- Key Laboratory for Structure and Environment Disaster Preventing of Agriculture Greenhouse, Weifang University of Science and Technology, Shouguang, People’s Republic of China
- College of General Education, Weifang University of Science and Technology, Shouguang, People’s Republic of China
| | - Yinghua Feng
- College of General Education, Weifang University of Science and Technology, Shouguang, People’s Republic of China
| | - Wentao Li
- College of General Education, Weifang University of Science and Technology, Shouguang, People’s Republic of China
| |
Collapse
|
20
|
Zhang Y, Xu J. An analytic global potential energy surface of the CsH 2 system and the dynamic calculations of the H + CsH reaction. Mol Phys 2022. [DOI: 10.1080/00268976.2022.2151524] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/02/2022]
Affiliation(s)
- Yong Zhang
- NEST Lab, Department of Chemistry, Department of Physics, College of Science, Shanghai University, Shanghai, People’s Republic of China
- Department of Physics, Tonghua Normal University, Tonghua, People’s Republic of China
| | - Jiaqiang Xu
- NEST Lab, Department of Chemistry, Department of Physics, College of Science, Shanghai University, Shanghai, People’s Republic of China
| |
Collapse
|
21
|
State-to-state dynamical studies of the C+ + H2(v = 0, j = 0) → CH+ + H reaction using time-dependent wave packet method. Chem Phys 2022. [DOI: 10.1016/j.chemphys.2022.111670] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
22
|
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
|
23
|
Zhang Y, Xu J, Yang H, Xu J. The global potential energy surface of the RbH2 system and dynamics studies of the H + RbH → Rb + H2 reaction. Chem Phys 2022. [DOI: 10.1016/j.chemphys.2022.111586] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/03/2022]
|
24
|
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
|
25
|
Mao Y, Buren B, Yang Z, Chen M. Time-dependent wave packet dynamics study of the resonances in the H + LiH +( v = 0, j = 0) → Li + + H 2 reaction at low collision energies. Phys Chem Chem Phys 2022; 24:15532-15539. [PMID: 35713276 DOI: 10.1039/d1cp05601h] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The depletion process of LiH+ by H collision plays an important role in the evolution of the early universe and astrophysical processes, including the eventual charge-states, abundances of atomic and molecular species and ensuing astrochemistry. Here, a quantum dynamics study on the H + LiH+(v = 0, j = 0) → Li+ + H2 reaction is performed at the low collision energy range from 0.1 meV to 10 meV using the time-dependent wave packet method. A Feshbach resonance peak is observed near 0.8 meV collision energy on the total reaction probability curves. This resonance originates from the coupling with the v = 0, j = 1 energy level of the reactant LiH+, and it is dominated by the contributions of J = 0-4 partial waves. Another partial wave resonance is also found on the total integral cross section at 1.2 meV, which is closely connected to the opening of the J = 7 partial wave. The opening of the J = 7 partial wave generates a notable forward scattering peak, and the Feshbach resonance can promote both the forward and backward scatterings. Moreover, the total and product vibrational state-resolved rate coefficients for the temperature range of 1-100 K are also reported.
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, P. R. 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, 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
|
26
|
Zhang Z, Zhu Z, Li W. The quenching processes of the Na(3 p) + H 2 → Na(3 s) + H 2 reaction studied using the time-dependent wave packet method. Mol Phys 2022. [DOI: 10.1080/00268976.2022.2083713] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
Affiliation(s)
- Zhijun Zhang
- College of Mechanical and Electrical Technology, Weifang University of Science and Technology, Shouguang, People’s Republic of China
| | - Ziliang Zhu
- College of Mechanical and Electrical Technology, Weifang University of Science and Technology, Shouguang, People’s Republic of China
| | - Wentao Li
- College of Mechanical and Electrical Technology, Weifang University of Science and Technology, Shouguang, People’s Republic of China
| |
Collapse
|
27
|
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
|
28
|
Dynamic study of the D + DAu reaction based on a new ground potential energy surface. COMPUT THEOR CHEM 2022. [DOI: 10.1016/j.comptc.2021.113540] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
|
29
|
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]
|
30
|
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
|
31
|
Yang Z, Chen H, Chen M. Representing Globally Accurate Reactive Potential Energy Surfaces with Complex Topography by Combining Gaussian Process Regression and Neural Network. Phys Chem Chem Phys 2022; 24:12827-12836. [DOI: 10.1039/d2cp00719c] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
There has been increasing attention in using machine learning technologies, such as neural network (NN) and Gaussian process regression (GPR), to model multidimensional potential energy surfaces (PESs). NN PES features...
Collapse
|
32
|
Buren B, Chen M, Sun Z, Guo H. Quantum Wave Packet Treatment of Cold Nonadiabatic Reactive Scattering at the State-To-State Level. J Phys Chem A 2021; 125:10111-10120. [PMID: 34767377 DOI: 10.1021/acs.jpca.1c08105] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Cold and ultracold collisions are dominated by quantum effects, such as resonances, tunneling, and nonadiabatic transitions between different electronic states. Due to the extremely long de Broglie wavelength in such processes, quantum reactive scattering is most conveniently characterized using the time-independent close-coupling (TICC) methods. However, the TICC approach is difficult for systems with a large number of channels because of its steep numerical scaling laws. Here, a recently proposed quantum wave packet (WP) approach for solving adiabatic reactive scattering problems at low collision energies is extended to include nonadiabatic transitions. To impose the outgoing boundary conditions, the total scattering wavefunction is split into three parts, the interaction, the asymptotic, and the long-range regions. Each region is associated with a different set of basis functions, which could be optimized separately. In this way, an extremely long grid can be used to accommodate the characteristic long de Broglie wavelengths in the scattering coordinate. The better numerical scaling laws of the WP approach have the potential for handling larger nonadiabatic reactive systems at low temperatures in the future.
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, China.,Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| | - 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, China
| | - Zhigang Sun
- 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, China
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| |
Collapse
|
33
|
Ghosh S, Sharma R, Adhikari S, Varandas AJC. Dynamical calculations of O( 3P) + OH( 2Π) reaction on the CHIPR potential energy surface using the fully coupled time-dependent wave-packet approach in hyperspherical coordinates. Phys Chem Chem Phys 2021; 23:21784-21796. [PMID: 34550126 DOI: 10.1039/d1cp02488d] [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/21/2022]
Abstract
We have carried out quantum dynamics calculations for the O + OH → H + O2 reaction on the CHIPR [A. J. C. Varandas, J. Chem. Phys., 2013, 138, 134117] potential energy surface (PES) for ground state HO2 using the fully coupled 3D time-dependent wavepacket formalism [S. Adhikari and A. J. C. Varandas, Comput. Phys. Commun., 2013, 184, 270] in hyperspherical coordinates. Reaction probabilities for J > 0 are calculated for different initial rotational states of the OH radical (v = 0; j = 0, 1). State-to-state as well as total integral cross sections and rate-coefficients are evaluated and compared with previous theoretical calculations and available experimental studies. Using the rate constant for the forward (hereinafter considered to be H + O2 → O + OH) and backward (O + OH → H + O2) reactions of this reactive system, the equilibrium constant of the reversible process [H + O2 ⇌ O + OH] is calculated as a function of temperature and compared with previous experimental measurements.
Collapse
Affiliation(s)
- Sandip Ghosh
- School of Chemical Sciences, Indian Association for the Cultivation of Science, 2A & 2B Raja S. C. Mullick Road, Jadavpur, Kolkata-700032, West Bengal, India.
| | - Rahul Sharma
- Department of Chemistry, St. Xaviers' College, Kolkata-700016, West Bengal, India
| | - Satrajit Adhikari
- School of Chemical Sciences, Indian Association for the Cultivation of Science, 2A & 2B Raja S. C. Mullick Road, Jadavpur, Kolkata-700032, West Bengal, India.
| | - António J C Varandas
- School of Physics and Physical Engineering, Qufu Normal University, 273165 Qufu, China.,Department of Physics, Universidade Federal do Espírito Santo, 29075-910 Victória, Brazil.,Departamento de Qumica, and Centro de Qumica, Universidade de Coimbra, 3004-535 Coimbra, Portugal.
| |
Collapse
|
34
|
|
35
|
Li W, Liang Y, Xing W. Quantum calculations for the abstraction and exchange channels of the H + SH+(v0 = 0, j0 = 0) reaction. Chem Phys 2021. [DOI: 10.1016/j.chemphys.2021.111306] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
36
|
A new global potential energy surface of X2A′ state of Li2H system and quantum dynamics calculation of the H + Li2 → Li + LiH reaction. Chem Phys 2021. [DOI: 10.1016/j.chemphys.2021.111201] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
|
37
|
Zhang J. The quantum and quasi-classical dynamics calculations of the Au ++HD → AuH +/AuD ++D/H reaction. Mol Phys 2021. [DOI: 10.1080/00268976.2021.1913251] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
Affiliation(s)
- Jiali Zhang
- Weifang University of Science and Technology, Shouguang, People’s Republic of China
| |
Collapse
|
38
|
Zhao H, Xie D, Sun Z. Interaction-Asymptotic Region Decomposition Method for a Triatomic Reactive Scattering with Symmetry Adoption. J Phys Chem A 2021; 125:2460-2471. [DOI: 10.1021/acs.jpca.0c11438] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Hailin Zhao
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210093, China
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China
- Center for Advanced Chemical Physics and 2011 Frontier Centre for Quantum Science and Technology, University of Science and Technology of China, 96 Jinzhai Road, Hefei 230026, China
| | - Daiqian Xie
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210093, China
| | - Zhigang Sun
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China
- Center for Advanced Chemical Physics and 2011 Frontier Centre for Quantum Science and Technology, University of Science and Technology of China, 96 Jinzhai Road, Hefei 230026, China
| |
Collapse
|
39
|
Zhu Z, Zhang A, He D, Li W. A new global potential energy surface of the SH 2+(X 4A'') system and quantum calculations for the S + + H 2( v = 0-3, j = 0) reaction. Phys Chem Chem Phys 2021; 23:4757-4767. [PMID: 33599223 DOI: 10.1039/d0cp06335e] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
A new global potential energy surface (PES) for the ground state of the SH2+(X4A'') system is constructed using a permutation invariant polynomial neural network method. In ab initio calculations, the MRCI-F12 method with the AVTZ basis set is used. Furthermore, the dynamics calculations of the S+ + H2(v = 0-3, j = 0) → SH+ + H reaction are carried out based on the new PES. The reaction probabilities and integral cross sections are compared with available theoretical calculations. Present values are in general good agreement with the previous theoretical studies. However, some discrepancies can still be found due to different PESs used in the calculation. Furthermore, the vibrational energy of the reactant molecule can significantly enhance the reactivity compared to the translational energy. The differential cross sections indicated that the reaction mechanism is changed from the "head-on" rebound mechanism to the tripping mechanism with the increasing number of initial vibrational excitation state.
Collapse
Affiliation(s)
- Ziliang Zhu
- Shandong Peninsula Engineering Research Center of Comprehensive Brine Utilization and Key Laboratory for Structure and Environment Disaster Preventing of Agriculture Greenhouse, Weifang University of Science and Technology, Shouguang, 262700, People's Republic of China.
| | - Aijie Zhang
- Key Laboratory of Magnetic Molecules & Magnetic Information Materials Ministry of Education, School of Chemistry and Material Science, Shanxi Normal University, Linfen 041004, China
| | - Di He
- School of Physics and Optoelectronics Engineering, Ludong University, Yantai 264025, China
| | - Wentao Li
- Shandong Peninsula Engineering Research Center of Comprehensive Brine Utilization and Key Laboratory for Structure and Environment Disaster Preventing of Agriculture Greenhouse, Weifang University of Science and Technology, Shouguang, 262700, People's Republic of China.
| |
Collapse
|
40
|
Zhao H, Xie D, Sun Z. Interaction-Asymptotic Region Decomposition Method for an Insertion Reaction: Application to the S( 1D) + H 2 Reaction. J Phys Chem A 2021; 125:2007-2018. [PMID: 33625216 DOI: 10.1021/acs.jpca.1c00140] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
With adjusting principal axes hyperspherical (APH) coordinate in the interaction region, and the Jacobi coordinates in the asymptotic regions, an efficient multidomain interaction-asymptotic region decomposition (IARD) method has been developed to solve the "coordinate problem" in a product-state-resolved reactive scattering calculation using the quantum wave packet method. Although the APH coordinate treats with all three channels equally, and is efficient for describing the interaction region for some direct reactions, it is inefficient for describing the insertion-type reaction due to the singularity problem, such as the S(1D) + H2 reaction. To deal with this issue, in this work, the channel-dependent Delves hyperspherical (DH) coordinate is proposed to describe the interaction region using the IARD method. The proposed DH-IARD method was applied to calculate the product-state-resolved reaction probabilities of the H + HD reaction, and the differential and integral cross sections of the typical insertion reaction S(1D) + H2. It is found that the new DH-IARD method is much more efficient than the previous APH-IARD method for dealing with insertion reactions. The partial wave resonance structures were observed in the integral cross section. It is found that at a low collision energy, the position of the initial wave packet has to be put far away. Otherwise, the partial wave resonance structures could not be correctly reproduced due to the reef well arising with a large total angular momentum J.
Collapse
Affiliation(s)
- Hailin Zhao
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210093, China.,State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China.,Center for Advanced Chemical Physics and 2011 Frontier Centre for Quantum Science and Technology, University of Science and Technology of China, 96 Jinzhai Road, Hefei 230026, China
| | - Daiqian Xie
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210093, China
| | - Zhigang Sun
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China.,Center for Advanced Chemical Physics and 2011 Frontier Centre for Quantum Science and Technology, University of Science and Technology of China, 96 Jinzhai Road, Hefei 230026, China
| |
Collapse
|
41
|
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]
|
42
|
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
|
43
|
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
|
44
|
Goswami S, Sahoo J, Paul SK, Rao TR, Mahapatra S. Effect of Reagent Vibration and Rotation on the State-to-State Dynamics of the Hydrogen Exchange Reaction, H + H 2 → H 2 + H. J Phys Chem A 2020; 124:9343-9359. [PMID: 33124827 DOI: 10.1021/acs.jpca.0c06707] [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/30/2022]
Abstract
State-to-state dynamics of the benchmark hydrogen exchange reaction H + H2 (v = 0-4, j = 0-3) → H2 (v', j') + H is investigated with the aid of the real wave packet approach of Gray and Balint-Kurti (J. Chem. Phys. 1998, 108, 950-962) and electronic ground BKMP2 potential energy surface of Boothroyd et al. (J. Chem. Phys. 1996, 104, 7139-7152). Initial state-selected and product state-resolved reaction probabilities, integral cross section, and product diatom vibrational and rotational level populations at a few collision energies are reported to elucidate the energy disposal mechanism. State-specific thermal rate constants are also calculated and compared with the available literature results. Coriolis coupling terms of the nuclear Hamiltonian are included, and calculations are parallelized over the helicity quantum number, Ω'. Attempts are made, in particular, to study the effect of reagent vibrational and rotational excitations on the dynamical attributes. It is found that the calculations become computationally expensive with reagent vibrational and rotational excitation. Reagent vibrational excitation is found to enhance the reactivity and has significant impact on the energy disposal to the vibrational and rotational degrees of freedom of the product. The interplay of reagent translational and vibrational energy on the product vibrational distribution unfolds an important aspect of the energy disposal mechanism. The effect of reagent rotation on the state-to-state dynamics is found not to be very significant, and the weak effect turns out to be specific to v'.
Collapse
Affiliation(s)
- Sugata Goswami
- School of Chemistry, University of Hyderabad, Hyderabad 500046, India
| | - Jayakrushna Sahoo
- School of Chemistry, University of Hyderabad, Hyderabad 500046, India
| | - Suranjan K Paul
- School of Chemistry, University of Hyderabad, Hyderabad 500046, India
| | - T Rajagopala Rao
- School of Chemistry, University of Hyderabad, Hyderabad 500046, India
| | - S Mahapatra
- School of Chemistry, University of Hyderabad, Hyderabad 500046, India
| |
Collapse
|
45
|
Zhao B, Manthe U. Eight-Dimensional Wave Packet Dynamics Within the Quantum Transition-State Framework: State-to-State Reactive Scattering for H2 + CH3 ⇆ H + CH4. J Phys Chem A 2020; 124:9400-9412. [DOI: 10.1021/acs.jpca.0c08461] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Affiliation(s)
- Bin Zhao
- Theoretische Chemie, Fakultät für Chemie, Universität Bielefeld, Universitätsstraße 25, D-33615 Bielefeld, Germany
| | - Uwe Manthe
- Theoretische Chemie, Fakultät für Chemie, Universität Bielefeld, Universitätsstraße 25, D-33615 Bielefeld, Germany
| |
Collapse
|
46
|
|
47
|
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
|
48
|
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
|
49
|
Zhang AJ, Ma XQ, Jia JF, Wu HS, Yuan JC. A time-dependent quantum dynamical study of H + AuH reaction. Chem Phys Lett 2020. [DOI: 10.1016/j.cplett.2020.137295] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
|
50
|
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]
|