1
|
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
|
2
|
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
|
3
|
He D, Yuan J, Li H, Chen M. Global diabatic potential energy surfaces and quantum dynamical studies for the Li(2p) + H2(X(1)Σ(+)g) → LiH(X(1)Σ(+)) + H reaction. Sci Rep 2016; 6:25083. [PMID: 27125781 PMCID: PMC4850413 DOI: 10.1038/srep25083] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/04/2015] [Accepted: 04/11/2016] [Indexed: 11/09/2022] Open
Abstract
The global diabatic potential energy surfaces which are correlated with the ground state 1A' and the excited state 2A' of the Li(2p) + H2 reaction are presented in this study. The multi-reference configuration interaction method and large basis sets (aug-cc-pVQZ for H atom and cc-pwCVQZ for Li atom) were employed in the ab initio single-point energy calculations. The diabatic potential energies were generated by the diabatization scheme based on transition dipole moment operators. The neural network method was utilized to fit the matrix elements of the diabatic energy surfaces, and the root mean square errors were extremely small (3.69 meV for , 5.34 meV for and 5.06 meV for ). The topographical features of the diabatic potential energy surfaces were characterized and the surfaces were found to be sufficiently smooth for the dynamical calculation. The crossing seam of the conical intersections between the and surfaces were pinpointed. Based on this new analytical diabatic potential energy surfaces, time-dependent wave packet calculation were conducted to investigate the mechanism of the title reaction. At low collision energies, the product LiH molecule tends to forward scattering, while at high collision energies, the forward and backward scatterings 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, PR 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, PR 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, 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
|
4
|
Yuan J, He D, Chen M. A new potential energy surface for the H2S system and dynamics study on the S((1)D) + H2(X(1)Σg(+)) reaction. Sci Rep 2015; 5:14594. [PMID: 26435516 PMCID: PMC4592959 DOI: 10.1038/srep14594] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/08/2015] [Accepted: 09/01/2015] [Indexed: 11/09/2022] Open
Abstract
We constructed a new global potential energy surface (PES) for the electronic ground state ((1)A') of H2S based on 21,300 accurate ab initio energy points over a large configuration space. The ab initio energies are obtained from multireference configuration interaction calculations with a Davidson correction using basis sets of quadruple zeta quality. The neural network method is applied to fit the PES, and the root mean square error of fitting is small (1.68 meV). Time-dependent wave packet studies for the S((1)D) + H2(X(1)Σg(+)) → H((2)S) + SH(X(2)Π) reaction on the new PES are conducted to study the reaction dynamics. The calculated integral cross sections decrease with increasing collision energy and remain fairly constant within the high collision energy range. Both forward and backward scatterings can be observed as expected for a barrierless reaction with a deep well on the PES. The calculated integral cross sections and differential cross sections are in good agreement with the experimental results.
Collapse
Affiliation(s)
- 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, 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
| | - 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
|
5
|
Cheng D, Yuan J, Chen M. Non-Born–Oppenheimer stereodynamics study for the D + + H 2 ( v, j = 0) reaction using coherent switching with decay of mixing method. CAN J CHEM 2015. [DOI: 10.1139/cjc-2015-0031] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
Abstract
A non-Born–Oppenheimer (BO) investigation using coherent switching with decay of mixing method is carried out for both reactive non-charge transfer and charge transfer channels of reaction D+ + H2 (v, j = 0). For the non-charge transfer channel, it is found that the non-BO effect has a small influence on the differential cross section, but it has a large impact on the rotational polarization of the HD product, which is reflected by the joint distribution P(θr, ϕr). And it is also found that the polarization of the charge transfer products shows a quite different form. The rotational polarizations are also investigated at various collision energies and initial vibrational states for both reactive channels.
Collapse
Affiliation(s)
- Dahai Cheng
- 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
- 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
- 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
- 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
|
6
|
Yuan J, He D, Chen M. A new potential energy surface for the ground electronic state of the LiH2 system, and dynamics studies on the H((2)S) + LiH(X(1)Σ(+)) → Li((2)S) + H2(X(1)Σg(+)) reaction. Phys Chem Chem Phys 2015; 17:11732-9. [PMID: 25870863 DOI: 10.1039/c4cp05352d] [Citation(s) in RCA: 40] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
A new global potential energy surface (PES) is obtained for the ground electronic state of the LiH2 system based on high-level energies. The energy points are calculated at the multireference configuration interaction level with aug-cc-pVXZ (X = Q, 5) basis sets, and these energies are extrapolated to the complete basis set limit. The neural network method and hierarchical construction scheme are applied in the fitting process and the root mean square error of the fitting result is very small (0.004 eV). The dissociation energies and equilibrium distances for LiH(X(1)Σ(+)) and H2(X(1)Σg(+)) obtained from the new PES are in good agreement with the experimental data. On the new PES, time-dependent wave packet studies for the H((2)S) + LiH(X(1)Σ(+)) → Li((2)S) + H2(X(1)Σg(+)) reaction have been carried out. In this reaction, no threshold is found due to the absence of an energy barrier on the minimum energy path. The calculated integral cross sections are high at low collision energy and will decrease with the increase of the collision energy. The product molecule H2 tends to be forward scattering due to direct reactive collisions, which becomes more evident at higher collision energies.
Collapse
Affiliation(s)
- 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.
| | | | | |
Collapse
|
7
|
Yuan J, Cheng D, Chen M. Time-dependent wave packet and quasiclassical trajectory studies of the Au + HD reaction: competition between the reactive channels. RSC Adv 2014. [DOI: 10.1039/c4ra06297c] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022] Open
|