1
|
Shu Y, Varga Z, Kanchanakungwankul S, Zhang L, Truhlar DG. Diabatic States of Molecules. J Phys Chem A 2022; 126:992-1018. [PMID: 35138102 DOI: 10.1021/acs.jpca.1c10583] [Citation(s) in RCA: 31] [Impact Index Per Article: 15.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Quantitative simulations of electronically nonadiabatic molecular processes require both accurate dynamics algorithms and accurate electronic structure information. Direct semiclassical nonadiabatic dynamics is expensive due to the high cost of electronic structure calculations, and hence it is limited to small systems, limited ensemble averaging, ultrafast processes, and/or electronic structure methods that are only semiquantitatively accurate. The cost of dynamics calculations can be made manageable if analytic fits are made to the electronic structure data, and such fits are most conveniently carried out in a diabatic representation because the surfaces are smooth and the couplings between states are smooth scalar functions. Diabatic representations, unlike the adiabatic ones produced by most electronic structure methods, are not unique, and finding suitable diabatic representations often involves time-consuming nonsystematic diabatization steps. The biggest drawback of using diabatic bases is that it can require large amounts of effort to perform a globally consistent diabatization, and one of our goals has been to develop methods to do this efficiently and automatically. In this Feature Article, we introduce the mathematical framework of diabatic representations, and we discuss diabatization methods, including adiabatic-to-diabatic transformations and recent progress toward the goal of automatization.
Collapse
Affiliation(s)
- Yinan Shu
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Zoltan Varga
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Siriluk Kanchanakungwankul
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Linyao Zhang
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States.,School of Energy Science and Engineering, Harbin Institute of Technology, Harbin 150001, P. R. China
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|
2
|
Zhang L, Shu Y, Sun S, Truhlar DG. Direct coherent switching with decay of mixing for intersystem crossing dynamics of thioformaldehyde: The effect of decoherence. J Chem Phys 2021; 154:094310. [PMID: 33685154 DOI: 10.1063/5.0037878] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023] Open
Abstract
We evaluate the effect of electronic decoherence on intersystem crossing in the photodynamics of thioformaldehyde. First, we show that the state-averaged complete-active-space self-consistent field electronic structure calculations with a properly chosen active space of 12 active electrons in 10 active orbitals can predict the potential energy surfaces and the singlet-triplet spin-orbit couplings quite well for CH2S, and we use this method for direct dynamics by coherent switching with decay of mixing (CSDM). We obtain similar dynamical results with CSDM or by adding energy-based decoherence to trajectory surface hopping, with the population of triplet states tending to a small steady-state value over 500 fs. Without decoherence, the state populations calculated by the conventional trajectory surface hopping method or the semiclassical Ehrenfest method gradually increase. This difference shows that decoherence changes the nature of the results not just quantitatively but qualitatively.
Collapse
Affiliation(s)
- Linyao Zhang
- School of Energy Science and Engineering, Harbin Institute of Technology, Harbin 150001, People's Republic of China
| | - Yinan Shu
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, USA
| | - Shaozeng Sun
- School of Energy Science and Engineering, Harbin Institute of Technology, Harbin 150001, People's Republic of China
| | - Donald G Truhlar
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, USA
| |
Collapse
|
3
|
Talotta F, Morisset S, Rougeau N, Lauvergnat D, Agostini F. Internal Conversion and Intersystem Crossing with the Exact Factorization. J Chem Theory Comput 2020; 16:4833-4848. [PMID: 32633509 DOI: 10.1021/acs.jctc.0c00493] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We present a detailed derivation of the generalized coupled-trajectory mixed quantum-classical (G-CT-MQC) algorithm based on the exact-factorization equations. The ultimate goal is to propose an algorithm that can be employed for molecular dynamics simulations of nonradiative phenomena, as the spin-allowed internal conversions and the spin-forbidden intersystem crossings. Internal conversions are nonadiabatic processes driven by the kinetic coupling between electronic states, whereas intersystem crossings are mediated by the spin-orbit coupling. In this paper, we discuss computational issues related to the suitable representation for electronic dynamics and the different natures of kinetic and spin-orbit coupling. Numerical studies on model systems allow us to test the performance of the G-CT-MQC algorithm in different situations.
Collapse
Affiliation(s)
- Francesco Talotta
- Institut de Chimie Physique, UMR8000, CNRS, Université Paris-Saclay, 91405 Orsay, France.,Institut des Sciences Moléculaires d'Orsay, CNRS, Université Paris-Saclay, 91405 Orsay, France
| | - Sabine Morisset
- Institut des Sciences Moléculaires d'Orsay, CNRS, Université Paris-Saclay, 91405 Orsay, France
| | - Nathalie Rougeau
- Institut des Sciences Moléculaires d'Orsay, CNRS, Université Paris-Saclay, 91405 Orsay, France
| | - David Lauvergnat
- Institut de Chimie Physique, UMR8000, CNRS, Université Paris-Saclay, 91405 Orsay, France
| | - Federica Agostini
- Institut de Chimie Physique, UMR8000, CNRS, Université Paris-Saclay, 91405 Orsay, France
| |
Collapse
|
4
|
Shu Y, Zhang L, Mai S, Sun S, González L, Truhlar DG. Implementation of Coherent Switching with Decay of Mixing into the SHARC Program. J Chem Theory Comput 2020; 16:3464-3475. [DOI: 10.1021/acs.jctc.0c00112] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Yinan Shu
- Department of Chemistry and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Linyao Zhang
- Department of Chemistry and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
- School of Energy Science and Engineering, Harbin Institute of Technology, Harbin 150001, People’s Republic of China
| | - Sebastian Mai
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Str. 17, 1090, Vienna, Austria
| | - Shaozeng Sun
- School of Energy Science and Engineering, Harbin Institute of Technology, Harbin 150001, People’s Republic of China
| | - Leticia González
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Str. 17, 1090, Vienna, Austria
| | - Donald G. Truhlar
- Department of Chemistry and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|
5
|
Talotta F, Morisset S, Rougeau N, Lauvergnat D, Agostini F. Spin-Orbit Interactions in Ultrafast Molecular Processes. PHYSICAL REVIEW LETTERS 2020; 124:033001. [PMID: 32031839 DOI: 10.1103/physrevlett.124.033001] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/27/2019] [Indexed: 06/10/2023]
Abstract
We investigate spin-orbit interactions in ultrafast molecular processes employing the exact factorization of the electron-nuclear wave function. We revisit the original derivation by including spin-orbit coupling, and show how the dynamics driven by the time-dependent potential energy surface alleviates inconsistencies arising from different electronic representations. We propose a novel trajectory-based scheme to simulate spin-forbidden non-radiative processes, and we show its performance in the treatment of excited-state dynamics where spin-orbit effects couple different spin multiplets.
Collapse
Affiliation(s)
- Francesco Talotta
- Laboratoire de Chimie Physique, UMR 8000 CNRS/University Paris-Sud, University Paris-Saclay, 91405 Orsay, France
- Institut de Sciences Moléculaires d'Orsay, UMR 8214 CNRS/University Paris-Sud, University Paris-Saclay, 91405 Orsay, France
| | - Sabine Morisset
- Institut de Sciences Moléculaires d'Orsay, UMR 8214 CNRS/University Paris-Sud, University Paris-Saclay, 91405 Orsay, France
| | - Nathalie Rougeau
- Institut de Sciences Moléculaires d'Orsay, UMR 8214 CNRS/University Paris-Sud, University Paris-Saclay, 91405 Orsay, France
| | - David Lauvergnat
- Laboratoire de Chimie Physique, UMR 8000 CNRS/University Paris-Sud, University Paris-Saclay, 91405 Orsay, France
| | - Federica Agostini
- Laboratoire de Chimie Physique, UMR 8000 CNRS/University Paris-Sud, University Paris-Saclay, 91405 Orsay, France
| |
Collapse
|
6
|
Wittenbrink N, Eisfeld W. Extension of the effective relativistic coupling by asymptotic representation (ERCAR) approach to multi-dimensional potential energy surfaces: 3D model for CH3I. J Chem Phys 2018. [DOI: 10.1063/1.5011757] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Nils Wittenbrink
- Theoretische Chemie, Universität Bielefeld, Postfach 100131, D-33501 Bielefeld, Germany
| | - Wolfgang Eisfeld
- Theoretische Chemie, Universität Bielefeld, Postfach 100131, D-33501 Bielefeld, Germany
| |
Collapse
|
7
|
Takayanagi T, Nakatomi T. Automated reaction path searches for spin-forbidden reactions. J Comput Chem 2018; 39:1319-1326. [PMID: 29504140 DOI: 10.1002/jcc.25202] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/10/2018] [Revised: 02/02/2018] [Accepted: 02/14/2018] [Indexed: 01/09/2023]
Abstract
Many catalytic and biomolecular reactions containing transition metals involve changes in the electronic spin state. These processes are referred to as "spin-forbidden" reactions within nonrelativistic quantum mechanics framework. To understand detailed reaction mechanisms of spin-forbidden reactions, one must characterize reaction pathways on potential energy surfaces with different spin states and then identify crossing points. Here we propose a practical computational scheme, where only the lowest mixed-spin eigenstate obtained from the diagonalization of the spin-coupled Hamiltonian matrix is used in reaction path search calculations. We applied this method to the 6,4 FeO+ + H2 → 6,4 Fe+ + H2 O, 6,4 FeO+ + CH4 → 6,4 Fe+ + CH3 OH, and 7 Mn+ + OCS → 5 MnS+ + CO reactions, for which crossings between the different spin states are known to play essential roles in the overall reaction kinetics. © 2018 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Toshiyuki Takayanagi
- Department of Chemistry, Saitama University, Shimo-Okubo 255, Sakura-Ku, Saitama City, Saitama, 338-8570, Japan
| | - Taiki Nakatomi
- Department of Chemistry, Saitama University, Shimo-Okubo 255, Sakura-Ku, Saitama City, Saitama, 338-8570, Japan
| |
Collapse
|
8
|
Bie L, Liu F, Li Y, Dong T, Gao J, Du L, Yuan Q. Spin crossover dynamics studies on the thermally activated molecular oxygen binding mechanism on a model copper complex. Phys Chem Chem Phys 2018; 20:15852-15862. [DOI: 10.1039/c8cp02482k] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
The theoretical description of the primary dioxygen (O2) binding and activation step in many copper or iron enzymes, suffers from the instrinsically electronic non-adiabaticity of the spin flip events of the triplet dioxygen molecule (3O2), mediated by spin–orbit couplings.
Collapse
Affiliation(s)
- Lihua Bie
- Hubei Key Laboratory of Agricultural Bioinformatics
- College of Informatics
- Huazhong Agricultural University
- Wuhan
- P. R. China
| | - Fang Liu
- Hubei Key Laboratory of Agricultural Bioinformatics
- College of Informatics
- Huazhong Agricultural University
- Wuhan
- P. R. China
| | - Yanwei Li
- Environment Research Institute
- Shandong University
- Jinan
- P. R. China
| | - Tiange Dong
- Hubei Key Laboratory of Agricultural Bioinformatics
- College of Informatics
- Huazhong Agricultural University
- Wuhan
- P. R. China
| | - Jun Gao
- Hubei Key Laboratory of Agricultural Bioinformatics
- College of Informatics
- Huazhong Agricultural University
- Wuhan
- P. R. China
| | - Likai Du
- Hubei Key Laboratory of Agricultural Bioinformatics
- College of Informatics
- Huazhong Agricultural University
- Wuhan
- P. R. China
| | - Qiaoxia Yuan
- College of Engineering
- Huazhong Agricultural University
- Wuhan 430070
- P. R. China
| |
Collapse
|
9
|
Affiliation(s)
- Tao Zeng
- Department of Chemistry, Carleton University, Ottawa, Ontario K1S5B6, Canada
| |
Collapse
|
10
|
Wittenbrink N, Venghaus F, Williams D, Eisfeld W. A new approach for the development of diabatic potential energy surfaces: Hybrid block-diagonalization and diabatization by ansatz. J Chem Phys 2016; 145:184108. [DOI: 10.1063/1.4967258] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Nils Wittenbrink
- Theoretische Chemie, Universität Bielefeld, Postfach 100131, D-33501 Bielefeld, Germany
| | - Florian Venghaus
- Theoretische Chemie, Universität Bielefeld, Postfach 100131, D-33501 Bielefeld, Germany
| | - David Williams
- Theoretische Chemie, Universität Bielefeld, Postfach 100131, D-33501 Bielefeld, Germany
| | - Wolfgang Eisfeld
- Theoretische Chemie, Universität Bielefeld, Postfach 100131, D-33501 Bielefeld, Germany
| |
Collapse
|
11
|
Hoyer CE, Parker K, Gagliardi L, Truhlar DG. The DQ and DQΦ electronic structure diabatization methods: Validation for general applications. J Chem Phys 2016; 144:194101. [DOI: 10.1063/1.4948728] [Citation(s) in RCA: 35] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Chad E. Hoyer
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, 207 Pleasant St. SE, Minneapolis, Minnesota 55455-0431, USA
| | - Kelsey Parker
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, 207 Pleasant St. SE, Minneapolis, Minnesota 55455-0431, USA
| | - Laura Gagliardi
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, 207 Pleasant St. SE, Minneapolis, Minnesota 55455-0431, USA
| | - Donald G. Truhlar
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, 207 Pleasant St. SE, Minneapolis, Minnesota 55455-0431, USA
| |
Collapse
|
12
|
Hoyer CE, Li Manni G, Truhlar DG, Gagliardi L. Comment on "Fe2: As simple as a Herculean labour. Neutral (Fe2), cationic (Fe2(+)), and anionic (Fe2(-)) species" [J. Chem. Phys. 142, 244304 (2015)]. J Chem Phys 2016; 144:027101. [PMID: 26772588 DOI: 10.1063/1.4939241] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A recent paper on Fe2 [A. Kalemos, J. Chem. Phys. 142, 244304 (2015)] critiqued our previous work on the system [Hoyer et al., J. Chem. Phys. 141, 204309 (2014)]. In this comment, we explain the nature of our previously reported potential energy curve for Fe2 and we discuss our computed properties for Fe2. Additionally, we fix a labeling error that was present in our previous work, although this error is unrelated to the main point of discussion.
Collapse
Affiliation(s)
- Chad E Hoyer
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, 207 Pleasant St. SE, Minneapolis, Minnesota 55455-0431, USA
| | - Giovanni Li Manni
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, 207 Pleasant St. SE, Minneapolis, Minnesota 55455-0431, USA
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, 207 Pleasant St. SE, Minneapolis, Minnesota 55455-0431, USA
| | - Laura Gagliardi
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, 207 Pleasant St. SE, Minneapolis, Minnesota 55455-0431, USA
| |
Collapse
|
13
|
Cao J, Xie ZZ. Internal conversion and intersystem crossing in α,β-enones: a combination of electronic structure calculations and dynamics simulations. Phys Chem Chem Phys 2016; 18:6931-45. [DOI: 10.1039/c5cp06833a] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The geometrical constraint of the ring gives rise to a smaller spin–orbital coupling in the singlet–triplet crossing region, resulting in a lower intersystem crossing rate.
Collapse
Affiliation(s)
- Jun Cao
- Guizhou Provincial Key Laboratory of Computational Nano-material Science
- Guizhou Education University
- Guiyang
- China
- Guizhou Synergetic Innovation Center of Scientific Big Data for Advanced Manufacturing Technology
| | - Zhi-Zhong Xie
- Department of Chemistry
- School of Chemistry
- Chemical Engineering and Life Sciences
- Wuhan University of Technology
- Wuhan 430070
| |
Collapse
|
14
|
Cui G, Thiel W. Generalized trajectory surface-hopping method for internal conversion and intersystem crossing. J Chem Phys 2015; 141:124101. [PMID: 25273406 DOI: 10.1063/1.4894849] [Citation(s) in RCA: 117] [Impact Index Per Article: 13.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023] Open
Abstract
Trajectory-based fewest-switches surface-hopping (FSSH) dynamics simulations have become a popular and reliable theoretical tool to simulate nonadiabatic photophysical and photochemical processes. Most available FSSH methods model internal conversion. We present a generalized trajectory surface-hopping (GTSH) method for simulating both internal conversion and intersystem crossing processes on an equal footing. We consider hops between adiabatic eigenstates of the non-relativistic electronic Hamiltonian (pure spin states), which is appropriate for sufficiently small spin-orbit coupling. This choice allows us to make maximum use of existing electronic structure programs and to minimize the changes to available implementations of the traditional FSSH method. The GTSH method is formulated within the quantum mechanics (QM)/molecular mechanics framework, but can of course also be applied at the pure QM level. The algorithm implemented in the GTSH code is specified step by step. As an initial GTSH application, we report simulations of the nonadiabatic processes in the lowest four electronic states (S0, S1, T1, and T2) of acrolein both in vacuo and in acetonitrile solution, in which the acrolein molecule is treated at the ab initio complete-active-space self-consistent-field level. These dynamics simulations provide detailed mechanistic insight by identifying and characterizing two nonadiabatic routes to the lowest triplet state, namely, direct S1 → T1 hopping as major pathway and sequential S1 → T2 → T1 hopping as minor pathway, with the T2 state acting as a relay state. They illustrate the potential of the GTSH approach to explore photoinduced processes in complex systems, in which intersystem crossing plays an important role.
Collapse
Affiliation(s)
- Ganglong Cui
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| | - Walter Thiel
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| |
Collapse
|
15
|
Li SL, Truhlar DG, Schmidt MW, Gordon MS. Model space diabatization for quantum photochemistry. J Chem Phys 2015; 142:064106. [DOI: 10.1063/1.4907038] [Citation(s) in RCA: 41] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Shaohong L. Li
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, USA
| | - Donald G. Truhlar
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, USA
| | | | - Mark S. Gordon
- Department of Chemistry, Iowa State University, Ames, Iowa 50011, USA
| |
Collapse
|
16
|
Hoyer CE, Xu X, Ma D, Gagliardi L, Truhlar DG. Diabatization based on the dipole and quadrupole: The DQ method. J Chem Phys 2014; 141:114104. [DOI: 10.1063/1.4894472] [Citation(s) in RCA: 52] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Chad E. Hoyer
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, 207 Pleasant St. SE, Minneapolis, Minnesota 55455-0431, USA
| | - Xuefei Xu
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, 207 Pleasant St. SE, Minneapolis, Minnesota 55455-0431, USA
| | - Dongxia Ma
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, 207 Pleasant St. SE, Minneapolis, Minnesota 55455-0431, USA
| | - Laura Gagliardi
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, 207 Pleasant St. SE, Minneapolis, Minnesota 55455-0431, USA
| | - Donald G. Truhlar
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, 207 Pleasant St. SE, Minneapolis, Minnesota 55455-0431, USA
| |
Collapse
|
17
|
Martínez-Fernández L, Corral I, Granucci G, Persico M. Competing ultrafast intersystem crossing and internal conversion: a time resolved picture for the deactivation of 6-thioguanine. Chem Sci 2014. [DOI: 10.1039/c3sc52856a] [Citation(s) in RCA: 114] [Impact Index Per Article: 11.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023] Open
|
18
|
Westermann T, Eisfeld W, Manthe U. Coupled potential energy surface for the F(2P) + CH4 → HF + CH3 entrance channel and quantum dynamics of the CH4 · F− photodetachment. J Chem Phys 2013; 139:014309. [DOI: 10.1063/1.4812251] [Citation(s) in RCA: 36] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
19
|
Wittenbrink N, Ndome H, Eisfeld W. Toward spin-orbit coupled diabatic potential energy surfaces for methyl iodide using effective relativistic coupling by asymptotic representation. J Phys Chem A 2013; 117:7408-20. [PMID: 23590710 DOI: 10.1021/jp401438x] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The theoretical treatment of state-state interactions and the development of coupled multidimensional potential energy surfaces (PESs) is of fundamental importance for the theoretical investigation of nonadiabatic processes. Usually, only derivative or vibronic coupling is considered, but the presence of heavy atoms in a system can render spin-orbit (SO) coupling important as well. In the present study, we apply a new method recently developed by us (J. Chem. Phys. 2012, 136, 034103, and J. Chem. Phys. 2012, 137, 064101) to generate SO coupled diabatic PESs along the C-I dissociation coordinate for methyl iodide (CH3I). This is the first and mandatory step toward the development of fully coupled full-dimensional PESs to describe the multistate photodynamics of this benchmark system. The method we use here is based on the diabatic asymptotic representation of the molecular fine structure states and an effective relativistic coupling operator. It therefore is called effective relativistic coupling by asymptotic representation (ERCAR). This approach allows the efficient and accurate generation of fully coupled PESs including derivative and SO coupling based on high-level ab initio calculations. In this study we develop a specific ERCAR model for CH3I that so far accounts only for the C-I bond cleavage. Details of the diabatization and the accuracy of the results are investigated in comparison to reference ab initio calculations and experiments. The energies of the adiabatic fine structure states are reproduced in excellent agreement with ab initio SO-CI data. The model is also compared to available literature data, and its performance is evaluated critically. This shows that the new method is very promising for the construction of fully coupled full-dimensional PESs for CH3I to be used in future quantum dynamics studies.
Collapse
Affiliation(s)
- Nils Wittenbrink
- Theoretische Chemie, Universität Bielefeld, Postfach 100131, D-33501 Bielefeld, Germany
| | | | | |
Collapse
|
20
|
Valero R, Truhlar DG. Photochemistry in a dense manifold of electronic states: Photodissociation of CH2ClBr. J Chem Phys 2012; 137:22A539. [DOI: 10.1063/1.4747704] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022] Open
|
21
|
Granucci G, Persico M, Spighi G. Surface hopping trajectory simulations with spin-orbit and dynamical couplings. J Chem Phys 2012; 137:22A501. [DOI: 10.1063/1.4707737] [Citation(s) in RCA: 105] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
22
|
Ndome H, Eisfeld W. Spin-orbit coupled potential energy surfaces and properties using effective relativistic coupling by asymptotic representation. J Chem Phys 2012; 137:064101. [DOI: 10.1063/1.4740248] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
|
23
|
Ndome H, Welsch R, Eisfeld W. A new method to generate spin-orbit coupled potential energy surfaces: Effective relativistic coupling by asymptotic representation. J Chem Phys 2012; 136:034103. [DOI: 10.1063/1.3675846] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
24
|
Valero R, Song L, Gao J, Truhlar DG. Perspective on Diabatic Models of Chemical Reactivity as Illustrated by the Gas-Phase S(N)2 Reaction of Acetate Ion with 1,2-Dichloroethane. J Chem Theory Comput 2009; 5:1-22. [PMID: 20047005 PMCID: PMC2658610 DOI: 10.1021/ct800318h] [Citation(s) in RCA: 36] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Diabatic models are widely employed for studying chemical reactivity in condensed phases and enzymes, but there has been little discussion of the pros and cons of various diabatic representations for this purpose. Here we discuss and contrast six different schemes for computing diabatic potentials for a charge rearrangement reaction. They include (i) the variational diabatic configurations (VDC) constructed by variationally optimizing individual valence bond structures and (ii) the consistent diabatic configurations (CDC) obtained by variationally optimizing the ground-state adiabatic energy, both in the nonorthogonal molecular orbital valence bond (MOVB) method, along with the orthogonalized (iii) VDC-MOVB and (iv) CDC-MOVB models. In addition, we consider (v) the fourfold way (based on diabatic molecular orbitals and configuration uniformity), and (vi) empirical valence bond (EVB) theory. To make the considerations concrete, we calculate diabatic electronic states and diabatic potential energies along the reaction path that connects the reactant and the product ion-molecule complexes of the gas-phase bimolecular nucleophilic substitution (S(N)2) reaction of 1,2-dichloethane (DCE) with acetate ion, which is a model reaction corresponding to the reaction catalyzed by haloalkane dehalogenase. We utilize ab initio block-localized molecular orbital theory to construct the MOVB diabatic states and ab initio multi-configuration quasidegenerate perturbation theory to construct the fourfold-way diabatic states; the latter are calculated at reaction path geometries obtained with the M06-2X density functional. The EVB diabatic states are computed with parameters taken from the literature. The MOVB and fourfold-way adiabatic and diabatic potential energy profiles along the reaction path are in qualitative but not quantitative agreement with each other. In order to validate that these wave-function-based diabatic states are qualitatively correct, we show that the reaction energy and barrier for the adiabatic ground state, obtained with these methods, agree reasonably well with the results of high-level calculations using the composite G3SX and G3SX(MP3) methods and the BMC-CCSD multi-coefficient correlation method. However, a comparison of the EVB gas-phase adiabatic ground-state reaction path with those obtained from MOVB and with the fourfold way reveals that the EVB reaction path geometries show a systematic shift towards the products region, and that the EVB lowest-energy path has a much lower barrier. The free energies of solvation and activation energy in water reported from dynamical calculations based on EVB also imply a low activation barrier in the gas phase. In addition, calculations of the free energy of solvation using the recently proposed SM8 continuum solvation model with CM4M partial atomic charges lead to an activation barrier in reasonable agreement with experiment only when the geometries and the gas-phase barrier are those obtained from electronic structure calculations, i.e., methods i-v. These comparisons show the danger of basing the diabatic states on molecular mechanics without the explicit calculation of electronic wave functions. Furthermore, comparison of schemes i-v with one another shows that significantly different quantitative results can be obtained by using different methods for extracting diabatic states from wave function calculations, and it is important for each user to justify the choice of diabatization method in the context of its intended use.
Collapse
Affiliation(s)
- Rosendo Valero
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, MN 55455-0431
| | - Lingchun Song
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, MN 55455-0431
| | - Jiali Gao
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, MN 55455-0431
| | - Donald G. Truhlar
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, MN 55455-0431
| |
Collapse
|
25
|
Valero R, Truhlar DG, Jasper AW. Adiabatic States Derived from a Spin-Coupled Diabatic Transformation: Semiclassical Trajectory Study of Photodissociation of HBr and the Construction of Potential Curves for LiBr+. J Phys Chem A 2008; 112:5756-69. [DOI: 10.1021/jp800738b] [Citation(s) in RCA: 28] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Rosendo Valero
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431
| | - Donald G. Truhlar
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431
| | - Ahren W. Jasper
- Combustion Research Facility, Sandia National Laboratories, P. O. Box 969, Livermore, California 94551-0969
| |
Collapse
|
26
|
Rozgonyi T, González L. A Two-Dimensional Wavepacket Study of the Nonadiabatic Dynamics of CH2BrCl. J Phys Chem A 2008; 112:5573-81. [DOI: 10.1021/jp8011427] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Affiliation(s)
- Tamás Rozgonyi
- Institute of Structural Chemistry, Chemical Research Center, Hungarian Academy of Sciences, 1025 Budapest, Pusztaszeri út 59-67, Hungary, Institut für Physikalische Chemie, Friedrich-Schiller Universität, Jena Helmholtzweg 4, 07743 Jena, Germany, and Institut für Chemie und Biochemie, Freie Universität Berlin, Takustrasse 3, 14195 Berlin, Germany
| | - Leticia González
- Institute of Structural Chemistry, Chemical Research Center, Hungarian Academy of Sciences, 1025 Budapest, Pusztaszeri út 59-67, Hungary, Institut für Physikalische Chemie, Friedrich-Schiller Universität, Jena Helmholtzweg 4, 07743 Jena, Germany, and Institut für Chemie und Biochemie, Freie Universität Berlin, Takustrasse 3, 14195 Berlin, Germany
| |
Collapse
|