Optimizing conical intersections without derivative coupling vectors: application to multistate multireference second-order perturbation theory (MS-CASPT2)

Electron-driven proton transfer enables nonradiative photodeactivation in microhydrated 2-aminoimidazole

Prebiotically credible activator of non-enzymatic RNA template-copying, 2-aminoimidazole, is protected from destructive photochemistry by photoacidity.

How Inter- and Intramolecular Processes Dictate Aggregation-Induced Emission in Crystals Undergoing Excited-State Proton Transfer

Aggregation-induced emission (AIE) offers a route for the development of luminescent technologies with high quantum efficiencies. Excited-state intramolecular proton transfer (ESIPT) coupled to AIE can produce devices with emission across the visible spectrum. We use a combination of theoretical models to determine the factors that mediate fluorescence in molecular crystals undergoing ESIPT. Using two materials based on 2'-hydroxychalcone as exemplar cases, we analyze how inter- and intramolecular processes determine the emissive properties in the crystal environment. This systematic investigation extends the current interpretation of AIE to polar chromophores with multiple decay pathways. We find that population of nonradiative pathways is dictated by the electronic effects of the substituents and the degree of distortion allowed in the crystal environment. Localization of the electron density is crucial to maximize fluorescence via ESIPT. Our conclusions offer design strategies for the development of luminescent molecular crystals.

Resolving the Notorious Case of Conical Intersections for Coupled Cluster Dynamics

The motion of electrons and nuclei in photochemical events often involves conical intersections, or degeneracies between electronic states. They serve as funnels in nuclear relaxation processes where the electrons and nuclei couple nonadiabatically. Accurate ab initio quantum chemical models are essential for interpreting experimental measurements of such phenomena. In this Letter, we resolve a long-standing problem in coupled cluster theory, presenting the first formulation of the theory that correctly describes conical intersections between excited electronic states of the same symmetry. This new development demonstrates that the highly accurate coupled cluster theory can be applied to describe dynamics on excited electronic states involving conical intersections.

α-CASSCF: An Efficient, Empirical Correction for SA-CASSCF To Closely Approximate MS-CASPT2 Potential Energy Surfaces

Because of its computational efficiency, the state-averaged complete active-space self-consistent field (SA-CASSCF) method is commonly employed in nonadiabatic ab initio molecular dynamics. However, SA-CASSCF does not effectively recover dynamical correlation. As a result, there can be qualitative differences between SA-CASSCF potential energy surfaces (PESs) and more accurate reference surfaces computed using multistate complete active space second-order perturbation theory (MS-CASPT2). Here we introduce an empirical correction to SA-CASSCF that scales the splitting between individual states and the state-averaged energy. We call this the α-CASSCF method, and we show here that it significantly improves the accuracy of relative energies and PESs compared with MS-CASPT2 for the chromophores of green fluorescent and photoactive yellow proteins. As such, this method may prove to be quite valuable for nonadiabatic dynamics.

Dual-Functional Tamm-Dancoff Approximation: A Convenient Density Functional Method that Correctly Describes S1/S0 Conical Intersections

Time-dependent Kohn-Sham density functional theory has been used successfully to compute vertical excitation energies, especially for large molecular systems. However, the lack of double excitation character in the excited amplitudes produced by linear response in the adiabatic approximation holds it back from broader applications in photochemistry; for example, it shows (3N - 7)-dimensional conical intersection seams (where N is the number of atoms) between ground and excited states, although the correct dimensionality is 3N - 8. In this letter, we present a new, conceptually simple, easy-to-implement, and easy-to-use way to employ time-dependent Kohn-Sham density functional theory that has global accuracy comparable with the conventional single-functional version and that recovers the double cone topology of the potential energy surfaces at S₁/S₀ conical intersection seams. The new method is called the dual-functional Tamm-Dancoff approximation (DF-TDA).

Insights into the deactivation of 5-bromouracil after ultraviolet excitation

5-Bromouracil is a nucleobase analogue that can replace thymine in DNA strands and acts as a strong radiosensitizer, with potential applications in molecular biology and cancer therapy. Here, the deactivation of 5-bromouracil after ultraviolet irradiation is investigated in the singlet and triplet manifold by accurate quantum chemistry calculations and non-adiabatic dynamics simulations. It is found that, after irradiation to the bright ππ* state, three main relaxation pathways are, in principle, possible: relaxation back to the ground state, intersystem crossing (ISC) and C-Br photodissociation. Based on accurate MS-CASPT2 optimizations, we propose that ground-state relaxation should be the predominant deactivation pathway in the gas phase. We then employ different electronic structure methods to assess their suitability to carry out excited-state dynamics simulations. MRCIS (multi-reference configuration interaction including single excitations) was used in surface hopping simulations to compute the ultrafast ISC dynamics, which mostly involves the 1nOπ* and 3ππ* states.This article is part of the themed issue 'Theoretical and computational studies of non-equilibrium and non-statistical dynamics in the gas phase, in the condensed phase and at interfaces'.

