1
|
Wang J, Xie C, Hu X, Guo H, Xie D. Impact of Geometric Phase on Dynamics of Complex-Forming Reactions: H + O 2 → OH + O. J Phys Chem Lett 2024; 15:4237-4243. [PMID: 38602563 DOI: 10.1021/acs.jpclett.4c00789] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/12/2024]
Abstract
Reaction dynamics on the ground electronic state might be significantly influenced by conical intersections (CIs) via the geometric phase (GP), as demonstrated for activated reactions (i.e., the H + H2 exchange reaction). However, there have been few investigations of GP effects in complex-forming reactions. Here, we report a full quantum dynamical study of an important reaction in combustion (H + O2 → OH + O), which serves as a proving ground for studying GP effects therein. The results reveal significant differences in reaction probabilities and differential cross sections (DCSs) obtained with and without GP, underscoring its strong impact. However, the GP effects are less pronounced for the reaction integral cross sections, apparently due to the integral of the DCS over the scattering angle. Further analysis indicated that the cross section has roughly the same contributions from the two topologically distinct paths around the CI, namely, the direct and looping paths.
Collapse
Affiliation(s)
- Junyan Wang
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, Jiangsu 210023, China
| | - Changjian Xie
- Institute of Modern Physics, Shaanxi Key Laboratory for Theoretical Physics Frontiers, Northwest University, Xi'an, Shaanxi 710127, China
| | - Xixi Hu
- Kuang Yaming Honors School, Institute for Brain Sciences, Jiangsu Key Laboratory of Vehicle Emissions Control, Center of Modern Analysis, Nanjing University, Nanjing, Jiangsu 210023, China
- Hefei National Laboratory, Hefei, Anhui 230088, China
| | - Hua Guo
- Department of Chemistry and Chemical Biology, Center for Computational Chemistry, University of New Mexico, Albuquerque, New Mexico 87131, United States
| | - Daiqian Xie
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, Jiangsu 210023, China
- Hefei National Laboratory, Hefei, Anhui 230088, China
| |
Collapse
|
2
|
Ibele LM, Agostini F. Exploring Exact-Factorization-Based Trajectories for Low-Energy Dynamics near a Conical Intersection. J Phys Chem A 2024. [PMID: 38660710 DOI: 10.1021/acs.jpca.4c00555] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/26/2024]
Abstract
We study low-energy dynamics generated by a two-dimensional two-state Jahn-Teller Hamiltonian in the vicinity of a conical intersection using quantum wave packet and trajectory dynamics. Recently, these dynamics were studied by comparing the adiabatic representation and the exact factorization, with the purpose to highlight the different nature of topological-phase and geometric-phase effects arising in the two theoretical representations of the same problem. Here, we employ the exact factorization to understand how to accurately model low-energy dynamics in the vicinity of a conical intersection using an approximate description of the nuclear motion that uses trajectories. We find that since nonadiabatic effects are weak but non-negligible, the trajectory-based description that invokes the classical approximation struggles to capture the correct behavior.
Collapse
Affiliation(s)
- Lea M Ibele
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, 91405 Orsay, France
| | - Federica Agostini
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, 91405 Orsay, France
| |
Collapse
|
3
|
Guan Y, Chen Q, Varandas AJC. Accurate diabatization based on combined-hyperbolic-inverse-power-representation: 1,2 2A' states of BeH2. J Chem Phys 2024; 160:154105. [PMID: 38624109 DOI: 10.1063/5.0200732] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/28/2024] [Accepted: 03/28/2024] [Indexed: 04/17/2024] Open
Abstract
A diabatic potential energy matrix (DPEM) for the two lowest states of BeH2+ has been constructed using the combined-hyperbolic-inverse-power-representation (CHIPR) method. By imposing symmetry constraints on the coefficients of polynomials, the complete nuclear permutation inversion symmetry is correctly preserved in the CHIPR functional form. The symmetrized CHIPR functional form is then used in the diabatization by ansatz procedure. The ab initio energies are reproduced with satisfactory accuracy. In addition, the CHIPR-based DPEM also reproduces the local topology of a conical intersection. Future work will focus on a complete four-state diabatic representation with emphasis on the long-range interactions and spin-orbit couplings, which will enable accurate quantum scattering calculations for the Be+(2P) + H2 → BeH+(X1Σ+) + H(2S) reaction.
Collapse
Affiliation(s)
- Yafu Guan
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, People's Republic of China
| | - Qun Chen
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, People's Republic of China
| | - António J C Varandas
- School of Physics and Physical Engineering, Qufu Normal University, 273165 Qufu, People's Republic of China
- Department of Physics, Universidade Federal do Espírito Santo, 29075-910 Vitória, Brazil
- Department of Chemistry and Chemistry Centre, University of Coimbra, 3004-535 Coimbra, Portugal
| |
Collapse
|
4
|
Wang Y, Mazziotti DA. Quantum simulation of conical intersections. Phys Chem Chem Phys 2024; 26:11491-11497. [PMID: 38587679 DOI: 10.1039/d4cp00391h] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/09/2024]
Abstract
We explore the simulation of conical intersections (CIs) on quantum devices, setting the groundwork for potential applications in nonadiabatic quantum dynamics within molecular systems. The intersecting potential energy surfaces of H3+ are computed from a variance-based contracted quantum eigensolver. We show how the CIs can be correctly described on quantum devices using wavefunctions generated by the anti-Hermitian contracted Schrödinger equation ansatz, which is a unitary transformation of wavefunctions that preserves the topography of CIs. A hybrid quantum-classical procedure is used to locate the seam of CIs. Additionally, we discuss the quantum implementation of the adiabatic to diabatic transformation and its relation to the geometric phase effect. Results on noisy intermediate-scale quantum devices showcase the potential of quantum computers in dealing with problems in nonadiabatic chemistry.
Collapse
Affiliation(s)
- Yuchen Wang
- Department of Chemistry and The James Franck Institute, The University of Chicago, Chicago, Illinois 60637, USA.
| | - David A Mazziotti
- Department of Chemistry and The James Franck Institute, The University of Chicago, Chicago, Illinois 60637, USA.
| |
Collapse
|
5
|
Gu B. Nonadiabatic Conical Intersection Dynamics in the Local Diabatic Representation with Strang Splitting and Fourier Basis. J Chem Theory Comput 2024; 20:2711-2718. [PMID: 38536965 DOI: 10.1021/acs.jctc.3c01317] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/10/2024]
Abstract
We develop and implement an exact conical intersection nonadiabatic wave packet dynamics method that combines the local diabatic representation, Strang splitting for the total molecular propagator, and discrete variable representation with uniform grids. By employing the local diabatic representation, this method captures all nonadiabatic effects, including nonadiabatic transitions, electronic coherences, and geometric phase. Moreover, it is free of singularities in the first and second derivative couplings and does not require the electronic wave function to be continuous with respect to the nuclear coordinates. We further show that in contrast to the adiabatic representation, the split-operator method can be directly applied to the full molecular propagator with the locally diabatic ansatz. The Fourier series, employed as the primitive nuclear basis functions, is universal and can be applied to all types of reactive coordinates. The combination of local diabatic representation, Strang splitting, and Fourier basis allows numerically exact modeling of conical intersection quantum dynamics directly with adiabatic electronic states that can be obtained from standard electronic structure computations.
Collapse
Affiliation(s)
- Bing Gu
- Department of Chemistry and Department of Physics, Westlake University, Hangzhou, Zhejiang 310030, China
- Institute of Natural Sciences, Westlake Institute for Advanced Study, Hangzhou, Zhejiang 310024, China
| |
Collapse
|
6
|
Morreale D, Persico M. Topology of Conical Intersection Seams and the Geometric Phase. J Phys Chem A 2024; 128:1707-1714. [PMID: 38408203 DOI: 10.1021/acs.jpca.3c07042] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/28/2024]
Abstract
In this paper, we demonstrate two topological properties of crossing seams, that is, the sets of points in the N-dimensional space of nuclear coordinates where two electronic eigenstates are degenerate. We shall examine the typical case of states of the same spin with accidental degeneracies, whereby the crossing seam is of dimension N - 2. The first property we demonstrate is that a crossing seam has no boundary, therefore, it must either extend asymptotically to infinite values of one or more coordinates or wrap on itself. The second property is that two (or more) crossing seams can intersect each other but in such a way that neither of them ends at the intersection. When N = 3, the crossing seam is a line in a 3D space; this is so in triatomic molecules but also in reduced dimensionality treatments of larger polyatomics. The above-mentioned rules then mean that the crossing seam is a line of infinite length or a closed loop and can split into three branches but not in two. The example of the first two excited 1A' states of H2Cl+ illustrates these rules and shows their usefulness for computational search and characterization of crossing seams.
Collapse
Affiliation(s)
- Daniel Morreale
- Department of Chemistry and Industrial Chemistry, University of Pisa, v. G. Moruzzi 13, Pisa 56126, Italy
| | - Maurizio Persico
- Department of Chemistry and Industrial Chemistry, University of Pisa, v. G. Moruzzi 13, Pisa 56126, Italy
| |
Collapse
|
7
|
Han S, Xie C, Hu X, Yarkony DR, Guo H, Xie D. Quantum Dynamics of Photodissociation: Recent Advances and Challenges. J Phys Chem Lett 2023; 14:10517-10530. [PMID: 37970789 DOI: 10.1021/acs.jpclett.3c02735] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2023]
Abstract
Recent advances in constructing accurate potential energy surfaces and nonadiabatic couplings from high-level ab initio data have revealed detailed potential landscapes in not only the ground electronic state but also excited ones. They enabled quantitatively accurate characterization of photoexcited reactive systems using quantum mechanical methods. In this Perspective, we survey the recent progress in quantum mechanical studies of adiabatic and nonadiabatic photodissociation dynamics, focusing on initial state control and product energy disposal. These new insights helped to understand quantum effects in small prototypical systems, and the results serve as benchmarks for developing more approximate theoretical methods.
Collapse
Affiliation(s)
- Shanyu Han
- International Center for Isotope Effects Research, School of Earth Sciences and Engineering, Nanjing University, Nanjing 210023, China
| | - Changjian Xie
- Institute of Modern Physics, Shaanxi Key Laboratory for Theoretical Physics Frontiers, Northwest University, Xi'an 710127, China
| | - Xixi Hu
- Kuang Yaming Honors School, Institute for Brain Sciences, Nanjing University, Nanjing 210023, China
| | - David R Yarkony
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Hua Guo
- Department of Chemistry and Chemical Biology, Center for Computational Chemistry, University of New Mexico, Albuquerque, New Mexico 87131, United States
| | - Daiqian Xie
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| |
Collapse
|
8
|
Fang W, Heller ER, Richardson JO. Competing quantum effects in heavy-atom tunnelling through conical intersections. Chem Sci 2023; 14:10777-10785. [PMID: 37829019 PMCID: PMC10566476 DOI: 10.1039/d3sc03706a] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/18/2023] [Accepted: 09/07/2023] [Indexed: 10/14/2023] Open
Abstract
Thermally activated chemical reactions are typically understood in terms of overcoming potential-energy barriers. However, standard rate theories break down in the presence of a conical intersection (CI) because these processes are inherently nonadiabatic, invalidating the Born-Oppenheimer approximation. Moreover, CIs give rise to intricate nuclear quantum effects such as tunnelling and the geometric phase, which are neglected by standard trajectory-based simulations and remain largely unexplored in complex molecular systems. We present new semiclassical transition-state theories based on an extension of golden-rule instanton theory to describe nonadiabatic tunnelling through CIs and thus provide an intuitive picture for the reaction mechanism. We apply the method in conjunction with first-principles electronic-structure calculations to the electron transfer in the bis(methylene)-adamantyl cation. Our study reveals a strong competition between heavy-atom tunnelling and geometric-phase effects.
Collapse
Affiliation(s)
- Wei Fang
- Department of Chemistry, Fudan University Shanghai 200438 P. R. China
- Department of Chemistry and Applied Biosciences, ETH Zürich 8093 Zürich Switzerland
| | - Eric R Heller
- Department of Chemistry and Applied Biosciences, ETH Zürich 8093 Zürich Switzerland
| | - Jeremy O Richardson
- Department of Chemistry and Applied Biosciences, ETH Zürich 8093 Zürich Switzerland
| |
Collapse
|
9
|
Davis JP, Neisser RW, Kidwell NM. Infrared Activated Signatures and Jahn-Teller Dynamics of NO-CH 4 Collision Complexes. J Phys Chem A 2023. [PMID: 37285367 DOI: 10.1021/acs.jpca.3c01410] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
Bimolecular collision outcomes sensitively depend on the chemical functionality and relative orientations of the colliding partners that define the accessible reactive and nonreactive pathways. Accurate predictions from multidimensional potential energy surfaces demand a full characterization of the available mechanisms. Therefore, there is a need for experimental benchmarks to control and characterize the collision conditions with spectroscopic accuracy to accelerate the predictive modeling of chemical reactivity. To this end, the bimolecular collision outcomes can be investigated systematically by preparing reactants in the entrance channel prior to reaction. Herein, we investigate the vibrational spectroscopy and infrared-driven dynamics of the bimolecular collision complex between nitric oxide and methane (NO-CH4). We recorded the vibrational spectroscopy of NO-CH4 in the CH4 asymmetric stretching region using resonant ion-depletion infrared spectroscopy and infrared action spectroscopy, thus revealing a significantly broad spectrum centered at 3030 cm-1 that extends over 50 cm-1. The asymmetric CH stretch feature of NO-CH4 is explained by CH4 internal rotation and attributed to transitions involving three different nuclear spin isomers of CH4. The vibrational spectra also show extensive homogeneous broadening due to the ultrafast vibrational predissociation of NO-CH4. Additionally, we combine infrared activation of NO-CH4 with velocity map imaging of NO (X2Π, ν″ = 0, J″, Fn, Λ) products to develop a molecular-level understanding of the nonreactive collisions of NO with CH4. The anisotropy of the ion image features is largely determined by the probed rotational quantum number of NO (J″) products. For a subset of NO fragments, the ion images and total kinetic energy release (TKER) distributions show an anisotropic component at low relative translation (∼225 cm-1) indicating a prompt dissociation mechanism. However, for other detected NO products, the ion images and TKER distributions are bimodal, in which the anisotropic component is accompanied by an isotropic feature at high relative translation (∼1400 cm-1) signifying a slow dissociation pathway. In addition to the predissociation dynamics following vibrational excitation, the Jahn-Teller dynamics prior to infrared activation need to be considered to fully describe the product spin-orbit distributions. Therefore, we correlate the Jahn-Teller mechanisms of NO-CH4 to the symmetry-restricted NO (X2Π, ν″ = 0, J″, Fn, Λ) + CH4 (ν″) product outcomes.
Collapse
Affiliation(s)
- John P Davis
- Department of Chemistry, The College of William & Mary, Williamsburg, Virginia 23187-8795, United States
| | - Ruby W Neisser
- Department of Chemistry, The College of William & Mary, Williamsburg, Virginia 23187-8795, United States
| | - Nathanael M Kidwell
- Department of Chemistry, The College of William & Mary, Williamsburg, Virginia 23187-8795, United States
| |
Collapse
|
10
|
Li C, Hou S, Xie C. Constructing Diabatic Potential Energy Matrices with Neural Networks Based on Adiabatic Energies and Physical Considerations: Toward Quantum Dynamic Accuracy. J Chem Theory Comput 2023. [PMID: 37216273 DOI: 10.1021/acs.jctc.2c01074] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/24/2023]
Abstract
A permutation invariant polynomial-neural network (PIP-NN) approach for constructing the global diabatic potential energy matrices (PEMs) of the coupled states of molecules is proposed. Specifically, the diabatization scheme is based merely on the adiabatic energy data of the system, which is ideally a most convenient way due to not requiring additional ab initio calculations for the data of the derivative coupling or any other physical properties of the molecule. Considering the permutation and coupling characteristics of the system, particularly in the presence of conical intersections, some vital treatments for the off-diagonal terms in diabatic PEM are essentially needed. Taking the photodissociation of H2O(X~/B~)/NH3(X~/A~) and nonadiabatic reaction Na(3p) + H2 → NaH(Σ+) + H for example, this PIP-NN method is shown to build up the global diabatic PEMs effectively and accurately. The root-mean-square errors of the adiabatic potential energies in the fitting for three different systems are all small (<10 meV). Further quantum dynamic calculations show that the absorption spectra and product branching ratios in both H2O(X~/B~) and NH3(X~/A~) nonadiabatic photodissociation are well reproduced on the new diabatic PEMs, and the nonadiabatic reaction probability of Na(3p) + H2 → NaH(Σ+) + H obtained on the new diabatic PEMs of the 12A1 and 12B2 states is in reasonably good agreement with previous theoretical result as well, validating this new PIP-NN method.
Collapse
Affiliation(s)
- Chaofan Li
- Institute of Modern Physics, Shaanxi Key Laboratory for Theoretical Physics Frontiers, Northwest University, Xi'an 710127, China
| | - Siting Hou
- Institute of Modern Physics, Shaanxi Key Laboratory for Theoretical Physics Frontiers, Northwest University, Xi'an 710127, China
| | - Changjian Xie
- Institute of Modern Physics, Shaanxi Key Laboratory for Theoretical Physics Frontiers, Northwest University, Xi'an 710127, China
| |
Collapse
|
11
|
Wang J, An F, Chen J, Hu X, Guo H, Xie D. Accurate Full-Dimensional Global Diabatic Potential Energy Matrix for the Two Lowest-Lying Electronic States of the H + O 2 ↔ HO + O Reaction. J Chem Theory Comput 2023; 19:2929-2938. [PMID: 37161259 DOI: 10.1021/acs.jctc.3c00291] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/11/2023]
Abstract
A new and more accurate diabatic potential energy matrix (DPEM) is developed for the two lowest-lying electronic states of HO2, covering both the strong interaction region and reaction asymptotes. The ab initio calculations were performed at the Davidson corrected multireference configuration interaction level with the augmented correlation-consistent polarized valence quintuple-zeta basis set (MRCI+Q/AV5Z). The accuracy of the electronic structure calculations is validated by excellent agreement with the experimental HO2 equilibrium geometry, fundamental vibrational frequencies, and H + O2 ↔ OH + O reaction energy. Through the combination of an electronic angular momentum-method and a configuration interaction vector-based method, the mixing angle between the first two 2A″ states of HO2 was successfully determined. Elements of the 2×2 DPEM were fit to neural networks with a proper account of the complete nuclear permutation inversion symmetry of HO2. The DPEM correctly predicted the properties of conical intersection seams at linear and T-shape geometries, thus providing a reliable platform for studying both the spectroscopy of HO2 and the nonadiabatic dynamics for the H + O2 ↔ OH + O reaction.
Collapse
Affiliation(s)
- Junyan Wang
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Feng An
- Research Center for Graph Computing, Zhejiang Lab, Hangzhou 311121, China
| | - Junjie Chen
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Xixi Hu
- Kuang Yaming Honors School, Jiangsu Key Laboratory of Vehicle Emissions Control, Center of Modern Analysis, Nanjing University, Nanjing 210023, China
- Hefei National Laboratory, Hefei 230088, China
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| | - 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
|
12
|
Shu Y, Zhang L, Wu D, Chen X, Sun S, Truhlar DG. New Gradient Correction Scheme for Electronically Nonadiabatic Dynamics Involving Multiple Spin States. J Chem Theory Comput 2023; 19:2419-2429. [PMID: 37079755 DOI: 10.1021/acs.jctc.2c01173] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/22/2023]
Abstract
It has been recommended that the best representation to use for trajectory surface hopping (TSH) calculations is the fully adiabatic basis in which the Hamiltonian is diagonal. Simulations of intersystem crossing processes with conventional TSH methods require an explicit computation of nonadiabatic coupling vectors (NACs) in the molecular-Coulomb-Hamiltonian (MCH) basis, also called the spin-orbit-free basis, in order to compute the gradient in the fully adiabatic basis (also called the diagonal representation). This explicit requirement destroys some of the advantages of the overlap-based algorithms and curvature-driven algorithms that can be used for the most efficient TSH calculations. Therefore, although these algorithms allow one to perform NAC-free simulations for internal conversion processes, one still requires NACs for intersystem crossing. Here, we show that how the NAC requirement is circumvented by a new computation scheme called the time-derivative-matrix scheme.
Collapse
Affiliation(s)
- Yinan Shu
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Linyao Zhang
- School of Energy Science and Engineering, Harbin Institute of Technology, Harbin 150001, P. R. China
| | - Dihua Wu
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Xiye Chen
- School of Energy Science and Engineering, Harbin Institute of Technology, Harbin 150001, P. R. China
| | - Shaozeng Sun
- School of Energy Science and Engineering, Harbin Institute of Technology, Harbin 150001, P. R. China
| | - Donald G Truhlar
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|
13
|
Han S, Schröder M, Gatti F, Meyer HD, Lauvergnat D, Yarkony DR, Guo H. Representation of Diabatic Potential Energy Matrices for Multiconfiguration Time-Dependent Hartree Treatments of High-Dimensional Nonadiabatic Photodissociation Dynamics. J Chem Theory Comput 2022; 18:4627-4638. [PMID: 35839299 DOI: 10.1021/acs.jctc.2c00370] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Conventional quantum mechanical characterization of photodissociation dynamics is restricted by steep scaling laws with respect to the dimensionality of the system. In this work, we examine the applicability of the multi-configurational time-dependent Hartree (MCTDH) method in treating nonadiabatic photodissociation dynamics in two prototypical systems, taking advantage of its favorable scaling laws. To conform to the sum-of-product form, elements of the ab initio diabatic potential energy matrix (DPEM) are re-expressed using the recently proposed Monte Carlo canonical polyadic decomposition method, with enforcement of proper symmetry. The MCTDH absorption spectra and product branching ratios are shown to compare well with those calculated using conventional grid-based methods, demonstrating its promise for treating high-dimensional nonadiabatic photodissociation problems.
Collapse
Affiliation(s)
- Shanyu Han
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| | - Markus Schröder
- Theoretische Chemie, Physikalisch Chemisches Institut, Ruprecht-Karls Universität Heidelberg, D-69120 Heidelberg, Germany
| | - Fabien Gatti
- ISMO, Institut des Sciences Moléculaires d'Orsay─UMR 8214 CNRS/Université Paris-Saclay, F-91405 Orsay, France
| | - Hans-Dieter Meyer
- Theoretische Chemie, Physikalisch Chemisches Institut, Ruprecht-Karls Universität Heidelberg, D-69120 Heidelberg, Germany
| | - David Lauvergnat
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, Orsay 91405, France
| | - David R Yarkony
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| |
Collapse
|
14
|
Maskri R, Joubert-Doriol L. The moving crude adiabatic alternative to the adiabatic representation in excited state dynamics. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2022; 380:20200379. [PMID: 35341311 DOI: 10.1098/rsta.2020.0379] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/06/2021] [Accepted: 12/09/2021] [Indexed: 06/14/2023]
Abstract
The choice of the electronic representation in on-the-fly quantum dynamics is crucial. The adiabatic representation is appealing since adiabatic states are readily available from quantum chemistry packages. The nuclear wavepackets are then expanded in a basis of Gaussian functions, which follow trajectories to explore the potential energy surfaces and approximate the potential using a local expansion of the adiabatic quantities. Nevertheless, the adiabatic representation is plagued with severe limitations when conical intersections are involved: the diagonal Born-Oppenheimer corrections (DBOCs) are non-integrable, and the geometric phase effect on the nuclear wavepackets cannot be accounted for unless a model is available. To circumvent these difficulties, the moving crude adiabatic (MCA) representation was proposed and successfully tested in low energy dynamics where the wavepacket skirts the conical intersection. We assess the MCA representation in the case of non-adiabatic transitions through conical intersections. First, we show that using a Gaussian basis in the adiabatic representation indeed exhibits the aforementioned difficulties with a special emphasis on the possibility to regularize the DBOC terms. Then, we show that MCA is indeed able to properly model non-adiabatic transitions. Tests are done on linear vibronic coupling models for the bis(methylene) adamantyl cation and the butatriene cation. This article is part of the theme issue 'Chemistry without the Born-Oppenheimer approximation'.
Collapse
Affiliation(s)
- Rosa Maskri
- Univ Gustave Eiffel, Univ Paris Est Creteil, CNRS, UMR 8208, MSME, F-77454 Marne-la-Vallée, France
| | - Loïc Joubert-Doriol
- Univ Gustave Eiffel, Univ Paris Est Creteil, CNRS, UMR 8208, MSME, F-77454 Marne-la-Vallée, France
| |
Collapse
|
15
|
Quantum-classical simulations of rhodopsin reveal excited-state population splitting and its effects on quantum efficiency. Nat Chem 2022; 14:441-449. [PMID: 35241801 PMCID: PMC8983576 DOI: 10.1038/s41557-022-00892-6] [Citation(s) in RCA: 13] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/28/2020] [Accepted: 01/12/2022] [Indexed: 01/12/2023]
Abstract
The activation of rhodopsin, the light-sensitive G-protein coupled receptor responsible for dim-light vision in vertebrates, is driven by an ultrafast excited-state double-bond isomerization with a quantum efficiency of almost 70%. The origin of such light sensitivity is not understood and a key question is whether in-phase nuclear motion controls the quantum efficiency value. Here, we use hundreds of quantum-classical trajectories to show that, 15 femtoseconds after light absorption, a degeneracy between the reactive excited state and a neighboring state causes the splitting of the rhodopsin population into subpopulations. These subpopulations propagate with different velocities and lead to distinct contributions to the quantum efficiency. We also show that such splitting is modulated by protein electrostatics, thus linking amino-acid sequence variations to quantum efficiency modulation. Finally, we discuss how such a linkage that in principle could be exploited to achieve higher quantum efficiencies, would simultaneously increase the receptor thermal noise leading to a trade-off that may have played a role in rhodopsin evolution.
Collapse
|
16
|
Kim J, Woo KC, Kim KK, Kang M, Kim SK. Tunneling dynamics dictated by the multidimensional conical intersection seam in the πσ*‐mediated photochemistry of heteroaromatic molecules. B KOREAN CHEM SOC 2021. [DOI: 10.1002/bkcs.12453] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Affiliation(s)
- Junggil Kim
- Department of Chemistry, KAIST Daejeon Republic of Korea
| | - Kyung Chul Woo
- Department of Chemistry, KAIST Daejeon Republic of Korea
- Division of Chemistry and Biological Chemistry, School of Physical and Mathematical Sciences Nanyang Technological University Singapore Singapore
| | - Kuk Ki Kim
- Department of Chemistry, KAIST Daejeon Republic of Korea
| | - Minseok Kang
- Department of Chemistry, KAIST Daejeon Republic of Korea
| | - Sang Kyu Kim
- Department of Chemistry, KAIST Daejeon Republic of Korea
| |
Collapse
|
17
|
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: 26] [Impact Index Per Article: 8.7] [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
|
18
|
Westermayr J, Marquetand P. Machine Learning for Electronically Excited States of Molecules. Chem Rev 2021; 121:9873-9926. [PMID: 33211478 PMCID: PMC8391943 DOI: 10.1021/acs.chemrev.0c00749] [Citation(s) in RCA: 162] [Impact Index Per Article: 54.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/17/2020] [Indexed: 12/11/2022]
Abstract
Electronically excited states of molecules are at the heart of photochemistry, photophysics, as well as photobiology and also play a role in material science. Their theoretical description requires highly accurate quantum chemical calculations, which are computationally expensive. In this review, we focus on not only how machine learning is employed to speed up such excited-state simulations but also how this branch of artificial intelligence can be used to advance this exciting research field in all its aspects. Discussed applications of machine learning for excited states include excited-state dynamics simulations, static calculations of absorption spectra, as well as many others. In order to put these studies into context, we discuss the promises and pitfalls of the involved machine learning techniques. Since the latter are mostly based on quantum chemistry calculations, we also provide a short introduction into excited-state electronic structure methods and approaches for nonadiabatic dynamics simulations and describe tricks and problems when using them in machine learning for excited states of molecules.
Collapse
Affiliation(s)
- Julia Westermayr
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
| | - Philipp Marquetand
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Vienna
Research Platform on Accelerating Photoreaction Discovery, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Data
Science @ Uni Vienna, University of Vienna, Währinger Strasse 29, 1090 Vienna, Austria
| |
Collapse
|
19
|
Abstract
Electronically excited states of molecules are at the heart of photochemistry, photophysics, as well as photobiology and also play a role in material science. Their theoretical description requires highly accurate quantum chemical calculations, which are computationally expensive. In this review, we focus on not only how machine learning is employed to speed up such excited-state simulations but also how this branch of artificial intelligence can be used to advance this exciting research field in all its aspects. Discussed applications of machine learning for excited states include excited-state dynamics simulations, static calculations of absorption spectra, as well as many others. In order to put these studies into context, we discuss the promises and pitfalls of the involved machine learning techniques. Since the latter are mostly based on quantum chemistry calculations, we also provide a short introduction into excited-state electronic structure methods and approaches for nonadiabatic dynamics simulations and describe tricks and problems when using them in machine learning for excited states of molecules.
Collapse
Affiliation(s)
- Julia Westermayr
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
| | - Philipp Marquetand
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Vienna Research Platform on Accelerating Photoreaction Discovery, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Data Science @ Uni Vienna, University of Vienna, Währinger Strasse 29, 1090 Vienna, Austria
| |
Collapse
|
20
|
Farag MH, Mandal A, Huo P. Polariton induced conical intersection and berry phase. Phys Chem Chem Phys 2021; 23:16868-16879. [PMID: 34328152 DOI: 10.1039/d1cp00943e] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
Abstract
We investigate the Polariton induced conical intersection (PICI) created from coupling a diatomic molecule with the quantized photon mode inside an optical cavity, and the corresponding Berry Phase effects. We use the rigorous Pauli-Fierz Hamiltonian to describe the quantum light-matter interactions between a LiF molecule and the cavity, and use the exact quantum propagation to investigate the polariton quantum dynamics. The molecular rotations relative to the cavity polarization direction play a role as the tuning mode of the PICI, resulting in an effective CI even within a diatomic molecule. To clearly demonstrate the dynamical effects of the Berry phase, we construct two additional models that have the same Born-Oppenheimer surface, but the effects of the geometric phase are removed. We find that when the initial wavefunction is placed in the lower polaritonic surface, the Berry phase causes a π phase-shift in the wavefunction after the encirclement around the CI, indicated from the nuclear probability distribution. On the other hand, when the initial wavefunction is placed in the upper polaritonic surface, the geometric phase significantly influences the couplings between polaritonic states and therefore, the population dynamics between them. These BP effects are further demonstrated through the photo-fragment angular distribution. PICI created from the quantized radiation field has the promise to open up new possibilities to modulate photochemical reactivities.
Collapse
Affiliation(s)
- Marwa H Farag
- Department of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, USA.
| | | | | |
Collapse
|
21
|
Choi S, Vaníček J. How important are the residual nonadiabatic couplings for an accurate simulation of nonadiabatic quantum dynamics in a quasidiabatic representation? J Chem Phys 2021; 154:124119. [PMID: 33810696 DOI: 10.1063/5.0046067] [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/15/2022] Open
Abstract
Diabatization of the molecular Hamiltonian is a standard approach to remove the singularities of nonadiabatic couplings at conical intersections of adiabatic potential energy surfaces. In general, it is impossible to eliminate the nonadiabatic couplings entirely-the resulting "quasidiabatic" states are still coupled by smaller but nonvanishing residual nonadiabatic couplings, which are typically neglected. Here, we propose a general method for assessing the validity of this potentially drastic approximation by comparing quantum dynamics simulated either with or without the residual couplings. To make the numerical errors negligible to the errors due to neglecting the residual couplings, we use the highly accurate and general eighth-order composition of the implicit midpoint method. The usefulness of the proposed method is demonstrated on nonadiabatic simulations in the cubic Jahn-Teller model of nitrogen trioxide and in the induced Renner-Teller model of hydrogen cyanide. We find that, depending on the system, initial state, and employed quasidiabatization scheme, neglecting the residual couplings can result in wrong dynamics. In contrast, simulations with the exact quasidiabatic Hamiltonian, which contains the residual couplings, always yield accurate results.
Collapse
Affiliation(s)
- Seonghoon Choi
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Jiří Vaníček
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| |
Collapse
|
22
|
Fábri C, Lasorne B, Halász GJ, Cederbaum LS, Vibók Á. Quantum light-induced nonadiabatic phenomena in the absorption spectrum of formaldehyde: Full- and reduced-dimensionality studies. J Chem Phys 2020; 153:234302. [DOI: 10.1063/5.0035870] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Csaba Fábri
- Laboratory of Molecular Structure and Dynamics, Institute of Chemistry, Eötvös Loránd University, Pázmány Péter Sétány 1/A, H-1117 Budapest, Hungary
- MTA-ELTE Complex Chemical Systems Research Group, P.O. Box 32, H-1518 Budapest 112, Hungary
| | - Benjamin Lasorne
- Institut Charles Gerhardt Montpellier (ICGM), Université de Montpellier, CNRS, ENSCM, F-34095 Montpellier, France
| | - Gábor J. Halász
- Department of Information Technology, University of Debrecen, P.O. Box 400, H-4002 Debrecen, Hungary
| | - Lorenz S. Cederbaum
- Theoretische Chemie, Physikalisch-Chemisches Institut, Universität Heidelberg, Im Neuenheimer Feld 229, 69120 Heidelberg, Germany
| | - Ágnes Vibók
- Department of Theoretical Physics, University of Debrecen, P.O. Box 400, H-4002 Debrecen, Hungary
- ELI-ALPS, ELI-HU Non-Profit Ltd., Dugonics tér 13, H-6720 Szeged, Hungary
| |
Collapse
|
23
|
Choi S, Vaníček J. Which form of the molecular Hamiltonian is the most suitable for simulating the nonadiabatic quantum dynamics at a conical intersection? J Chem Phys 2020; 153:211101. [DOI: 10.1063/5.0033410] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Seonghoon Choi
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Jiří Vaníček
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| |
Collapse
|
24
|
Fábri C, Halász GJ, Cederbaum LS, Vibók Á. Born-Oppenheimer approximation in optical cavities: from success to breakdown. Chem Sci 2020; 12:1251-1258. [PMID: 34163887 PMCID: PMC8179040 DOI: 10.1039/d0sc05164k] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022] Open
Abstract
The coupling of a molecule and a cavity induces nonadiabaticity in the molecule which makes the description of its dynamics complicated. For polyatomic molecules, reduced-dimensional models and the use of the Born-Oppenheimer approximation (BOA) may remedy the situation. It is demonstrated that contrary to expectation, BOA may even fail in a one-dimensional model and is generally expected to fail in two- or more-dimensional models due to the appearance of conical intersections induced by the cavity.
Collapse
Affiliation(s)
- Csaba Fábri
- Laboratory of Molecular Structure and Dynamics, Institute of Chemistry, Eötvös Loránd University Pázmány Péter sétány 1/A H-1117 Budapest Hungary .,MTA-ELTE Complex Chemical Systems Research Group P.O. Box 32 H-1518 Budapest Hungary
| | - Gábor J Halász
- Department of Information Technology, University of Debrecen P.O. Box 400 H-4002 Debrecen Hungary
| | - Lorenz S Cederbaum
- Theoretische Chemie, Physikalisch-Chemisches Institut, Universität Heidelberg Im Neuenheimer Feld 229 69120 Heidelberg Germany
| | - Ágnes Vibók
- Department of Theoretical Physics, University of Debrecen PO Box 400 H-4002 Debrecen Hungary .,ELI-ALPS, ELI-HU Non-Profit Ltd Dugonics tér 13 H-6720 Szeged Hungary
| |
Collapse
|
25
|
Yan Y, Sharma K, Miller TA, Liu J. Rotational and fine structure of open-shell molecules in nearly degenerate electronic states. II. Interpretation of experimentally determined interstate coupling parameters of alkoxy radicals. J Chem Phys 2020; 153:174306. [PMID: 33167646 DOI: 10.1063/5.0026212] [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/14/2022] Open
Abstract
Rotationally and fine-structure resolved B̃←X̃ laser-induced fluorescence (LIF) spectra of alkoxy radicals have been simulated with a "coupled two-states model" [J. Liu, J. Chem. Phys. 148, 124112 (2018)], in which the nearly degenerate X̃ and à states are considered together. These two electronic states are separated by the "difference potential" and coupled by the spin-orbit (SO) interaction and the Coriolis interaction. Molecular constants determined in fitting the LIF spectra using the coupled two-states model provide quantitative insight into the SO and Coriolis interactions, as well as other intramolecular dynamics, including the pseudo-Jahn-Teller effect. The spectroscopic model also allows semi-quantitative prediction of effective spin-rotation constants using molecular geometry and SO constants, which can be calculated ab initio with considerable accuracy. The dependence of fit values of molecular constants on the size and conformation of alkoxy radicals is discussed.
Collapse
Affiliation(s)
- Yi Yan
- Department of Chemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Ketan Sharma
- Department of Chemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Terry A Miller
- Department of Chemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Jinjun Liu
- Department of Chemistry, University of Louisville, Louisville, Kentucky 40292, USA
| |
Collapse
|
26
|
Li J, Zhao B, Xie D, Guo H. Advances and New Challenges to Bimolecular Reaction Dynamics Theory. J Phys Chem Lett 2020; 11:8844-8860. [PMID: 32970441 DOI: 10.1021/acs.jpclett.0c02501] [Citation(s) in RCA: 32] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
Dynamics of bimolecular reactions in the gas phase are of foundational importance in combustion, atmospheric chemistry, interstellar chemistry, and plasma chemistry. These collision-induced chemical transformations are a sensitive probe of the underlying potential energy surface(s). Despite tremendous progress in past decades, our understanding is still not complete. In this Perspective, we survey the recent advances in theoretical characterization of bimolecular reaction dynamics, stimulated by new experimental observations, and identify key new challenges.
Collapse
Affiliation(s)
- Jun Li
- School of Chemistry and Chemical Engineering & Chongqing Key Laboratory of Theoretical and Computational Chemistry, Chongqing University, Chongqing 401331, China
| | - Bin Zhao
- Theoretische Chemie, Fakultät für Chemie, Universität Bielefeld, Universitätsstraße 25, D-33615 Bielefeld, Germany
| | - Daiqian Xie
- Institute of Theoretical and Computational Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| |
Collapse
|
27
|
Han S, Wang Y, Guan Y, Yarkony DR, Guo H. Impact of Diabolical Singular Points on Nonadiabatic Dynamics and a Remedy: Photodissociation of Ammonia in the First Band. J Chem Theory Comput 2020; 16:6776-6784. [DOI: 10.1021/acs.jctc.0c00811] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Shanyu Han
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| | - Yucheng Wang
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Yafu Guan
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - David R. Yarkony
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| |
Collapse
|
28
|
Yuan D, Huang Y, Chen W, Zhao H, Yu S, Luo C, Tan Y, Wang S, Wang X, Sun Z, Yang X. Observation of the geometric phase effect in the H+HD→H 2+D reaction below the conical intersection. Nat Commun 2020; 11:3640. [PMID: 32686682 PMCID: PMC7371868 DOI: 10.1038/s41467-020-17381-4] [Citation(s) in RCA: 25] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/09/2020] [Accepted: 06/26/2020] [Indexed: 11/24/2022] Open
Abstract
It has long been known that there is a conical intersection (CI) between the ground and first excited electronic state in the H3 system. Its associated geometric phase (GP) effect has been theoretically predicted to exist below the CI since a long time. However, the experimental evidence has not been established yet and its dynamical origin is waiting to be elucidated. Here we report a combined crossed molecular beam and quantum reactive scattering dynamics study of the H+HD → H2+D reaction at 2.28 eV, which is well below the CI. The GP effect is clearly identified by the observation of distinct oscillations in the differential cross section around the forward direction. Quantum dynamics theory reveals that the GP effect arises from the phase alteration of a small part of the wave function, which corresponds to an unusual roaming-like abstraction pathway, as revealed by quasi-classical trajectory calculations. The geometric phase effect associated with a conical intersection between the ground and first excited electronic state has been predicted in the H3 system below the conical intersection energy. The authors, by a crossed molecular beam technique and quantum dynamic calculations, provide experimental evidence and insight into its origin.
Collapse
Affiliation(s)
- Daofu Yuan
- Hefei National Laboratory for Physical Sciences at the Microscale and Department of Chemical Physics, University of Science and Technology of China, Hefei, 230026, China
| | - Yin Huang
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian, 116023, China
| | - Wentao Chen
- Hefei National Laboratory for Physical Sciences at the Microscale and Department of Chemical Physics, University of Science and Technology of China, Hefei, 230026, China
| | - Hailin Zhao
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian, 116023, China
| | - Shengrui Yu
- Hangzhou Institute of Advanced Studies, Zhejiang Normal University, Hangzhou, 311231, China
| | - Chang Luo
- Hefei National Laboratory for Physical Sciences at the Microscale and Department of Chemical Physics, University of Science and Technology of China, Hefei, 230026, China
| | - Yuxin Tan
- Hefei National Laboratory for Physical Sciences at the Microscale and Department of Chemical Physics, University of Science and Technology of China, Hefei, 230026, China
| | - Siwen Wang
- Hefei National Laboratory for Physical Sciences at the Microscale and Department of Chemical Physics, University of Science and Technology of China, Hefei, 230026, China
| | - Xingan Wang
- Hefei National Laboratory for Physical Sciences at the Microscale and Department of Chemical Physics, University of Science and Technology of China, Hefei, 230026, China.
| | - Zhigang Sun
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian, 116023, China.
| | - Xueming Yang
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian, 116023, China.,Department of Chemistry, College of Science, Southern University of Science and Technology, Shenzhen, 518055, China
| |
Collapse
|
29
|
Fábri C, Lasorne B, Halász GJ, Cederbaum LS, Vibók Á. Striking Generic Impact of Light-Induced Non-Adiabaticity in Polyatomic Molecules. J Phys Chem Lett 2020; 11:5324-5329. [PMID: 32530631 DOI: 10.1021/acs.jpclett.0c01574] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
Non-adiabaticity, i.e., the effect of mixing electronic states by nuclear motion, is a central phenomenon in molecular science. The strongest nonadiabatic effects arise due to the presence of conical intersections of electronic energy surfaces. These intersections are abundant in polyatomic molecules. Laser light can induce in a controlled manner new conical intersections, called light-induced conical intersections, which lead to strong nonadiabatic effects similar to those of the natural conical intersections. These effects are, however, controllable and may even compete with those of the natural intersections. In this work we show that the standard low-energy vibrational spectrum of the electronic ground state can change dramatically by inducing non-adiabaticity via a light-induced conical intersection. This generic effect is demonstrated for an explicit example by full-dimensional high-level quantum calculations using a pump-probe scheme with a moderate-intensity pump laser and a weak probe laser.
Collapse
Affiliation(s)
- Csaba Fábri
- Laboratory of Molecular Structure and Dynamics, Institute of Chemistry, Eötvös Loránd University, Pázmány Péter sétány 1/A, H-1117 Budapest, Hungary
- MTA-ELTE Complex Chemical Systems Research Group, P.O. Box 32, H-1518 Budapest 112, Hungary
| | - Benjamin Lasorne
- Institut Charles Gerhardt Montpellier (ICGM), Université de Montpellier, CNRS, ENSCM, F-34095 Montpellier, France
| | - Gábor J Halász
- Department of Information Technology, University of Debrecen, P.O. Box 400, H-4002 Debrecen, Hungary
| | - Lorenz S Cederbaum
- Theoretische Chemie, Physikalisch-Chemisches Institut, Universität Heidelberg, Im Neuenheimer Feld 229, 69120 Heidelberg, Germany
| | - Ágnes Vibók
- Department of Theoretical Physics, University of Debrecen, P.O. Box 400, H-4002 Debrecen, Hungary
- ELI-ALPS, ELI-HU Non-Profit Ltd., Dugonics tér 13, H-6720 Szeged, Hungary
| |
Collapse
|
30
|
Prlj A, Begušić T, Zhang ZT, Fish GC, Wehrle M, Zimmermann T, Choi S, Roulet J, Moser JE, Vaníček J. Semiclassical Approach to Photophysics Beyond Kasha's Rule and Vibronic Spectroscopy Beyond the Condon Approximation. The Case of Azulene. J Chem Theory Comput 2020; 16:2617-2626. [PMID: 32119547 DOI: 10.1021/acs.jctc.0c00079] [Citation(s) in RCA: 21] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Azulene is a prototypical molecule with an anomalous fluorescence from the second excited electronic state, thus violating Kasha's rule, and with an emission spectrum that cannot be understood within the Condon approximation. To better understand the photophysics and spectroscopy of azulene and other nonconventional molecules, we developed a systematic, general, and efficient computational approach combining the semiclassical dynamics of nuclei with ab initio electronic structure. First, to analyze the nonadiabatic effects, we complement the standard population dynamics by a rigorous measure of adiabaticity, estimated with the multiple-surface dephasing representation. Second, we propose a new semiclassical method for simulating non-Condon spectra, which combines the extended thawed Gaussian approximation with the efficient single-Hessian approach. S1 ← S0 and S2 ← S0 absorption and S2 → S0 emission spectra of azulene, recorded in a new set of experiments, agree very well with our calculations. We find that accuracy of the evaluated spectra requires the treatment of anharmonicity, Herzberg-Teller, and mode-mixing effects.
Collapse
Affiliation(s)
- Antonio Prlj
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Tomislav Begušić
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Zhan Tong Zhang
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - George Cameron Fish
- Photochemical Dynamics Group, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Marius Wehrle
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Tomáš Zimmermann
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Seonghoon Choi
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Julien Roulet
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Jacques-Edouard Moser
- Photochemical Dynamics Group, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Jiří Vaníček
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| |
Collapse
|
31
|
Xie C, Zhao B, Malbon CL, Yarkony DR, Xie D, Guo H. Insights into the Mechanism of Nonadiabatic Photodissociation from Product Vibrational Distributions. The Remarkable Case of Phenol. J Phys Chem Lett 2020; 11:191-198. [PMID: 31821757 DOI: 10.1021/acs.jpclett.9b03407] [Citation(s) in RCA: 20] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
The fate of a photoexcited molecule is often strongly influenced by electronic degeneracies, such as conical intersections, which break the Born-Oppenheimer separation of electronic and nuclear motion. Detailed information concerning internal energy redistribution in a nonadiabatic process can be extracted from the product state distribution of a photofragment in photodissociation. Here, we focus on the nonadiabatic photodissociation of phenol and discuss the internal excitation of the phenoxyl fragment using both symmetry analysis and wave packet dynamics. It is shown that unique and general selection rules exist, which can be attributed to the geometric phase in the adiabatic representation. Further, our results provide a reinterpretation of the experimental data, shedding light on the impact of conical intersections on the product state distribution.
Collapse
Affiliation(s)
- Changjian Xie
- Institute of Modern Physics, Shaanxi Key Laboratory for Theoretical Physics Frontiers , Northwest University , Xian , Shaanxi 710127 , China
- Department of Chemistry and Chemical Biology , University of New Mexico , Albuquerque , New Mexico 87131 , United States
| | - Bin Zhao
- Department of Chemistry and Chemical Biology , University of New Mexico , Albuquerque , New Mexico 87131 , United States
| | - Christopher L Malbon
- Department of Chemistry , Johns Hopkins University , Baltimore , Maryland 21218 , United States
| | - David R Yarkony
- Department of Chemistry , Johns Hopkins University , Baltimore , Maryland 21218 , United States
| | - Daiqian Xie
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering , Nanjing University , Nanjing 210023 , China
| | - Hua Guo
- Department of Chemistry and Chemical Biology , University of New Mexico , Albuquerque , New Mexico 87131 , United States
| |
Collapse
|
32
|
Wang Y, Guan Y, Yarkony DR. On the Impact of Singularities in the Two-State Adiabatic to Diabatic State Transformation: A Global Treatment. J Phys Chem A 2019; 123:9874-9880. [DOI: 10.1021/acs.jpca.9b08519] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Yuchen Wang
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Yafu Guan
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - David R. Yarkony
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States
| |
Collapse
|
33
|
Schäpers D, Manthe U. Vibronic coupling in the F·CH4 prereactive complex. J Chem Phys 2019; 151:104106. [DOI: 10.1063/1.5110246] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Daniela Schäpers
- Theoretische Chemie, Fakultät für Chemie, Universität Bielefeld, Universitätsstr. 25, D-33615 Bielefeld, Germany
| | - Uwe Manthe
- Theoretische Chemie, Fakultät für Chemie, Universität Bielefeld, Universitätsstr. 25, D-33615 Bielefeld, Germany
| |
Collapse
|
34
|
Wang Y, Xie C, Guo H, Yarkony DR. A Quasi-Diabatic Representation of the 1,21A States of Methylamine. J Phys Chem A 2019; 123:5231-5241. [DOI: 10.1021/acs.jpca.9b03801] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
Affiliation(s)
- Yuchen Wang
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Changjian Xie
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| | - 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
|
35
|
Yarkony DR, Xie C, Zhu X, Wang Y, Malbon CL, Guo H. Diabatic and adiabatic representations: Electronic structure caveats. COMPUT THEOR CHEM 2019. [DOI: 10.1016/j.comptc.2019.01.020] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
|
36
|
Xie C, Malbon CL, Xie D, Yarkony DR, Guo H. Nonadiabatic Dynamics in Photodissociation of Hydroxymethyl in the 32A(3px) Rydberg State: A Nine-Dimensional Quantum Study. J Phys Chem A 2019; 123:1937-1944. [DOI: 10.1021/acs.jpca.8b12184] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Changjian Xie
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| | - Christopher L. Malbon
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Daiqian Xie
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - David R. Yarkony
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| |
Collapse
|