1
|
Richardson JO. Nonadiabatic Tunneling in Chemical Reactions. J Phys Chem Lett 2024; 15:7387-7397. [PMID: 38995660 DOI: 10.1021/acs.jpclett.4c01098] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/13/2024]
Abstract
Quantum tunneling can have a dramatic effect on chemical reaction rates. In nonadiabatic reactions such as electron transfers or spin crossovers, nuclear tunneling effects can be even stronger than for adiabatic proton transfers. Ring-polymer instanton theory enables molecular simulations of tunneling in full dimensionality and has been shown to be far more reliable than commonly used separable approximations. First-principles instanton calculations predict significant nonadiabatic tunneling of heavy atoms even at room temperature and give excellent agreement with experimental measurements for the intersystem crossing of two nitrenes in cryogenic matrix isolation, the spin-forbidden relaxation of photoexcited thiophosgene in the gas phase, and singlet oxygen deactivation in water at ambient conditions. Finally, an outlook of further theoretical developments is discussed.
Collapse
Affiliation(s)
- Jeremy O Richardson
- Department of Chemistry and Applied Biosciences, ETH Zurich, 8093 Zurich, Switzerland
| |
Collapse
|
2
|
Hsiao IY, Teranishi Y, Nakamura H. Classically forbidden nonadiabatic transitions in multidimensional chemical dynamics. Phys Chem Chem Phys 2024; 26:3795-3799. [PMID: 38251768 DOI: 10.1039/d3cp04794f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2024]
Abstract
An accurate method is proposed to deal with such nonadiabatic transitions as those energetically inaccessible, namely, classically forbidden transitions. This is formulated by using the corresponding Zhu-Nakamura formulas and finding the optimal paths in the classically forbidden tunneling regions that maximize the overall transition probabilities. This can be done for both the nonadiabatic tunneling type (so-called normal case in electron transfer) in which two diabatic potentials have opposite signs of slopes and the Landau-Zener type (inverted case) in which two diabatic potentials have the same sign of slopes. The method is numerically demonstrated to be useful for clarifying chemical and biological dynamics.
Collapse
Affiliation(s)
- I-Yun Hsiao
- Institute of Physics, National Yang Ming Chiao Tung University, 1001 Ta Hsueh Rd., Hsinchu, 30010, Taiwan.
| | - Yoshiaki Teranishi
- Institute of Physics, National Yang Ming Chiao Tung University, 1001 Ta Hsueh Rd., Hsinchu, 30010, Taiwan.
| | - Hiroki Nakamura
- Institute for Molecular Science, National Institutes of Natural Sciences, Okazaki, Japan.
| |
Collapse
|
3
|
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
|
4
|
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
|
5
|
Dergachev ID, Dergachev VD, Rooein M, Mirzanejad A, Varganov SA. Predicting Kinetics and Dynamics of Spin-Dependent Processes. Acc Chem Res 2023; 56:856-866. [PMID: 36926853 DOI: 10.1021/acs.accounts.2c00843] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/18/2023]
Abstract
ConspectusPredicting mechanisms and rates of nonadiabatic spin-dependent processes including photoinduced intersystem crossings, thermally activated spin-forbidden reactions, and spin crossovers in metal centers is a very active field of research. These processes play critical roles in transition-metal-based and metalloenzymatic catalysis, molecular magnets, light-harvesting materials, organic light-emitting diodes, photosensitizers for photodynamic therapy, and many other applications. Therefore, accurate modeling of spin-dependent processes in complex systems and on different time scales is important for many problems in chemistry, biochemistry, and materials sciences.Nonadiabatic statistical theory (NAST) and nonadiabatic molecular dynamics (NAMD) are two complementary approaches to modeling the kinetics and dynamics of spin-dependent processes. NAST predicts the probabilities and rate constants of nonradiative transitions between electronic states with different spin multiplicities using molecular properties at only few critical points on the potential energy surfaces (PESs), including the reactant minimum and the minimum energy crossing point (MECP) between two spin states. This makes it possible to obtain molecular properties for NAST calculations using accurate but often computationally expensive electronic structure methods, which is critical for predicting the rate constants of spin-dependent processes. Alternatively, NAST can be used to study spin-dependent processes in very large complex molecular systems using less computationally expensive electronic structure methods. The nuclear quantum effects, such as zero-point vibrational energy, tunneling, and interference between reaction paths can be easily incorporated. However, the statistical and local nature of NAST makes it more suitable for large systems and slow kinetics. In contrast, NAMD explores entire PESs of interacting electronic states, making it ideal for modeling fast barrierless spin-dependent processes. Because the knowledge of large portions of PESs is often needed, the simulations require a very large number of electronic structure calculations, which limits the NAMD applicability to relatively small molecular systems and ultrafast kinetics.In this Account, we discuss our contribution to the development of the NAST and NAMD approaches for predicting the rates and mechanism of spin-dependent processes. First, we briefly describe our NAST and NAMD implementations. The NAST implementation is an extension of the transition state theory to the processes involving two crossing potential energy surfaces of different spin multiplicities. The NAMD approach includes the trajectory surface hopping (TSH) and ab initio multiple spawning (AIMS) methods. Second, we discuss several applications of NAST and NAMD to model spin-dependent processes in different systems. The NAST applicability to large complex systems is demonstrated by the studies of the spin-forbidden isomerization of the active sites of metal-sulfur proteins. Our implementation of the MECP search algorithm within the fully ab initio fragment molecular orbital method allows applying NAST to systems with thousands of atoms, such as the solvated protein rubredoxin. Applications of NAMD to ultrafast spin-dependent processes are represented by the generalized AIMS simulations utilizing the fast GPU-based TeraChem electronic structure program to gain insight into the complex photoexcited state relaxation in 2-cyclopentenone.
Collapse
Affiliation(s)
- Ilya D Dergachev
- Department of Chemistry, University of Nevada, Reno, 1664 N. Virginia Street, Reno, Nevada 89557-0216, United States
| | - Vsevolod D Dergachev
- Department of Chemistry, University of Nevada, Reno, 1664 N. Virginia Street, Reno, Nevada 89557-0216, United States
| | - Mitra Rooein
- Department of Chemistry, University of Nevada, Reno, 1664 N. Virginia Street, Reno, Nevada 89557-0216, United States
| | - Amir Mirzanejad
- Department of Chemistry, University of Nevada, Reno, 1664 N. Virginia Street, Reno, Nevada 89557-0216, United States
| | - Sergey A Varganov
- Department of Chemistry, University of Nevada, Reno, 1664 N. Virginia Street, Reno, Nevada 89557-0216, United States
| |
Collapse
|
6
|
Rooein M, Varganov SA. How to calculate the rate constants for nonradiative transitions between the MS components of spin multiplets? Mol Phys 2022. [DOI: 10.1080/00268976.2022.2116364] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/14/2022]
Affiliation(s)
- Mitra Rooein
- Department of Chemistry, University of Nevada, Reno, NV, USA
| | | |
Collapse
|
7
|
Mirzanejad A, Varganov SA. The role of the intermediate triplet state in iron-catalyzed multi-state C-H activation. Phys Chem Chem Phys 2022; 24:20721-20727. [PMID: 36018581 DOI: 10.1039/d2cp02733j] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Efficient activation and functionalization of the C-H bond under mild conditions are of a great interest in chemical synthesis. We investigate the previously proposed spin-accelerated activation of the C(sp2)-H bond by a Fe(II)-based catalyst to clarify the role of the intermediate triplet state in the reaction mechanism. High-level electronic structure calculations on a small model of a catalytic system utilizing the coupled cluster with the single, double, and perturbative triple excitations [CCSD(T)] are used to select the density functional for the full-size model. Our analysis indicates that the previously proposed two-state quintet-singlet reaction pathway is unlikely to be efficient due to a very weak spin-orbit coupling between these two spin states. We propose a more favorable multi-state quintet-triplet-singlet reaction pathway and discuss the importance of the intermediate triplet state. This triplet state facilitates a spin-accelerated reaction mechanism by strongly coupling to both quintet and singlet states. Our calculations show that the C-H bond activation through the proposed quintet-triplet-singlet reaction pathway is more thermodynamically favorable than the single-state quintet and two-state singlet-quintet mechanisms.
Collapse
Affiliation(s)
- Amir Mirzanejad
- Department of Chemistry, University of Nevada, Reno, 1664 N. Virginia Street, Reno, NV 89557-0216, USA.
| | - Sergey A Varganov
- Department of Chemistry, University of Nevada, Reno, 1664 N. Virginia Street, Reno, NV 89557-0216, USA.
| |
Collapse
|
8
|
Heller ER, Richardson JO. Heavy‐Atom Quantum Tunnelling in Spin Crossovers of Nitrenes**. Angew Chem Int Ed Engl 2022; 61:e202206314. [PMID: 35698730 PMCID: PMC9540336 DOI: 10.1002/anie.202206314] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/03/2022] [Indexed: 01/01/2023]
Abstract
We simulate two recent matrix‐isolation experiments at cryogenic temperatures, in which a nitrene undergoes spin crossover from its triplet state to a singlet state via quantum tunnelling. We detail the failure of the commonly applied weak‐coupling method (based on a linear approximation of the potentials) in describing these deep‐tunnelling reactions. The more rigorous approach of semiclassical golden‐rule instanton theory in conjunction with double‐hybrid density‐functional theory and multireference perturbation theory does, however, provide rate constants and kinetic isotope effects in good agreement with experiment. In addition, these calculations locate the optimal tunnelling pathways, which provide a molecular picture of the reaction mechanism. The reactions involve substantial heavy‐atom quantum tunnelling of carbon, nitrogen and oxygen atoms, which unexpectedly even continues to play a role at room temperature.
Collapse
Affiliation(s)
- Eric R. Heller
- Laboratory of Physical Chemistry ETH Zürich 8093 Zürich Switzerland
| | | |
Collapse
|
9
|
Heller ER, Richardson JO. Heavy‐Atom Quantum Tunnelling in Spin Crossovers of Nitrenes. Angew Chem Int Ed Engl 2022. [DOI: 10.1002/ange.202206314] [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)
- Eric R Heller
- Eidgenossische Technische Hochschule Zurich Lab. Physical Chemistry SWITZERLAND
| | - Jeremy O Richardson
- Eidgenössische Technische Hochschule Zürich Lab. Physical Chemistry Vladimir-Prelog-Weg 2 8093 Zurich SWITZERLAND
| |
Collapse
|
10
|
Ansari IM, Heller ER, Trenins G, Richardson JO. Instanton theory for Fermi's golden rule and beyond. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2022; 380:20200378. [PMID: 35341312 PMCID: PMC8958279 DOI: 10.1098/rsta.2020.0378] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/16/2021] [Accepted: 09/22/2021] [Indexed: 05/29/2023]
Abstract
Instanton theory provides a semiclassical approximation for computing quantum tunnelling effects in complex molecular systems. It is typically applied to proton-transfer reactions for which the Born-Oppenheimer approximation is valid. However, many processes in physics, chemistry and biology, such as electron transfers, are non-adiabatic and are correctly described instead using Fermi's golden rule. In this work, we discuss how instanton theory can be generalized to treat these reactions in the golden-rule limit. We then extend the theory to treat fourth-order processes such as bridge-mediated electron transfer and apply the method to simulate an electron moving through a model system of three coupled quantum dots. By comparison with benchmark quantum calculations, we demonstrate that the instanton results are much more reliable than alternative approximations based on superexchange-mediated effective coupling or a classical sequential mechanism. This article is part of the theme issue 'Chemistry without the Born-Oppenheimer approximation'.
Collapse
Affiliation(s)
| | - Eric R. Heller
- Laboratory of Physical Chemistry, ETH, Zürich, Switzerland
| | - George Trenins
- Laboratory of Physical Chemistry, ETH, Zürich, Switzerland
| | | |
Collapse
|
11
|
NAST: Nonadiabatic Statistical Theory Package for Predicting Kinetics of Spin-Dependent Processes. Top Curr Chem (Cham) 2022; 380:15. [PMID: 35201520 DOI: 10.1007/s41061-022-00366-w] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/03/2021] [Accepted: 01/15/2022] [Indexed: 10/19/2022]
Abstract
We present a nonadiabatic statistical theory (NAST) package for predicting kinetics of spin-dependent processes, such as intersystem crossings, spin-forbidden unimolecular reactions, and spin crossovers. The NAST package can calculate the probabilities and rates of transitions between the electronic states of different spin multiplicities. Both the microcanonical (energy-dependent) and canonical (temperature-dependent) rate constants can be obtained. Quantum effects, including tunneling, zero-point vibrational energy, and reaction path interference, can be accounted for. In the limit of an adiabatic unimolecular reaction proceeding on a single electronic state, NAST reduces to the traditional transition state theory. Because NAST requires molecular properties at only a few points on potential energy surfaces, it can be applied to large molecular systems, used with accurate high-level electronic structure methods, and employed to study slow nonadiabatic processes. The essential NAST input data include the nuclear Hessian at the reactant minimum, as well as the nuclear Hessians, energy gradients, and spin-orbit coupling at the minimum energy crossing point (MECP) between two states. The additional computational tools included in the NAST package can be used to extract the required input data from the output files of electronic structure packages, calculate the effective Hessian at the MECP, and fit the reaction coordinate for more advanced NAST calculations. We describe the theory, its implementation, and three examples of application to different molecular systems.
Collapse
|
12
|
Heller ER, Richardson JO. Spin Crossover of Thiophosgene via Multidimensional Heavy-Atom Quantum Tunneling. J Am Chem Soc 2021; 143:20952-20961. [PMID: 34846871 DOI: 10.1021/jacs.1c10088] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
The spin-crossover reaction of thiophosgene has drawn broad attention from both experimenters and theoreticians as a prime example of radiationless intramolecular decay via intersystem crossing. Despite multiple attempts over 20 years, theoretical predictions have typically been orders of magnitude in error relative to the experimentally measured triplet lifetime. We address the T1 → S0 transition by the first application of semiclassical golden-rule instanton theory in conjunction with on-the-fly electronic-structure calculations based on multireference perturbation theory. Our first-principles approach provides excellent agreement with the experimental rates. This was only possible because instanton theory goes beyond previous methods by locating the optimal tunneling pathway in full dimensionality and thus captures "corner cutting" effects. Since the reaction is situated in the Marcus inverted regime, the tunneling mechanism can be interpreted in terms of two classical trajectories, one traveling forward and one backward in imaginary time, which are connected by particle-antiparticle creation and annihilation events. The calculated mechanism indicates that the spin crossover is sped up by many orders of magnitude due to multidimensional quantum tunneling of the carbon atom even at room temperature.
Collapse
Affiliation(s)
- Eric R Heller
- Laboratory of Physical Chemistry, ETH Zürich, 8093 Zürich, Switzerland
| | | |
Collapse
|
13
|
Hermsmeier R, Kłos J, Kotochigova S, Tscherbul TV. Quantum Spin State Selectivity and Magnetic Tuning of Ultracold Chemical Reactions of Triplet Alkali-Metal Dimers with Alkali-Metal Atoms. PHYSICAL REVIEW LETTERS 2021; 127:103402. [PMID: 34533330 DOI: 10.1103/physrevlett.127.103402] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/13/2021] [Revised: 06/29/2021] [Accepted: 07/28/2021] [Indexed: 06/13/2023]
Abstract
We demonstrate that it is possible to efficiently control ultracold chemical reactions of alkali-metal atoms colliding with open-shell alkali-metal dimers in their metastable triplet states by choosing the internal hyperfine and rovibrational states of the reactants as well as by inducing magnetic Feshbach resonances with an external magnetic field. We base these conclusions on coupled-channel statistical calculations that include the effects of hyperfine contact and magnetic-field-induced Zeeman interactions on ultracold chemical reactions of hyperfine-resolved ground-state Na and the triplet NaLi(a^{3}Σ^{+}) producing singlet Na_{2}(^{1}Σ_{g}^{+}) and a Li atom. We find that the reaction rates are sensitive to the initial hyperfine states of the reactants. The chemical reaction of fully spin-polarized, high-spin states of rotationless NaLi(a^{3}Σ^{+},v=0,N=0) molecules with fully spin-polarized Na is suppressed by a factor of 10-100 compared to that of unpolarized reactants. We interpret these findings within the adiabatic state model, which treats the reaction as a sequence of nonadiabatic transitions between the initial nonreactive high-spin state and the final low-spin states of the reaction complex. In addition, we show that magnetic Feshbach resonances can similarly change reaction rate coefficients by several orders of magnitude. Some of these resonances are due to resonant trimer bound states dissociating to the N=2 rotational state of NaLi(a^{3}Σ^{+},v=0) and would thus exist in systems without hyperfine interactions.
Collapse
Affiliation(s)
| | - Jacek Kłos
- Department of Physics, Joint Quantum Institute, University of Maryland College Park, College Park, Maryland 20742, USA
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, USA
| | | | - Timur V Tscherbul
- Department of Physics, University of Nevada, Reno, Nevada 89557, USA
| |
Collapse
|
14
|
Shchepanovska D, Shannon RJ, Curchod BFE, Glowacki DR. Nonadiabatic Kinetics in the Intermediate Coupling Regime: Comparing Molecular Dynamics to an Energy-Grained Master Equation. J Phys Chem A 2021; 125:3473-3488. [PMID: 33880919 DOI: 10.1021/acs.jpca.1c01260] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
We propose and test an extension of the energy-grained master equation (EGME) for treating nonadiabatic (NA) hopping between different potential energy surfaces, which enables us to model the competition between stepwise collisional relaxation and kinetic processes which transfer population between different electronic states of the same spin symmetry. By incorporating Zhu-Nakamura theory into the EGME, we are able to treat NA passages beyond the simple Landau-Zener approximation, along with the corresponding treatments of zero-point energy and tunneling probability. To evaluate the performance of this NA-EGME approach, we carried out detailed studies of the UV photodynamics of the volatile organic compound C6-hydroperoxy aldehyde (C6-HPALD) using on-the-fly ab initio molecular dynamics and trajectory surface hopping. For this multichromophore molecule, we show that the EGME is able to capture important aspects of the dynamics, including kinetic timescales, and diabatic trapping. Such an approach provides a promising and efficient strategy for treating the long-time dynamics of photoexcited molecules in regimes which are difficult to capture using atomistic on-the-fly molecular dynamics.
Collapse
Affiliation(s)
| | - Robin J Shannon
- School of Chemistry, University of Bristol, Bristol BS8 1TS, U.K
| | | | - David R Glowacki
- School of Chemistry, University of Bristol, Bristol BS8 1TS, U.K.,Intangible Realities Laboratory, University of Bristol, Bristol BS8 1UB, U.K.,Department of Computer Science, University of Bristol, Bristol BS8 1UB, U.K
| |
Collapse
|
15
|
Lykhin AO, Truhlar DG, Gagliardi L. Role of Triplet States in the Photodynamics of Aniline. J Am Chem Soc 2021; 143:5878-5889. [PMID: 33843225 DOI: 10.1021/jacs.1c00989] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The dynamics of excited heteroaromatic molecules is a key to understanding the photoprotective properties of many biologically relevant chromophores that dissipate their excitation energy nonreactively and thereby prevent the detrimental effects of ultraviolet radiation. Despite their structural variability, most substituted aromatic compounds share a common feature of a repulsive 1πσ* potential energy surface. This surface can lead to photoproducts, and it can also facilitate the population transfer back to the ground electronic state by means of a 1πσ*/S0 conical intersection. Here, we explore a hidden relaxation route involving the triplet electronic state of aniline, which has recently been discovered by means of time-selected photofragment translational spectroscopy [J. Chem. Phys. 2019, 151, 141101]. By using the recently available analytical gradients for multiconfiguration pair-density functional theory, it is now possible to locate the minimum-energy crossing points between states of different spin and therefore compute the intersystem crossing rates with a multireference method, rather than with the less reliable single-reference methods. Using such calculations, we demonstrate that the population loss of aniline in the T1(3ππ*) state is dominated by C6H5NH2 → C6H5NH· + H· dissociation, and we explain the long nonradiative lifetimes of the T1(3ππ*) state at the excitation wavelengths of 294-264 nm.
Collapse
Affiliation(s)
- Aleksandr O Lykhin
- Department of Chemistry, Pritzker School of Molecular Engineering, The James Franck Institute and Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, United States
| | - Laura Gagliardi
- Department of Chemistry, Pritzker School of Molecular Engineering, The James Franck Institute and Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
| |
Collapse
|
16
|
Abstract
In this article, we review nonadiabatic molecular dynamics (NAMD) methods for modeling spin-crossover transitions. First, we discuss different representations of electronic states employed in the grid-based and direct NAMD simulations. The nature of interstate couplings in different representations is highlighted, with the main focus on nonadiabatic and spin-orbit couplings. Second, we describe three NAMD methods that have been used to simulate spin-crossover dynamics, including trajectory surface hopping, ab initio multiple spawning, and multiconfiguration time-dependent Hartree. Some aspects of employing different electronic structure methods to obtain information about potential energy surfaces and interstate couplings for NAMD simulations are also discussed. Third, representative applications of NAMD to spin crossovers in molecular systems of different sizes and complexities are highlighted. Finally, we pose several fundamental questions related to spin-dependent processes. These questions should be possible to address with future methodological developments in NAMD.
Collapse
Affiliation(s)
- Saikat Mukherjee
- Institut de Chimie Radicalaire, CNRS 7273, Aix-Marseille University, 13013 Marseille, France;
| | - Dmitry A Fedorov
- Oak Ridge Associated Universities, Oak Ridge, Tennessee 37830, USA;
| | - Sergey A Varganov
- Department of Chemistry, University of Nevada, Reno, Nevada 89557-0216, USA;
| |
Collapse
|
17
|
Farcaş AA, Bende A. Theoretical modeling of the singlet-triplet spin transition in different Ni(II)-diketo-pyrphyrin-based metal-ligand octahedral complexes. Phys Chem Chem Phys 2021; 23:4784-4795. [PMID: 33599640 DOI: 10.1039/d0cp05366j] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The structural stability, charge transfer effects and strength of the spin-orbit couplings in different Ni(ii)-ligand complexes have been studied at the DFT (B3LYP and CAM-B3LYP) and coupled cluster (DLPNO-CCSD(T)) levels of theory. Accordingly, two different, porphyrin- and diketo-pyrphyrin-based four-coordination macrocycles as planar ligands as well as pyridine (or pyrrole) and mesylate anion molecular groups as vertical ligands were considered in order to build metal-organic complexes with octahedral coordination configurations. For each molecular system, the identification of equilibrium geometries and the intersystem crossing (the minimum energy crossing) points between the potential energy surfaces of the singlet and triplet spin states is followed by computing the spin-orbit couplings between the two spin states. Structures, based on the diketo-pyrphyrin macrocycle as the planar ligand, show stronger six-coordination metal-organic complexes due to the extra electrostatic interaction between the positively charged central metal cation and the negatively charged vertical ligands. The results also show that the magnitude of the spin-orbit coupling is influenced by the atomic positions of deprotonations of the ligands, and implicitly the direction of the charge transfer between the ligand and the central metal ion.
Collapse
Affiliation(s)
- Alex-Adrian Farcaş
- Faculty of Physics, "Babeş-Bolyai" University, Mihail Kogalniceanu Street No. 1, Ro-400084 Cluj-Napoca, Romania
| | - Attila Bende
- Molecular and Biomolecular Physics Department, National Institute for Research and Development of Isotopic and Molecular Technologies, Donat Street, No. 67-103, Ro-400293 Cluj-Napoca, Romania.
| |
Collapse
|
18
|
Viegas LP, M. Nunes C, Fausto R. Spin-forbidden heavy-atom tunneling in the ring-closure of triplet cyclopentane-1,3-diyl. Phys Chem Chem Phys 2021; 23:5797-5803. [DOI: 10.1039/d1cp00076d] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The putative spin-forbidden heavy-atom tunneling process for the ring closure of cyclopentane-1,3-diyl at cryogenic temperatures is confirmed with calculations employing the weak-coupling formulation of nonadiabatic transition state theory.
Collapse
Affiliation(s)
- Luís P. Viegas
- University of Coimbra
- CQC
- Department of Chemistry
- 3004-535 Coimbra
- Portugal
| | - Cláudio M. Nunes
- University of Coimbra
- CQC
- Department of Chemistry
- 3004-535 Coimbra
- Portugal
| | - Rui Fausto
- University of Coimbra
- CQC
- Department of Chemistry
- 3004-535 Coimbra
- Portugal
| |
Collapse
|
19
|
Ortega P, Zanchet A, Sanz-Sanz C, Gómez-Carrasco S, González-Sánchez L, Jambrina PG. DpgC-Catalyzed Peroxidation of 3,5-Dihydroxyphenylacetyl-CoA (DPA-CoA): Insights into the Spin-Forbidden Transition and Charge Transfer Mechanisms*. Chemistry 2020; 27:1700-1712. [PMID: 32975323 DOI: 10.1002/chem.202002993] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/22/2020] [Indexed: 11/06/2022]
Abstract
Despite being a very strong oxidizing agent, most organic molecules are not oxidized in the presence of O2 at room temperature because O2 is a diradical whereas most organic molecules are closed-shell. Oxidation then requires a change in the spin state of the system, which is forbidden according to non-relativistic quantum theory. To overcome this limitation, oxygenases usually rely on metal or redox cofactors to catalyze the incorporation of, at least, one oxygen atom into an organic substrate. However, some oxygenases do not require any cofactor, and the detailed mechanism followed by these enzymes remains elusive. To fill this gap, here the mechanism for the enzymatic cofactor-independent oxidation of 3,5-dihydroxyphenylacetyl-CoA (DPA-CoA) is studied by combining multireference calculations on a model system with QM/MM calculations. Our results reveal that intersystem crossing takes place without requiring the previous protonation of molecular oxygen. The characterization of the electronic states reveals that electron transfer is concomitant with the triplet-singlet transition. The enzyme plays a passive role in promoting the intersystem crossing, although spontaneous reorganization of the water wire connecting the active site with the bulk presets the substrate for subsequent chemical transformations. The results show that the stabilization of the singlet radical-pair between dioxygen and enolate is enough to promote spin-forbidden reaction without the need for neither metal cofactors nor basic residues in the active site.
Collapse
Affiliation(s)
- Pablo Ortega
- Departamento de Química Física, University of Salamanca, Salamanca, 37008, Spain
| | - Alexandre Zanchet
- Departamento de Química Física, University of Salamanca, Salamanca, 37008, Spain.,Instituto de Física Fundamental (CSIC), Madrid, 28006, Spain
| | - Cristina Sanz-Sanz
- Departamento de Química Física Aplicada, University Autónoma de Madrid, Madrid, 28049, Spain
| | | | | | - Pablo G Jambrina
- Departamento de Química Física, University of Salamanca, Salamanca, 37008, Spain
| |
Collapse
|
20
|
Pohlman AJ, Kaliakin DS, Varganov SA, Casey SM. Spin controlled surface chemistry: alkyl desorption from Si(100)-2×1 by nonadiabatic hydrogen elimination. Phys Chem Chem Phys 2020; 22:16641-16647. [PMID: 32661543 DOI: 10.1039/d0cp01913e] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
An understanding of the role that spin states play in semiconductor surface chemical reactions is currently limited. Herein, we provide evidence of a nonadiabatic reaction involving a localized singlet to triplet thermal excitation of the Si(100) surface dimer dangling bond. By comparing the β-hydrogen elimination kinetics of ethyl adsorbates probed by thermal desorption experiments to electronic structure calculation results, we determined that a coverage-dependent change in mechanism occurs. At low coverage, a nonadiabatic, inter-dimer mechanism is dominant, while adiabatic mechanisms become dominant at higher coverage. Computational results indicate that the spin crossover is rapid near room temperature and the nonadiabatic path is accelerated by a barrier that is 40 kJ mol-1 less than the adiabatic path. Simulated thermal desorption reactions using nonadiabatic transition state theory (NA-TST) for the surface dimer intersystem crossing are in close agreement with experimental observations.
Collapse
Affiliation(s)
- Andrew J Pohlman
- Department of Chemistry and Chemical Physics Program, University of Nevada, 1664 N. Virginia Street, Reno, Nevada 89557-0216, USA.
| | | | | | | |
Collapse
|
21
|
Gacka E, Burdzinski G, Marciniak B, Kubas A, Lewandowska-Andralojc A. Interaction of light with a non-covalent zinc porphyrin-graphene oxide nanohybrid. Phys Chem Chem Phys 2020; 22:13456-13466. [PMID: 32519730 DOI: 10.1039/d0cp02545c] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The present study explores the influence of graphene oxide (GO) on deactivation pathways of the excited states of zinc 5,10,15,20-tetrakis(4-(hydroxyphenyl))porphyrin (ZnTPPH). The interaction of light with free ZnTPPH molecules and with ZnTPPH molecules adsorbed on graphene oxide sheets was probed via UV-vis spectroscopy, fluorescence spectroscopy, femtosecond pump-probe technique and nanosecond flash photolysis. Formation of the ground-state ZnTPPH-GO complex in solution was monitored by the red-shift of the porphyrin Soret absorption band. It was found that Stern-Volmer fluorescence quenching can be described in terms of two different quenching regimes depending on the GO concentration. In addition, our comprehensive analysis of the steady-state and time-resolved emission experiments led to the conclusion that the observed quenching was entirely attributable to a static mechanism. Laser flash photolysis showed that the triplet lifetime of the ZnTPPH increased in the presence of GO from 174 μs to 292 μs, which is related to the decrease in the rate constant of a radiationless decay mechanism involving rotation of the peripheral hydroxyphenyl rings of the porphyrin. Femtosecond transient absorption spectroscopy demonstrated the presence of a fast photoinduced electron transfer from the singlet excited state of ZnTPPH to the GO sheets, as indicated by the formation of a porphyrin radical cation. Quantum chemical calculations were used to gain deeper insights into the nature of the electronically excited states in the free ZnTPPH as well as in the ZnTPPH-GO complex.
Collapse
Affiliation(s)
- Ewelina Gacka
- Faculty of Chemistry, Adam Mickiewicz University, Uniwersytetu Poznanskiego 8, 61-614 Poznan, Poland.
| | - Gotard Burdzinski
- Faculty of Physics, Adam Mickiewicz University, Uniwersytetu Poznanskiego 2, 61-614 Poznan, Poland
| | - Bronisław Marciniak
- Faculty of Chemistry, Adam Mickiewicz University, Uniwersytetu Poznanskiego 8, 61-614 Poznan, Poland. and Center for Advanced Technology, Adam Mickiewicz University, Uniwersytetu Poznanskiego 10, 61-614 Poznan, Poland
| | - Adam Kubas
- Institute of Physical Chemistry, Polish Academy of Sciences, Kasprzaka 44/52, 01-224 Warsaw, Poland.
| | - Anna Lewandowska-Andralojc
- Faculty of Chemistry, Adam Mickiewicz University, Uniwersytetu Poznanskiego 8, 61-614 Poznan, Poland. and Center for Advanced Technology, Adam Mickiewicz University, Uniwersytetu Poznanskiego 10, 61-614 Poznan, Poland
| |
Collapse
|