Nonradiative Relaxation Mechanisms of UV Excited Phenylalanine Residues: A Comparative Computational Study

The present work is directed toward understanding the mechanisms of excited state deactivation in three neutral model peptides containing the phenylalanine residue. The excited state dynamics of theγL(g+)folded form of N-acetylphenylalaninylamide (NAPA B) and its amide-N-methylated derivative (NAPMA B) is reviewed and compared to the dynamics of the monohydrated structure of NAPA (NAPAH). The goal is to unravel how the environment, and in particular solvation, impacts the photodynamics of peptides. The systems are investigated using reaction path calculations and surface hopping nonadiabatic dynamics based on the coupled cluster doubles (CC2) method and time-dependent density functional theory. The work emphasizes the role that excitation transfer from the phenylππ*to amidenπ*state plays in the deactivation of the three systems and shows how the ease of out-of-plane distortions of the amide group determines the rate of population transfer between the two electronic states. The subsequent dynamics on thenπ*state is barrierless along several pathways and leads to fast deactivation to the ground electronic state.

QM/MM Study of Substituent and Solvent Effects on the Excited State Dynamics of the Photoactive Yellow Protein Chromophore

Substituent and solvent effects on the excited state dynamics of the Photoactive Yellow Protein chromophore are studied using the average solvent electrostatic potential from molecular dynamics (ASEP/MD) method. Four molecular models were considered: the ester and thioester derivatives of the p-coumaric acid anion and their methylated derivatives. We found that the solvent produces dramatic modifications on the free energy profile of the S1 state: 1) Two twisted structures that are minima in the gas phase could not be located in aqueous solution. 2) Conical intersections (CIs) associated with the rotation of the single bond adjacent to the phenyl group are found for the four derivatives in water solution but only for thio derivatives in the gas phase. 3) The relative stability of minima and CIs is reverted with respect to the gas phase values, affecting the prevalent de-excitation paths. As a consequence of these changes, three competitive de-excitation channels are open in aqueous solution: the fluorescence emission from a planar minimum on S1, the trans-cis photoisomerization through a CI that involves the rotation of the vinyl double bond, and the nonradiative, nonreactive, de-excitation through the CI associated with the rotation of the single bond adjacent to the phenyl group. In the gas phase, the minima are the structures with the lower energy, while in solution these are the conical intersections. In solution, the de-excitation prevalent path seems to be the photoisomerization for oxo compounds, while thio compounds return to the initial trans ground state without emission.

On the decay of the triplet state of thionucleobases

The double-well triplet state of thionucleobases allows for a two-step mechanistic control of their triplet decay lifetime.

Mechanism of ultrafast non-reactive deactivation of the retinal chromophore in non-polar solvents

Counterion sensitive photodynamics of the retinal chromophore in solution.

Ab initio molecular dynamics relaxation and intersystem crossing mechanisms of 5-azacytosine

Nonadiabatic SHARC dynamics simulations reveal the molecular deformations involved in the photodeactivation pathways of 5-azacytosine.

Estimating the Derivative Coupling Vector Using Gradients

The difference gradient and derivative coupling vectors span the branching planes of conical intersections between electronic states. While gradients are commonly available in many electronic structure methods, the derivative coupling vectors are not always implemented and ready for use in characterizing conical intersections. This Letter shows how the derivative coupling vectors can be computed to high accuracy (direction and magnitude) using energy and gradient information. The new method is based on the combination of a linear-coupling two-state Hamiltonian and a finite-difference Davidson approach for computing the branching plane. Benchmark cases are provided showing these vectors can be efficiently computed near conical intersections.

New Insights into the State Trapping of UV-Excited Thymine

After UV excitation, gas phase thymine returns to a ground state in 5 to 7 ps, showing multiple time constants. There is no consensus on the assignment of these processes, with a dispute between models claiming that thymine is trapped either in the first (S₁) or in the second (S₂) excited states. In the present study, a nonadiabatic dynamics simulation of thymine is performed on the basis of ADC(2) surfaces, to understand the role of dynamic electron correlation on the deactivation pathways. The results show that trapping in S₂ is strongly reduced in comparison to previous simulations considering only non-dynamic electron correlation on CASSCF surfaces. The reason for the difference is traced back to the energetic cost for formation of a CO π bond in S₂.

