1
|
Chen J, Zhang H, Zhou L, Hu X, Xie D. New accurate diabatic potential energy surfaces for the two lowest 1A'' states of H 2S and photodissociation dynamics in its first absorption band. Phys Chem Chem Phys 2023; 25:26032-26042. [PMID: 37750311 DOI: 10.1039/d3cp03026a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/27/2023]
Abstract
In this work, state-to-state photodissociation dynamics of H2S in its first absorption band has been studied quantum mechanically with a new set of coupled potential energy surfaces (PESs) for the first two 1A'' excited states, which were developed at the explicitly correlated internally contracted multi-reference configuration interaction level with the cc-pVQZ-F12 basis set and a large active space. The calculated absorption spectrum, product state distributions, and angular distributions are in excellent agreement with available experimental data, validating the accuracy of the PESs and the non-adiabatic couplings. Detailed analysis of the dynamics reveals that there are strong non-adiabatic couplings between the bound 11B1 and dissociative 11A2 states around the Franck-Condon region, leading to very fast predissociation to ro-vibrationally cold SH(X̃) fragments, during which marginal angular anisotropy of the PESs is involved. This study provides quantitatively accurate characterization of the electronic structure and detailed fragmentation dynamics of this prototypical photodissociation system, which is desirable for improving astrochemical modelling.
Collapse
Affiliation(s)
- Junjie Chen
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Hanzi Zhang
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Linsen Zhou
- Institute of Materials, China Academy of Engineering Physics, Mianyang 621907, China.
| | - Xixi Hu
- Kuang Yaming Honors School, Institute for Brain Sciences, Nanjing University, Nanjing 210023, China.
- Hefei National Laboratory, Hefei 230088, China
| | - Daiqian Xie
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
- Hefei National Laboratory, Hefei 230088, China
| |
Collapse
|
2
|
Yang D, Guo H, Xie D. Recent advances in quantum theory on ro-vibrationally inelastic scattering. Phys Chem Chem Phys 2023; 25:3577-3594. [PMID: 36602236 DOI: 10.1039/d2cp05069b] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
Molecular collisions are of fundamental importance in understanding intermolecular interaction and dynamics. Its importance is accentuated in cold and ultra-cold collisions because of the dominant quantum mechanical nature of the scattering. We review recent advances in the time-independent approach to quantum mechanical characterization of non-reactive scattering in tetratomic systems, which is ideally suited for large collisional de Broglie wavelengths characteristic in cold and ultracold conditions. We discuss quantum scattering algorithms between two diatoms and between a triatom and an atom and their implementation, as well as various approximate schemes. They not only enable the characterization of collision dynamics in realistic systems but also serve as benchmarks for developing more approximate methods.
Collapse
Affiliation(s)
- Dongzheng Yang
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, USA.
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, USA.
| | - Daiqian Xie
- Institute of Theoretical and Computational Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China. .,Hefei National Laboratory, Hefei 230088, China
| |
Collapse
|
3
|
Wang Y, Guo H, Yarkony DR. Internal conversion and intersystem crossing dynamics based on coupled potential energy surfaces with full geometry-dependent spin-orbit and derivative couplings. Nonadiabatic photodissociation dynamics of NH 3(A) leading to the NH(X 3Σ -, a 1Δ) + H 2 channel. Phys Chem Chem Phys 2022; 24:15060-15067. [PMID: 35696936 DOI: 10.1039/d2cp01271e] [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
We simulate the photodissociation of NH3 originating from its first excited singlet state S1 into the NH2 + H (radical) and NH + H2 (molecular) channels. The states considered are the ground singlet state S0, the first excited singlet state S1 and the lowest-lying triplet state T1, which permit for the first time a uniform treatment of the internal conversion and intersystem crossing. The simulations are based on a diabatic potential energy matrix (DPEM) of S0, S1 coupled by a conical intersection seam, as well as a potential energy surface (PES) for T1 coupled by spin-orbit coupling (SOC) to the two singlet states. The DPEM and PES are fitted to ab initio electronic structure data (ESD) including energies, energy gradients, and derivative couplings. The DPEM also defines an adiabatic to diabatic state (AtD) transformation, which is used to transform the singular adiabatic SOC into a smooth function of the nuclear coordinates in the diabatic representation, allowing the diabatic SOC to be fit to an analytical functional form. ESD and SOC data obtained from these surfaces can serve as input for either quantum or semi-classical characterization of the nonadiabatic dynamics. Using the SHARC suite of programs, nonadiabatic simulations based on over 40 000 semi-classical trajectories assess the convergence of our results. The production of NH + H2 is not direct, but is only achieved through a quasi-statistical dissociation mechanism after internal conversion to the ground electronic state. This leads to a much lower yield comparing with the main NH2 + H channel. The NH(X3Σ_) radical produced through the intersystem crossing from S0 to T1 is rare (∼0.2%) compared to NH(a1Δ) due to the process being spin forbidden.
Collapse
Affiliation(s)
- Yuchen Wang
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, USA.
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, USA
| | - David R Yarkony
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, USA.
| |
Collapse
|
4
|
Li Y, Liu J, Li J, Zhai Y, Yang J, Qu Z, Li H. A new permutation-symmetry-adapted machine learning diabatization procedure and its application in MgH 2 system. J Chem Phys 2021; 155:214102. [PMID: 34879675 DOI: 10.1063/5.0072004] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
This work introduces a new permutation-symmetry-adapted machine learning diabatization procedure, termed the diabatization by equivariant neural network (DENN). In this approach, the permutation symmetric and anti-symmetric elements in diabatic potential energy metrics (DPEMs) were simultaneously simulated by the equivariant neural network. The diabatization by deep neural network scheme was adopted for machine learning diabatization, and non-zero diabatic coupling was included to increase accuracy in the near degenerate region. Based on DENN, the global DPEMs for 11A' and 21A' states of MgH2 have been constructed. To the best of our knowledge, these are the first global DPEMs for the MgH2 system. The root-mean-square-errors (RMSEs) for diagonal elements (H11 and H22) and the off-diagonal element (H12) around the conical intersection region were 5.824, 5.307, and 5.796 meV, respectively. The RMSEs of global adiabatic energies for two adiabatic states were 4.613 and 12.755 meV, respectively. The spectroscopic calculations of the 11A' state in the linear HMgH region are in good agreement with the experiment and previous theoretical results. The differences between calculated frequencies and corresponding experiment values are 1.38 and 1.08 cm-1 for anti-symmetric stretching fundamental vibrational frequency and first overtone. The global DPEMs obtained in this work should be useful for further quantum mechanics dynamic simulations on the MgH2 system.
Collapse
Affiliation(s)
- You Li
- Institute of Theoretical Chemistry, College of Chemistry, Jilin University, 2519 Jiefang Road, Changchun 130023, People's Republic of China
| | - Jingmin Liu
- Institute of Theoretical Chemistry, College of Chemistry, Jilin University, 2519 Jiefang Road, Changchun 130023, People's Republic of China
| | - Jiarui Li
- Institute of Theoretical Chemistry, College of Chemistry, Jilin University, 2519 Jiefang Road, Changchun 130023, People's Republic of China
| | - Yu Zhai
- Institute of Theoretical Chemistry, College of Chemistry, Jilin University, 2519 Jiefang Road, Changchun 130023, People's Republic of China
| | - Jitai Yang
- Institute of Theoretical Chemistry, College of Chemistry, Jilin University, 2519 Jiefang Road, Changchun 130023, People's Republic of China
| | - Zexing Qu
- Institute of Theoretical Chemistry, College of Chemistry, Jilin University, 2519 Jiefang Road, Changchun 130023, People's Republic of China
| | - Hui Li
- Institute of Theoretical Chemistry, College of Chemistry, Jilin University, 2519 Jiefang Road, Changchun 130023, People's Republic of China
| |
Collapse
|
5
|
Guan Y, Xie C, Yarkony DR, Guo H. High-fidelity first principles nonadiabaticity: diabatization, analytic representation of global diabatic potential energy matrices, and quantum dynamics. Phys Chem Chem Phys 2021; 23:24962-24983. [PMID: 34473156 DOI: 10.1039/d1cp03008f] [Citation(s) in RCA: 28] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
Nonadiabatic dynamics, which goes beyond the Born-Oppenheimer approximation, has increasingly been shown to play an important role in chemical processes, particularly those involving electronically excited states. Understanding multistate dynamics requires rigorous quantum characterization of both electronic and nuclear motion. However, such first principles treatments of multi-dimensional systems have so far been rather limited due to the lack of accurate coupled potential energy surfaces and difficulties associated with quantum dynamics. In this Perspective, we review recent advances in developing high-fidelity analytical diabatic potential energy matrices for quantum dynamical investigations of polyatomic uni- and bi-molecular nonadiabatic processes, by machine learning of high-level ab initio data. Special attention is paid to methods of diabatization, high fidelity construction of multi-state coupled potential energy surfaces and property surfaces, as well as quantum mechanical characterization of nonadiabatic nuclear dynamics. To illustrate the tremendous progress made by these new developments, several examples are discussed, in which direct comparison with quantum state resolved measurements led to either confirmation of the observation or sometimes reinterpretation of the experimental data. The insights gained in these prototypical systems greatly advance our understanding of nonadiabatic dynamics in chemical systems.
Collapse
Affiliation(s)
- Yafu Guan
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, USA.
| | - Changjian Xie
- Institute of Modern Physics, Northwest University, Xi'an, Shaanxi 710069, China.
| | - David R Yarkony
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, USA.
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico, 87131, USA.
| |
Collapse
|
6
|
Wang Y, Yarkony DR. Conical intersection seams in spin-orbit coupled systems with an even number of electrons: A numerical study based on neural network fit surfaces. J Chem Phys 2021; 155:174115. [PMID: 34742185 DOI: 10.1063/5.0067660] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In this work, we consider the existence and topography of seams of conical intersections (CIs) for two key singlet-triplet systems, including a uniformly scaled spin-orbit interaction. The basic one triplet and one singlet state system denoted as (S0,T1) and the two singlets and one triplet system denoted as (S0,S1,T1) are treated. Essential to this analysis are realistic electronic structure data taken from a recently reported neural network fit for the 1,21A and 13A states of NH3, including Hsf (spin-free) and Hso (spin-orbit) surfaces derived from high quality ab initio wavefunctions. Three types of seams for the (S0,S1,T1) system are reported, which depend on the choice of the electronic Hamiltonian, He. The nonrelativistic CI seam [He = Hsf, (S0,S1)], the energy minimized nonrelativistic singlet-triplet intersection seam [He = Hsf, (S0,T1)], and the fully relativistic seam in the spin-diabatic representation (He = Htot = Hsf + Hso) are reported as functions of R(N-H). The derivative couplings are computed using He = Htot and Hsf from the fit data. The line integral of the derivative coupling is employed to juxtapose the geometric phase in the relativistic, He = Htot, and nonrelativistic, He = Hsf, cases. It is found for the (S0,T1) system that there is no CI in the spin-adiabatic representation, while for the (S0,S1,T1) system, CI can only be formed for two pairs of spin-adiabatic electronic states. The geometric phase effect thus needs to be handled with care when it comes to spin-nonconserving dynamics simulations.
Collapse
Affiliation(s)
- Yuchen Wang
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, USA
| | - David R Yarkony
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, USA
| |
Collapse
|
7
|
Pershin AA, Torbin AP, Mikheyev PA, Kaiser RI, Mebel AM, Azyazov VN. Ozone destruction due to the recombination of oxygen atoms. J Chem Phys 2021; 155:164307. [PMID: 34717353 DOI: 10.1063/5.0064361] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Kinetics of ozone destruction due to the recombination of oxygen atoms produced by pulsed 266 nm laser photolysis of O3/M (M = CO2 and/or N2) mixtures was studied using the absorption and emission spectroscopy to follow time evolutions of O3 and electronically excited molecules O2* formed in the recombination process 2O(3P) + M → O2* + M. An unexpected high ozone destruction rate was observed when O2* was present in the system. The kinetic model developed for the oxygen nightglow on the terrestrial planets was adapted to interpret the detected temporal profiles of the ozone number density and the O2* emission intensities. It was deduced that the vibrationally excited singlet delta oxygen molecule O2(a1Δ, υ) formed in the secondary processes reacts efficiently with ozone in the process O2(a1Δ, υ ≥ 3) + O3 → 2O2 + O, and the rate constant of this process was estimated to be 3 × 10-11 cm3 s-1. Ab initio calculations at the CASPT2(14, 12)/cc-pVTZ/UωB97XD/cc-pVTZ level of theory were applied to find the reaction pathway from the reactants to products on the O5 potential energy surface. These calculations revealed that the O2(a1Δ) + O3 reaction is likely to proceed via singlet-triplet intersystem crossing exhibiting an energy barrier of 9.6 kcal/mol, which lies between two and three quanta of vibrational excitation of O2(a1Δ), and hence, O2(a1Δ, υ) with υ ≥ 3 could rapidly react with ozone.
Collapse
|
8
|
Zhao B, Han S, Malbon CL, Manthe U, Yarkony DR, Guo H. Full-dimensional quantum stereodynamics of the non-adiabatic quenching of OH(A 2Σ +) by H 2. Nat Chem 2021; 13:909-915. [PMID: 34373597 PMCID: PMC8440216 DOI: 10.1038/s41557-021-00730-1] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/16/2020] [Accepted: 05/11/2021] [Indexed: 11/16/2022]
Abstract
The Born–Oppenheimer approximation, assuming separable nuclear and electronic motion, is widely adopted for characterizing chemical reactions in a single electronic state. However, the breakdown of the Born–Oppenheimer approximation is omnipresent in chemistry, and a detailed understanding of the non-adiabatic dynamics is still incomplete. Here we investigate the non-adiabatic quenching of electronically excited OH(A2Σ+) molecules by H2 molecules using full-dimensional quantum dynamics calculations for zero total nuclear angular momentum using a high-quality diabatic-potential-energy matrix. Good agreement with experimental observations is found for the OH(X2Π) ro-vibrational distribution, and the non-adiabatic dynamics are shown to be controlled by stereodynamics, namely the relative orientation of the two reactants. The uncovering of a major (in)elastic channel, neglected in a previous analysis but confirmed by a recent experiment, resolves a long-standing experiment–theory disagreement concerning the branching ratio of the two electronic quenching channels. ![]()
The breakdown of the Born–Oppenheimer approximation is omnipresent in chemistry and detailed understanding of non-adiabatic dynamics is still incomplete. Now, the non-adiabatic quenching of electronically excited OH(A2Σ+) molecules by H2 has been investigated using full-dimensional quantum dynamics calculations and a high-quality diabatic-potential-energy matrix, providing insight into the branching ratio of the two electronic quenching channels.
Collapse
Affiliation(s)
- Bin Zhao
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, NM, USA. .,Theoretische Chemie, Fakultät für Chemie, Universität Bielefeld, Bielefeld, Germany.
| | - Shanyu Han
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, NM, USA
| | | | - Uwe Manthe
- Theoretische Chemie, Fakultät für Chemie, Universität Bielefeld, Bielefeld, Germany.
| | - David R Yarkony
- Department of Chemistry, Johns Hopkins University, Baltimore, MD, USA.
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, NM, USA.
| |
Collapse
|
9
|
Guan Y, Xie C, Guo H, Yarkony DR. Enabling a Unified Description of Both Internal Conversion and Intersystem Crossing in Formaldehyde: A Global Coupled Quasi-Diabatic Hamiltonian for Its S 0, S 1, and T 1 States. J Chem Theory Comput 2021; 17:4157-4168. [PMID: 34132545 DOI: 10.1021/acs.jctc.1c00370] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
In our recent work, a diabatic Hamiltonian that couples the S0 and S1 states of formaldehyde was constructed using a robust fitting-and-diabatizing procedure with artificial neural networks, which is capable of representing adiabatic energies, energy gradients, and derivative couplings over a wide range of geometries including seams of conical intersection. In this work, based on the diabatization of S0 and S1, the spin-orbit couplings between singlet states (S0, S1) and triplet state T1 are also determined in the same diabatic representation. The diabatized spin-orbit couplings are then fit with a symmetrized neural-network functional form. The ab initio spin-orbit couplings are well reproduced in large configuration space. Together with the neural-network-based potential energy surface for T1, the full quasi-diabatic Hamiltonian for the S0, S1, and T1 states is completed, enabling a unified description of both internal conversion and intersystem crossing in formaldehyde. The vibrational levels on the three adiabatic states are found to be in good agreement with known experimental band origins.
Collapse
Affiliation(s)
- Yafu Guan
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Changjian Xie
- Institute of Modern Physics, Northwest University, Xi'an, Shaanxi 710069, People's Republic of China
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| | - David R Yarkony
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States
| |
Collapse
|
10
|
Guan Y, Xie C, Guo H, Yarkony DR. Neural Network Based Quasi-diabatic Representation for S0 and S1 States of Formaldehyde. J Phys Chem A 2020; 124:10132-10142. [DOI: 10.1021/acs.jpca.0c08948] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Yafu Guan
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Changjian Xie
- Institute of Modern Physics, Northwest University, Xi’an, Shaanxi 710069, People’s Republic of China
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| | - David R. Yarkony
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States
| |
Collapse
|