1
|
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
|
2
|
Jara-Toro RA, Roncero O, Lique F. The temperature variation of the CH + + H reaction rate coefficients: a puzzle finally understood? Phys Chem Chem Phys 2024; 26:21370-21378. [PMID: 39099553 PMCID: PMC11323935 DOI: 10.1039/d4cp01902d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/07/2024] [Accepted: 07/19/2024] [Indexed: 08/06/2024]
Abstract
CH+ was the first molecular ion identified in the interstellar medium and is found to be ubiquitous in interstellar clouds. However, its formation and destruction paths are not well understood, especially at low temperatures. A new theoretical approach based on the canonical variational transition state theory was used to study the H + CH+ reactive collisions. Rate coefficients for formation of C+ ions are calculated as a function of temperature. We considered the participation of a direct path and an indirect path in which the reactants should overcome an entropic barrier to form a van der Waals complex or pass through a CH2+ intermediate complex, respectively. We show that the contribution of both pathways to the formation of C+ has to be taken into account. The new reactive rate coefficients for the title reaction, complemented by reactive data for CH+/CH2+ in the H/H2/He mixture, have been used to simulate the corresponding kinetics experimentally measured using an Atomic Beam 22 Pole Trap apparatus at low temperature. A good agreement with the experimental findings was found at 50 K. At a lower temperature, the model overestimates the formation of C+. This shows that secondary reactions are not responsible for the weak C+ production in the experiments at such temperature. Then, we discuss the possible impact of non-adiabatic effects in the study of the H + CH+ reactive collisions and we found that such effects can be responsible for the decrease of the H + CH+ rate coefficients at low temperature. This study offers an explanation for the disagreement between H + CH+ theoretical and experimental rate coefficients which has been going on for 20 years and highlights the need for performing non-adiabatic studies for this simple chemical reaction.
Collapse
Affiliation(s)
- Rafael A Jara-Toro
- Univ Rennes, CNRS, IPR (Institut de Physique de Rennes), UMR 6251, F-35000 Rennes, France.
| | - Octavio Roncero
- Instituto de Física Fundamental (IFF-CSIC), C.S.I.C., Serrano 123, 28006 Madrid, Spain
| | - François Lique
- Univ Rennes, CNRS, IPR (Institut de Physique de Rennes), UMR 6251, F-35000 Rennes, France.
| |
Collapse
|
3
|
Nakashima H, Nakatsuji H. Potential Energy Curves of the Low-Lying Five 1Σ + and 1Π States of a CH + Molecule Based on the Free Complement - Local Schrödinger Equation Theory and the Chemical Formula Theory. J Chem Theory Comput 2023; 19:6733-6744. [PMID: 37706317 DOI: 10.1021/acs.jctc.3c00645] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/15/2023]
Abstract
The potential energy curves (PECs) of the low-lying five 1Σ+ and 1Π states (X1Σ+, C1Σ+, 31Σ+, A1Π, and D1Π states) of a CH+ molecule, an important interstellar molecule, were calculated by the free complement (FC) - local Schrödinger equation (LSE) theory with the direct local sampling scheme. The FC wave functions were constructed based on the chemical formula theory (CFT), whose local characters correspond to the covalent dissociations: C+(2P°(s2p))) + H(2S) of the X1Σ+ and A1Π states and the ionic dissociations: C(1D(s2p2)) + H+ of the C1Σ+ and D1Π states. All the calculated PECs were obtained with satisfying the chemical accuracy, i.e., error less than 1 kcal/mol, as absolute total energy of the Schrödinger equation without any energy shift. The spectroscopic data calculated from the PECs agreed well with both experimental and other accurate theoretical references. We also analyzed the wave functions using the inverse overlap weights proposed by Gallup et al. with the CFT configurations. For the X1Σ+ and A1Π states, the covalent C+(sp2) and C+(p3) configurations played important roles for bond formation. In the small internuclear distances of the C1Σ+, D1Π, and 31Σ+ states, the covalent character was also dominant as a result of the electron charge transfer from C to H+. Thus, the present FC-LSE results not only are accurate but also can provide chemical understanding according to the CFT.
Collapse
Affiliation(s)
- Hiroyuki Nakashima
- Quantum Chemistry Research Institute, Kyoto Technoscience Center 16, 14 Yoshida Kawaramachi, Sakyo-ku, Kyoto 606-8305, Japan
| | - Hiroshi Nakatsuji
- Quantum Chemistry Research Institute, Kyoto Technoscience Center 16, 14 Yoshida Kawaramachi, Sakyo-ku, Kyoto 606-8305, Japan
| |
Collapse
|
4
|
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]
|
5
|
Lü Y, Zhang C, Wang H, Guo Q, Li Y. An accurate many-body expansion potential energy surface for AlH 2 (2 2A') and quantum dynamics in Al( 3P) + H 2 ( v0 = 0-3, j0 = 0, 2, 4, 6) collisions. Phys Chem Chem Phys 2022; 24:16637-16646. [PMID: 35766326 DOI: 10.1039/d2cp01802k] [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
An accurate potential energy surface is constructed for the excited state of AlH2 by fitting extensive ab initio points calculated at the multi-reference configuration interaction level based on aug-cc-pV(Q+d)Z and aug-cc-pV(5+d)Z basis sets. All the calculated energies are corrected via the many-body expansion method and extrapolated to the complete basis set limit. The various topographic features of the new potential energy surface are investigated to demonstrate the correct behavior of Al(3P) + H2(X1Σg+) and AlH(a3Π) + H(2S) dissociation limits. By employing the time-dependent wave packet approach, the integral scattering cross-sections obtained from the Coriolis coupling calculation and the centrifugal sudden approximation, respectively, are compared in detail and show that the former has a higher effect on the reaction. Moreover, the thermal rate constants for Al(3P) + H2 (v0 = 0-3, j0 = 0, 2, 4, 6) in the temperature range of 0-5000 K are calculated, thereby providing insights into the influence of ro-vibrational quantum numbers on the thermal rate constants.
Collapse
Affiliation(s)
- Yanling Lü
- Department of Physics, Liaoning University, Shenyang 110036, China.
| | - ChengYuan Zhang
- Department of Physics, Liaoning University, Shenyang 110036, China.
| | - Hainan Wang
- Department of Physics, Liaoning University, Shenyang 110036, China.
| | - Qiang Guo
- Department of Physics, Liaoning University, Shenyang 110036, China. .,Lvyuan Institute of Energy & Environmental Science and Technology, Liaoning University, Shenyang, 110036, China
| | - Yongqing Li
- Department of Physics, Liaoning University, Shenyang 110036, China. .,Lvyuan Institute of Energy & Environmental Science and Technology, Liaoning University, Shenyang, 110036, China
| |
Collapse
|
6
|
Li X, Qin Z, Li J, Liu L. An accurate NH 2(X 2A′′) CHIPR potential energy surface via extrapolation to the complete basis set limit and dynamics of the N( 2D) + H 2(X 1Σ+g) reaction. Phys Chem Chem Phys 2022; 24:26564-26574. [DOI: 10.1039/d2cp01961b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
An accurate CHIPR potential energy surface for NH2(X2A′′) is structured to study the N(2D) + H2(X1Σ+g) reaction using the time-dependent wave packet and quasi-classical trajectory method.
Collapse
Affiliation(s)
- Ximing Li
- School of Energy and Power Engineering, Shandong University, 250061, Jinan, China
| | - Zhi Qin
- School of Energy and Power Engineering, Shandong University, 250061, Jinan, China
- Optics and Thermal Radiation Research Center, Institute of Frontier and Interdisciplinary Science, Shandong University, 266237, Qingdao, China
| | - Jing Li
- School of Physics and Physical Engineering, Qufu Normal University, 273165, Qufu, China
| | - Linhua Liu
- School of Energy and Power Engineering, Shandong University, 250061, Jinan, China
- Optics and Thermal Radiation Research Center, Institute of Frontier and Interdisciplinary Science, Shandong University, 266237, Qingdao, China
- School of Energy Science and Engineering, Harbin Institute of Technology, 150001, Harbin, China
| |
Collapse
|
7
|
Konings M, Desrousseaux B, Lique F, Loreau J. Benchmarking an improved statistical adiabatic channel model for competing inelastic and reactive processes. J Chem Phys 2021; 155:104302. [PMID: 34525820 DOI: 10.1063/5.0062388] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Inelastic collisions and elementary chemical reactions proceeding through the formation and subsequent decay of an intermediate collision complex, with an associated deep well on the potential energy surface, pose a challenge for accurate fully quantum mechanical approaches, such as the close-coupling method. In this study, we report on the theoretical prediction of temperature-dependent state-to-state rate coefficients for these complex-mode processes, using a statistical quantum method. This statistical adiabatic channel model is benchmarked by a direct comparison using accurate rate coefficients from the literature for a number of systems (H2 + H+, HD + H+, SH+ + H, and CH+ + H) of interest in astrochemistry and astrophysics. For all of the systems considered, an error of less than factor 2 was found, at least for the dominant transitions and at low temperatures, which is sufficiently accurate for applications in the above mentioned disciplines.
Collapse
Affiliation(s)
- Maarten Konings
- KU Leuven, Division of Quantum Chemistry and Physical Chemistry, Department of Chemistry, Celestijnenlaan 200F, 3001 Leuven, Belgium
| | - Benjamin Desrousseaux
- Université de Rennes 1, CNRS, IPR (Institut de Physique de Rennes) - UMR 6251, F-35000 Rennes, France
| | - François Lique
- Université de Rennes 1, CNRS, IPR (Institut de Physique de Rennes) - UMR 6251, F-35000 Rennes, France
| | - Jérôme Loreau
- KU Leuven, Division of Quantum Chemistry and Physical Chemistry, Department of Chemistry, Celestijnenlaan 200F, 3001 Leuven, Belgium
| |
Collapse
|
8
|
Ma H, Zhang C, Song Y, Ma F, Li Y. Accurate High-Level Ab Initio-Based Global Potential Energy Surface and Quantum Dynamics Calculation for the First Excited State of CH 2. J Phys Chem A 2021; 125:5490-5498. [PMID: 34137628 DOI: 10.1021/acs.jpca.1c02413] [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/29/2022]
Abstract
A full three-dimensional global potential energy surface (PES), covering the whole configuration space, is reported first for the title system by fitting high-level ab initio energies at the multireference configuration interaction level with the aug-cc-pV6Z basis set. In this work, the many-body expansion method is invoked to fit the innate character of the CH2+(12A″) PES. The topographical features are examined in detail based on the new global PES and in accordance with the other calculations from the ab initio energies, which show the correct behavior at the C+(2P) + H2(X1Σg+) and CH+(a3Π) + H(2S) dissociation limits. Using a time-dependent wave packet method, we provide insights into the dynamics behavior for reaction of C+(2P) + H2(X1Σg+) → CH+(a3Π) + H(2S). The integral cross sections and reaction probabilities increase monotonically in terms of the collision energy.
Collapse
Affiliation(s)
- Hongyu Ma
- Department of Physics, Liaoning University, Shenyang 110036, China
| | - Chengyuan Zhang
- Department of Physics, Liaoning University, Shenyang 110036, China
| | - Yuzhi Song
- Collaborative Innovation Center of Light Manipulations and Applications, Shandong Normal University, Jinan 250358, China
| | - Fengcai Ma
- Department of Physics, Liaoning University, Shenyang 110036, China
| | - Yongqing Li
- Department of Physics, Liaoning University, Shenyang 110036, China
| |
Collapse
|
9
|
Hillenbrand PM, Bowen KP, Dayou F, Miller KA, de Ruette N, Urbain X, Savin DW. Experimental study of the proton-transfer reaction C + H 2+ → CH + + H and its isotopic variant (D 2+). Phys Chem Chem Phys 2020; 22:27364-27384. [PMID: 33231243 DOI: 10.1039/d0cp04810k] [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
We report absolute integral cross section (ICS) measurements using a dual-source merged-fast-beams apparatus to study the titular reactions over the relative translational energy range of Er ∼ 0.01-10 eV. We used photodetachment of C- to produce a pure beam of atomic C in the ground electronic 3P term, with statistically populated fine-structure levels. The H2+ and D2+ were formed in an electron impact ionization source, with well known vibrational and rotational distributions. The experimental work is complemented by a theoretical study of the CH2+ electronic system in the reactant and product channels, which helps to clarify the possible reaction mechanisms underlying the ICS measurements. Our measurements provide evidence that the reactions are barrierless and exoergic. They also indicate the apparent absence of an intermolecular isotope effect, to within the total experimental uncertainties. Capture models, taking into account either the charge-induced dipole interaction potential or the combined charge-quadrupole and charge-induced dipole interaction potentials, produce reaction cross sections that lie a factor of ∼4 above the experimental results. Based on our theoretical study, we hypothesize that the reaction is most likely to proceed adiabatically through the 14A' and 14A'' states of CH2+via the reaction C(3P) + H2+(2Σ+g) → CH+(3Π) + H(2S). We also hypothesize that at low collision energies only H2+(v ≤ 2) and D2+(v ≤ 3) contribute to the titular reactions, due to the onset of dissociative charge transfer for higher vibrational v levels. Incorporating these assumptions into the capture models brings them into better agreement with the experimental results. Still, for energies ⪅0.1 eV where capture models are most relevant, the modified charge-induced dipole model yields reaction cross sections with an incorrect energy dependence and lying ∼10% below the experimental results. The capture cross section obtained from the combined charge-quadrupole and charge-induced dipole model better matches the measured energy dependence but lies ∼30-50% above the experimental results. These findings provide important guidance for future quasiclassical trajectory and quantum mechanical treatments of this reaction.
Collapse
|
10
|
Chai S, Wang H, Lü Y, Guo Q, Li Y. Accurate global adiabatic potential energy surfaces for three low-lying electronic states of AlH 2 free radicals. Phys Chem Chem Phys 2020; 22:26544-26551. [PMID: 33200164 DOI: 10.1039/d0cp02939d] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
In order to obtain the all-round molecular properties of the AlH2 system and the corresponding dynamical characteristics of the Al + H2 (v = 0, j = 0) → H + AlH reaction, three significant global adiabatic potential energy surfaces of AlH2 (X2A1, 2B1, and 2B2) free radicals were constructed for the first time. Ab initio energies were calculated under the multi-reference configuration interaction method and the aug-cc-pV(T,Q)Z basis sets; then the ab initio energies were extrapolated to the complete basis sets limit. The three adiabatic potential energy surfaces were constructed by the many-body expansion theory. The maximum root-mean square error was just 50 cm-1, which was small enough to ensure that the potential energy surfaces were accurate. The concerned T-type insertion topographical features, dissociation schemes, C2v geometry reaction mechanisms, and minimum energy curve paths were investigated and are discussed in detail. Several differences from previous studies are also pointed out. Eventually, the integral cross-sections of Al + H2 (v = 0, j = 0) → H + AlH reaction as calculated by quasi-classical trajectory method were employed to predict the dynamical properties of AlH2, providing the most reliable theoretical reference of the dynamical characteristics known thus far for such a reaction. These new potential energy surfaces can be treated as a reliable basis for investigation of the dynamics and as a component for constructing larger aluminum-/hydrogen-containing systems.
Collapse
Affiliation(s)
- Shijie Chai
- Department of Physics, Liaoning University, Shenyang 110036, China.
| | | | | | | | | |
Collapse
|
11
|
Wu H, Duan Z, Chen G. State-resolved quantum mechanical study of the intramolecular isotope effect in the C+ + HD reaction. Chem Phys Lett 2020. [DOI: 10.1016/j.cplett.2020.137783] [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]
|
12
|
Liu D, Zhao J, Wang L, Song Y, Meng Q, Zhang L. Exploring reaction mechanism and vibrational excitation effect in H + CH( v,j = 0) reaction. Chem Phys Lett 2020. [DOI: 10.1016/j.cplett.2020.137398] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
13
|
Chai SJ, Ma HY, Lü YL, Liu JY, Li YQ. Accurate global adiabatic potential energy surface for the ground state of AlH 2+by extrapolation to the complete basis set limit. Mol Phys 2020. [DOI: 10.1080/00268976.2019.1655599] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
Affiliation(s)
- S. J. Chai
- Department of Physics, Liaoning University, Shenyang, People’s Republic of China
| | - H. Y. Ma
- Department of Physics, Liaoning University, Shenyang, People’s Republic of China
| | - Y. L. Lü
- Department of Physics, Liaoning University, Shenyang, People’s Republic of China
| | - J. Y. Liu
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian, People’s Republic of China
| | - Y. Q. Li
- Department of Physics, Liaoning University, Shenyang, People’s Republic of China
| |
Collapse
|
14
|
Effects of rovibrational excitation of LiH on the LiH depletion and H exchange channels for the reaction H ( 2S) + LiH (X 1Σ+) on a new potential energy surface. Chem Phys Lett 2020. [DOI: 10.1016/j.cplett.2019.137043] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
15
|
Zhang L, Yue D, Zhao J, Song Y, Meng Q. Potential energy curves, spectroscopic constants, and vibrational energy levels of CS +( X2Σ +/ A2Π). Mol Phys 2020. [DOI: 10.1080/00268976.2019.1597199] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
Affiliation(s)
- Lulu Zhang
- School of Science, Shandong Jiaotong University, Jinan, China
| | - Daguang Yue
- School of Science, Shandong Jiaotong University, Jinan, China
| | - Juan Zhao
- School of Science, Shandong Jiaotong University, Jinan, China
| | - Yuzhi Song
- Shandong Province Key Laboratory of Medical Physics and Image Processing Technology, School of Physics and Electronics, Shandong Normal University, Jinan, China
| | - Qingtian Meng
- Shandong Province Key Laboratory of Medical Physics and Image Processing Technology, School of Physics and Electronics, Shandong Normal University, Jinan, China
| |
Collapse
|
16
|
Lü YL, Chai SJ, Ma HY, Gao SB, Li YQ. Globally accurate potential energy surface for PH2+ (11 A′) by using the switching function formalism. Mol Phys 2019. [DOI: 10.1080/00268976.2019.1688875] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
Affiliation(s)
- Y. L. Lü
- School of Physics, Liaoning University, Shenyang, People’s Republic of China
| | - S. J. Chai
- School of Physics, Liaoning University, Shenyang, People’s Republic of China
| | - H. Y. Ma
- School of Physics, Liaoning University, Shenyang, People’s Republic of China
| | - S. B. Gao
- School of Physics and Electronics, Shandong Normal University, Jinan, People’s Republic of China
| | - Y. Q. Li
- School of Physics, Liaoning University, Shenyang, People’s Republic of China
| |
Collapse
|
17
|
Gao F, Zhang L, Zhao W, Meng Q, Song Y. Accurate global potential energy surface for SiH 2 +(X 2A 1) and quantum dynamics of related reaction H( 2S) + SiH +(X 1Σ +). J Chem Phys 2019; 150:224304. [PMID: 31202249 DOI: 10.1063/1.5088637] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
With the many-body expansion method, an accurate global potential energy surface (PES) is constructed for SiH2 +(X2A1) by mapping 4762 ab initio energy points calculated on the multireference configuration interaction level including Davidson corrections with aug-cc-pV6Z Dunning's basis set. The dissociation energies and equilibrium geometries of SiH+(X1Σ+) and H2(X1Σg +) agree well with the experimental results. The topographical characteristics of all stationary points for the SiH2 +(X2A1) PES are discussed in detail and compared with other theoretical and experimental results. In order to verify the validity and usability of the present PES, the dynamics calculations based on the Chebyshev quantum wave packet method are performed for the H(S2)+SiH+(X1Σ+)→Si+(P2)+H2(X1Σg +) reaction. The probabilities, the total integral cross sections, and the rate constants are computed, and the analogies with the corresponding ones of reaction H(S2) + CH+(X1Σ+)→C+(P2) + H2(X1Σg +) are also made. The reasonable dynamical behavior throughout the entire configuration space indicates that the PES is suitable for relevant dynamics investigations and serves as a building block for constructing the PES of larger molecular systems containing Si+/H.
Collapse
Affiliation(s)
- Feng Gao
- School of Physics and Electronics, Shandong Normal University, Jinan 250358, People's Republic of China
| | - Lulu Zhang
- School of Science, Shandong Jiaotong University, Jinan 250357, People's Republic of China
| | - Wenli Zhao
- School of Information Science and Engineering, Shandong Agricultural University, Taian 271018, People's Republic of China
| | - Qingtian Meng
- School of Physics and Electronics, Shandong Normal University, Jinan 250358, People's Republic of China
| | - Yuzhi Song
- School of Physics and Electronics, Shandong Normal University, Jinan 250358, People's Republic of China
| |
Collapse
|
18
|
Zhang L, Gao S, Meng Q, Pan J, Song Y. Accurate potential energy surface of H 2S +( X 2 A″) via extrapolation to the complete basis set limit and its use in dynamics study of S + ( D 2 ) + H 2 ( X 1 Σ g + ) reaction. J Chem Phys 2018; 149:154303. [PMID: 30342440 DOI: 10.1063/1.5046315] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The single-sheeted potential energy surface (PES) of H 2 S + ( X 2 A ' ' ) is developed based on the ab initio energies calculated by the multi-reference configuration interaction method including the Davidson correction. All the ab initio energies are first calculated using aug-cc-pVQdZ and aug-cc-pV5dZ basis sets, which are then extrapolated to the complete basis set (CBS) limit. A switching function is developed to model the transition of S + D 2 to S + S 4 . The many-body expansion formalism is employed to obtain the H 2 S + ( X 2 A ' ' ) PES by fitting such CBS energies and the root-mean square derivation is 0.0367 eV. The topographical features of the present PES are examined in detail, which are well consistent with previous studies. The quasiclassical trajectory method is subsequently utilized to study the S + D 2 + H 2 ( X 1 Σ g + ) → S H + ( X 3 Σ - ) + H ( S 2 ) reaction. The capture time, integral cross sections, and rovibrational distributions are calculated. By examining the capture time, it can be concluded that the title reaction is mainly controlled by the indirect mechanism for lower collision energies, while the direct and indirect mechanisms coexist and the latter plays a dominant role for higher collision energies.
Collapse
Affiliation(s)
- Lulu Zhang
- Shandong Province Key Laboratory of Medical Physics and Image Processing Technology, School of Physics and Electronics, Shandong Normal University, 250358 Jinan, China
| | - Shoubao Gao
- Shandong Province Key Laboratory of Medical Physics and Image Processing Technology, School of Physics and Electronics, Shandong Normal University, 250358 Jinan, China
| | - Qingtian Meng
- Shandong Province Key Laboratory of Medical Physics and Image Processing Technology, School of Physics and Electronics, Shandong Normal University, 250358 Jinan, China
| | - Jie Pan
- Shandong Province Key Laboratory of Medical Physics and Image Processing Technology, School of Physics and Electronics, Shandong Normal University, 250358 Jinan, China
| | - Yuzhi Song
- Shandong Province Key Laboratory of Medical Physics and Image Processing Technology, School of Physics and Electronics, Shandong Normal University, 250358 Jinan, China
| |
Collapse
|
19
|
Yang Z, Yuan J, Wang S, Chen M. Global diabatic potential energy surfaces for the BeH 2 + system and dynamics studies on the Be +( 2P) + H 2(X 1Σ g +) → BeH +(X 1Σ +) + H( 2S) reaction. RSC Adv 2018; 8:22823-22834. [PMID: 35539737 PMCID: PMC9081383 DOI: 10.1039/c8ra04305a] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/21/2018] [Revised: 07/13/2018] [Accepted: 06/14/2018] [Indexed: 11/21/2022] Open
Abstract
The Be+(2P) + H2(X1Σg +) → BeH+(X1Σ+) + H(2S) reaction has great significance for studying diabatic processes and ultracold chemistry. The first global diabatic potential energy surfaces (PESs) which are correlated with the lowest two adiabatic states 12A' and 22A' of the BeH2 + system are constructed by using the neural network method. Ab initio energy points are calculated using the multi-reference configuration interaction method with the Davidson correction and AVQZ basis set. The diabatic energies are obtained from the transformation of ab initio data based on the dipole moment operators. The topographical characteristics of the diabatic PESs are described in detail, and the positions of crossing between the V d 11 and V d 22 are pinpointed. On new diabatic PESs, the time-dependent quantum wave packet method is carried out to study the mechanism of the title reaction. The results of dynamics calculations indicate the reaction has no threshold and the product BeH+ is excited to high vibrational states easily. In addition, the product BeH+ tends to backward scattering at most collision energies.
Collapse
Affiliation(s)
- Zijiang Yang
- Key Laboratory of Materials Modification by Laser, Electron, and Ion Beams (Ministry of Education), School of Physics, Dalian University of Technology Dalian 116024 P. R. China
| | - Jiuchuang Yuan
- Key Laboratory of Materials Modification by Laser, Electron, and Ion Beams (Ministry of Education), School of Physics, Dalian University of Technology Dalian 116024 P. R. China
| | - Shufen Wang
- Key Laboratory of Materials Modification by Laser, Electron, and Ion Beams (Ministry of Education), School of Physics, Dalian University of Technology Dalian 116024 P. R. China
| | - Maodu Chen
- Key Laboratory of Materials Modification by Laser, Electron, and Ion Beams (Ministry of Education), School of Physics, Dalian University of Technology Dalian 116024 P. R. China
| |
Collapse
|
20
|
Sundaram P, Padmanaban R. Dynamical resonances of the deuterated CH 2+ complex in the electronic ground state: A quantum wavepacket study. J Chem Phys 2018; 148:164306. [PMID: 29716227 DOI: 10.1063/1.5008783] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
We here investigate the effects of isotopic substituents on the vibrational energy levels of the CH2+ complex, supported by the electronic ground (1 2A') potential energy surface (PES) of the H + CH+ reaction. We calculate the transition state spectrum by Fourier transforming the time-autocorrelation function of the initial wavepacket (WP) chosen in the interaction region of the PES. Using the time-dependent WP approach, the dynamical resonances are identified as bound and quasibound in nature, and they are characterized in terms of the eigenfunctions and lifetimes. The present work on the isotopic variants [CHD+(CDH+) and CD2+] is compared with our earlier work [P. Sundaram et al., Phys. Chem. Chem. Phys. 19, 20172 (2017)] on the parent CH2+ species. The isotopic variants reveal a large number of peaks in the spectra and the eigenfunctions exhibit the systematic nodal progressions and periodic orbits, the same as in CH2+. While the CD2+ complex exactly mimics the resonance behaviors (local and hyperspherical modes) of the bound and quasibound CH2+ complex, the CHD+(CDH+) complex reveals only the local mode behaviors at low energies and significantly less number of resonance structures at high energies. Lifetime analysis of the isotopic variants implies that the CD2+ complex survives much longer than the CHD+(CDH+) complex and concludes the work by noting the following order in the decay profile of the deuterated CH2+ resonances as CH2+>CHD+(CDH+) >CD2+.
Collapse
Affiliation(s)
- P Sundaram
- Department of Chemistry, Pondicherry University, Puducherry 605 014, India
| | - R Padmanaban
- Department of Chemistry, Pondicherry University, Puducherry 605 014, India
| |
Collapse
|
21
|
Guo L, Ma H, Zhang L, Song Y, Li Y. Accurate global potential energy surface for the ground state of CH 2 + by extrapolation to the complete basis set limit. RSC Adv 2018; 8:13635-13642. [PMID: 35539314 PMCID: PMC9079815 DOI: 10.1039/c8ra02228c] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/14/2018] [Accepted: 03/25/2018] [Indexed: 11/21/2022] Open
Abstract
A full three-dimensional global potential energy surface is reported for the ground state of CH2 + by fitting accurate multireference configuration interaction energies calculated using aug-cc-pVQZ and aug-cc-pV5Z basis sets with extrapolation of the electron correlation energy to the complete basis set limit. The topographical characteristics have been compared in detail with a potential energy surface of the same type recently reported [J. Chem. Phys., 2015, 142, 124302] based on a least-squares fit to accurate high level ab initio MRCI(Q) energies, calculated using AV6Z basis set. The new three-dimensional global potential energy surface is then used in quasiclassical trajectory calculations for H(2S) + CH+(X 1Σ+) → C+(2P) + H2(X 1Σg +) reaction. The integral cross sections, differential cross sections and the rate coefficients have been computed. A comparison shows that our potential energy surface can be applied to any type of dynamic study.
Collapse
Affiliation(s)
- Lu Guo
- Department of Physics, Liaoning University Shenyang 110036 China
| | - Hongyu Ma
- Department of Physics, Liaoning University Shenyang 110036 China
| | - Lulu Zhang
- School of Physics and Electronics, Shandong Normal University Jinan 250014 China
| | - Yuzhi Song
- School of Physics and Electronics, Shandong Normal University Jinan 250014 China
| | - Yongqing Li
- Department of Physics, Liaoning University Shenyang 110036 China
| |
Collapse
|
22
|
Guo J, Zhang AJ, Zhou Y, Liu JY, Jia JF, Wu HS. A time-dependent quantum dynamical study of the C+(2P) + H2(D2,HD) → CH+(CD+) + H(D) reaction. Chem Phys Lett 2017. [DOI: 10.1016/j.cplett.2017.09.056] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
23
|
Sundaram P, Manivannan V, Padmanaban R. Dynamics and resonances of the H( 2S) + CH +(X 1Σ +) reaction in the electronic ground state: a detailed quantum wavepacket study. Phys Chem Chem Phys 2017; 19:20172-20187. [PMID: 28726890 DOI: 10.1039/c7cp03110f] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Initial state-selected and energy resolved channel-specific reaction probabilities, integral cross sections and thermal rate constants of the H(2S) + CH+(X1Σ+) reaction are calculated within the coupled states approximation by a time-dependent wave packet propagation method. The new ab initio global potential energy surface (PES) of the electronic ground state (1 2A') of the system, recently reported by Li et al. [J. Chem. Phys., 2015, 142, 124302], is employed for this purpose. All partial wave contributions up to the total angular momentum J = 60 are considered to obtain the converged integral reaction cross section up to a collision energy of 1.0 eV. Thermal rate constants are calculated by averaging the reaction cross sections over the Boltzmann distribution of energies and compared with the available theoretical and experimental results for the temperature range 10-1000 K. Investigation of the channel-specific reaction attributes shows that the H abstraction (CH+ destruction) channel is highly favored over the H exchange channel. The effect of rotational and vibrational excitations of the CH+ reagent on the dynamics is also studied. The resonances formed during the course of the reaction are also identified by calculating the transition state spectrum and characterized in terms of the eigenfunctions and lifetimes. More than 260 vibrational levels are obtained and their eigenfunctions are calculated, which are represented in terms of the nodal assignments and the eigenenergies. They reveal both the local and hyperspherical behavior for the bound and quasibound states of the CH2+ complex in the ground 1 2A' surface. The lifetime analysis of the quasibound states indicates that the CH2+ resonances survive for as long as ∼400 fs at high energies (E ∼ 2.0 eV) and are expected to decay faster with further increasing energy. Finally, the type of mechanism for the formation of the product (C+ + H2) is elucidated.
Collapse
Affiliation(s)
- P Sundaram
- Department of Chemistry, Pondicherry University, Puducherry - 605 014, India.
| | | | | |
Collapse
|
24
|
Faure A, Halvick P, Stoecklin T, Honvault P, Epée Epée MD, Mezei JZ, Motapon O, Schneider IF, Tennyson J, Roncero O, Bulut N, Zanchet A. State-to-state chemistry and rotational excitation of CH + in photon-dominated regions. MONTHLY NOTICES OF THE ROYAL ASTRONOMICAL SOCIETY 2017; 469:612-620. [PMID: 28690343 PMCID: PMC5500105 DOI: 10.1093/mnras/stx892] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
We present a detailed theoretical study of the rotational excitation of CH+ due to reactive and nonreactive collisions involving C+(2P), H2, CH+, H and free electrons. Specifically, the formation of CH+ proceeds through the reaction between C+(2P) and H2(νH2 = 1, 2), while the collisional (de)excitation and destruction of CH+ is due to collisions with hydrogen atoms and free electrons. State-to-state and initial-state-specific rate coefficients are computed in the kinetic temperature range 10-3000 K for the inelastic, exchange, abstraction and dissociative recombination processes using accurate potential energy surfaces and the best scattering methods. Good agreement, within a factor of 2, is found between the experimental and theoretical thermal rate coefficients, except for the reaction of CH+ with H atoms at kinetic temperatures below 50 K. The full set of collisional and chemical data are then implemented in a radiative transfer model. Our Non-LTE calculations confirm that the formation pumping due to vibrationally excited H2 has a substantial effect on the excitation of CH+ in photon-dominated regions. In addition, we are able to reproduce, within error bars, the far-infrared observations of CH+ toward the Orion Bar and the planetary nebula NGC 7027. Our results further suggest that the population of νH2 = 2 might be significant in the photon-dominated region of NGC 7027.
Collapse
Affiliation(s)
- A. Faure
- Univ. Grenoble Alpes, CNRS, IPAG, F-38000 Grenoble, France
| | - P. Halvick
- Univ. Bordeaux, CNRS, ISM, F-33400 Talence, France
| | - T. Stoecklin
- Univ. Bordeaux, CNRS, ISM, F-33400 Talence, France
| | - P. Honvault
- Univ. Bourgogne Franche-Comté, Laboratoire ICB, F-21078 Dijon, France
| | - M. D. Epée Epée
- UFD Mathématiques, Informatique Appliqúee et Physique Fondamentale, University of Douala, P. O. Box 24157, Douala, Cameroon
| | - J. Zs. Mezei
- Univ. Normandie, CNRS, LOMC, F-76058 Le Havre, France
- LSPM, Univ. Paris 13, 99 avenue Jean-Baptiste Clément, F-93430 Villetaneuse, France
- Univ. Paris-Sud, CNRS, Laboratoire Aimé Cotton, F-91405 Orsay, France
- Institute of Nuclear Research of the Hungarian Academy of Sciences, P.O. Box 51, Debrecen H-4001, Hungary
| | - O. Motapon
- UFD Mathématiques, Informatique Appliqúee et Physique Fondamentale, University of Douala, P. O. Box 24157, Douala, Cameroon
- University of Maroua, Faculty of Science, P. O. Box 814 Maroua, Cameroon
| | - I. F. Schneider
- Univ. Normandie, CNRS, LOMC, F-76058 Le Havre, France
- Univ. Paris-Sud, CNRS, Laboratoire Aimé Cotton, F-91405 Orsay, France
| | - J. Tennyson
- Department of Physics and Astronomy, University College, London, Gower St., London WC1E 6BT, UK
| | - O. Roncero
- Instituto de Física Fundamental, CSIC, C/ Serrano, 123, E-28006 Madrid, Spain
| | - N. Bulut
- Firat University, Department of Physics, 23169 Elazig̃, Turkey
| | - A. Zanchet
- Instituto de Física Fundamental, CSIC, C/ Serrano, 123, E-28006 Madrid, Spain
| |
Collapse
|
25
|
Song YZ, Zhang LL, Gao SB, Meng QT. Globally accurate potential energy surface for the ground-state HCS(X 2A') and its use in reaction dynamics. Sci Rep 2016; 6:37734. [PMID: 27898106 PMCID: PMC5127192 DOI: 10.1038/srep37734] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/19/2016] [Accepted: 10/31/2016] [Indexed: 11/09/2022] Open
Abstract
A globally accurate many-body expansion potential energy surface is reported for HCS(X2A') by fitting a wealth of accurate ab initio energies calculated at the multireference configuration interaction level using aug-cc-pVQZ and aug-cc-pV5Z basis sets via extrapolation to the complete basis set limit. The topographical features of the present potential energy surface are examined in detail and is in good agreement with the raw ab initio results, as well as other theoretical results available in literatures. By utilizing the potential energy surface of HCS(X2A'), the dynamic studies of the C(3P) + SH(X2Π) → H(2S) + CS(X1∑+) reaction has been carried out using quasi-classical trajectory method.
Collapse
Affiliation(s)
- Yu-Zhi Song
- School of Physics and Electronics, Shandong Normal University, Jinan 250014, China
| | - Lu-Lu Zhang
- School of Physics and Electronics, Shandong Normal University, Jinan 250014, China
| | - Shou-Bao Gao
- School of Physics and Electronics, Shandong Normal University, Jinan 250014, China
| | - Qing-Tian Meng
- School of Physics and Electronics, Shandong Normal University, Jinan 250014, China
| |
Collapse
|
26
|
He X, Lv S, Hayat T, Han K. Potential Energy Surfaces for the First Two Lowest-Lying Electronic States of the LiH2+ System, and Dynamics of the H+ + LiH ⇌ H2+ + Li + Reactions. J Phys Chem A 2016; 120:2459-70. [PMID: 27022663 DOI: 10.1021/acs.jpca.6b02007] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Xiaohu He
- School
of Chemical Engineering, Dalian University of Technology, Dalian 116023, China
- State
Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of
Chemical Physics, Chinese Academy of Science, Dalian 116023, China
| | - Shuangjiang Lv
- State
Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of
Chemical Physics, Chinese Academy of Science, Dalian 116023, China
| | - Tasawar Hayat
- Department
of Mathematics, Quaid-I-Azam University, Islamabad 44000, Pakistan
- Nonlinear
Analysis and Applied Mathematics (NAAM) Research Group, Department
of Mathematics, Faculty of Science, King Abdulaziz University, Jeddah 21589, Saudi Arabia
| | - Keli Han
- State
Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of
Chemical Physics, Chinese Academy of Science, Dalian 116023, China
| |
Collapse
|
27
|
Werfelli G, Halvick P, Honvault P, Kerkeni B, Stoecklin T. Low temperature rate coefficients of the H + CH(+) → C(+) + H2 reaction: New potential energy surface and time-independent quantum scattering. J Chem Phys 2015; 143:114304. [PMID: 26395702 DOI: 10.1063/1.4931103] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The observed abundances of the methylidyne cation, CH(+), in diffuse molecular clouds can be two orders of magnitude higher than the prediction of the standard gas-phase models which, in turn, predict rather well the abundances of neutral CH. It is therefore necessary to investigate all the possible formation and destruction processes of CH(+) in the interstellar medium with the most abundant species H, H2, and e(-). In this work, we address the destruction process of CH(+) by hydrogen abstraction. We report a new calculation of the low temperature rate coefficients for the abstraction reaction, using accurate time-independent quantum scattering and a new high-level ab initio global potential energy surface including a realistic model of the long-range interaction between the reactants H and CH(+). The calculated thermal rate coefficient is in good agreement with the experimental data in the range 50 K-800 K. However, at lower temperatures, the experimental rate coefficient takes exceedingly small values which are not reproduced by the calculated rate coefficient. Instead, the latter rate coefficient is close to the one given by the Langevin capture model, as expected for a reaction involving an ion and a neutral species. Several recent theoretical works have reported a seemingly good agreement with the experiment below 50 K, but an analysis of these works show that they are based on potential energy surfaces with incorrect long-range behavior. The experimental results were explained by a loss of reactivity of the lowest rotational states of the reactant; however, the quantum scattering calculations show the opposite, namely, a reactivity enhancement with rotational excitation.
Collapse
Affiliation(s)
- Ghofran Werfelli
- Université de Bordeaux, Institut des Sciences Moléculaires, CNRS UMR 5255, 33405 Talence Cedex, France
| | - Philippe Halvick
- Université de Bordeaux, Institut des Sciences Moléculaires, CNRS UMR 5255, 33405 Talence Cedex, France
| | - Pascal Honvault
- Université de Bourgogne Franche-Comté, Laboratoire Interdisciplinaire Carnot de Bourgogne, CNRS UMR 6303, 21078 Dijon Cedex, France
| | - Boutheïna Kerkeni
- Faculté des Sciences de Tunis, Université de Tunis El-Manar, Laboratoire de Physique de la Matière Condensée, 2092 Tunis, Tunisia
| | - Thierry Stoecklin
- Université de Bordeaux, Institut des Sciences Moléculaires, CNRS UMR 5255, 33405 Talence Cedex, France
| |
Collapse
|