DOCSLIB.ORG

Femtosecond Intersystem Crossing in the DNA Nucleobase Cytosine

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
Femtosecond Intersystem Crossing in the DNA Nucleobase Cytosine Femtosecond Intersystem Crossing in the DNA Nucleobase Cytosine Martin Richter,y Philipp Marquetand,z Jes´usGonz´alez-Vazquez,{ Ignacio Sola,{ and Leticia Gonz´alez∗,z yInstitute of Physical Chemistry, Friedrich Schiller University Jena, Helmholtzweg 4, 07743 Jena, Germany zInstitute of Theoretical Chemistry, University of Vienna, W¨ahringerStr. 17, 1090 Vienna, Austria {Departamento de Qu´ımica F´ısica I, Universidad Complutense, 28040 Madrid, Spain E-mail: [email protected] Keywords: DNA photostability, excited Graphical TOC Entry state dynamics, intersystem crossing, spin-orbit coupling, conical intersection Abstract Ab initio molecular dynamics including non- adiabatic and spin-orbit couplings on equal footing is used to unravel the deactivation of cytosine after UV light absorption. Intersys- tem crossing (ISC) is found to compete directly with internal conversion in tens of femtosec- onds, thus making cytosine the organic com- pound with the fastest triplet population cal- culated so far. It is found that close degen- eracy between singlet and triplet states can more than compensate for very small spin-orbit couplings, leading to efficient ISC. The fem- tosecond nature of the intersystem crossing pro- cess highlights its importance in photochem- arXiv:2103.13070v1 [physics.chem-ph] 24 Mar 2021 istry and challenges the conventional view that large singlet-triplet couplings are required for an efficient population flow into triplet states. These findings are important to understand DNA photostability and the photochemistry and dynamics of organic molecules in general. 1 The interaction of DNA and RNA with ra- goes beyond what can be learned from quan- diation, from mobile-phone emissions1 to UV tum chemical calculations alone. wavelengths,2 has enthralled the scientific com- Cytosine presents three tautomers: enol-, munity for years due to its implications in pho- keto- and keto-imine-cytosine. Keto-cytosine todamage.3 Of particular interest is to under- is the biological relevant tautomer found in stand photostability, i.e. the relaxation mech- the DNA's nucleotides linked to the deoxyri- anisms that bring DNA4{6 to the ground state bose sugar moiety and the only one for which before any other photoreaction can occur. This a crystalline structure exists.40 Therefore, here means that, instead of fluorescence or phospho- we focus on keto-cytosine. Several station- rescence, the electronic energy provided upon ary points3,41,42 of keto-cytosine, including two- photoexcitation in DNA is transferred to the 43{48 and three-state49{51 conical intersections nuclear degrees of freedom of the molecular sys- involved in the process of IC have been calcu- tem in different ways. It is precisely the atom- lated with ab initio methods. Time-dependent istic description of these different relaxation calculations have indicated that the dynami- pathways that is still heavily discussed in the cal behavior of cytosine after photoexcitation literature. In the last years it has been clearly is one of the most complicated among nucle- established that excited states of isolated DNA obases, involving delocalization of the excited nucleobases undergo ultrafast internal conver- wave packet and relaxation through multiple sion (IC) allowing for an efficient radiationless competing pathways in the singlet excited state decay towards lower-lying electronic states.7{15 manifold.6,52{54 The possible triplet state for- The role of intersystem crossing (ISC) in the mation via ISC along the internal conversion process of photostability is, however, much less pathway of excited singlet keto-cytosine has discussed,15{19 probably because it is thought been discussed by Merch´anet al.55,56 in the to be a much slower process in comparison to light of quantum chemical calculations. IC20 and also because the quantum yields of Our ab initio molecular dynamic simula- triplet states population in DNA and RNA nu- tions are performed on seven states simultane- cleobases are generally very small and thus dif- ously: four singlets and three triplets. Ener- ficult to access from the experimental point of gies, energy gradients, non-adiabatic and spin- view.10,15,21 We note, however, that ultrafast orbit couplings are computed on-the-fly using time scales for ISC in other organic molecules the state-average Complete Active Space Self- have been experimentally reported or predicted Consistent Field (CASSCF) method.57,58 Fur- before.22{36 ther details are given in the Supporting Infor- In this work we present the first excited state mation (SI). The first singlet excited state, S1, dynamical study of a DNA nucleobase includ- has ππ∗ character at the equilibrium geometry ing singlet and triplet states. Such simulations and it is bright while the states higher in energy, ∗ are done using the newly developed surface- S2 and S3, correspond to dark nπ excitations, hopping method SHARC.37 SHARC stands for i.e. they have negligible oscillator strengths surface hopping including arbitrary couplings. when vertically excited. The order of states at With SHARC one can treat non-adiabatic and equilibrium geometry is not altered when go- spin-orbit couplings (which mediate IC and ing to higher levels of theory that include dy- ISC, respectively) on equal footing. The appli- namical correlation (see Table S2 of SI). How- cability of SHARC to include spin-orbit as well ever, the on-the-fly approach used in this work as dipole couplings is documented in Refs.37{39 prohibits the use of higher level methods such Here, SHARC is employed to investigate the role as CASPT2 and therefore we employ CASSCF. of the triplet states in the deactivation of cy- The Franck-Condon region {from which exci- tosine within the framework of nonadiabatic tations take place{ does not only comprise the molecular dynamics based on ab initio multi- equilibrium geometry but also slightly distorted configurational methods. Such study is nec- geometries. These distortions are due to the essary to provide a mechanistic insight that different vibrations included within the zero- 2 point energy of the system. Because in cyto- that in some cases numbers can be higher than sine rather small deviations of the equilibrium the initial population. Percentages not adding geometry lead to a different ordering of the to 100% are due to minor pathways not indi- state character, the S2 and S3 states can also cated. The decay times are obtained fitting be bright states and contribute to the absorp- the net amount of hops between two particu- tion spectrum (see Refs.6,53 and Fig. S1 of SI). lar states (see Fig. S3 and Table S1 of SI) to an The character of the lowest three excited triplet exponential function. The branching efficien- ∗ states at equilibrium is ππ for T1 and T2, and cies given in % are also graphically indicated by ∗ nπ for T3. A comprehensive report of vertical arrows of different thickness according to their excitation energies at different levels of theory importance. After photo-excitation, which cor- can be found in the SI. responds to time zero in our simulations, the In order to obtain a global picture of the re- population of the ππ∗ is distributed as 13% in laxation mechanisms of keto-cytosine, we have S3, 47% in S2, and 40% in S1, as dictated by the first used initial conditions spanning the whole weight of the oscillator strength of each state. first absorbtion band of the UV spectrum, i.e. Since the character of a state can adiabatically covering excitation energies from ca 4 to 7 eV. change during the simulation, hereafter we shall As explained above, this requires launching tra- refer to the states by its energetic order rather jectories from the first three excited states, S1, than by their character. S2 and S3. Most time-resolved spectroscopic ex- Analyzing the 13% population of S3, 10% re- 9,12,48,59,60 periments in cytosine use a pumping laxes non-adiabatically to the S2 and from there wavelength of 267 nm (4.64 eV), just below the to the S1 within 25 fs. After ca 155 fs, the sys- center of the first absorption band located at tem populates the electronic ground state S0. 260 nm (4.77 eV). In order to narrow the ini- The remaining 3% of the population of S3 deac- tial conditions to the energy range correspond- tivates directly to S1 via a threefold degeneracy 51 ing to the experimental one, we have also ana- S3/S2/S1, as proposed in Ref. This process is lyzed the results (Figure S2 of SI) limited to the calculated to be slower than the previous one, bandwidth 4.75±0.25 eV, just below our theo- with a time constant of 110 fs. Most of the pop- retically predicted first absorption band maxi- ulation in S2 transfers preferably to the lower- mum. Also in this energy range, states S1 to lying electronic states within less than 100 fs. S3 are excited. The results from both energy Also in this case the process of IC is possible via ranges qualitatively agree with each other so a cascade of subsequent S2/S1 and S1/S0 conical that we will limit the discussion to the more intersections, or directly via three-state conical 49,51 general broad range. intersections S2/S1/S0, as proposed in Refs. Figure 1 displays the time evolution of all Both pathways to the ground state are relevant, the state populations and Fig. 2 summarizes in agreement with the time-dependent simula- the most important deactivation paths found in tions of Ref.53 As deduced from the time con- keto-cytosine with SHARC, including decay times stants, the three-state conical intersection path- and branching efficiencies. One should note way is faster (25 fs) than the two-step pathway that because the calculations are done at lev- (25 fs / 155 fs).
Load more

Recommended publications
  • Multidimensional Dynamics Involving a Conical Intersection: Wavepacket Calculations Using the MCTDH Method
    Multidimensional dynamics involving a conical intersection: Wavepacket calculations using the MCTDH method G. A. Worth†, H.-D. Meyer∗, and L. S. Cederbaum∗ ∗ Theoretical Chemistry, University of Heidelberg, INF 229, 69120 Heidelberg, Germany † Dept. of Chemistry, King’s College London, The Strand, London, WC2R 2LS, U.K. This article will appear in Conical Intersections, W. Domcke, D. R. Yarkony, and H. K¨oppel, Eds., World Scientific Co., (Singapore, 2004), pp. 573–607. June 1, 2004 1 Introduction The time-evolution of a molecular system through a conical intersection is multi- dimensional in nature, i.e. its description cannot be reduced to a single reaction coordinate. In addition to this, such processes are inherently quantum mechan- ical due to the mixing of electronic states as the intersection is traversed, and a powerful quantum dynamics method is required to treat these, often large, systems accurately. The most intuitive way of obtaining a numerically exact solution of the time- dependent Schr¨odinger equation is to use standard wavepacket propagation meth- ods. These represent the evolving system wavefunction using basis set expansion techniques or grid representations [1–3]. Unfortunately, the standard methods are restricted by an exponential increase in the computational resources required with the number of degrees of freedom treated, usually no more than 4-6 nuclear modes can be included. Systems in which vibronic coupling plays a role are often much larger than this. To treat these systems using the standard methods it is then necessary to turn to reduced dimensionality models in which only a few degrees of freedom are explicitly included, and the remainder either ignored or included in a phenomenological manner.
    [Show full text]
  • European Conference on the Dynamics of Molecular Systems Book of Abstracts
    European conference on the dynamics of molecular systems Book of abstracts Department of Chemistry Jesus College University of Oxford 1 2 Plenary talks 3 PLENARY TALK P1 Imaging studies of inelastic collisions David H. Parker Dept of Molecular and Laser Physics, IMM, Radboud University, Nijmegen, the Netherlands. E‐mail address: [email protected] State-to-state imaging studies of the interaction of small molecules including O2,H2O, OH, CO and ammonia with H2 and with He will be described. These molecules are prominent components of molecular interstellar matter. In order to extract quantitative information on these molecules from telescope spectra, models are employed that depend critically on the rates of rotational energy exchange due to collisions with molecular hydrogen and helium. Collision rates are currently determined by theory from the multidimensional Potential Energy Surface (PES) describing the interaction of H2O and H2 [1] or He [2] and OH with H2 [3], or He. Our velocity map imaging [4] measurements of state-to-state differential and relative integral cross sections of rotational inelastic collisions, also as a function of collision energy, are used to test these PESs. For the H2O/H2, He system experiment is compared with state-of-the-art theoretical calculations by the group of L. Wiesenfeld (Grenoble) [5]. Our studies of astrochemistry relevant small molecules with the collision partners He and H2 at collision energy relevant to that of the interstellar media should place the theoretically determined PESs and the collision rates extracted from the PES on a firmer basis. References 1. P. Valiron M. Wernli, A.
    [Show full text]
  • The Dissipative Photochemical Origin of Life: UVC Abiogenesis of Adenine
    Preprints (www.preprints.org) | NOT PEER-REVIEWED | Posted: 25 January 2021 doi:10.20944/preprints202101.0500.v1 Article The Dissipative Photochemical Origin of Life: UVC Abiogenesis of Adenine Karo Michaelian Department of Nuclear Physics and Application of Radiation, Instituto de Física, Universidad Nacional Autónoma de México, Circuito Interior de la Investigación Científica, Cuidad Universitaria, Cuidad de México, C.P. 04510.; karo@fisica.unam.mx Academic Editor: Karo Michaelian Version January 21, 2021 submitted to Entropy 1 Abstract: I describe the non-equilibrium thermodynamics and the photochemical mechanisms which 2 may have been involved in the dissipative structuring, proliferation and evolution of the fundamental 3 molecules at the origin of life from simpler and more common precursor molecules under the 4 impressed UVC photon flux of the Archean. Dissipative structuring of the fundamental molecules is 5 evidenced by their strong and broad wavelength absorption bands and rapid radiationless dexcitation 6 in this wavelength region. Molecular configurations with large photon dissipative efficacy become 7 dominant through dissipative selection. Proliferation arises from the auto- and cross-catalytic nature 8 of the intermediate products. This inherent non-linearity gives rise to numerous stationary states. 9 Amplification of a molecular concentration fluctuation near a bifurcation allows evolution of the 10 concentration profile towards states of generally greater photon disspative efficacy. An example is 11 given of photochemical dissipative abiogenesis of adenine from the precursors HCN and H2O within 12 a fatty acid vesicle on a hot ocean surface, driven far from equilibrium by the impressed UVC light. 13 The kinetic equations for the photochemical reactions with diffusion are resolved under different 14 environmental conditions and the results analyzed within the framework of Classical Irreversible 15 Thermodynamic theory.
    [Show full text]
  • The Dissipative Photochemical Origin of Life: UVC Abiogenesis of Adenine
    entropy Article The Dissipative Photochemical Origin of Life: UVC Abiogenesis of Adenine Karo Michaelian Department of Nuclear Physics and Applications of Radiation, Instituto de Física, Universidad Nacional Autónoma de México, Circuito Interior de la Investigación Científica, Cuidad Universitaria, Mexico City, C.P. 04510, Mexico; karo@fisica.unam.mx Abstract: The non-equilibrium thermodynamics and the photochemical reaction mechanisms are described which may have been involved in the dissipative structuring, proliferation and complex- ation of the fundamental molecules of life from simpler and more common precursors under the UVC photon flux prevalent at the Earth’s surface at the origin of life. Dissipative structuring of the fundamental molecules is evidenced by their strong and broad wavelength absorption bands in the UVC and rapid radiationless deexcitation. Proliferation arises from the auto- and cross-catalytic nature of the intermediate products. Inherent non-linearity gives rise to numerous stationary states permitting the system to evolve, on amplification of a fluctuation, towards concentration profiles providing generally greater photon dissipation through a thermodynamic selection of dissipative efficacy. An example is given of photochemical dissipative abiogenesis of adenine from the precursor HCN in water solvent within a fatty acid vesicle floating on a hot ocean surface and driven far from equilibrium by the incident UVC light. The kinetic equations for the photochemical reactions with diffusion are resolved under different environmental conditions and the results analyzed within the framework of non-linear Classical Irreversible Thermodynamic theory. Keywords: origin of life; dissipative structuring; prebiotic chemistry; abiogenesis; adenine; organic molecules; non-equilibrium thermodynamics; photochemical reactions Citation: Michaelian, K. The Dissipative Photochemical Origin of MSC: 92-10; 92C05; 92C15; 92C40; 92C45; 80Axx; 82Cxx Life: UVC Abiogenesis of Adenine.
    [Show full text]
  • Hydrogen-Adduction to Open-Shell Graphene Fragments: Spectroscopy, Thermochemistry and Cite This: Chem
    Chemical Science View Article Online EDGE ARTICLE View Journal | View Issue Hydrogen-adduction to open-shell graphene fragments: spectroscopy, thermochemistry and Cite this: Chem. Sci.,2017,8,1186 astrochemistry Gerard D. O'Connor,a Bun Chan,bc Julian A. Sanelli,d Katie M. Cergol,b Viktoras Dryza,d Richard J. Payne,b Evan J. Bieske,d Leo Radomb and Timothy W. Schmidt*a We apply a combination of state-of-the-art experimental and quantum-chemical methods to elucidate the electronic and chemical energetics of hydrogen adduction to a model open-shell graphene fragment. The lowest-energy adduct, 1H-phenalene, is determined to have a bond dissociation energy of 258.1 kJ molÀ1, while other isomers exhibit reduced or in some cases negative bond dissociation energies, the metastable species being bound by the emergence of a conical intersection along the high-symmetry dissociation coordinate. The gas-phase excitation spectrum of 1H-phenalene and its radical cation are recorded Received 24th August 2016 using laser spectroscopy coupled to mass-spectrometry. Several electronically excited states of both Accepted 24th September 2016 Creative Commons Attribution 3.0 Unported Licence. species are observed, allowing the determination of the excited-state bond dissociation energy. The DOI: 10.1039/c6sc03787a ionization energy of 1H-phenalene is determined to be 7.449(17) eV, consistent with high-level W1X-2 www.rsc.org/chemicalscience calculations. 1 Introduction Using density functional theory it has been calculated that semihydrogenated graphene (graphone) becomes a ferromag- Graphene is a material consisting of a single layer of carbon netic semiconductor with a small indirect gap.7 It has also been 1,2 atoms bonded in a hexagon lattice.
    [Show full text]
  • Metal-Organic Charge Transfer Can Produce Biradical States and Is Mediated by Conical Intersections
    Metal-organic charge transfer can produce biradical INAUGURAL ARTICLE states and is mediated by conical intersections Oksana Tishchenkoa,1, Ruifang Lia,b, and Donald G. Truhlara,1 aDepartment of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, MN 55455-0431; and bDepartment of Chemistry, Nankai University, Tianjin 300071, People’s Republic of China This article is part of the special series of Inaugural Articles by members of the National Academy of Sciences elected in 2009. Edited by Nicholas J. Turro, Columbia University, New York, NY, and approved September 20, 2010 (received for review July 14, 2010) The present paper illustrates key features of charge transfer be- theory revealed the presence of two types of equilibrium struc- tween calcium atoms and prototype conjugated hydrocarbons tures: one at a distance R between Ca and the closest carbon (ethylene, benzene, and coronene) as elucidated by electronic of about 3.8 Å with slightly negative charge on the calcium atom, structure calculations. One- and two-electron charge transfer is and another state with a smaller R and with considerable shift of controlled by two sequential conical intersections. The two lowest electrons from Ca to coronene; the latter is metastable. To get a electronic states that undergo a conical intersection have closed- better understanding of the structure of potential energy surfaces shell and open-shell dominant configurations correlating with that govern the charge transfer process and of the nature of the 4s2 and 4s13d1 states of Ca, respectively. Unlike the neutral- electronic states involved, we performed a detailed electronic ionic state crossing in, for example, hydrogen halides or alkali structure study using multiconfigurational wave functions for halides, the path from separated reactants to the conical intersec- the interaction of Ca with the smaller prototype molecules ben- tion region is uphill and the charge-transferred state is a biradical.
    [Show full text]
  • Photoprotective Properties of Eumelanin: Computational Insights Into the Photophysics of a Catechol:Quinone Heterodimer Model System
    Article Photoprotective Properties of Eumelanin: Computational Insights into the Photophysics of a Catechol:Quinone Heterodimer Model System Victoria C. Frederick 1, Thomas A. Ashy 1, Barbara Marchetti 1, Michael N. R. Ashfold 2 and Tolga N. V. Karsili 1,* 1 Department of Chemistry, University of Louisiana at Lafayette, Lafayette, LA 70503, USA; [email protected] (V.C.F.); [email protected] (T.A.A.); [email protected] (B.M.) 2 School of Chemistry, University of Bristol, Bristol BS8 1TS, UK; [email protected] * Correspondence: [email protected] Abstract: Melanins are skin-centered molecular structures that block harmful UV radiation from the sun and help protect chromosomal DNA from UV damage. Understanding the photodynamics of the chromophores that make up eumelanin is therefore paramount. This manuscript presents a multi-reference computational study of the mechanisms responsible for the experimentally observed photostability of a melanin-relevant model heterodimer comprising a catechol (C)–benzoquinone (Q) pair. The present results validate a recently proposed photoinduced intermolecular transfer of an H atom from an OH moiety of C to a carbonyl-oxygen atom of the Q. Photoexcitation of the ground state C:Q heterodimer (which has a π-stacked “sandwich” structure) results in population of a locally excited ππ* state (on Q), which develops increasing charge-transfer (biradical) character as it evolves Citation: Frederick, V.C.; Ashy, T.A.; to a “hinged” minimum energy geometry and drives proton transfer (i.e., net H atom transfer) from Marchetti, B.; Ashfold, M.N.R.; C to Q. The study provides further insights into excited state decay mechanisms that could contribute Karsili, T.N.V.
    [Show full text]
  • Sequential Electron Transfer Governs the UV-Induced Self-Repair of DNA Photolesions
    Edinburgh Research Explorer Sequential electron transfer governs the UV-induced self-repair of DNA photolesions Citation for published version: Szabla, R, Kruse, H, Stadlbauer, P, Šponer, J & Sobolewski, AL 2018, 'Sequential electron transfer governs the UV-induced self-repair of DNA photolesions', Chemical Science, vol. 9, no. 12, pp. 3131-3140. https://doi.org/10.1039/c8sc00024g Digital Object Identifier (DOI): 10.1039/c8sc00024g Link: Link to publication record in Edinburgh Research Explorer Document Version: Publisher's PDF, also known as Version of record Published In: Chemical Science General rights Copyright for the publications made accessible via the Edinburgh Research Explorer is retained by the author(s) and / or other copyright owners and it is a condition of accessing these publications that users recognise and abide by the legal requirements associated with these rights. Take down policy The University of Edinburgh has made every reasonable effort to ensure that Edinburgh Research Explorer content complies with UK legislation. If you believe that the public display of this file breaches copyright please contact [email protected] providing details, and we will remove access to the work immediately and investigate your claim. Download date: 28. Sep. 2021 Chemical Science View Article Online EDGE ARTICLE View Journal | View Issue Sequential electron transfer governs the UV- induced self-repair of DNA photolesions† Cite this: Chem. Sci.,2018,9,3131 ab b bc b Rafał Szabla, * Holger Kruse, Petr Stadlbauer, Jiˇr´ı Sponerˇ and Andrzej L. Sobolewskia Cyclobutane pyrimidine dimers (CpDs) are among the most common DNA lesions occurring due to the interaction with ultraviolet light.
    [Show full text]
  • Evidence for Conical Intersection Dynamics Mediating Ultrafast Singlet Exciton Fission
    View metadata, citation and similar papers at core.ac.uk brought to you by CORE provided by Apollo Evidence for Conical Intersection Dynamics Mediating Ultrafast Singlet Exciton Fission Andrew J. Musser1*, Matz Liebel2*, Christoph Schnedermann2, Torsten Wende2, Tom B. 1 1 2 Kehoe , Akshay Rao , Philipp Kukura 1 Cavendish Laboratory, University of Cambridge J J Thomson Avenue, Cambridge, UK 2 Physical and Theoretical Chemistry Laboratory, University of Oxford South Parks Road, Oxford, UK Email: [email protected], [email protected] *These authors contributed equally to this work. Introductory Paragraph Singlet exciton fission is the process in organic semiconductors through which a spin-singlet exciton converts into a pair of spin-triplet excitons residing on different chromophores, entangled in an overall spin-zero state. For some systems, singlet fission has been shown to occur on the 100 fs timescale and with a 200% yield, but the mechanism of this process remains uncertain. Here we study a model singlet fission system, TIPS-pentacene, using ultrafast vibronic spectroscopy. We observe that vibrational coherence in the initially photogenerated singlet state is transferred to the triplet state and show that this behaviour is effectively identical to that observed in ultrafast internal conversion for polyenes in solution. This similarity in vibronic dynamics suggest that both multi-molecular singlet fission and single-molecular internal conversion are mediated by the same underlying relaxation processes, based on strong coupling between nuclear and electronic degrees of freedom. In its most efficient form this leads to a conical intersection between the coupled electronic states. 1 Main Text Singlet exciton fission has long commanded interest as an exceptionally fast channel to generate triplet excitons in organic materials.
    [Show full text]
  • Ultrafast Deactivation of an Excited Cytosine-Guanine Base Pair in DNA Gerrit Groenhof,† Lars V
    Published on Web 05/08/2007 Ultrafast Deactivation of an Excited Cytosine-Guanine Base Pair in DNA Gerrit Groenhof,† Lars V. Scha¨fer,† Martial Boggio-Pasqua,‡ Maik Goette,† Helmut Grubmu¨ller,† and Michael A. Robb*,‡ Contribution from the Department of Theoretical and Computational Biophysics, Max-Planck-Institute for Biophysical Chemistry, Am Fassberg 11, D-37077 Go¨ttingen, Germany, and Department of Chemistry, Imperial College, London SW 2AZ, U.K. Received December 31, 2006; E-mail: [email protected] Abstract: Multiconfigurational ab initio calculations and QM/MM molecular dynamics simulations of a photoexcited cytosine-guanine base pair in both gas phase and embedded in the DNA provide detailed structural and dynamical insights into the ultrafast radiationless deactivation mechanism. Photon absorption promotes transfer of a proton from the guanine to the cytosine. This proton transfer is followed by an efficient radiationless decay of the excited state via an extended conical intersection seam. The optimization of the conical intersection revealed that it has an unusual topology, in that there is only one degeneracy- lifting coordinate. This is the central mechanistic feature for the decay both in vacuo and in the DNA. Radiationless decay occurs along an extended hyperline nearly parallel to the proton-transfer coordinate, indicating the proton transfer itself is not directly responsible for the deactivation. The seam is displaced from the minimum energy proton-transfer path along a skeletal deformation of the bases. Decay can thus occur anywhere along the single proton-transfer coordinate, accounting for the remarkably short excited- state lifetime of the Watson-Crick base pair. In vacuo, decay occurs after a complete proton transfer, whereas in DNA, decay can also occur much earlier.
    [Show full text]
  • Ultrafast 25-Fs Relaxation in Highly Excited States of Methyl Azide Mediated by Strong Nonadiabatic Coupling
    Ultrafast 25-fs relaxation in highly excited states of PNAS PLUS methyl azide mediated by strong nonadiabatic coupling William K. Petersa,b,1, David E. Coucha,b, Benoit Mignoletc, Xuetao Shid, Quynh L. Nguyena,b, Ryan C. Fortenberrye, H. Bernhard Schlegeld, Françoise Remaclec, Henry C. Kapteyna,b, Margaret M. Murnanea,b,1, and Wen Lid aJILA, University of Colorado, Boulder, CO 80309; bDepartment of Physics, University of Colorado, Boulder, CO 80309; cTheoretical Physical Chemistry, Unité de Recherche Molecular Systems (UR MOLSYS) University of Liège, B4000 Liège, Belgium; dDepartment of Chemistry, Wayne State University, Detroit, MI 48202; and eDepartment of Chemistry and Biochemistry, Georgia Southern University, Statesboro, GA 30460 Contributed by Margaret M. Murnane, October 14, 2017 (sent for review July 16, 2017; reviewed by Elliot Bernstein and Arthur Suits) Highly excited electronic states are challenging to explore exper- When molecules are excited by vacuum-UV (VUV) light in imentally and theoretically—due to the large density of states and the 7- to 10-eV range, just below the ionization potential, one or the fact that small structural changes lead to large changes in several excited states can be accessed depending on the band- electronic character with associated strong nonadiabatic dynamics. width of the pulse. In this energy range, there is typically a dense They can play a key role in astrophysical and ionospheric chemis- manifold of valence and Rydberg states that are often strongly try, as well as the detonation chemistry of high-energy density coupled, each with an associated vibrational structure (9). As a materials. Here, we implement ultrafast vacuum-UV (VUV)-driven result, VUV molecular absorption spectra typically contain electron–ion coincidence imaging spectroscopy to directly probe quasicontinuous but structured bands with overlapping sharp the reaction pathways of highly excited states of energetic mole- peaks, indicating the wide variety of timescales and states in- — cules in this case, methyl azide.
    [Show full text]
  • Ultrafast Radiationless Decay Mechanisms Through Conical Intersections in Cytosine. a New Semi-Planar Conical Intersection
    Indian Journal of Chemistry Vol. 53A, February 2014, pp. 143-151 Ultrafast radiationless decay mechanisms through conical intersections in cytosine. A new semi-planar conical intersection Saadullah G Aziza, Walid I Hassanb, Shaaban K Alroubya, c, Abdulrahman Alyoubia & Rifaat Hilala, b, * aChemistry Department, Faculty of Science, King Abdulaziz University, Jeddah, SA bChemistry Department, Faculty of Science Cairo University, Giza, Egypt cChemistry Department, Faculty of Science Beni-Suief University, Beni-Suief, Egypt Email: [email protected] Received 26 November 2013; revised and accepted 27 January 2014 Geometry, energetics and dipole moment of all possible conformers of cytosine in the ground state are calculated using density functional theory B3LYP method and the 6-311++G(3df,3pd) basis set. The most stable conformer is keto-amino conformer. The amino group hydrogen atoms are slightly out of plane by about 6.3° and 9.9° each. Ultrafast radiationless decay mechanism has been theoretically investigated using Complete Active Space Multiconfiguration SCF calculations. Effective pathways of ultrafast radiationless transitions from the optically allowed π π* state to the ground state S0 of cytosine are explored. The nπ*, nσ*, and the ππ* states have been taken into account as states involved in the radiationless process. Optimized geometry and conical intersections have been searched in the full dimensional space for the vibrational degrees of freedom. Three competing direct decay mechanisms through three possible conical intersections have been found to exist. The first pathway is through the bending of molecule in a sofa-like structure leading to conical intersection with ground state at 4.23 eV. The second pathway occurs through a twisted structure that has hydrogen twisted and the cytosine ring slightly deformed leading to conical intersection at 4.08 eV.
    [Show full text]