The origin of efficient triplet state population in sulfur-substituted nucleobases

Elucidating the photophysical mechanisms in sulfur-substituted nucleobases (thiobases) is essential for designing prospective drugs for photo- and chemotherapeutic applications. Although it has long been established that the phototherapeutic activity of thiobases is intimately linked to efficient intersystem crossing into reactive triplet states, the molecular factors underlying this efficiency are poorly understood. Herein we combine femtosecond transient absorption experiments with quantum chemistry and nonadiabatic dynamics simulations to investigate 2-thiocytosine as a necessary step to unravel the electronic and structural elements that lead to ultrafast and near-unity triplet-state population in thiobases in general. We show that different parts of the potential energy surfaces are stabilized to different extents via thionation, quenching the intrinsic photostability of canonical DNA and RNA nucleobases. These findings satisfactorily explain why thiobases exhibit the fastest intersystem crossing lifetimes measured to date among bio-organic molecules and have near-unity triplet yields, whereas the triplet yields of canonical nucleobases are nearly zero.

Electronic structure of the para-benzoquinone radical anion revisited

Photoinduced processes in para-benzoquinone anion are studied with multistate multireference perturbation theory: an interplay between autodetachment and internal conversion.

Solvent effects on de-excitation channels in the p-coumaric acid methyl ester anion, an analogue of the photoactive yellow protein (PYP) chromophore

Environmental effects on the deactivation channels of the PYP chromophore in the gas phase and water solution are compared at the CASPT2//CASSCF/cc-pVDZ level.

Photorelaxation of imidazole and adenine via electron-driven proton transfer along H2O wires

Photochemically created πσ* states were classified among the most prominent factors determining the ultrafast radiationless deactivation and photostability of many biomolecular building blocks. In the past two decades, the gas phase photochemistry of πσ* excitations was extensively investigated and was attributed to N–H and O–H bond fission processes. However, complete understanding of the complex photorelaxation pathways of πσ* states in the aqueous environment was very challenging, owing to the direct participation of solvent molecules in the excited-state deactivation. Here, we present non-adiabatic molecular dynamics simulations and potential energy surface calculations of the photoexcited imidazole–(H₂O)₅ cluster using the algebraic diagrammatic construction method to the second-order [ADC(2)]. We show that electron driven proton transfer (EDPT) along a wire of at least two water molecules may lead to the formation of a πσ*/S₀ state crossing, similarly to what we suggested for 2-aminooxazole. We expand on our previous findings by direct comparison of the imidazole–(H₂O)₅ cluster to non-adiabatic molecular dynamics simulations of imidazole in the gas phase, which reveal that the presence of water molecules extends the overall excited-state lifetime of the chromophore. To embed the results in a biological context, we provide calculations of potential energy surface cuts for the analogous photorelaxation mechanism present in adenine, which contains an imidazole ring in its structure.

Accurate nonadiabatic dynamics

This Perspective addresses the use of coupled diabatic potential energy surfaces (PESs) together with rigorous quantum dynamics in full or reduced dimensional coordinate spaces to obtain accurate solutions to problems in nonadiabatic dynamics.

A DFT and multi-configurational perturbation theory study on O2 binding to a model heme compound via the spin-change barrier

Potential energy surface of the O₂ binding to a model heme compound via minimum energy intersystem crossing point (MEISCP) was investigated with DFT and MS-CASPT2. The porphyrin's symmetric shrinking vibration mode contributes to reach the MESICP from the T₁ minimum.

Excited-state deactivation in 8-oxo-deoxyguanosine: comparison between anionic and neutral forms

Ab initio explorations of excited-state potential-energy surfaces show that a radiationless deactivation mechanism via intramolecular excited-state proton transfer is available in neutral 8-oxo-deoxyguanosine, whereas it is not available in the anionic form.

Assessment of Approximate Coupled-Cluster and Algebraic-Diagrammatic-Construction Methods for Ground- and Excited-State Reaction Paths and the Conical-Intersection Seam of a Retinal-Chromophore Model

