1
|
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
|
2
|
Janoš J, Slavíček P. What Controls the Quality of Photodynamical Simulations? Electronic Structure Versus Nonadiabatic Algorithm. J Chem Theory Comput 2023; 19:8273-8284. [PMID: 37939301 PMCID: PMC10688183 DOI: 10.1021/acs.jctc.3c00908] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/18/2023] [Revised: 10/18/2023] [Accepted: 10/19/2023] [Indexed: 11/10/2023]
Abstract
The field of nonadiabatic dynamics has matured over the last decade with a range of algorithms and electronic structure methods available at the moment. While the community currently focuses more on developing and benchmarking new nonadiabatic dynamics algorithms, the underlying electronic structure controls the outcome of nonadiabatic simulations. Yet, the electronic-structure sensitivity analysis is typically neglected. In this work, we present a sensitivity analysis of the nonadiabatic dynamics of cyclopropanone to electronic structure methods and nonadiabatic dynamics algorithms. In particular, we compare wave function-based CASSCF, FOMO-CASCI, MS- and XMS-CASPT2, density-functional REKS, and semiempirical MRCI-OM3 electronic structure methods with the Landau-Zener surface hopping, fewest switches surface hopping, and ab initio multiple spawning with informed stochastic selection algorithms. The results clearly demonstrate that the electronic structure choice significantly influences the accuracy of nonadiabatic dynamics for cyclopropanone even when the potential energy surfaces exhibit qualitative and quantitative similarities. Thus, selecting the electronic structure solely on the basis of the mapping of potential energy surfaces can be misleading. Conversely, we observe no discernible differences in the performance of the nonadiabatic dynamics algorithms across the various methods. Based on the above results, we discuss the present-day practice in computational photodynamics.
Collapse
Affiliation(s)
- Jiří Janoš
- Department of Physical Chemistry, University of Chemistry and Technology, Technická 5, 16628 Prague 6, Czech Republic
| | - Petr Slavíček
- Department of Physical Chemistry, University of Chemistry and Technology, Technická 5, 16628 Prague 6, Czech Republic
| |
Collapse
|
3
|
Sen S, Sisodiya DS, Chattopadhyay A. The mechanism of photoconversion of cyclic dinitrone to oxaziridine and dioxaziridine: A computational investigation of an experimentally reported photochemical reaction. J PHYS ORG CHEM 2021. [DOI: 10.1002/poc.4310] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Affiliation(s)
- Sindhuja Sen
- Department of Chemistry Birla Institute of Technology and Science (BITS), Pilani – K.K. Birla Goa Campus Zuarinagar Goa India
| | - Dilawar Singh Sisodiya
- Department of Chemistry Birla Institute of Technology and Science (BITS), Pilani – K.K. Birla Goa Campus Zuarinagar Goa India
| | - Anjan Chattopadhyay
- Department of Chemistry Birla Institute of Technology and Science (BITS), Pilani – K.K. Birla Goa Campus Zuarinagar Goa India
| |
Collapse
|
4
|
Segarra‐Martí J, Bearpark MJ. Modelling Photoionisation in Isocytosine: Potential Formation of Longer-Lived Excited State Cations in its Keto Form. Chemphyschem 2021; 22:2172-2181. [PMID: 34370368 PMCID: PMC8597144 DOI: 10.1002/cphc.202100402] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/25/2021] [Revised: 08/09/2021] [Indexed: 11/24/2022]
Abstract
Studying the effects of UV and VUV radiation on non-canonical DNA/RNA nucleobases allows us to compare how they release excess energy following absorption with respect to their canonical counterparts. This has attracted much research attention in recent years because of its likely influence on the origin of our genetic lexicon in prebiotic times. Here we present a CASSCF and XMS-CASPT2 theoretical study of the photoionisation of non-canonical pyrimidine nucleobase isocytosine in both its keto and enol tautomeric forms. We analyse their lowest energy cationic excited states including 2 π + , 2 n O + and 2 n N + and compare these to the corresponding electronic states in cytosine. Investigating lower-energy decay pathways we find - unexpectedly - that keto-isocytosine+ presents a sizeable energy barrier potentially inhibiting decay to its cationic ground state, whereas enol-isocytosine+ features a barrierless and consequently ultrafast pathway analogous to the one previously found for the canonical (keto) form of cytosine+ . Dynamic electron correlation reduces the energy barrier in the keto form substantially (by ∼1 eV) but it is nevertheless still present. We additionally compute the UV/Vis absorption signals of the structures encountered along these decay channels to provide spectroscopic fingerprints to assist future experiments in monitoring these intricate photo-processes.
Collapse
Affiliation(s)
- Javier Segarra‐Martí
- Department of ChemistryMolecular Sciences Research HubImperial College LondonWhite City Campus, 82 Wood LaneLondonW12 0BZUK
- Present address: Instituto de Ciencia MolecularUniversitat de ValenciaP.O. Box 22085ES-46071ValenciaSpain
| | - Michael J. Bearpark
- Department of ChemistryMolecular Sciences Research HubImperial College LondonWhite City Campus, 82 Wood LaneLondonW12 0BZUK
| |
Collapse
|
5
|
Kochman MA, Palczewski K, Kubas A. Theoretical Study of the Photoisomerization Mechanism of All- Trans-Retinyl Acetate. J Phys Chem A 2021; 125:8358-8372. [PMID: 34546761 PMCID: PMC8488936 DOI: 10.1021/acs.jpca.1c05533] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/04/2022]
Abstract
![]()
The compound 9-cis-retinyl acetate (9-cis-RAc) is a precursor
to 9-cis-retinal,
which has potential application in the treatment of some hereditary
diseases of the retina. An attractive synthetic route to 9-cis-RAc is based on the photoisomerization reaction of the
readily available all-trans-RAc. In the present study,
we examine the mechanism of the photoisomerization reaction with the
use of state-of-the-art electronic structure calculations for two
polyenic model compounds: tEtEt-octatetraene and tEtEtEc-2,6-dimethyl-1,3,5,7,9-decapentaene. The occurrence
of photoisomerization is attributed to a chain-kinking mechanism,
whereby a series of S1/S0 conical intersections
associated with kinking deformations at different positions along
the polyenic chain mediate internal conversion to the S0 state, and subsequent isomerization around one of the double bonds.
Two other possible photoisomerization mechanisms are taken into account,
but they are rejected as incompatible with simulation results and/or
the available spectroscopic data.
Collapse
Affiliation(s)
- Michał Andrzej Kochman
- Institute of Physical Chemistry, Polish Academy of Sciences, Ul. Marcina Kasprzaka 44/52, 01-224 Warszawa, Poland
| | - Krzysztof Palczewski
- Department of Ophthalmology, Gavin Herbert Eye Institute, University of California, Irvine, California 92697, United States.,Department of Physiology and Biophysics, University of California, Irvine, California 92697, United States.,Department of Chemistry, University of California, Irvine, California 92697, United States
| | - Adam Kubas
- Institute of Physical Chemistry, Polish Academy of Sciences, Ul. Marcina Kasprzaka 44/52, 01-224 Warszawa, Poland
| |
Collapse
|
6
|
Werner HJ, Knowles PJ, Manby FR, Black JA, Doll K, Heßelmann A, Kats D, Köhn A, Korona T, Kreplin DA, Ma Q, Miller TF, Mitrushchenkov A, Peterson KA, Polyak I, Rauhut G, Sibaev M. The Molpro quantum chemistry package. J Chem Phys 2020; 152:144107. [PMID: 32295355 DOI: 10.1063/5.0005081] [Citation(s) in RCA: 475] [Impact Index Per Article: 118.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Abstract
Molpro is a general purpose quantum chemistry software package with a long development history. It was originally focused on accurate wavefunction calculations for small molecules but now has many additional distinctive capabilities that include, inter alia, local correlation approximations combined with explicit correlation, highly efficient implementations of single-reference correlation methods, robust and efficient multireference methods for large molecules, projection embedding, and anharmonic vibrational spectra. In addition to conventional input-file specification of calculations, Molpro calculations can now be specified and analyzed via a new graphical user interface and through a Python framework.
Collapse
Affiliation(s)
- Hans-Joachim Werner
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Peter J Knowles
- School of Chemistry, Cardiff University, Main Building, Park Place, Cardiff CF10 3AT, United Kingdom
| | - Frederick R Manby
- School of Chemistry, University of Bristol, Cantock's Close, Bristol BS8 1TS, United Kingdom
| | - Joshua A Black
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Klaus Doll
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Andreas Heßelmann
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Daniel Kats
- Max-Planck Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | - Andreas Köhn
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Tatiana Korona
- Faculty of Chemistry, University of Warsaw, L. Pasteura 1 St., 02-093 Warsaw, Poland
| | - David A Kreplin
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Qianli Ma
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Thomas F Miller
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | | | - Kirk A Peterson
- Department of Chemistry, Washington State University, Pullman, Washington 99164-4630, USA
| | - Iakov Polyak
- School of Chemistry, Cardiff University, Main Building, Park Place, Cardiff CF10 3AT, United Kingdom
| | - Guntram Rauhut
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Marat Sibaev
- School of Chemistry, Cardiff University, Main Building, Park Place, Cardiff CF10 3AT, United Kingdom
| |
Collapse
|
7
|
Park JW, Al-Saadon R, MacLeod MK, Shiozaki T, Vlaisavljevich B. Multireference Electron Correlation Methods: Journeys along Potential Energy Surfaces. Chem Rev 2020; 120:5878-5909. [PMID: 32239929 DOI: 10.1021/acs.chemrev.9b00496] [Citation(s) in RCA: 54] [Impact Index Per Article: 13.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Multireference electron correlation methods describe static and dynamical electron correlation in a balanced way and, therefore, can yield accurate and predictive results even when single-reference methods or multiconfigurational self-consistent field theory fails. One of their most prominent applications in quantum chemistry is the exploration of potential energy surfaces. This includes the optimization of molecular geometries, such as equilibrium geometries and conical intersections and on-the-fly photodynamics simulations, both of which depend heavily on the ability of the method to properly explore the potential energy surface. Because such applications require nuclear gradients and derivative couplings, the availability of analytical nuclear gradients greatly enhances the scope of quantum chemical methods. This review focuses on the developments and advances made in the past two decades. A detailed account of the analytical nuclear gradient and derivative coupling theories is presented. Emphasis is given to the software infrastructure that allows one to make use of these methods. Notable applications of multireference electron correlation methods to chemistry, including geometry optimizations and on-the-fly dynamics, are summarized at the end followed by a discussion of future prospects.
Collapse
Affiliation(s)
- Jae Woo Park
- Department of Chemistry, Chungbuk National University, Chungdae-ro 1, Cheongju 28644, Korea
| | - Rachael Al-Saadon
- Department of Chemistry, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208, United States
| | - Matthew K MacLeod
- Workday, 4900 Pearl Circle East, Suite 100, Boulder, Colorado 80301, United States
| | - Toru Shiozaki
- Department of Chemistry, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208, United States.,Quantum Simulation Technologies, Inc., 625 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
| | - Bess Vlaisavljevich
- Department of Chemistry, University of South Dakota, 414 East Clark Street, Vermillion, South Dakota 57069, United States
| |
Collapse
|
8
|
Sen S, Chattopadhyay A. A Computational Study on the Photochemistry of 2,4,4‐Trimethyl‐1‐pyrroline 1‐Oxide and Investigation on the Reaction Paths of Its Photoproduct Oxaziridine. ChemistrySelect 2019. [DOI: 10.1002/slct.201903230] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Affiliation(s)
- Sindhuja Sen
- Department of ChemistryBirla Institute of Technology and Science (BITS)Pilani –K.K. Birla Goa Campus Goa 403 726 India
| | - Anjan Chattopadhyay
- Department of ChemistryBirla Institute of Technology and Science (BITS)Pilani –K.K. Birla Goa Campus Goa 403 726 India
| |
Collapse
|
9
|
Ibele LM, Nicolson A, Curchod BFE. Excited-state dynamics of molecules with classically driven trajectories and Gaussians. Mol Phys 2019. [DOI: 10.1080/00268976.2019.1665199] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
Affiliation(s)
- Lea M. Ibele
- Department of Chemistry, Durham University, Durham, UK
| | | | | |
Collapse
|
10
|
Segarra‐Martí J, Tran T, Bearpark MJ. Computing the Ultrafast and Radiationless Electronic Excited State Decay of Cytosine and 5‐methyl‐cytosine Cations: Uncovering the Role of Dynamic Electron Correlation. CHEMPHOTOCHEM 2019. [DOI: 10.1002/cptc.201900105] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Affiliation(s)
- Javier Segarra‐Martí
- Department of Chemistry, Molecular Sciences Research HubImperial College London White City Campus, 80 Wood Lane W12 0BZ London UK
| | - Thierry Tran
- Department of Chemistry, Molecular Sciences Research HubImperial College London White City Campus, 80 Wood Lane W12 0BZ London UK
| | - Michael J. Bearpark
- Department of Chemistry, Molecular Sciences Research HubImperial College London White City Campus, 80 Wood Lane W12 0BZ London UK
| |
Collapse
|
11
|
Borin VA, Wiebeler C, Schapiro I. A QM/MM study of the initial excited state dynamics of green-absorbing proteorhodopsin. Faraday Discuss 2019; 207:137-152. [PMID: 29393940 DOI: 10.1039/c7fd00198c] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Abstract
The primary photochemical reaction of the green-absorbing proteorhodopsin is studied by means of a hybrid quantum mechanics/molecular mechanics (QM/MM) approach. The simulations are based on a homology model derived from the blue-absorbing proteorhodopsin crystal structure. The geometry of retinal and the surrounding sidechains in the protein binding pocket were optimized using the QM/MM method. Starting from this geometry the isomerization was studied with a relaxed scan along the C13[double bond, length as m-dash]C14 dihedral. It revealed an "aborted bicycle pedal" mechanism of isomerization that was originally proposed by Warshel for bovine rhodopsin and bacteriorhodopsin. However, the isomerization involved the concerted rotation about C13[double bond, length as m-dash]C14 and C15[double bond, length as m-dash]N, with the latter being highly twisted but not isomerized. Further, the simulation showed an increased steric interaction between the hydrogen at the C14 of the isomerizing bond and the hydroxyl group at the neighbouring tyrosine 200. In addition, we have simulated a nonadiabatic trajectory which showed the timing of the isomerization. In the first 20 fs upon excitation the order of the conjugated double and single bonds is inverted, consecutively the C13[double bond, length as m-dash]C14 rotation is activated for 200 fs until the S1-S0 transition is detected. However, the isomerization is reverted due to the specific interaction with the tyrosine as observed along the relaxed scan calculation. Our simulations indicate that the retinal - tyrosine 200 interaction plays an important role in the outcome of the photoisomerization.
Collapse
Affiliation(s)
- Veniamin A Borin
- Fritz Haber Center for Molecular Dynamics, Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem 91904, Israel.
| | | | | |
Collapse
|
12
|
Kats D, Werner HJ. Multi-state local complete active space second-order perturbation theory using pair natural orbitals (PNO-MS-CASPT2). J Chem Phys 2019; 150:214107. [DOI: 10.1063/1.5097644] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Daniel Kats
- Max-Planck Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | - Hans-Joachim Werner
- Institute for Theoretical Chemistry, University of Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| |
Collapse
|
13
|
Segarra-Martí J, Tran T, Bearpark MJ. Ultrafast and radiationless electronic excited state decay of uracil and thymine cations: computing the effects of dynamic electron correlation. Phys Chem Chem Phys 2019; 21:14322-14330. [PMID: 30698175 DOI: 10.1039/c8cp07189f] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
Abstract
In this article we characterise the radiationless decay of the first few electronic excited states of the cations of DNA/RNA nucleobases uracil and thymine, including the effects of dynamic electron correlation on energies and geometries (optimised with XMS-CASPT2). In both systems, we find that one state of 2n and another two of 2π+ character can be populated following photoionisation, and their different minima and interstate crossings are located. We find strong similarities between uracil and thymine cations: with accessible conical intersections suggesting that depopulation of their electronic excited states takes place on ultrafast timescales in both systems, suggesting that they are photostable in agreement with previous theoretical (uracil+) evidence. We find that dynamic electron correlation separates the energy levels of the "3-state" conical intersection (D2/D1/D0)CI previously located with CASSCF for uracil+, which will therefore have a different geometry and higher energy. Simulating the electronic and vibrational absorptions allows us to characterise spectral fingerprints that could be used to monitor these cation photo-processes experimentally.
Collapse
Affiliation(s)
- Javier Segarra-Martí
- Department of Chemistry, Molecular Sciences Research Hub, Imperial College London, White City Campus, 80 Wood Lane, W12 0BZ, London, UK.
| | - Thierry Tran
- Department of Chemistry, Molecular Sciences Research Hub, Imperial College London, White City Campus, 80 Wood Lane, W12 0BZ, London, UK.
| | - Michael J Bearpark
- Department of Chemistry, Molecular Sciences Research Hub, Imperial College London, White City Campus, 80 Wood Lane, W12 0BZ, London, UK.
| |
Collapse
|
14
|
Sen S, Oruganti Y, Chattopadhyay A. Investigation on the Unexplored Photochemistry of 5,5-Dimethyl-1-pyrroline 1-Oxide (DMPO). J Phys Chem A 2019; 123:163-170. [DOI: 10.1021/acs.jpca.8b08181] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Sindhuja Sen
- Department of Chemistry, Birla Institute of Technology and Science (BITS), Pilani−K.K. Birla Goa Campus, Goa 403 726, India
| | - Yasaswini Oruganti
- Department of Chemistry, Birla Institute of Technology and Science (BITS), Pilani−K.K. Birla Goa Campus, Goa 403 726, India
| | - Anjan Chattopadhyay
- Department of Chemistry, Birla Institute of Technology and Science (BITS), Pilani−K.K. Birla Goa Campus, Goa 403 726, India
| |
Collapse
|
15
|
Ray J, Ramesh SG. Conical intersections involving the lowest 1πσ∗ state in aniline: Role of the NH2 group. Chem Phys 2018. [DOI: 10.1016/j.chemphys.2018.03.015] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
16
|
Sen S, Schapiro I. A comprehensive benchmark of the XMS-CASPT2 method for the photochemistry of a retinal chromophore model. Mol Phys 2018. [DOI: 10.1080/00268976.2018.1501112] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023]
Affiliation(s)
- Saumik Sen
- Fritz Haber Center for Molecular Dynamics, Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem, Israel
| | - Igor Schapiro
- Fritz Haber Center for Molecular Dynamics, Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem, Israel
| |
Collapse
|
17
|
Chattopadhyay A, Saini P, Hakim R, Komainda A, Köppel H. A computational investigation on the photo-isomerization of 2,4,6-octatriene and its UV–visible spectrum. Chem Phys 2018. [DOI: 10.1016/j.chemphys.2017.12.006] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
18
|
Affiliation(s)
- Basile F. E. Curchod
- Department of Chemistry, Durham University, South Road, Durham DH1 3LE, United Kingdom
| | - Todd J. Martínez
- Department of Chemistry and PULSE Institute, Stanford University, Stanford, California 94305, United States
- SLAC National Accelerator Laboratory, Menlo Park, California 94025, United States
| |
Collapse
|
19
|
Szefczyk B, Grabarek D, Walczak E, Andruniów T. Excited-state minima and emission energies of retinal chromophore analogues: Performance of CASSCF and CC2 methods as compared with CASPT2. J Comput Chem 2017; 38:1799-1810. [PMID: 28512740 DOI: 10.1002/jcc.24821] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/06/2017] [Revised: 04/10/2017] [Accepted: 04/13/2017] [Indexed: 11/08/2022]
Abstract
This study provides gas-phase S1 excited-state geometries along with emission and adiabatic energies for methylated/demethylated and ring-locked analogues of protonated Schiff base retinal models comprising system of five conjugated double bonds (PSB5), using second order multiconfiguration perturbation theory (CASPT2). CASPT2 results serve as reference data to assess the performance of CC2 (second-order approximate coupled cluster singles and doubles) and a commonly used CASSCF/CASPT2 protocol, that is, complete active space self-consistent field (CASSCF) geometry optimization followed by CASPT2 energy calculation. We find that the CASSCF methodology fails to locate planar S1 minimum energy structures for four out of five investigated planar models in contrast to CC2 and CASPT2 methods. However, for those which were found: one planar and two twisted minima, there is an excellent agreement between CASSCF and CASPT2 results in terms of geometrical parameters, one-electron properties, as well as emission and adiabatic energies. CC2 performs well for in-plane S1 minima and their spectroscopic and electronic properties. However, this picture deteriorates for twisted minima. As expected, the CC2 description of the S2 electronic state, with strong multireference and significant double excitation character, is very poor, exhibiting errors in transition energies exceeding 1 eV. They may be substantially diminished by recalculating transition energies with CASPT2 method. Our work shows that CASSCF/CASPT2 and CC2 shortcomings may influence gas-phase retinal analogues' excited state description in a dramatic way. © 2017 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Borys Szefczyk
- Advanced Materials Engineering and Modelling Group, Department of Chemistry, Wroclaw University of Science and Technology, Wyb. Wyspianskiego 27, Wroclaw, 50-370, Poland
| | - Dawid Grabarek
- Advanced Materials Engineering and Modelling Group, Department of Chemistry, Wroclaw University of Science and Technology, Wyb. Wyspianskiego 27, Wroclaw, 50-370, Poland
| | - Elżbieta Walczak
- Advanced Materials Engineering and Modelling Group, Department of Chemistry, Wroclaw University of Science and Technology, Wyb. Wyspianskiego 27, Wroclaw, 50-370, Poland
| | - Tadeusz Andruniów
- Advanced Materials Engineering and Modelling Group, Department of Chemistry, Wroclaw University of Science and Technology, Wyb. Wyspianskiego 27, Wroclaw, 50-370, Poland
| |
Collapse
|
20
|
Segarra-Martí J, Francés-Monerris A, Roca-Sanjuán D, Merchán M. Assessment of the Potential Energy Hypersurfaces in Thymine within Multiconfigurational Theory: CASSCF vs. CASPT2. Molecules 2016; 21:molecules21121666. [PMID: 27918489 PMCID: PMC6274573 DOI: 10.3390/molecules21121666] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/15/2016] [Revised: 11/30/2016] [Accepted: 12/01/2016] [Indexed: 11/16/2022] Open
Abstract
The present study provides new insights into the topography of the potential energy hypersurfaces (PEHs) of the thymine nucleobase in order to rationalize its main ultrafast photochemical decay paths by employing two methodologies based on the complete active space self-consistent field (CASSCF) and the complete active space second-order perturbation theory (CASPT2) methods: (i) CASSCF optimized structures and energies corrected with the CASPT2 method at the CASSCF geometries and (ii) CASPT2 optimized geometries and energies. A direct comparison between these strategies is drawn, yielding qualitatively similar results within a static framework. A number of analyses are performed to assess the accuracy of these different computational strategies under study based on a variety of numerical thresholds and optimization methods. Several basis sets and active spaces have also been calibrated to understand to what extent they can influence the resulting geometries and subsequent interpretation of the photochemical decay channels. The study shows small discrepancies between CASSCF and CASPT2 PEHs, displaying a shallow planar or twisted 1(ππ*) minimum, respectively, and thus featuring a qualitatively similar scenario for supporting the ultrafast bi-exponential deactivation registered in thymine upon UV-light exposure. A deeper knowledge of the PEHs at different levels of theory provides useful insight into its correct characterization and subsequent interpretation of the experimental observations. The discrepancies displayed by the different methods studied here are then discussed and framed within their potential consequences in on-the-fly non-adiabatic molecular dynamics simulations, where qualitatively diverse outcomes are expected.
Collapse
Affiliation(s)
- Javier Segarra-Martí
- Instituto de Ciencia Molecular, Universitat de València, P. O. Box 22085, ES-46071 Valencia, Spain.
- Present Address: Laboratoire de Chimie UMR 5182, École Normale Supérieure de Lyon, CNRS, Université de Lyon, 46 Allée d'Italie, F-69364 Lyon Cedex 07, France.
| | - Antonio Francés-Monerris
- Instituto de Ciencia Molecular, Universitat de València, P. O. Box 22085, ES-46071 Valencia, Spain.
| | - Daniel Roca-Sanjuán
- Instituto de Ciencia Molecular, Universitat de València, P. O. Box 22085, ES-46071 Valencia, Spain.
| | - Manuela Merchán
- Instituto de Ciencia Molecular, Universitat de València, P. O. Box 22085, ES-46071 Valencia, Spain.
| |
Collapse
|
21
|
Vlaisavljevich B, Shiozaki T. Nuclear Energy Gradients for Internally Contracted Complete Active Space Second-Order Perturbation Theory: Multistate Extensions. J Chem Theory Comput 2016; 12:3781-7. [DOI: 10.1021/acs.jctc.6b00572] [Citation(s) in RCA: 76] [Impact Index Per Article: 9.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Bess Vlaisavljevich
- Department of Chemistry, Northwestern University, 2145 Sheridan
Road, Evanston, Illinois 60208, United States
| | - Toru Shiozaki
- Department of Chemistry, Northwestern University, 2145 Sheridan
Road, Evanston, Illinois 60208, United States
| |
Collapse
|
22
|
Saini P, Banerjee M, Chattopadhyay A. Computational Investigation of the Photochemical Reaction Path of Some Synthesized and Experimentally Analyzed Small-Chain Conjugated Nitrones. J Phys Chem A 2016; 120:396-406. [PMID: 26697926 DOI: 10.1021/acs.jpca.5b11069] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Praveen Saini
- Department
of Chemistry, Birla Institute of Technology and Science (BITS), Pilani−K.K.
Birla Goa Campus, Zuarinagar, Goa 403 726, India
| | - Mainak Banerjee
- Department
of Chemistry, Birla Institute of Technology and Science (BITS), Pilani−K.K.
Birla Goa Campus, Zuarinagar, Goa 403 726, India
| | - Anjan Chattopadhyay
- Department
of Chemistry, Birla Institute of Technology and Science (BITS), Pilani−K.K.
Birla Goa Campus, Zuarinagar, Goa 403 726, India
| |
Collapse
|
23
|
Chattopadhyay A, Saini P, Pandharkar R. Exploring the isomerization paths of push–pull hexatrienes. RSC Adv 2016. [DOI: 10.1039/c6ra16812d] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022] Open
Abstract
The first photo-excited singlet state of all-trans donor (amine), acceptor (cyano)-substituted hexatriene forms donor side-twisted minima and its passage towards the cis–trans–trans isomer is characterized by an S0/S1 conical intersection.
Collapse
Affiliation(s)
- Anjan Chattopadhyay
- Department of Chemistry
- Birla Institute of Technology and Science (BITS)
- Pilani – K.K. Birla Goa Campus
- India
| | - Praveen Saini
- Department of Chemistry
- Birla Institute of Technology and Science (BITS)
- Pilani – K.K. Birla Goa Campus
- India
| | - Riddhish Pandharkar
- Department of Chemistry
- Birla Institute of Technology and Science (BITS)
- Pilani – K.K. Birla Goa Campus
- India
| |
Collapse
|
24
|
Zimmermann T, Vaníček J. Efficient on-the-fly ab initio semiclassical method for computing time-resolved nonadiabatic electronic spectra with surface hopping or Ehrenfest dynamics. J Chem Phys 2015; 141:134102. [PMID: 25296779 DOI: 10.1063/1.4896735] [Citation(s) in RCA: 34] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We derive a somewhat crude, yet very efficient semiclassical approximation for computing nonadiabatic spectra. The resulting method, which is a generalization of the multiple-surface dephasing representation, includes quantum effects through interference of mixed quantum-classical trajectories and through quantum treatment of the collective electronic degree of freedom. The method requires very little computational effort beyond the fewest-switches surface hopping or Ehrenfest locally mean-field dynamics and is very easy to implement. The proposed approximation is tested by computing the absorption and time-resolved stimulated emission spectra of pyrazine using the four-dimensional three-surface model which allows for comparison with the numerically exact quantum spectra. As expected, the multiple-surface dephasing representation is not suitable for high-resolution linear spectra, yet it seems to capture all the important features of pump-probe spectra. Finally, the method is combined with on-the-fly ab initio evaluation of the electronic structure (i.e., energies, forces, electric-dipole, and nonadiabatic couplings) in order to compute fully dimensional nonadiabatic spectra of pyrazine without approximations inherent to analytical, including vibronic-coupling models. The Appendix provides derivations of perturbative expressions for linear and pump-probe spectra of arbitrary mixed states and for arbitrary laser pulse shapes.
Collapse
Affiliation(s)
- Tomáš Zimmermann
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Jiří Vaníček
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| |
Collapse
|
25
|
Sala M, Lasorne B, Gatti F, Guérin S. The role of the low-lying dark nπ* states in the photophysics of pyrazine: a quantum dynamics study. Phys Chem Chem Phys 2015; 16:15957-67. [PMID: 24964033 DOI: 10.1039/c4cp02165g] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The excited state dynamics of pyrazine has attracted considerable attention in the last three decades. It has long been recognized that after UV excitation, the dynamics of the molecule is impacted by strong non-adiabatic effects due to the existence of a conical intersection between the B2u(ππ*) and B3u(nπ*) electronic states. However, a recent study based on trajectory surface hopping dynamics simulations suggested the participation of the Au(nπ*) and B2g(nπ*) low-lying dark electronic states in the ultrafast radiationless decay of the molecule after excitation to the B2u(ππ*) state. The purpose of this work was to pursue the investigation of the role of the Au(nπ*) and B2g(nπ*) states in the photophysics of pyrazine. A linear vibronic coupling model hamiltonian including the four lowest excited electronic states and the sixteen most relevant vibrational degrees of freedom was constructed using high level XMCQDPT2 electronic structure calculations. Wavepacket propagations using the MCTDH method were then performed and used to simulate the absorption spectrum and the electronic state population dynamics of the system. Our results show that the Au(nπ*) state plays an important role in the photophysics of pyrazine.
Collapse
Affiliation(s)
- Matthieu Sala
- Laboratoire Interdisciplinaire Carnot de Bourgogne UMR 6303 CNRS, Université de Bourgogne, BP 47870, F-21078 Dijon, France.
| | | | | | | |
Collapse
|
26
|
Sala M, Guérin S, Gatti F. Quantum dynamics of the photostability of pyrazine. Phys Chem Chem Phys 2015; 17:29518-30. [DOI: 10.1039/c5cp04605j] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
We propose a new mechanism for the radiationless decay of photoexcited pyrazine to its ground electronic state involving a conical intersection between the dark Au(nπ) state and the ground state.
Collapse
Affiliation(s)
- Matthieu Sala
- Laboratoire Interdisciplinaire Carnot de Bourgogne UMR 6303 CNRS
- Université de Bourgogne Franche-Comté
- F-21078 Dijon
- France
| | - Stéphane Guérin
- Laboratoire Interdisciplinaire Carnot de Bourgogne UMR 6303 CNRS
- Université de Bourgogne Franche-Comté
- F-21078 Dijon
- France
| | - Fabien Gatti
- CTMM
- Institut Charles Gerhardt UMR 5253 CNRS
- F-34095 Montpellier
- France
| |
Collapse
|
27
|
Saab M, Sala M, Lasorne B, Gatti F, Guérin S. Full-dimensional control of the radiationless decay in pyrazine using the dynamic Stark effect. J Chem Phys 2014; 141:134114. [DOI: 10.1063/1.4896938] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Mohamad Saab
- CTMM, Institut Charles Gerhardt Montpellier (UMR5253), CC 15001, Université Montpellier 2, F-34095 Montpellier, Cedex 05, France
| | - Matthieu Sala
- Laboratoire Interdisciplinaire Carnot de Bourgogne (UMR 6303) CNRS, Université de Bourgogne, BP 47870, F-21078 Dijon, France
| | - Benjamin Lasorne
- CTMM, Institut Charles Gerhardt Montpellier (UMR5253), CC 15001, Université Montpellier 2, F-34095 Montpellier, Cedex 05, France
| | - Fabien Gatti
- CTMM, Institut Charles Gerhardt Montpellier (UMR5253), CC 15001, Université Montpellier 2, F-34095 Montpellier, Cedex 05, France
| | - Stéphane Guérin
- Laboratoire Interdisciplinaire Carnot de Bourgogne (UMR 6303) CNRS, Université de Bourgogne, BP 47870, F-21078 Dijon, France
| |
Collapse
|
28
|
Sala M, Saab M, Lasorne B, Gatti F, Guérin S. Laser control of the radiationless decay in pyrazine using the dynamic Stark effect. J Chem Phys 2014; 140:194309. [DOI: 10.1063/1.4875736] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
29
|
Győrffy W, Shiozaki T, Knizia G, Werner HJ. Analytical energy gradients for second-order multireference perturbation theory using density fitting. J Chem Phys 2013; 138:104104. [DOI: 10.1063/1.4793737] [Citation(s) in RCA: 86] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
30
|
|