1
|
Padula D. Discriminating Clockwise and Counterclockwise Photoisomerization Paths in Achiral Photoswitches by Excited-State Electronic Circular Dichroism. J Phys Chem B 2024. [PMID: 39171863 DOI: 10.1021/acs.jpcb.4c02384] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/23/2024]
Abstract
Despite the numerous investigations of photoisomerization reactions from both the computational and experimental points of view, even in complex environments, to date there is no direct demonstration of the direction of rotation of the retinal chromophore, initiating the vision process in several organisms, occurring upon light irradiation. In the literature, many proposals have been formulated to shed light on the details of this process, most of which are extracted from semiclassical simulations. Although high hopes are held in the development of time-resolved X-ray spectroscopy, I argue in this work that simpler but less known techniques can be used to unravel the details of this fascinating photochemical process. In fact, chiroptical spectroscopy would unambiguously prove the direction of the rotatory motion of the chromophore during the photoisomerization process by probing excited state chirality, a piece of information that, so far, has been exclusively extracted from atomistic simulations. I demonstrate this statement by computing the expected chiroptical response along photoisomerization pathways for several models of the retinal chromophores that are found in nature bound to rhodopsins, including nuclear ensemble spectra from semiclassical dynamics simulations, that can be compared with time-resolved experiments.
Collapse
Affiliation(s)
- Daniele Padula
- Dipartimento di Biotecnologie, Chimica e Farmacia, Università di Siena, Via A. Moro 2, Siena 53100, Italy
| |
Collapse
|
2
|
Kabir M, Ghosh P, Gozem S. Electronic Structure Methods for Simulating Flavin's Spectroscopy and Photophysics: Comparison of Multi-reference, TD-DFT, and Single-Reference Wave Function Methods. J Phys Chem B 2024; 128:7545-7557. [PMID: 39074870 PMCID: PMC11317985 DOI: 10.1021/acs.jpcb.4c03748] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/05/2024] [Revised: 07/14/2024] [Accepted: 07/22/2024] [Indexed: 07/31/2024]
Abstract
The use of flavins and flavoproteins in photocatalytic, sensing, and biotechnological applications has led to a growing interest in computationally modeling the excited-state electronic structure and photophysics of flavin. However, there is limited consensus regarding which computational methods are appropriate for modeling flavin's photophysics. We compare the energies of low-lying excited states of flavin computed with time-dependent density functional theory (TD-DFT), equation-of-motion coupled cluster (EOM-EE-CCSD), scaled opposite-spin configuration interaction [SOS-CIS(D)], multiconfiguration pair-density functional theory (MC-PDFT), and several multireference perturbation theory (MR-PT2) methods. In the first part, we focus on excitation energies of the first singlet excited state (S1) of five different redox and protonation states of flavin, with the goal of finding a suitable active space for MR-PT2 calculations. In the second part, we construct two sets of one-dimensional potential energy surfaces connecting the S0 and S1 equilibrium geometries (S0-S1 path) and the S1 (π,π*) and S2 (n,π*) equilibrium geometries (S1-S2 path). The first path therefore follows a Franck-Condon active mode of flavin while the second path maps crossings points between low-lying singlet and triplet states in flavin. We discuss the similarities and differences in the TD-DFT, EOM-EE-CCSD, SOS-CIS(D), MC-PDFT and MR-PT2 energy profiles along these paths. We find that (TD-)DFT methods are suitable for applications such as simulating the spectra of flavins but are inconsistent with several other methods when used for some geometry optimizations and when describing the energetics of dark (n,π*) states. MR-PT2 methods show promise for the simulation of flavin's low-lying excited states, but the selection of orbitals for the active space and the number of roots used for state averaging must be done carefully to avoid artifacts. Some properties, such as the intersystem crossing geometry and energy between the S1 (π,π*) and T2 (n,π*) states, may require additional benchmarking before they can be determined quantitatively.
Collapse
Affiliation(s)
- Mohammad
Pabel Kabir
- Department of Chemistry, Georgia State University, Atlanta, Georgia 30302, United States
| | - Paulami Ghosh
- Department of Chemistry, Georgia State University, Atlanta, Georgia 30302, United States
| | - Samer Gozem
- Department of Chemistry, Georgia State University, Atlanta, Georgia 30302, United States
| |
Collapse
|
3
|
Oliveira LMF, Valverde D, Costa GJ, Borin AC. The copious photochemistry of 2,6-diaminopurine: Luminescence, triplet population, and ground state recovery. Photochem Photobiol 2024; 100:323-338. [PMID: 37403286 DOI: 10.1111/php.13833] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/01/2023] [Revised: 06/15/2023] [Accepted: 06/19/2023] [Indexed: 07/06/2023]
Abstract
9H- and 7H-2,6-Diaminopurine (26DAP) photoinduced events in vacuum were studied at the MS-CASPT2/cc-pVDZ level of theory. The S1 1 (ππ* La ) state is initially populated evolving barrierless towards its minimum energy structure, from where two photochemical events can take place in both tautomers. The first is the return of the electronic population to the ground state via the C6 conical intersection (CI-C6). The second involves an internal conversion to the ground through the C2 conical intersection (CI-C2). According to our geodesic interpolated paths connecting the critical structures, the second route is less favorable in both tautomers, due to the presence of high energy barriers. Our calculations suggest a competition between fluorescence and ultrafast relaxation to the electronic ground state via internal conversion process. Based on our calculated potential energy surfaces and experimental excited state lifetimes from the literature, we can infer that the 7H- must have a greater fluorescence yield than the 9H-tautomer. We also explored the triplet state population mechanisms on the 7H-26DAP to understand their long-lived components observed experimentally.
Collapse
Affiliation(s)
- Leonardo M F Oliveira
- Department of Fundamental Chemistry, Institute of Chemistry, University of São Paulo, São Paulo, Brazil
| | - Danillo Valverde
- Unité de Chimie Physique Théorique et Structurale, Namur Institute of Structured Matter, Université de Namur, Namur, Belgium
| | - Gustavo J Costa
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas, USA
| | - Antonio Carlos Borin
- Department of Fundamental Chemistry, Institute of Chemistry, University of São Paulo, São Paulo, Brazil
| |
Collapse
|
4
|
Gate G, Williams A, Boldissar S, Šponer J, Szabla R, de Vries M. The tautomer-specific excited state dynamics of 2,6-diaminopurine using resonance-enhanced multiphoton ionization and quantum chemical calculations. Photochem Photobiol 2024; 100:404-418. [PMID: 38124372 DOI: 10.1111/php.13897] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2023] [Revised: 10/25/2023] [Accepted: 11/14/2023] [Indexed: 12/23/2023]
Abstract
2,6-Diaminopurine (2,6-dAP) is an alternative nucleobase that potentially played a role in prebiotic chemistry. We studied its excited state dynamics in the gas phase by REMPI, IR-UV hole burning, and ps pump-probe spectroscopy and performed quantum chemical calculations at the SCS-ADC(2) level of theory to interpret the experimental results. We found the 9H tautomer to have a small barrier to ultrafast relaxation via puckering of its 6-membered ring. The 7H tautomer has a larger barrier to reach a conical intersection and also has a sizable triplet yield. These results are discussed relative to other purines, for which 9H tautomerization appears to be more photostable than 7H and homosubstituted purines appear to be less photostable than heterosubstituted or singly substituted purines.
Collapse
Affiliation(s)
- Gregory Gate
- Department of Chemistry and Biochemistry, University of California, Santa Barbara, California, USA
| | - Ann Williams
- Department of Chemistry and Biochemistry, University of California, Santa Barbara, California, USA
| | - Samuel Boldissar
- Department of Chemistry and Biochemistry, University of California, Santa Barbara, California, USA
| | - Jiří Šponer
- Institute of Biophysics of the Czech Academy of Sciences, Brno, Czech Republic
- Regional Centre of Advanced Technologies and Materials, Czech Advanced Technology and Research Institute (CATRIN), Palacky University Olomouc, Olomouc-Holice, Czech Republic
| | - Rafal Szabla
- Institute of Advanced Materials, Faculty of Chemistry, Wrocław University of Science and Technology, Wrocław, Poland
| | - Mattanjah de Vries
- Department of Chemistry and Biochemistry, University of California, Santa Barbara, California, USA
| |
Collapse
|
5
|
Papineau TV, Jacquemin D, Vacher M. Which Electronic Structure Method to Choose in Trajectory Surface Hopping Dynamics Simulations? Azomethane as a Case Study. J Phys Chem Lett 2024; 15:636-643. [PMID: 38205955 DOI: 10.1021/acs.jpclett.3c03014] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2024]
Abstract
Nonadiabatic dynamics simulations have become a standard approach to explore photochemical reactions. Such simulations require underlying potential energy surfaces and couplings between them, calculated at a chosen level of theory, yet this aspect is rarely assessed. Here, in combination with the popular trajectory surface hopping dynamics method, we use a high-accuracy XMS-CASPT2 electronic structure level as a benchmark for assessing the performances of various post-Hartree-Fock methods (namely, CIS, ADC(2), CC2, and CASSCF) and exchange-correlation functionals (PBE, PBE0, and CAM-B3LYP) in a TD-DFT/TDA context, using the isomerization around a double bond as test case. Different relaxation pathways are identified, and the ability of the different methods to reproduce their relative importance and time scale is discussed. The results show that multireference electronic structure methods should be preferred, when studying nonadiabatic decay between excited and ground states. If not affordable, TD-DFT with TDA and hybrid functionals and ADC(2) are efficient alternatives but overestimate the nonradiative decay yield and thus may miss deexcitation pathways.
Collapse
Affiliation(s)
| | - Denis Jacquemin
- Nantes Université, CNRS, CEISAM UMR 6230, Nantes F-44000, France
- Institut Universitaire de France (IUF), Paris F-75005, France
| | - Morgane Vacher
- Nantes Université, CNRS, CEISAM UMR 6230, Nantes F-44000, France
| |
Collapse
|
6
|
Taylor JT, Tozer DJ, Curchod BFE. On the description of conical intersections between excited electronic states with LR-TDDFT and ADC(2). J Chem Phys 2023; 159:214115. [PMID: 38059547 DOI: 10.1063/5.0176140] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/11/2023] [Accepted: 11/14/2023] [Indexed: 12/08/2023] Open
Abstract
Conical intersections constitute the conceptual bedrock of our working understanding of ultrafast, nonadiabatic processes within photochemistry (and photophysics). Accurate calculation of potential energy surfaces within the vicinity of conical intersections, however, still poses a serious challenge to many popular electronic structure methods. Multiple works have reported on the deficiency of methods like linear-response time-dependent density functional theory within the adiabatic approximation (AA LR-TDDFT) or algebraic diagrammatic construction to second-order [ADC(2)]-approaches often used in excited-state molecular dynamics simulations-to describe conical intersections between the ground and excited electronic states. In the present study, we focus our attention on conical intersections between excited electronic states and probe the ability of AA LR-TDDFT and ADC(2) to describe their topology and topography, using protonated formaldimine and pyrazine as two exemplar molecules. We also take the opportunity to revisit the performance of these methods in describing conical intersections involving the ground electronic state in protonated formaldimine-highlighting in particular how the intersection ring exhibited by AA LR-TDDFT can be perceived either as a (near-to-linear) seam of intersection or two interpenetrating cones, depending on the magnitude of molecular distortions within the branching space.
Collapse
Affiliation(s)
- Jack T Taylor
- Department of Chemistry, Durham University, South Road, Durham DH1 3LE, United Kingdom
| | - David J Tozer
- Department of Chemistry, Durham University, South Road, Durham DH1 3LE, United Kingdom
| | - Basile F E Curchod
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Cantock's Close, Bristol BS8 1TS, United Kingdom
| |
Collapse
|
7
|
Barneschi L, Kaliakin D, Huix-Rotllant M, Ferré N, Filatov Gulak M, Olivucci M. Assessment of the Electron Correlation Treatment on the Quantum-Classical Dynamics of Retinal Protonated Schiff Base Models: XMS-CASPT2, RMS-CASPT2, and REKS Methods. J Chem Theory Comput 2023; 19:8189-8200. [PMID: 37937990 DOI: 10.1021/acs.jctc.3c00879] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2023]
Abstract
We compare the performance of three different multiconfigurational wave function-based electronic structure methods and two implementations of the spin-restricted ensemble-referenced Kohn-Sham (REKS) method. The study is characterized by three features: (i) it uses a small set of quantum-classical trajectories rather than potential energy surface mapping, (ii) it focuses, exclusively, on the photoisomerization of retinal protonated Schiff base models, and (iii) it probes the effect of both methyl substitution and the increase in length of the conjugate π-system. For each tested method, the corresponding analytical gradients are used to drive the quantum-classical (Tully's FSSH method) trajectory propagation, including the recent multistate XMS-CASPT2 and RMS-CASPT2 gradients. It is shown that while CASSCF, XMS-CASPT2, and RMS-CASPT2 yield consistent photoisomerization dynamics descriptions, REKS produces, in some of these systems, qualitatively different behavior that is attributed to a flatter and topographically different excited state potential energy surface. The origin of this behavior can be traced back to the effect of the employed density functional approximation. The above studies are further expanded by benchmarking, at the CASSCF and REKS levels, the electronic structure methods using a QM/MM model of the visual pigment rhodopsin.
Collapse
Affiliation(s)
- Leonardo Barneschi
- Dipartimento di Biotecnologie, Chimica e Farmacia, Università di Siena, I-53100 Siena, Italy
| | - Danil Kaliakin
- Department of Chemistry, Bowling Green State University, Bowling Green, Ohio 43403, United States
| | - Miquel Huix-Rotllant
- Aix-Marseille Université, CNRS, Institut Chimie Radicalaire, 13013 Marseille, France
| | - Nicolas Ferré
- Aix-Marseille Université, CNRS, Institut Chimie Radicalaire, 13013 Marseille, France
| | - Michael Filatov Gulak
- Department of Chemistry, Southern Methodist University, Dallas, Texas 75275, United States
| | - Massimo Olivucci
- Dipartimento di Biotecnologie, Chimica e Farmacia, Università di Siena, I-53100 Siena, Italy
- Department of Chemistry, Bowling Green State University, Bowling Green, Ohio 43403, United States
| |
Collapse
|
8
|
Wen J, Mai S, González L. Excited-State Dynamics Simulations of a Light-Driven Molecular Motor in Solution. J Phys Chem A 2023; 127:9520-9529. [PMID: 37917883 PMCID: PMC10658450 DOI: 10.1021/acs.jpca.3c05841] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2023] [Revised: 10/14/2023] [Accepted: 10/17/2023] [Indexed: 11/04/2023]
Abstract
Molecular motors, where light can be transformed into motion, are promising in the design of nanomechanical devices. For applications, however, finding relationships between molecular motion and the environment is important. Here, we report the study of excited-state dynamics of an overcrowded alkene in solution using a hybrid quantum mechanics/molecular mechanics (QM/MM) approach combined with excited-state molecular dynamics simulations. Using QM/MM surface-hopping trajectories, we calculated time-resolved emission and transient absorption spectra. These show the rise of a short-lived Franck-Condon state, followed by the formation of a dark state in the first 150 fs before the molecular motor relaxes to the ground state in about 1 ps. From the analysis of radial distribution functions, we infer that the orientation of the solvent with respect to the molecular motor in the electronic excited state is similar to that in the ground state during the photoisomerization.
Collapse
Affiliation(s)
- Jin Wen
- State
Key Laboratory for Modification of Chemical Fibers and Polymer Materials,
College of Materials Science and Engineering, Donghua University, Shanghai 201620, China
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Str. 17, Vienna 1090, Austria
| | - Sebastian Mai
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Str. 17, Vienna 1090, Austria
| | - Leticia González
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Str. 17, Vienna 1090, Austria
| |
Collapse
|
9
|
Komarov K, Park W, Lee S, Huix-Rotllant M, Choi CH. Doubly Tuned Exchange-Correlation Functionals for Mixed-Reference Spin-Flip Time-Dependent Density Functional Theory. J Chem Theory Comput 2023; 19:7671-7684. [PMID: 37844129 DOI: 10.1021/acs.jctc.3c00884] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2023]
Abstract
It is demonstrated that significant accuracy improvements in MRSF-TDDFT can be achieved by introducing two different exchange-correlation (XC) functionals for the reference Kohn-Sham DFT and the response part of the calculations, respectively. Accordingly, two new XC functionals of doubly tuned Coulomb attenuated method-vertical excitation energy (DTCAM-VEE) and DTCAM-AEE were developed on the basis of the "adaptive exact exchange (AEE)" concept in the framework of the Coulomb-attenuating XC functionals. The values by DTCAM-VEE are in excellent agreement with those of Thiel's set [mean absolute errors (MAEs) and the interquartile range (IQR) values of 0.218 and 0.327 eV, respectively]. On the other hand, DTCAM-AEE faithfully reproduced the qualitative aspects of conical intersections (CIs) of trans-butadiene and thymine and the nonadiabatic molecular dynamics (NAMD) simulations on thymine. The latter functional also remarkably exhibited the exact 1/R asymptotic behavior of the charge-transfer state of an ethylene-tetrafluoroethylene dimer and the accurate potential energy surfaces (PESs) along the two torsional angles of retinal protonated Schiff base model with six double bonds (rPSB6). Overall, DTCAM-AEE generally performs well, as its MAE (0.237) and IQR (0.41 eV) are much improved as compared to BH&HLYP. The current idea can also be applied to other XC functionals as well as other variants of linear response theories, opening a new way of developing XC functionals.
Collapse
Affiliation(s)
- Konstantin Komarov
- Center for Quantum Dynamics, Pohang University of Science and Technology, Pohang 37673, South Korea
| | - Woojin Park
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| | - Seunghoon Lee
- Department of Chemistry, Seoul National University, Seoul, 151-747, South Korea
| | | | - Cheol Ho Choi
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| |
Collapse
|
10
|
Park W, Komarov K, Lee S, Choi CH. Mixed-Reference Spin-Flip Time-Dependent Density Functional Theory: Multireference Advantages with the Practicality of Linear Response Theory. J Phys Chem Lett 2023; 14:8896-8908. [PMID: 37767969 PMCID: PMC10561896 DOI: 10.1021/acs.jpclett.3c02296] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/16/2023] [Accepted: 09/21/2023] [Indexed: 09/29/2023]
Abstract
The density functional theory (DFT) and linear response (LR) time-dependent (TD)-DFT are of the utmost importance for routine computations. However, the single reference formulation of DFT suffers in the description of open-shell singlet systems such as diradicals and bond-breaking. LR-TDDFT, on the other hand, finds difficulties in the modeling of conical intersections, doubly excited states, and core-level excitations. In this Perspective, we demonstrate that many of these limitations can be overcome by recently developed mixed-reference (MR) spin-flip (SF)-TDDFT, providing an alternative yet accurate route for such challenging situations. Empowered by the practicality of the LR formalism, it is anticipated that MRSF-TDDFT can become one of the major workhorses for general routine tasks.
Collapse
Affiliation(s)
- Woojin Park
- Department
of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| | - Konstantin Komarov
- Center
for Quantum Dynamics, Pohang University
of Science and Technology, Pohang 37673, South Korea
| | - Seunghoon Lee
- Division
of Chemistry and Chemical Engineering, California
Institute of Technology, Pasadena, California 91125, United States
| | - Cheol Ho Choi
- Department
of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| |
Collapse
|
11
|
Mu D, Li QS. A theoretical study on the photochemical generation of phenylborylene from phenyldiazidoborane. Phys Chem Chem Phys 2023; 25:8074-8081. [PMID: 36876656 DOI: 10.1039/d2cp05349g] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/17/2023]
Abstract
Organic borylenes are a kind of highly reactive species, which play important roles in a lot of reactions as vigorous intermediates. In this work, we investigated the photochemical generation mechanisms of phenylborylene (PhB) together with the side product N-phenylnitrenoiminoborane (PhNBN) from phenyldiazidoborane (PhBN6) by extrusion of dinitrogen in the two lowest electronic singlet states (S0 and S1) based on the complete active space self-consistent field (CASSCF) and its second-order perturbation (CASPT2) methods combined with time-dependent density functional theory (TD-DFT) calculations. Our results show that the reaction PhBN6 → PhB + 3N2 involves stepwise N2 extrusion three times and the azido region rearrangement. Moreover, we found that the studied photo-induced processes are kinetically feasible because the highest energy barrier is only 0.36 eV and excitation with light of wavelength 254 nm can provide enough excess energy to overcome these energy barriers. Importantly, we revealed that several conical intersections between S1 and S0 states participate and facilitate the studied photochemical processes. Our results not only clarify the experimental observations, (H. F. Bettinger, J. Am. Chem. Soc. 2006, 128, 2534), but also provide valuable insights into borylene chemistry.
Collapse
Affiliation(s)
- Di Mu
- Key Laboratory of Cluster Science of Ministry of Education, Beijing Key Laboratory of Photoelectronic/Electrophotonic Conversion Materials, School of Chemistry and Chemical Engineering, Beijing Institute of Technology, Beijing, 100081, China.
| | - Quan-Song Li
- Key Laboratory of Cluster Science of Ministry of Education, Beijing Key Laboratory of Photoelectronic/Electrophotonic Conversion Materials, School of Chemistry and Chemical Engineering, Beijing Institute of Technology, Beijing, 100081, China.
| |
Collapse
|
12
|
Hutton L, Curchod BFE. Photodynamics of Gas‐Phase Pyruvic Acid Following Light Absorption in the Actinic Region. CHEMPHOTOCHEM 2022. [DOI: 10.1002/cptc.202200151] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Affiliation(s)
- Lewis Hutton
- Centre for Computational Chemistry School of Chemistry University of Bristol Bristol BS8 1TS UK
| | - Basile F. E. Curchod
- Centre for Computational Chemistry School of Chemistry University of Bristol Bristol BS8 1TS UK
| |
Collapse
|
13
|
Valverde D, Mai S, Canuto S, Borin AC, González L. Ultrafast Intersystem Crossing Dynamics of 6-Selenoguanine in Water. JACS AU 2022; 2:1699-1711. [PMID: 35911449 PMCID: PMC9327080 DOI: 10.1021/jacsau.2c00250] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/22/2023]
Abstract
Rationalizing the photochemistry of nucleobases where an oxygen is replaced by a heavier atom is essential for applications that exploit near-unity triplet quantum yields. Herein, we report on the ultrafast excited-state deactivation mechanism of 6-selenoguanine (6SeGua) in water by combining nonadiabatic trajectory surface-hopping dynamics with an electrostatic embedding quantum mechanics/molecular mechanics (QM/MM) scheme. We find that the predominant relaxation mechanism after irradiation starts on the bright singlet S2 state that converts internally to the dark S1 state, from which the population is transferred to the triplet T2 state via intersystem crossing and finally to the lowest T1 state. This S2 → S1 → T2 → T1 deactivation pathway is similar to that observed for the lighter 6-thioguanine (6tGua) analogue, but counterintuitively, the T1 lifetime of the heavier 6SeGua is shorter than that of 6tGua. This fact is explained by the smaller activation barrier to reach the T1/S0 crossing point and the larger spin-orbit couplings of 6SeGua compared to 6tGua. From the dynamical simulations, we also calculate transient absorption spectra (TAS), which provide two time constants (τ1 = 131 fs and τ2 = 191 fs) that are in excellent agreement with the experimentally reported value (τexp = 130 ± 50 fs) (Farrel et al. J. Am. Chem. Soc. 2018, 140, 11214). Intersystem crossing itself is calculated to occur with a time scale of 452 ± 38 fs, highlighting that the TAS is the result of a complex average of signals coming from different nonradiative processes and not intersystem crossing alone.
Collapse
Affiliation(s)
- Danillo Valverde
- Department
of Fundamental Chemistry, Institute of Chemistry, University of São Paulo, Avenida Professor Lineu Prestes, 748, São Paulo, São Paulo CEP 05508-000, Brazil
- Institute
of Physics, University of São Paulo, Rua do Matão 1371, São Paulo, São Paulo CEP 05508-090, Brazil
| | - Sebastian Mai
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, Vienna 1090, Austria
| | - Sylvio Canuto
- Institute
of Physics, University of São Paulo, Rua do Matão 1371, São Paulo, São Paulo CEP 05508-090, Brazil
| | - Antonio Carlos Borin
- Department
of Fundamental Chemistry, Institute of Chemistry, University of São Paulo, Avenida Professor Lineu Prestes, 748, São Paulo, São Paulo CEP 05508-000, Brazil
| | - Leticia González
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, Vienna 1090, Austria
| |
Collapse
|
14
|
Milovanović B, Novak J, Etinski M, Domcke W, Došlić N. On the propensity of formation of cyclobutane dimers in face-to-face and face-to-back uracil stacks in solution. Phys Chem Chem Phys 2022; 24:14836-14845. [PMID: 35697028 DOI: 10.1039/d2cp00495j] [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
UV irradiation of RNA leads to the formation of intra- and inter-strand crosslinks of cyclobutane type. Despite the importance of this reaction, relatively little is known about how the mutual orientation of the two bases affects the outcome of the reaction. Here we report a comparative nonadiabatic molecular dynamics study of face-to-back (F2B) and face-to-face (F2F) stacked uracil-water clusters. The computations were performed using the second-order algebraic-diagrammatic-construction (ADC(2)) method. We found that F2B stacked uracil-water clusters either relax non-reactively to the ground state by an ethylenic twist around the CC bond or remain in the lowest nπ* state in which the two bases gradually move away from each other. This finding is consistent with the low propensity for the formation of intra-strand cyclobutane dimers between adjacent RNA bases. On the contrary, in F2F stacked uracil-water clusters, in addition to non-reactive deactivation, we found a pro-reactive deactivation pathway, which may lead to the formation of cyclobutane uracil dimers in the electronic ground state. On a qualitative level, the observed photodynamics of F2F stacked uracil-water clusters explains the greater propensity of RNA to form inter-strand cyclobutane-type crosslinks.
Collapse
Affiliation(s)
| | - Jurica Novak
- Department of Biotechnology, University of Rijeka, HR-51000 Rijeka, Croatia.,Scientific and Educational Center "Biomedical Technologies" School of Medical Biology, South Ural State University, RU-454080, Chelyabinsk, Russia.,Center for Artificial Intelligence and Cybersecurity, University of Rijeka, 51000 Rijeka, Croatia
| | - Mihajlo Etinski
- University of Belgrade, Faculty of Physical Chemistry, Belgrade, Serbia
| | - Wolfgang Domcke
- Department of Chemistry, Technical University of Munich, Garching, Germany
| | - Nađa Došlić
- Department of Physical Chemistry, Ruđer Bošković Institute, Zagreb, Croatia.
| |
Collapse
|
15
|
Huang X, Domcke W. Ab initio trajectory surface-hopping dynamics studies of excited-state proton-coupled electron transfer reactions in trianisoleheptazine-phenol complexes. Phys Chem Chem Phys 2022; 24:15925-15936. [PMID: 35726762 DOI: 10.1039/d2cp01262f] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The excited-state proton-coupled electron-transfer (PCET) reaction in hydrogen-bonded complexes of trianisoleheptazine (TAHz), a chromophore related to polymeric carbon nitrides widely used in hydrogen-evolution photocatalysis, with several phenol derivatives were recently studied by Schlenker and coworkers with time-resolved photoluminescence quenching and pump-probe experiments. A pronounced dependence of the PCET reactivity on the electron-donating/electron-withdrawing character of the substituents on phenol was found, with indications of a barrierless or nearly barrierless PCET reaction for the most strongly electron-donating substituent, methoxy. In the present work, the excited-state PCET dynamics was explored with first-principles nonadiabatic dynamics simulations using the TDDFT/ωB97X-D electronic-structure model for two selected complexes, TAHz-phenol and TAHz-methoxyphenol. The qualitative reliability of the TDDFT/ωB97X-D electronic-structure model was assessed by extensive benchmarking of excitation energies and potential-energy profiles against a wave-function-based ab initio method, the algebraic-diagrammatic construction of second order (ADC(2)). The nonadiabatic dynamics simulations provide temporally and structurally resolved insights into paradigmatic PCET reactions in TAHz-phenol complexes. The radiationless relaxation of the photoexcited bright 1ππ* state to the long-lived dark S1 state of TAHz occurs in less than 100 fs. The ensuing PCET reaction on the adiabatic S1 surface is faster in TAHz-methoxyphenol complexes than in TAHz-phenol complexes due to a lower H-atom-transfer barrier, as observed in the experiments. The relaxation of the complexes to the electronic ground state is found to occur exclusively via PCET within the 250 fs time window covered by the present simulations, confirming the essential role of the hydrogen bond for the fluorescence quenching process. The absolute values of the computed PCET time constants are significantly shorter than those extracted from time-resolved photoluminescence measurements for mixtures of TAHz with phenolic substrates in toluene. The possible origins of this discrepancy are discussed.
Collapse
Affiliation(s)
- Xiang Huang
- Department of Chemistry, Technical University of Munich, D-85747 Garching, Germany.
| | - Wolfgang Domcke
- Department of Chemistry, Technical University of Munich, D-85747 Garching, Germany.
| |
Collapse
|
16
|
Liu MK, Li J, Li QS, Li ZS. Theoretical insights into photo-induced isomerization mechanisms of phenylsulfinyl radical PhSO˙. Phys Chem Chem Phys 2022; 24:6266-6273. [PMID: 35229842 DOI: 10.1039/d1cp05331k] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Sulfinyl radicals (R-SO˙) play important roles in lots of reactions, while the isomer oxathiyl radicals (R-OS˙) and the isomerization between them are rarely observed due to the poor stability of R-OS˙. In this work, the complete active space self-consistent field (CASSCF) and its multi-state second order perturbation (MS-CASPT2) methods were employed to study the photo-induced reaction mechanisms of phenylsulfinyl radical PhSO˙ 1 and its isomer phenoxathiyl radical PhOS˙ 2. Our results show that 1 and 2 have similar singly occupied molecular orbitals in the ground state but different properties in the excited state, which determine their diverse behaviors after irradiation. Radical 1 can generate 2 by light irradiation, but 2 produces isomerization product 3 (2-hydroxyphenylthiyl radical) and ring-opening product 4 (acyclic thioketoketene radical) in two paths via S atom migration intermediate Int1 (2-carbonylcyclohexadienthiyl radical). The former path involves consequent hydrogen shift reactions with a strongly exothermic process while the latter path involves both ring-expansion and ring-opening processes with a high barrier, resulting in a structural and energetic preference for the former path. Moreover, we revealed several conical intersections that participate in the reactions and facilitate the photochemical processes. Our calculations not only remain consistent with and clarify the experimental observations (X. Zeng, et al., J. Am. Chem. Soc., 2018, 140(31), 9972-9978) but also enrich the knowledge of sulfinyl radicals and isomer oxathiyl radicals.
Collapse
Affiliation(s)
- Meng-Kai Liu
- Key Laboratory of Cluster Science of Ministry of Education, Beijing Key Laboratory of Photoelectronic/Electrophotonic Conversion Materials, School of Chemistry and Chemical Engineering, Beijing Institute of Technology, Beijing 100081, China.
| | - Jian Li
- Key Laboratory of Cluster Science of Ministry of Education, Beijing Key Laboratory of Photoelectronic/Electrophotonic Conversion Materials, School of Chemistry and Chemical Engineering, Beijing Institute of Technology, Beijing 100081, China.
| | - Quan-Song Li
- Key Laboratory of Cluster Science of Ministry of Education, Beijing Key Laboratory of Photoelectronic/Electrophotonic Conversion Materials, School of Chemistry and Chemical Engineering, Beijing Institute of Technology, Beijing 100081, China.
| | - Ze-Sheng Li
- Key Laboratory of Cluster Science of Ministry of Education, Beijing Key Laboratory of Photoelectronic/Electrophotonic Conversion Materials, School of Chemistry and Chemical Engineering, Beijing Institute of Technology, Beijing 100081, China.
| |
Collapse
|
17
|
Understanding the Photolysis of CH 3ONO 2 with the On-the-fly Nonadiabatic Dynamics Simulation at the ADC(2) Level. CHINESE J CHEM PHYS 2022. [DOI: 10.1063/1674-0068/cjcp2201013] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022]
|
18
|
Huang X, Xie W, Došlić N, Gelin MF, Domcke W. Ab Initio Quasiclassical Simulation of Femtosecond Time-Resolved Two-Dimensional Electronic Spectra of Pyrazine. J Phys Chem Lett 2021; 12:11736-11744. [PMID: 34851116 DOI: 10.1021/acs.jpclett.1c03589] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
Two-dimensional (2D) electronic spectroscopy is a powerful nonlinear technique which provides spectroscopic information on two frequency axes as well as dynamical information as a function of the so-called waiting time. Herein, an ab initio theoretical framework for the simulation of electronic 2D spectra has been developed. The method is based on the classical approximation to the doorway-window representation of three-pulse photon-echo signals and the description of nuclear motion by classical trajectories. Nonadiabatic effects are taken into account by a trajectory surface-hopping algorithm. 2D electronic spectra were simulated with ab initio on-the-fly trajectory calculations using the ADC(2) electronic-structure method for the pyrazine molecule, which is a benchmark system for ultrafast radiationless decay through conical intersections. It is demonstrated that 2D spectroscopy with subfemtosecond UV pulses can provide unprecedented detailed information on the ultrafast photodynamics of polyatomic molecules.
Collapse
Affiliation(s)
- Xiang Huang
- Department of Chemistry, Technical University of Munich, Garching, D-85747, Germany
| | - Weiwei Xie
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, Karlsruhe, 76131, Germany
| | - Nađa Došlić
- Department of Physical Chemistry, Ruder Boscovic Institute, Zagreb, HR-10000, Croatia
| | - Maxim F Gelin
- School of Sciences, Hangzhou Dianzi University, Hangzhou, 310018, China
| | - Wolfgang Domcke
- Department of Chemistry, Technical University of Munich, Garching, D-85747, Germany
| |
Collapse
|
19
|
de Castro Araujo Valente D, Borges I, Cardozo TM. Nonradiative relaxation mechanisms of the elusive silole molecule. Phys Chem Chem Phys 2021; 23:26561-26574. [PMID: 34811562 DOI: 10.1039/d1cp03803f] [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
Silole derivatives have been extensively employed for developing organic optoelectronics, but few studies focused on the photophysical properties of the silole molecule. In this work, we investigate these properties by computing the absorption spectra and performing nonadiabatic molecular dynamics of silole employing the algebraic diagrammatic construction [ADC(2)] and extended multi-state XMS-CASPT2 ab initio electronic structure methods. For vertical excitations and excited state optimizations, the equation of motion coupled-cluster singles and doubles (EOM-CCSD) was also used. The nuclear ensemble and the fewest-switches surface hopping molecular dynamics methods coupled with the first two high-level electronic structure methods were applied to probe the relaxation mechanisms of silole. We could reproduce the experimental first absorption maximum value and found an ultrafast relaxation process occurring exclusively through ring-puckering distortions without breaking ring bonds or hydrogen elimination. Minimum energy conical intersection optimizations were carried out and potential energy curves, including triplet states, were calculated to further elucidate the relaxation process of silole.
Collapse
Affiliation(s)
| | - Itamar Borges
- Departamento de Química, Instituto Militar de Engenharia, Rio de Janeiro, RJ, Brazil
| | - Thiago Messias Cardozo
- Instituto de Química, Universidade Federal do Rio de Janeiro, Rio de Janeiro, RJ, Brazil.
| |
Collapse
|
20
|
Wang J, Durbeej B. Thermal Fluctuations in Conjugation and their Effect on Calculated Excitation Energies: A Case Study on the Astaxanthin Carotenoid. CHEMPHOTOCHEM 2021. [DOI: 10.1002/cptc.202100178] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Affiliation(s)
- Jun Wang
- Jiangsu Key Laboratory for Chemistry of Low-Dimensional Materials Jiangsu Engineering Laboratory for Environment Functional Materials School of Chemistry and Chemical Engineering Huaiyin Normal University No. 111 West Changjiang Road 223300 Huaian Jiangsu Province China
| | - Bo Durbeej
- Division of Theoretical Chemistry IFM Linköping University 581 83 Linköping Sweden
| |
Collapse
|
21
|
Kochman MA, Durbeej B, Kubas A. Simulation and Analysis of the Transient Absorption Spectrum of 4-( N, N-Dimethylamino)benzonitrile (DMABN) in Acetonitrile. J Phys Chem A 2021; 125:8635-8648. [PMID: 34550700 PMCID: PMC8503879 DOI: 10.1021/acs.jpca.1c06166] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/10/2021] [Revised: 09/03/2021] [Indexed: 12/17/2022]
Abstract
4-(N,N-Dimethylamino)benzonitrile (DMABN) is a well-known model compound for dual fluorescence-in sufficiently polar solvents, it exhibits two distinct fluorescence emission bands. The interpretation of its transient absorption (TA) spectrum in the visible range is the subject of a long-standing controversy. In the present study, we resolve this issue by calculating the TA spectrum on the basis of nonadiabatic molecular dynamics simulations. An unambiguous assignment of spectral signals to specific excited-state structures is achieved by breaking down the calculated spectrum into contributions from twisted and nontwisted molecular geometries. In particular, the much-discussed excited-state absorption band near 1.7 eV (ca. 700 nm) is attributed to the near-planar locally excited (LE) minimum on the S1 state. On the technical side, our study demonstrates that the second-order approximate coupled cluster singles and doubles (CC2) method can be used successfully to calculate the TA spectra of moderately large organic molecules, provided that the system in question does not approach a crossing between the lowest excited state and the singlet ground state within the time frame of the simulation.
Collapse
Affiliation(s)
- Michał Andrzej Kochman
- Institute
of Physical Chemistry, Polish Academy of
Sciences, Ul. Marcina Kasprzaka 44/52, 01-224 Warszawa, Poland
| | - Bo Durbeej
- Division
of Theoretical Chemistry, Department of Physics, Chemistry and Biology
(IFM), Linköping University, 581 83 Linköping, Sweden
| | - Adam Kubas
- Institute
of Physical Chemistry, Polish Academy of
Sciences, Ul. Marcina Kasprzaka 44/52, 01-224 Warszawa, Poland
| |
Collapse
|
22
|
Jones CM, List NH, Martínez TJ. Resolving the ultrafast dynamics of the anionic green fluorescent protein chromophore in water. Chem Sci 2021; 12:11347-11363. [PMID: 34667545 PMCID: PMC8447926 DOI: 10.1039/d1sc02508b] [Citation(s) in RCA: 28] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/06/2021] [Accepted: 07/13/2021] [Indexed: 11/21/2022] Open
Abstract
The chromophore of the green fluorescent protein (GFP) is critical for probing environmental influences on fluorescent protein behavior. Using the aqueous system as a bridge between the unconfined vacuum system and a constricting protein scaffold, we investigate the steric and electronic effects of the environment on the photodynamical behavior of the chromophore. Specifically, we apply ab initio multiple spawning to simulate five picoseconds of nonadiabatic dynamics after photoexcitation, resolving the excited-state pathways responsible for internal conversion in the aqueous chromophore. We identify an ultrafast pathway that proceeds through a short-lived (sub-picosecond) imidazolinone-twisted (I-twisted) species and a slower (several picoseconds) channel that proceeds through a long-lived phenolate-twisted (P-twisted) intermediate. The molecule navigates the non-equilibrium energy landscape via an aborted hula-twist-like motion toward the one-bond-flip dominated conical intersection seams, as opposed to following the pure one-bond-flip paths proposed by the excited-state equilibrium picture. We interpret our simulations in the context of time-resolved fluorescence experiments, which use short- and long-time components to describe the fluorescence decay of the aqueous GFP chromophore. Our results suggest that the longer time component is caused by an energetically uphill approach to the P-twisted intersection seam rather than an excited-state barrier to reach the twisted intramolecular charge-transfer species. Irrespective of the location of the nonadiabatic population events, the twisted intersection seams are inefficient at facilitating isomerization in aqueous solution. The disordered and homogeneous nature of the aqueous solvent environment facilitates non-selective stabilization with respect to I- and P-twisted species, offering an important foundation for understanding the consequences of selective stabilization in heterogeneous and rigid protein environments.
Collapse
Affiliation(s)
- Chey M Jones
- Department of Chemistry and the PULSE Institute, Stanford University Stanford CA 94305 USA
- SLAC National Accelerator Laboratory 2575 Sand Hill Road Menlo Park CA 94025 USA
| | - Nanna H List
- Department of Chemistry and the PULSE Institute, Stanford University Stanford CA 94305 USA
- SLAC National Accelerator Laboratory 2575 Sand Hill Road Menlo Park CA 94025 USA
| | - Todd J Martínez
- Department of Chemistry and the PULSE Institute, Stanford University Stanford CA 94305 USA
- SLAC National Accelerator Laboratory 2575 Sand Hill Road Menlo Park CA 94025 USA
| |
Collapse
|
23
|
Jankowska J, Sobolewski AL. Modern Theoretical Approaches to Modeling the Excited-State Intramolecular Proton Transfer: An Overview. Molecules 2021; 26:molecules26175140. [PMID: 34500574 PMCID: PMC8434569 DOI: 10.3390/molecules26175140] [Citation(s) in RCA: 22] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/15/2021] [Revised: 08/20/2021] [Accepted: 08/23/2021] [Indexed: 02/02/2023] Open
Abstract
The excited-state intramolecular proton transfer (ESIPT) phenomenon is nowadays widely acknowledged to play a crucial role in many photobiological and photochemical processes. It is an extremely fast transformation, often taking place at sub-100 fs timescales. While its experimental characterization can be highly challenging, a rich manifold of theoretical approaches at different levels is nowadays available to support and guide experimental investigations. In this perspective, we summarize the state-of-the-art quantum-chemical methods, as well as molecular- and quantum-dynamics tools successfully applied in ESIPT process studies, focusing on a critical comparison of their specific properties.
Collapse
Affiliation(s)
- Joanna Jankowska
- Faculty of Chemistry, University of Warsaw, 02-093 Warsaw, Poland
- Correspondence:
| | | |
Collapse
|
24
|
Matsika S. Electronic Structure Methods for the Description of Nonadiabatic Effects and Conical Intersections. Chem Rev 2021; 121:9407-9449. [PMID: 34156838 DOI: 10.1021/acs.chemrev.1c00074] [Citation(s) in RCA: 59] [Impact Index Per Article: 19.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
Nonadiabatic effects are ubiquitous in photophysics and photochemistry, and therefore, many theoretical developments have been made to properly describe them. Conical intersections are central in nonadiabatic processes, as they promote efficient and ultrafast nonadiabatic transitions between electronic states. A proper theoretical description requires developments in electronic structure and specifically in methods that describe conical intersections between states and nonadiabatic coupling terms. This review focuses on the electronic structure aspects of nonadiabatic processes. We discuss the requirements of electronic structure methods to describe conical intersections and nonadiabatic couplings, how the most common excited state methods perform in describing these effects, and what the recent developments are in expanding the methodology and implementing nonadiabatic couplings.
Collapse
Affiliation(s)
- Spiridoula Matsika
- Department of Chemistry, Temple University, Philadelphia, Pennsylvania 19122, United States
| |
Collapse
|
25
|
Marsili E, Prlj A, Curchod BFE. Caveat when using ADC(2) for studying the photochemistry of carbonyl-containing molecules. Phys Chem Chem Phys 2021; 23:12945-12949. [PMID: 34085679 PMCID: PMC8207513 DOI: 10.1039/d1cp02185k] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2021] [Accepted: 05/28/2021] [Indexed: 12/21/2022]
Abstract
Several electronic-structure methods are available to study the photochemistry and photophysics of organic molecules. Among them, ADC(2) stands as a sweet spot between computational efficiency and accuracy. As a result, ADC(2) has recently seen its number of applications booming, in particular to unravel the deactivation pathways and photodynamics of organic molecules. Despite this growing success, we demonstrate here that care has to be taken when studying the nonradiative pathways of carbonyl-containing molecules, as ADC(2) appears to suffer from a systematic flaw.
Collapse
Affiliation(s)
| | - Antonio Prlj
- Department of Chemistry, Durham University, Durham DH1 3LE, UK.
| | | |
Collapse
|
26
|
Tapavicza E, von Rudorff GF, De Haan DO, Contin M, George C, Riva M, von Lilienfeld OA. Elucidating an Atmospheric Brown Carbon Species-Toward Supplanting Chemical Intuition with Exhaustive Enumeration and Machine Learning. ENVIRONMENTAL SCIENCE & TECHNOLOGY 2021; 55:8447-8457. [PMID: 34080853 DOI: 10.1021/acs.est.1c00885] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
Brown carbon (BrC) is involved in atmospheric light absorption and climate forcing and can cause adverse health effects. Understanding the formation mechanisms and molecular structure of BrC is of key importance in developing strategies to control its environment and health impact. Structure determination of BrC is challenging, due to the lack of experiments providing molecular fingerprints and the sheer number of molecular candidates with identical mass. Suggestions based on chemical intuition are prone to errors due to the inherent bias. We present an unbiased algorithm, using graph-based molecule generation and machine learning, which can identify all molecular structures of compounds involved in biomass burning and the composition of BrC. We apply this algorithm to C12H12O7, a light-absorbing "test case" molecule identified in chamber experiments on the aqueous photo-oxidation of syringol, a prevalent marker in wood smoke. Of the 260 million molecular graphs, the algorithm leaves only 36,518 (0.01%) as viable candidates matching the spectrum. Although no unique molecular structure is obtained from only a chemical formula and a UV/vis absorption spectrum, we discuss further reduction strategies and their efficacy. With additional data, the method can potentially more rapidly identify isomers extracted from lab and field aerosol particles without introducing human bias.
Collapse
Affiliation(s)
- Enrico Tapavicza
- Department of Chemistry and Biochemistry, California State University, Long Beach, 1250 Bellflower Boulevard, Long Beach, California 90840, United States
| | - Guido Falk von Rudorff
- Faculty of Physics, University of Vienna, Kolingasse 14-16, AT-1090 Wien, Austria
- Institute of Physical Chemistry and National Center for Computational Design and Discovery of Novel Materials (MARVEL), Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| | - David O De Haan
- Department of Chemistry and Biochemistry, University of San Diego, 5998 Alcala Park, San Diego, California 92110, United States
| | - Mario Contin
- Facultad de Farmacia y Bioquímica, Departamento de Química Analitica y Fisicoquímica, Universidad de Buenos Aires, Junín 956, Buenos Aires C1113AAD, Argentina
| | - Christian George
- Université Lyon, Université Claude Bernard Lyon 1, CNRS, IRCELYON, 69626 Villeurbanne, France
| | - Matthieu Riva
- Université Lyon, Université Claude Bernard Lyon 1, CNRS, IRCELYON, 69626 Villeurbanne, France
| | - O Anatole von Lilienfeld
- Faculty of Physics, University of Vienna, Kolingasse 14-16, AT-1090 Wien, Austria
- Institute of Physical Chemistry and National Center for Computational Design and Discovery of Novel Materials (MARVEL), Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| |
Collapse
|
27
|
Szabla R, Zdrowowicz M, Spisz P, Green NJ, Stadlbauer P, Kruse H, Šponer J, Rak J. 2,6-diaminopurine promotes repair of DNA lesions under prebiotic conditions. Nat Commun 2021; 12:3018. [PMID: 34021158 PMCID: PMC8139960 DOI: 10.1038/s41467-021-23300-y] [Citation(s) in RCA: 20] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/16/2020] [Accepted: 04/20/2021] [Indexed: 01/04/2023] Open
Abstract
High-yielding and selective prebiotic syntheses of RNA and DNA nucleotides involve UV irradiation to promote the key reaction steps and eradicate biologically irrelevant isomers. While these syntheses were likely enabled by UV-rich prebiotic environment, UV-induced formation of photodamages in polymeric nucleic acids, such as cyclobutane pyrimidine dimers (CPDs), remains the key unresolved issue for the origins of RNA and DNA on Earth. Here, we demonstrate that substitution of adenine with 2,6-diaminopurine enables repair of CPDs with yields reaching 92%. This substantial self-repairing activity originates from excellent electron donating properties of 2,6-diaminopurine in nucleic acid strands. We also show that the deoxyribonucleosides of 2,6-diaminopurine and adenine can be formed under the same prebiotic conditions. Considering that 2,6-diaminopurine was previously shown to increase the rate of nonenzymatic RNA replication, this nucleobase could have played critical roles in the formation of functional and photostable RNA/DNA oligomers in UV-rich prebiotic environments. UV-induced photodamage that likely occurred during the prebiotic synthesis of DNA and RNA is still an untackled issue for their origin on early Earth. Here, the authors show that substitution of 2,6-diaminopurine for adenine enables repair of cyclobutane pyrimidine dimers with high yields, and demonstrate that both 2,6-diaminopurine and adenine nucleosides can be formed under the same prebiotic conditions.
Collapse
Affiliation(s)
- Rafał Szabla
- EaStCHEM School of Chemistry, University of Edinburgh, Edinburgh, UK. .,Institute of Physics, Polish Academy of Sciences, Warsaw, Poland.
| | | | - Paulina Spisz
- Faculty of Chemistry, University of Gdańsk, Gdańsk, Poland
| | | | - Petr Stadlbauer
- Institute of Biophysics of the Czech Academy of Sciences, Brno, Czech Republic
| | - Holger Kruse
- Institute of Biophysics of the Czech Academy of Sciences, Brno, Czech Republic
| | - Jiří Šponer
- Institute of Biophysics of the Czech Academy of Sciences, Brno, Czech Republic
| | - Janusz Rak
- Faculty of Chemistry, University of Gdańsk, Gdańsk, Poland
| |
Collapse
|
28
|
Lee S, Horbatenko Y, Filatov M, Choi CH. Fast and Accurate Computation of Nonadiabatic Coupling Matrix Elements Using the Truncated Leibniz Formula and Mixed-Reference Spin-Flip Time-Dependent Density Functional Theory. J Phys Chem Lett 2021; 12:4722-4728. [PMID: 33983029 DOI: 10.1021/acs.jpclett.1c00932] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
We present a fast and accurate numerical algorithm for computing the first-order nonadiabatic coupling matrix element (NACME). The algorithm employs the truncated Leibniz formula (TLF) approximation within the finite-difference method, which makes it easily applicable in connection with any wave function-based methodology. In this work, we used the algorithm in connection with the recently developed mixed-reference spin-flip time-dependent density functional theory (MRSF-TDDFT, MRSF for brevity). The accuracy is assessed for NACME between the singlet electronic states of a dissociating hydrogen molecule. It is demonstrated that an intermediate approximation, TLF(1), affords a negligible numeric error on the order of ∼10-10 a.u. while enabling a fast computation of NACME. As the MRSF method yields the correct description of the dissociation curves of H2 for all the electronic states involved, the numeric TLF(1)/MRSF NACME values are in excellent agreement with the reference analytical values obtained by the full configuration interaction. For polyatomic molecules, the MRSF NAC vectors agree very closely with the MRCISD NAC vectors. Hence, the proposed protocol is a promising tool for the evaluation of NACMEs.
Collapse
Affiliation(s)
- Seunghoon Lee
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Yevhen Horbatenko
- Department of Chemistry, Kyungpook National University, Daegu 702-701, South Korea
| | - Michael Filatov
- Department of Chemistry, Kyungpook National University, Daegu 702-701, South Korea
| | - Cheol Ho Choi
- Department of Chemistry, Kyungpook National University, Daegu 702-701, South Korea
| |
Collapse
|
29
|
Park W, Lee S, Huix-Rotllant M, Filatov M, Choi CH. Impact of the Dynamic Electron Correlation on the Unusually Long Excited-State Lifetime of Thymine. J Phys Chem Lett 2021; 12:4339-4346. [PMID: 33929858 DOI: 10.1021/acs.jpclett.1c00712] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
Non-radiative relaxation of the photoexcited thymine in the gas phase shows an unusually long excited-state lifetime, and, over the years, a number of models, i.e., S1-trapping, S2-trapping, and S1&S2-trapping, have been put forward to explain its mechanism. Here, we investigate this mechanism using non-adiabatic molecular dynamics (NAMD) simulations in connection with the recently developed mixed-reference spin-flip time-dependent density functional theory (MRSF-TDDFT) method. We show that the previously predicted S2-trapping model was due to an artifact caused by an insufficient account of the dynamic electron correlation. The current work supports the S1-trapping mechanism with two lifetimes, τ1 = 30 ± 1 fs and τ2 = 6.1 ± 0.035 ps, quantitatively consistent with the recent time-resolved experiments. Upon excitation to the S2 (ππ*) state, thymine undergoes an ultrafast (ca. 30 fs) S2→S1 internal conversion and resides around the minimum on the S1 (nOπ*) surface, slowly decaying to the ground state (ca. 6.1 ps). While the S2→S1 internal conversion is mediated by fast bond length alternation distortion, the subsequent S1→S0 occurs through several conical intersections, involving a slow puckering motion.
Collapse
Affiliation(s)
- Woojin Park
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| | - Seunghoon Lee
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | | | - Michael Filatov
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| | - Cheol Ho Choi
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| |
Collapse
|
30
|
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
|
31
|
Rivera M, Stojanović L, Crespo-Otero R. Role of Conical Intersections on the Efficiency of Fluorescent Organic Molecular Crystals. J Phys Chem A 2021; 125:1012-1024. [PMID: 33492964 DOI: 10.1021/acs.jpca.0c11072] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Abstract
Organic molecular crystals are attractive materials for luminescent applications because of their promised tunability. However, the link between the chemical structure and emissive behavior is poorly understood because of the numerous interconnected factors which are at play in determining radiative and nonradiative behaviors at the solid-state level. In particular, the decay through conical intersection dominates the nonadiabatic regions of the potential energy surface, and thus, their accessibility is a telling indicator of the luminosity of the material. In this study, we investigate the radiative mechanism for five organic molecular crystals which display a solid-state emission, with a focus on the role of conical intersections in their photomechanisms. The objective is to situate the importance of the accessibility of conical intersections with regards to emissive behavior, taking into account other nonradiative decay channels, namely, vibrational decay, and exciton hopping. We begin by giving a brief overview of the structural patterns of the five systems within a larger pool of 13 crystals for a richer comparison. We observe that because of the prevalence of sheet like and herringbone packing in organic molecular crystals, the conformational diversity of crystal dimers is limited. Additionally, similarly spaced dimers have exciton coupling values of a similar order within a 50 meV interval. Next, we focus on three exemplary cases, where we disentangle the role of nonradiative decay mechanisms and show how rotational minimum energy conical intersections in vacuum lead to puckered ones in the crystal, increasing their instability upon crystallization in typical packing motifs. In contrast, molecules with puckered conical intersections in vacuum tend to conserve this trait upon crystallization, and therefore, their quantum yield of fluorescence is determined predominantly by other nonradiative decay mechanisms.
Collapse
Affiliation(s)
- Miguel Rivera
- School of Biological and Chemical Sciences, Queen Mary University of London, Mile End Road, London E1 4NS, U.K
| | - Ljiljana Stojanović
- School of Biological and Chemical Sciences, Queen Mary University of London, Mile End Road, London E1 4NS, U.K
| | - Rachel Crespo-Otero
- School of Biological and Chemical Sciences, Queen Mary University of London, Mile End Road, London E1 4NS, U.K
| |
Collapse
|
32
|
Ha JK, Kim K, Min SK. Machine Learning-Assisted Excited State Molecular Dynamics with the State-Interaction State-Averaged Spin-Restricted Ensemble-Referenced Kohn-Sham Approach. J Chem Theory Comput 2021; 17:694-702. [PMID: 33470100 DOI: 10.1021/acs.jctc.0c01261] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
We present a machine learning-assisted excited state molecular dynamics (ML-ESMD) based on the ensemble density functional theory framework. Since we represent a diabatic Hamiltonian in terms of generalized valence bond ansatz within the state-interaction state-averaged spin-restricted ensemble-referenced Kohn-Sham (SI-SA-REKS) method, we can avoid singularities near conical intersections, which are crucial in excited state molecular dynamics simulations. We train the diabatic Hamiltonian elements and their analytical gradients with the SchNet architecture to construct machine learning models, while the phase freedom of off-diagonal elements of the Hamiltonian is cured by introducing the phase-less loss function. Our machine learning models show reasonable accuracy with mean absolute errors of ∼0.1 kcal/mol and ∼0.5 kcal/mol/Å for the diabatic Hamiltonian elements and their gradients, respectively, for penta-2,4-dieniminium cation. Moreover, by exploiting the diabatic representation, our models can predict correct conical intersection structures and their topologies. In addition, our ML-ESMD simulations give almost identical result with a direct dynamics at the same level of theory.
Collapse
Affiliation(s)
- Jong-Kwon Ha
- Department of Chemistry, School of Natural Science, Ulsan National Institute of Science and Technology (UNIST), 50 UNIST-gil, Ulju-gun, Ulsan 44919, South Korea
| | - Kicheol Kim
- Department of Chemistry, School of Natural Science, Ulsan National Institute of Science and Technology (UNIST), 50 UNIST-gil, Ulju-gun, Ulsan 44919, South Korea
| | - Seung Kyu Min
- Department of Chemistry, School of Natural Science, Ulsan National Institute of Science and Technology (UNIST), 50 UNIST-gil, Ulju-gun, Ulsan 44919, South Korea
| |
Collapse
|
33
|
Li J, Liu MK, Li QS, Li ZS. Theoretical study on the photochemistry of furoylazides: Curtius rearrangement and subsequent reactions. Phys Chem Chem Phys 2020; 22:28317-28324. [PMID: 33300534 DOI: 10.1039/d0cp05539e] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023]
Abstract
Organic azides are an efficient source of nitrenes, which serve as vigorous intermediates in many useful organic reactions. In this work, the complete active space self-consistent field (CASSCF) and its second-order perturbation (CASPT2) methods were employed to study the photochemistry of 2-furoylazide 1 and 3-furoylazide 5, including the Curtius rearrangement to two furylisocyanates (3 and 7) and subsequent reactions to the final product cyanoacrolein 9. Our calculations show that the photoinduced Curtius rearrangement of the two furoylazides takes place through similar stepwise mechanisms via two bistable furoylnitrenes 2 and 6. However, the decarbonylation and ring-opening process of 7 to 9 prefers a stepwise mechanism involving the 3-furoylnitrene intermediate 8, while 3 to 9 goes in a concerted asynchronous way without the corresponding 2-furoylnitrene intermediate 4. Importantly, we revealed that several conical intersections play key roles in the photochemistry of furoylazides. Our results are not only consistent and also make clear the experimental observations (X. Zeng, et al., J. Am. Chem. Soc., 2018, 140, 10-13), but additionally provide important information on the chemistry of furoylazides and nitrenes.
Collapse
Affiliation(s)
- Jian Li
- Beijing Key Laboratory of Photoelectronic/Electrophotonic Conversion Materials, Key Laboratory of Cluster Science of Ministry of Education, School of Chemistry and Chemical Engineering, Beijing Institute of Technology, Beijing 100081, China.
| | | | | | | |
Collapse
|
34
|
Chrayteh A, Blondel A, Loos PF, Jacquemin D. Mountaineering Strategy to Excited States: Highly Accurate Oscillator Strengths and Dipole Moments of Small Molecules. J Chem Theory Comput 2020; 17:416-438. [DOI: 10.1021/acs.jctc.0c01111] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Affiliation(s)
- Amara Chrayteh
- Université de Nantes, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
| | - Aymeric Blondel
- Université de Nantes, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
| | - Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques, Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| | - Denis Jacquemin
- Université de Nantes, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
| |
Collapse
|
35
|
Lee W, Matsika S. Stabilization of the Triplet Biradical Intermediate of 5-Methylcytosine Enhances Cyclobutane Pyrimidine Dimer (CPD) Formation in DNA. Chemistry 2020; 26:14181-14186. [PMID: 32809239 DOI: 10.1002/chem.202002834] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/12/2020] [Revised: 08/01/2020] [Indexed: 11/08/2022]
Abstract
Cyclobutane pyrimidine dimer (CPD) is a photoproduct formed by two stacked pyrimidine bases through a cycloaddition reaction upon irradiation. Owing to its close association with skin cancer, the mechanism of CPD formation has been studied thoroughly. Among many aspects of CPD, its formation involving 5-methylcytosine (5mC) has been of special interest because the CPD yield is known to increase with C5-methylation of cytosine. In this work, high-level quantum mechanics/molecular mechanics (QM/MM) calculations are used to examine a previously experimentally detected pathway for CPD formation in hetero (thymine-cytosine and thymine-5mC) dipyrimidines, which is facilitated through intersystem crossing in thymine and formation of a triplet biradical intermediate. A DNA duplex model system containing a core sequence TmCG or TCG is used. The stabilization of a radical center in the biradical intermediate by the methyl group of 5mC can lead to increased CPD yield in TmCG compared with its non-methylated counterpart, TCG, thereby suggesting the existence of a new pathway of CPD formation enhanced by 5mC.
Collapse
Affiliation(s)
- Wook Lee
- Department of Biochemistry, Kangwon National University, Chuncheon, 24341, South Korea
| | - Spiridoula Matsika
- Department of Chemistry, Temple University, Philadelphia, PA, 19122, USA
| |
Collapse
|
36
|
Marsili E, Olivucci M, Lauvergnat D, Agostini F. Quantum and Quantum-Classical Studies of the Photoisomerization of a Retinal Chromophore Model. J Chem Theory Comput 2020; 16:6032-6048. [PMID: 32931266 DOI: 10.1021/acs.jctc.0c00679] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
We report an in-depth analysis of the photo-induced isomerization of the 2-cis-penta-2,4-dieniminium cation: a minimal model of the 11-cis retinal protonated Schiff base chromophore of the dim-light photoreceptor rhodopsin. Based on recently developed three-dimensional potentials parametrized on ab initio multi-state multi-configurational second-order perturbation theory data, we perform quantum-dynamical studies. In addition, simulations based on various quantum-classical methods, among which Tully surface hopping and the coupled-trajectory approach derived from the exact factorization, allow us to validate their performance against vibronic wavepacket propagation and, therefore, a purely quantum treatment. Quantum-dynamics results uncover qualitative differences with respect to the two-dimensional Hahn-Stock potentials, widely used as model potentials for the isomerization of the same chromophore, due to the increased dimensionality and three-mode correlation. Quantum-classical simulations show, instead, that three-dimensional model potentials are capable of capturing a number of features revealed by atomistic simulations and experimental observations. In particular, a recently reported vibrational phase relationship between double-bond torsion and hydrogen-out-of-plane modes critical for rhodopsin isomerization efficiency is correctly reproduced.
Collapse
Affiliation(s)
- Emanuele Marsili
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, Orsay 91405, France.,Department of Chemistry, Durham University, South Road, Durham DH1 3LE, United Kingdom
| | - Massimo Olivucci
- Department of Biotechnology, Chemistry and Pharmacy, Università degli Studi di Siena, Via A. Moro 2, I-53100 Siena, Italy.,Department of Chemistry, Bowling Green State University, Bowling Green, Ohio 43403, United States
| | - David Lauvergnat
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, Orsay 91405, France
| | - Federica Agostini
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, Orsay 91405, France
| |
Collapse
|
37
|
Mai S, González L. Molecular Photochemistry: Recent Developments in Theory. Angew Chem Int Ed Engl 2020; 59:16832-16846. [PMID: 32052547 PMCID: PMC7540682 DOI: 10.1002/anie.201916381] [Citation(s) in RCA: 77] [Impact Index Per Article: 19.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2019] [Revised: 02/12/2020] [Indexed: 12/16/2022]
Abstract
Photochemistry is a fascinating branch of chemistry that is concerned with molecules and light. However, the importance of simulating light-induced processes is reflected also in fields as diverse as biology, material science, and medicine. This Minireview highlights recent progress achieved in theoretical chemistry to calculate electronically excited states of molecules and simulate their photoinduced dynamics, with the aim of reaching experimental accuracy. We focus on emergent methods and give selected examples that illustrate the progress in recent years towards predicting complex electronic structures with strong correlation, calculations on large molecules, describing multichromophoric systems, and simulating non-adiabatic molecular dynamics over long time scales, for molecules in the gas phase or in complex biological environments.
Collapse
Affiliation(s)
- Sebastian Mai
- Photonics InstituteVienna University of TechnologyGusshausstrasse 27–291040ViennaAustria
| | - Leticia González
- Institute of Theoretical ChemistryFaculty of ChemistryUniversity of ViennaWähringer Strasse 171090ViennaAustria
| |
Collapse
|
38
|
Hegger R, Binder R, Burghardt I. First-Principles Quantum and Quantum-Classical Simulations of Exciton Diffusion in Semiconducting Polymer Chains at Finite Temperature. J Chem Theory Comput 2020; 16:5441-5455. [PMID: 32786907 DOI: 10.1021/acs.jctc.0c00351] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
We report on first-principles quantum-dynamical and quantum-classical simulations of photoinduced exciton dynamics in oligothiophene chain segments, representative of intrachain exciton migration in the poly(3-hexylthiophene) (P3HT) polymer. Following up on our recent study (Binder R.; Burghardt, I. Faraday Discuss. 2020, 221, 406), multilayer multiconfiguration time-dependent Hartree calculations for a short oligothiophene segment comprising 20 monomer units (OT-20) are carried out to obtain full quantum-dynamical simulations at finite temperature. These are employed to benchmark mean-field Ehrenfest calculations, which are shown to give qualitatively correct results for the present system. Periodic boundary conditions turn out to significantly improve earlier estimates of diffusion coefficients. Using the Ehrenfest approach, a series of calculations are subsequently carried out for larger lattices (OT-40 to OT-80), leading to estimates for temperature-dependent mean-squared displacements, which are found to exhibit a near-linear dependence as a function of time. The resulting diffusion coefficient estimates are an increasing function of temperature, whose detailed functional form depends on the degree of static disorder. With a realistic static disorder parameter (σs ≃ 0.06 eV), the diffusion coefficients decrease from D ∼ 1 × 10-2 cm2 s-1 to D ∼ 1 × 10-3 cm2 s-1, in qualitative agreement with experimental data for P3HT. The dynamical scenario obtained from our simulations shows that exciton migration in P3HT-type chains is a largely adiabatic process throughout the temperature regime we investigated (i.e., T = 50-300 K). The resulting picture of exciton migration is a coherent, but not bandlike, motion of an exciton-polaron driven by fluctuations induced by low-frequency modes. This process acquires partial hopping character if static disorder becomes prominent and Anderson localization sets in.
Collapse
Affiliation(s)
- Rainer Hegger
- Institute for Physical and Theoretical Chemistry, Goethe University, Max-von-Laue-Strasse 7, 60438 Frankfurt/Main, Germany
| | - Robert Binder
- Institute for Physical and Theoretical Chemistry, Goethe University, Max-von-Laue-Strasse 7, 60438 Frankfurt/Main, Germany
| | - Irene Burghardt
- Institute for Physical and Theoretical Chemistry, Goethe University, Max-von-Laue-Strasse 7, 60438 Frankfurt/Main, Germany
| |
Collapse
|
39
|
Song C, Martínez TJ. Reduced scaling extended multi-state CASPT2 (XMS-CASPT2) using supporting subspaces and tensor hyper-contraction. J Chem Phys 2020; 152:234113. [PMID: 32571032 DOI: 10.1063/5.0007417] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023] Open
Abstract
We present a reduced scaling formulation of the extended multi-state CASPT2 (XMS-CASPT2) method, which is based on our recently developed state-specific CASPT2 (SS-CASPT2) formulation using supporting subspaces and tensor hyper-contraction. By using these two techniques, the off-diagonal elements of the effective Hamiltonian can be computed with only O(N3) operations and O(N2) memory, where N is the number of basis functions. This limits the overall computational scaling to O(N4) operations and O(N2) memory. Thus, excited states can now be obtained at the same reduced (relative to previous algorithms) scaling we achieved for SS-CASPT2. In addition, we also investigate how the energy denominators can be factorized with the Laplace quadrature when some of the denominators are negative, which is critical for excited state calculations. An efficient implementation of the method has been developed using graphical processing units while also exploiting spatial sparsity in tensor operations. We benchmark the accuracy of the new method by comparison to non-THC formulated XMS-CASPT2 for the excited states of various molecules. In our tests, the THC approximation introduces negligible errors (≈0.01 eV) compared to the non-THC reference method. Scaling behavior and computational timings are presented to demonstrate performance. The new method is also interfaced with quantum mechanics/molecular mechanics (QM/MM). In an example study of green fluorescent protein, we show how the XMS-CASPT2 potential energy surfaces and excitation energies are affected by increasing the size of the QM region up to 278 QM atoms with more than 2300 basis functions.
Collapse
Affiliation(s)
- Chenchen Song
- Department of Chemistry and the PULSE Institute, Stanford University, Stanford, California 94305, USA
| | - Todd J Martínez
- Department of Chemistry and the PULSE Institute, Stanford University, Stanford, California 94305, USA
| |
Collapse
|
40
|
Mai S, González L. Molekulare Photochemie: Moderne Entwicklungen in der theoretischen Chemie. Angew Chem Int Ed Engl 2020. [DOI: 10.1002/ange.201916381] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Affiliation(s)
- Sebastian Mai
- Institut für Photonik Technische Universität Wien Gußhausstraße 27–29 1040 Wien Österreich
| | - Leticia González
- Institut für theoretische Chemie Fakultät für Chemie Universität Wien Währinger Straße 17 1090 Wien Österreich
| |
Collapse
|
41
|
Cuzzocrea A, Scemama A, Briels WJ, Moroni S, Filippi C. Variational Principles in Quantum Monte Carlo: The Troubled Story of Variance Minimization. J Chem Theory Comput 2020; 16:4203-4212. [PMID: 32419451 PMCID: PMC7365558 DOI: 10.1021/acs.jctc.0c00147] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Abstract
![]()
We
investigate the use of different variational principles in quantum
Monte Carlo, namely, energy and variance minimization, prompted by
the interest in the robust and accurate estimation of electronic excited
states. For two prototypical, challenging molecules, we readily reach
the accuracy of the best available reference excitation energies using
energy minimization in a state-specific or state-average fashion for
states of different or equal symmetry, respectively. On the other
hand, in variance minimization, where the use of suitable functionals
is expected to target specific states regardless of the symmetry,
we encounter severe problems for a variety of wave functions: as the
variance converges, the energy drifts away from that of the selected
state. This unexpected behavior is sometimes observed even when the
target is the ground state and generally prevents the robust estimation
of total and excitation energies. We analyze this problem using a
very simple wave function and infer that the optimization finds little
or no barrier to escape from a local minimum or local plateau, eventually
converging to a lower-variance state instead of the target state.
For the increasingly complex systems becoming in reach of quantum
Monte Carlo simulations, variance minimization with current functionals
appears to be an impractical route.
Collapse
Affiliation(s)
- Alice Cuzzocrea
- MESA+ Institute for Nanotechnology, University of Twente, P.O. Box 217, 7500 AE Enschede, The Netherlands
| | - Anthony Scemama
- Laboratoire de Chimie et Physique Quantiques, Université de Toulouse, CNRS, UPS, 118 route de Narbonne, 31062 Toulouse Cedex 09, France
| | - Wim J Briels
- MESA+ Institute for Nanotechnology, University of Twente, P.O. Box 217, 7500 AE Enschede, The Netherlands
| | - Saverio Moroni
- CNR-IOM DEMOCRITOS, Istituto Officina dei Materiali, Via Bonomea 265, I-34136 Trieste, Italy.,SISSA Scuola Internazionale Superiore di Studi Avanzati, Via Bonomea 265, I-34136 Trieste, Italy
| | - Claudia Filippi
- MESA+ Institute for Nanotechnology, University of Twente, P.O. Box 217, 7500 AE Enschede, The Netherlands
| |
Collapse
|
42
|
Siddique F, Barbatti M, Cui Z, Lischka H, Aquino AJA. Nonadiabatic Dynamics of Charge-Transfer States Using the Anthracene–Tetracyanoethylene Complex as a Prototype. J Phys Chem A 2020; 124:3347-3357. [DOI: 10.1021/acs.jpca.0c01900] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- Farhan Siddique
- School of Pharmaceutical Science and Technology, Tianjin University, Tianjin 300072, China
| | | | - Zhonghua Cui
- Institute of Atomic and Molecular Physics, Jilin University, Changchun 130400, China
- Beijing National Laboratory for Molecular Sciences, Beijing 100080, China
| | - Hans Lischka
- School of Pharmaceutical Science and Technology, Tianjin University, Tianjin 300072, China
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409-1061, United States
| | - Adelia J. A. Aquino
- School of Pharmaceutical Science and Technology, Tianjin University, Tianjin 300072, China
- Department of Mechanical Engineering, Texas Tech University, Lubbock, Texas 79409, United States
| |
Collapse
|
43
|
List NH, Dempwolff AL, Dreuw A, Norman P, Martínez TJ. Probing competing relaxation pathways in malonaldehyde with transient X-ray absorption spectroscopy. Chem Sci 2020; 11:4180-4193. [PMID: 34122881 PMCID: PMC8152795 DOI: 10.1039/d0sc00840k] [Citation(s) in RCA: 22] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022] Open
Abstract
Excited-state intramolecular hydrogen transfer (ESIHT) is a fundamental reaction relevant to chemistry and biology. Malonaldehyde is the simplest example of ESIHT, yet only little is known experimentally about its excited-state dynamics. Several competing relaxation pathways have been proposed, including internal conversion mediated by ESIHT and C[double bond, length as m-dash]C torsional motion as well as intersystem crossing. We perform an in silico transient X-ray absorption spectroscopy (TRXAS) experiment at the oxygen K-edge to investigate its potential to monitor the proposed ultrafast decay pathways in malonaldehyde upon photoexcitation to its bright S2(ππ*) state. We employ both restricted active space perturbation theory and algebraic-diagrammatic construction for the polarization propagator along interpolated reaction coordinates as well as representative trajectories from ab initio multiple spawning simulations to compute the TRXAS signals from the lowest valence states. Our study suggests that oxygen K-edge TRXAS can distinctly fingerprint the passage through the H-transfer intersection and the concomitant population transfer to the S1(nπ*) state. Potential intersystem crossing to T1(ππ*) is detectable from reappearance of the double pre-edge signature and reversed intensities. Moreover, the torsional deactivation pathway induces transient charge redistribution from the enol side towards the central C-atom and manifests itself as substantial shifts of the pre-edge features. Given the continuous advances in X-ray light sources, our study proposes an experimental route to disentangle ultrafast excited-state decay channels in this prototypical ESIHT system and provides a pathway-specific mapping of the TRXAS signal to facilitate the interpretation of future experiments.
Collapse
Affiliation(s)
- Nanna H List
- Department of Chemistry, The PULSE Institute, Stanford University Stanford CA 94305 USA .,SLAC National Accelerator Laboratory 2575 Sand Hill Road Menlo Park CA 94025 USA
| | - Adrian L Dempwolff
- Interdisciplinary Center for Scientific Computing, Heidelberg University Im Neuenheimer Feld 205 D-69120 Heidelberg Germany
| | - Andreas Dreuw
- Interdisciplinary Center for Scientific Computing, Heidelberg University Im Neuenheimer Feld 205 D-69120 Heidelberg Germany
| | - Patrick Norman
- School of Engineering Sciences in Chemistry, Biotechnology and Health (CBH), Department of Theoretical Chemistry and Biology, KTH Royal Institute of Technology Sweden
| | - Todd J Martínez
- Department of Chemistry, The PULSE Institute, Stanford University Stanford CA 94305 USA .,SLAC National Accelerator Laboratory 2575 Sand Hill Road Menlo Park CA 94025 USA
| |
Collapse
|
44
|
Rivera M, Dommett M, Sidat A, Rahim W, Crespo-Otero R. fromage: A library for the study of molecular crystal excited states at the aggregate scale. J Comput Chem 2020; 41:1045-1058. [PMID: 31909830 PMCID: PMC7079081 DOI: 10.1002/jcc.26144] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/12/2019] [Revised: 11/27/2019] [Accepted: 12/03/2019] [Indexed: 12/31/2022]
Abstract
The study of photoexcitations in molecular aggregates faces the twofold problem of the increased computational cost associated with excited states and the complexity of the interactions among the constituent monomers. A mechanistic investigation of these processes requires the analysis of the intermolecular interactions, the effect of the environment, and 3D arrangements or crystal packing on the excited states. A considerable number of techniques have been tailored to navigate these obstacles; however, they are usually restricted to in‐house codes and thus require a disproportionate effort to adopt by researchers approaching the field. Herein, we present the FRamewOrk for Molecular AGgregate Excitations (fromage), which implements a collection of such techniques in a Python library complemented with ready‐to‐use scripts. The program structure is presented and the principal features available to the user are described: geometrical analysis, exciton characterization, and a variety of ONIOM schemes. Each is illustrated by examples of diverse organic molecules in condensed phase settings. The program is available at https://github.com/Crespo-Otero-group/fromage.
Collapse
Affiliation(s)
- Miguel Rivera
- Department of Chemistry, School of Biological and Chemical Sciences, Queen Mary University of London, London, UK
| | - Michael Dommett
- Department of Chemistry, School of Biological and Chemical Sciences, Queen Mary University of London, London, UK
| | - Amir Sidat
- Department of Chemistry, School of Biological and Chemical Sciences, Queen Mary University of London, London, UK
| | - Warda Rahim
- Department of Chemistry, School of Biological and Chemical Sciences, Queen Mary University of London, London, UK
| | - Rachel Crespo-Otero
- Department of Chemistry, School of Biological and Chemical Sciences, Queen Mary University of London, London, UK
| |
Collapse
|
45
|
Gromov EV, Domratcheva T. Four resonance structures elucidate double-bond isomerisation of a biological chromophore. Phys Chem Chem Phys 2020; 22:8535-8544. [DOI: 10.1039/d0cp00814a] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Four resonance structures determining the electronic structure of the chromophore’s ground and first excited states. Changing the relative energies of the structures by hydrogen-bonding interactions tunes all chromophore’s photochemical properties.
Collapse
Affiliation(s)
- Evgeniy V. Gromov
- Max-Planck Institute for Medical Research
- Jahnstraße 29
- 69120 Heidelberg
- Germany
| | - Tatiana Domratcheva
- Max-Planck Institute for Medical Research
- Jahnstraße 29
- 69120 Heidelberg
- Germany
| |
Collapse
|
46
|
Dupuy MS, Gloaguen E, Tardivel B, Mons M, Brenner V. CC2 Benchmark for Models of Phenylalanine Protein Chains: 0–0 Transition Energies and IR Signatures of the ππ* Excited State. J Chem Theory Comput 2019; 16:601-611. [DOI: 10.1021/acs.jctc.9b00923] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Mi-Song Dupuy
- LIDYL, CEA, CNRS, Université Paris-Saclay, 91191 Gif-sur-Yvette, France
| | - Eric Gloaguen
- LIDYL, CEA, CNRS, Université Paris-Saclay, 91191 Gif-sur-Yvette, France
| | - Benjamin Tardivel
- LIDYL, CEA, CNRS, Université Paris-Saclay, 91191 Gif-sur-Yvette, France
| | - Michel Mons
- LIDYL, CEA, CNRS, Université Paris-Saclay, 91191 Gif-sur-Yvette, France
| | - Valérie Brenner
- LIDYL, CEA, CNRS, Université Paris-Saclay, 91191 Gif-sur-Yvette, France
| |
Collapse
|
47
|
Abedini F, Omidyan R, Salehi M. Theoretical insights on nonradiative deactivation mechanisms of protonated xanthine. J Photochem Photobiol A Chem 2019. [DOI: 10.1016/j.jphotochem.2019.112067] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023]
|
48
|
Kochman MA, Durbeej B. Theoretical Study of Ground- and Excited-State Charge Transfer in Fulvene-Based Donor–Acceptor Systems. J Phys Chem A 2019; 123:6660-6673. [DOI: 10.1021/acs.jpca.9b02962] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Michał Andrzej Kochman
- Division of Theoretical Chemistry, Department of Physics, Chemistry and Biology (IFM), Linköping University, 581 83 Linköping, Sweden
| | - Bo Durbeej
- Division of Theoretical Chemistry, Department of Physics, Chemistry and Biology (IFM), Linköping University, 581 83 Linköping, Sweden
| |
Collapse
|
49
|
Lee S, Shostak S, Filatov M, Choi CH. Conical Intersections in Organic Molecules: Benchmarking Mixed-Reference Spin-Flip Time-Dependent DFT (MRSF-TD-DFT) vs Spin-Flip TD-DFT. J Phys Chem A 2019; 123:6455-6462. [PMID: 31283235 DOI: 10.1021/acs.jpca.9b06142] [Citation(s) in RCA: 22] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The mixed-reference spin-flip time-dependent density functional theory (MRSF-TD-DFT) method eliminates the erroneous spin contamination of the SF-TD-DFT methodology, while retaining the conceptual and practical simplicity of the latter. The availability of the analytic gradient of the energy of the MRSF-TD-DFT response states enables automatic geometry optimization of the targeted states. Here, we apply the new method to optimize the geometry of several S1/S0 conical intersections occurring in typical organic molecules. We demonstrate that MRSF-TD-DFT is capable of producing the correct double-cone topology of the intersections and describing the geometry of the lowest-energy conical intersections and their relative energies with accuracy matching that of the best multireference wavefunction ab initio methods. In this regard, MRSF-TD-DFT differs from many popular single-reference methods, such as, e.g., the linear response TD-DFT method, which fail to produce the correct topology of the intersections. As the new methodology completely eliminates the ambiguity with the identification of the response states as proper singlets or triplets, which is plaguing the SF-TD-DFT calculations, it can be used for automatic geometry optimization and molecular dynamic simulations not requiring constant human intervention.
Collapse
Affiliation(s)
- Seunghoon Lee
- Department of Chemistry , Seoul National University , Seoul 151-747 , South Korea
| | - Svetlana Shostak
- Department of Chemistry , Kyungpook National University , Daegu 702-701 , South Korea
| | - Michael Filatov
- Department of Chemistry , Kyungpook National University , Daegu 702-701 , South Korea
| | - Cheol Ho Choi
- Department of Chemistry , Kyungpook National University , Daegu 702-701 , South Korea
| |
Collapse
|
50
|
Mai S, Atkins AJ, Plasser F, González L. The Influence of the Electronic Structure Method on Intersystem Crossing Dynamics. The Case of Thioformaldehyde. J Chem Theory Comput 2019; 15:3470-3480. [DOI: 10.1021/acs.jctc.9b00282] [Citation(s) in RCA: 22] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Affiliation(s)
- Sebastian Mai
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straß 17, 1090 Vienna, Austria
| | - Andrew J. Atkins
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straß 17, 1090 Vienna, Austria
| | - Felix Plasser
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straß 17, 1090 Vienna, Austria
| | - Leticia González
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straß 17, 1090 Vienna, Austria
| |
Collapse
|