As a minimal model of the chromophore of rhodopsin proteins, the penta-2,4-dieniminium cation (PSB3) poses a challenging test system for the assessment of electronic-structure methods for the exploration of ground- and excited-state potential-energy surfaces, the topography of conical intersections, and the dimensionality (topology) of the branching space. Herein, we report on the performance of the approximate linear-response coupled-cluster method of second order (CC2) and the algebraic-diagrammatic-construction scheme of the polarization propagator of second and third orders (ADC(2) and ADC(3)). For the ADC(2) method, we considered both the strict and extended variants (ADC(2)-s and ADC(2)-x). For both CC2 and ADC methods, we also tested the spin-component-scaled (SCS) and spin-opposite-scaled (SOS) variants. We have explored several ground- and excited-state reaction paths, a circular path centered around the S1/S0 surface crossing, and a 2D scan of the potential-energy surfaces along the branching space. We find that the CC2 and ADC methods yield a different dimensionality of the intersection space. While the ADC methods yield a linear intersection topology, we find a conical intersection topology for the CC2 method. We present computational evidence showing that the linear-response CC2 method yields a surface crossing between the reference state and the first response state featuring characteristics that are expected for a true conical intersection. Finally, we test the performance of these methods for the approximate geometry optimization of the S1/S0 minimum-energy conical intersection and compare the geometries with available data from multireference methods. The present study provides new insight into the performance of linear-response CC2 and polarization-propagator ADC methods for molecular electronic spectroscopy and applications in computational photochemistry.

Ultrafast Charge Transfer of a Valence Double Hole in Glycine Driven Exclusively by Nuclear Motion

We explore theoretically the ultrafast transfer of a double electron hole between the functional groups of glycine after K-shell ionization and subsequent Auger decay. Although a large energy gap of about 15 eV initially exists between the two electronic states involved and coherent electronic dynamics play no role in the hole transfer, we find that the double hole is transferred within 3 to 4 fs between both functional ends of the glycine molecule driven solely by specific nuclear displacements and non-Born-Oppenheimer effects. The nuclear displacements along specific vibrational modes are of the order of 15% of a typical chemical bond between carbon, oxygen, and nitrogen atoms and about 30% for bonds involving hydrogen atoms. The time required for the hole transfer corresponds to less than half a vibrational period of the involved nuclear modes. This finding challenges the common wisdom that nuclear dynamics of the molecular skeleton are unimportant for charge transfer processes at the few-femtosecond time scale and shows that they can even play a prominent role. It also indicates that in x-ray imaging experiments, in which ionization is unavoidable, valence electron redistribution caused by nuclear dynamics might be much faster than previously anticipated. Thus, non-Born-Oppenheimer effects may affect the apparent electron densities extracted from such measurements.

Energy dissipative photoprotective mechanism of carotenoid spheroidene from the photoreaction center of purple bacteria Rhodobacter sphaeroides

Carotenoid spheroidene (SPO) functions for photoprotection in the photosynthetic reaction centers (RCs) and effectively dissipates its triplet excitation energy. Sensitized cis-to-trans isomerization was proposed as a possible mechanism for a singlet-triplet energy crossing for the 15,15'-cis-SPO; however, it has been questioned recently. To understand the dissipative photoprotective mechanism of this important SPO and to overcome the existing controversies on this issue, we carried out a theoretical investigation using density functional theory on the possible triplet energy relaxation mechanism through the cis-to-trans isomerization. Together with the earlier experimental observations, the possible mechanism was discussed for the triplet energy relaxation of the 15,15'-cis-SPO. The result shows that complete cis-to-trans isomerization is not necessary. Twisting the C15-C15' bond leads to singlet-triplet energy crossing at ϕ(14,15,15',14') = 77° with an energy 32.5 kJ mol(-1) (7.7 kcal mol(-1)) higher than that of the T1 15,15'-cis minimum. Further exploration of the minimum-energy intersystem crossing (MEISC) point shows that triplet relaxation could occur at a less distorted structure (ϕ = 58.4°) with the energy height of 26.5 KJ mol(-1) (6.3 kcal mol(-1)). Another important reaction coordinate to reach the MEISC point is the bond-length alternation. The model truncation effect, solvent effect, and spin-orbit coupling were also investigated. The singlet-triplet crossing was also investigated for the 13,14-cis stereoisomer and locked-13,14-cis-SPO. We also discussed the origin of the natural selection of the cis over trans isomer in the RC.

Defect-Induced Conical Intersections Promote Nonradiative Recombination

We apply multireference electronic structure calculations to demonstrate the presence of conical intersections between the ground and the first excited electronic states of three silicon nanocrystals containing defects characteristic of the oxidized silicon surface. These intersections are accessible upon excitation at visible wavelengths and are predicted to facilitate nonradiative recombination with a rate that increases with decreasing particle size. This work illustrates a new framework for identifying defects responsible for nonradiative recombination.