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

The copious photochemistry of 2,6-diaminopurine: Luminescence, triplet population, and ground state recovery

9H- and 7H-2,6-Diaminopurine (26DAP) photoinduced events in vacuum were studied at the MS-CASPT2/cc-pVDZ level of theory. The S1 1 (ππ* La ) state is initially populated evolving barrierless towards its minimum energy structure, from where two photochemical events can take place in both tautomers. The first is the return of the electronic population to the ground state via the C6 conical intersection (CI-C6). The second involves an internal conversion to the ground through the C2 conical intersection (CI-C2). According to our geodesic interpolated paths connecting the critical structures, the second route is less favorable in both tautomers, due to the presence of high energy barriers. Our calculations suggest a competition between fluorescence and ultrafast relaxation to the electronic ground state via internal conversion process. Based on our calculated potential energy surfaces and experimental excited state lifetimes from the literature, we can infer that the 7H- must have a greater fluorescence yield than the 9H-tautomer. We also explored the triplet state population mechanisms on the 7H-26DAP to understand their long-lived components observed experimentally.

The tautomer-specific excited state dynamics of 2,6-diaminopurine using resonance-enhanced multiphoton ionization and quantum chemical calculations

2,6-Diaminopurine (2,6-dAP) is an alternative nucleobase that potentially played a role in prebiotic chemistry. We studied its excited state dynamics in the gas phase by REMPI, IR-UV hole burning, and ps pump-probe spectroscopy and performed quantum chemical calculations at the SCS-ADC(2) level of theory to interpret the experimental results. We found the 9H tautomer to have a small barrier to ultrafast relaxation via puckering of its 6-membered ring. The 7H tautomer has a larger barrier to reach a conical intersection and also has a sizable triplet yield. These results are discussed relative to other purines, for which 9H tautomerization appears to be more photostable than 7H and homosubstituted purines appear to be less photostable than heterosubstituted or singly substituted purines.

Which Electronic Structure Method to Choose in Trajectory Surface Hopping Dynamics Simulations? Azomethane as a Case Study

Nonadiabatic dynamics simulations have become a standard approach to explore photochemical reactions. Such simulations require underlying potential energy surfaces and couplings between them, calculated at a chosen level of theory, yet this aspect is rarely assessed. Here, in combination with the popular trajectory surface hopping dynamics method, we use a high-accuracy XMS-CASPT2 electronic structure level as a benchmark for assessing the performances of various post-Hartree-Fock methods (namely, CIS, ADC(2), CC2, and CASSCF) and exchange-correlation functionals (PBE, PBE0, and CAM-B3LYP) in a TD-DFT/TDA context, using the isomerization around a double bond as test case. Different relaxation pathways are identified, and the ability of the different methods to reproduce their relative importance and time scale is discussed. The results show that multireference electronic structure methods should be preferred, when studying nonadiabatic decay between excited and ground states. If not affordable, TD-DFT with TDA and hybrid functionals and ADC(2) are efficient alternatives but overestimate the nonradiative decay yield and thus may miss deexcitation pathways.

On the description of conical intersections between excited electronic states with LR-TDDFT and ADC(2)

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.

Assessment of the Electron Correlation Treatment on the Quantum-Classical Dynamics of Retinal Protonated Schiff Base Models: XMS-CASPT2, RMS-CASPT2, and REKS Methods

We compare the performance of three different multiconfigurational wave function-based electronic structure methods and two implementations of the spin-restricted ensemble-referenced Kohn-Sham (REKS) method. The study is characterized by three features: (i) it uses a small set of quantum-classical trajectories rather than potential energy surface mapping, (ii) it focuses, exclusively, on the photoisomerization of retinal protonated Schiff base models, and (iii) it probes the effect of both methyl substitution and the increase in length of the conjugate π-system. For each tested method, the corresponding analytical gradients are used to drive the quantum-classical (Tully's FSSH method) trajectory propagation, including the recent multistate XMS-CASPT2 and RMS-CASPT2 gradients. It is shown that while CASSCF, XMS-CASPT2, and RMS-CASPT2 yield consistent photoisomerization dynamics descriptions, REKS produces, in some of these systems, qualitatively different behavior that is attributed to a flatter and topographically different excited state potential energy surface. The origin of this behavior can be traced back to the effect of the employed density functional approximation. The above studies are further expanded by benchmarking, at the CASSCF and REKS levels, the electronic structure methods using a QM/MM model of the visual pigment rhodopsin.

Excited-State Dynamics Simulations of a Light-Driven Molecular Motor in Solution

Molecular motors, where light can be transformed into motion, are promising in the design of nanomechanical devices. For applications, however, finding relationships between molecular motion and the environment is important. Here, we report the study of excited-state dynamics of an overcrowded alkene in solution using a hybrid quantum mechanics/molecular mechanics (QM/MM) approach combined with excited-state molecular dynamics simulations. Using QM/MM surface-hopping trajectories, we calculated time-resolved emission and transient absorption spectra. These show the rise of a short-lived Franck-Condon state, followed by the formation of a dark state in the first 150 fs before the molecular motor relaxes to the ground state in about 1 ps. From the analysis of radial distribution functions, we infer that the orientation of the solvent with respect to the molecular motor in the electronic excited state is similar to that in the ground state during the photoisomerization.

Doubly Tuned Exchange-Correlation Functionals for Mixed-Reference Spin-Flip Time-Dependent Density Functional Theory

It is demonstrated that significant accuracy improvements in MRSF-TDDFT can be achieved by introducing two different exchange-correlation (XC) functionals for the reference Kohn-Sham DFT and the response part of the calculations, respectively. Accordingly, two new XC functionals of doubly tuned Coulomb attenuated method-vertical excitation energy (DTCAM-VEE) and DTCAM-AEE were developed on the basis of the "adaptive exact exchange (AEE)" concept in the framework of the Coulomb-attenuating XC functionals. The values by DTCAM-VEE are in excellent agreement with those of Thiel's set [mean absolute errors (MAEs) and the interquartile range (IQR) values of 0.218 and 0.327 eV, respectively]. On the other hand, DTCAM-AEE faithfully reproduced the qualitative aspects of conical intersections (CIs) of trans-butadiene and thymine and the nonadiabatic molecular dynamics (NAMD) simulations on thymine. The latter functional also remarkably exhibited the exact 1/R asymptotic behavior of the charge-transfer state of an ethylene-tetrafluoroethylene dimer and the accurate potential energy surfaces (PESs) along the two torsional angles of retinal protonated Schiff base model with six double bonds (rPSB6). Overall, DTCAM-AEE generally performs well, as its MAE (0.237) and IQR (0.41 eV) are much improved as compared to BH&HLYP. The current idea can also be applied to other XC functionals as well as other variants of linear response theories, opening a new way of developing XC functionals.

Mixed-Reference Spin-Flip Time-Dependent Density Functional Theory: Multireference Advantages with the Practicality of Linear Response Theory

The density functional theory (DFT) and linear response (LR) time-dependent (TD)-DFT are of the utmost importance for routine computations. However, the single reference formulation of DFT suffers in the description of open-shell singlet systems such as diradicals and bond-breaking. LR-TDDFT, on the other hand, finds difficulties in the modeling of conical intersections, doubly excited states, and core-level excitations. In this Perspective, we demonstrate that many of these limitations can be overcome by recently developed mixed-reference (MR) spin-flip (SF)-TDDFT, providing an alternative yet accurate route for such challenging situations. Empowered by the practicality of the LR formalism, it is anticipated that MRSF-TDDFT can become one of the major workhorses for general routine tasks.

A theoretical study on the photochemical generation of phenylborylene from phenyldiazidoborane

Organic borylenes are a kind of highly reactive species, which play important roles in a lot of reactions as vigorous intermediates. In this work, we investigated the photochemical generation mechanisms of phenylborylene (PhB) together with the side product N-phenylnitrenoiminoborane (PhNBN) from phenyldiazidoborane (PhBN₆) by extrusion of dinitrogen in the two lowest electronic singlet states (S₀ and S₁) based on the complete active space self-consistent field (CASSCF) and its second-order perturbation (CASPT2) methods combined with time-dependent density functional theory (TD-DFT) calculations. Our results show that the reaction PhBN₆ → PhB + 3N₂ involves stepwise N₂ extrusion three times and the azido region rearrangement. Moreover, we found that the studied photo-induced processes are kinetically feasible because the highest energy barrier is only 0.36 eV and excitation with light of wavelength 254 nm can provide enough excess energy to overcome these energy barriers. Importantly, we revealed that several conical intersections between S₁ and S₀ states participate and facilitate the studied photochemical processes. Our results not only clarify the experimental observations, (H. F. Bettinger, J. Am. Chem. Soc. 2006, 128, 2534), but also provide valuable insights into borylene chemistry.

Ultrafast Intersystem Crossing Dynamics of 6-Selenoguanine in Water

Rationalizing the photochemistry of nucleobases where an oxygen is replaced by a heavier atom is essential for applications that exploit near-unity triplet quantum yields. Herein, we report on the ultrafast excited-state deactivation mechanism of 6-selenoguanine (6SeGua) in water by combining nonadiabatic trajectory surface-hopping dynamics with an electrostatic embedding quantum mechanics/molecular mechanics (QM/MM) scheme. We find that the predominant relaxation mechanism after irradiation starts on the bright singlet S₂ state that converts internally to the dark S₁ state, from which the population is transferred to the triplet T₂ state via intersystem crossing and finally to the lowest T₁ state. This S₂ → S₁ → T₂ → T₁ deactivation pathway is similar to that observed for the lighter 6-thioguanine (6tGua) analogue, but counterintuitively, the T₁ lifetime of the heavier 6SeGua is shorter than that of 6tGua. This fact is explained by the smaller activation barrier to reach the T₁/S₀ crossing point and the larger spin-orbit couplings of 6SeGua compared to 6tGua. From the dynamical simulations, we also calculate transient absorption spectra (TAS), which provide two time constants (τ₁ = 131 fs and τ₂ = 191 fs) that are in excellent agreement with the experimentally reported value (τ_exp = 130 ± 50 fs) (Farrel et al. J. Am. Chem. Soc. 2018, 140, 11214). Intersystem crossing itself is calculated to occur with a time scale of 452 ± 38 fs, highlighting that the TAS is the result of a complex average of signals coming from different nonradiative processes and not intersystem crossing alone.

On the propensity of formation of cyclobutane dimers in face-to-face and face-to-back uracil stacks in solution

UV irradiation of RNA leads to the formation of intra- and inter-strand crosslinks of cyclobutane type. Despite the importance of this reaction, relatively little is known about how the mutual orientation of the two bases affects the outcome of the reaction. Here we report a comparative nonadiabatic molecular dynamics study of face-to-back (F2B) and face-to-face (F2F) stacked uracil-water clusters. The computations were performed using the second-order algebraic-diagrammatic-construction (ADC(2)) method. We found that F2B stacked uracil-water clusters either relax non-reactively to the ground state by an ethylenic twist around the CC bond or remain in the lowest nπ* state in which the two bases gradually move away from each other. This finding is consistent with the low propensity for the formation of intra-strand cyclobutane dimers between adjacent RNA bases. On the contrary, in F2F stacked uracil-water clusters, in addition to non-reactive deactivation, we found a pro-reactive deactivation pathway, which may lead to the formation of cyclobutane uracil dimers in the electronic ground state. On a qualitative level, the observed photodynamics of F2F stacked uracil-water clusters explains the greater propensity of RNA to form inter-strand cyclobutane-type crosslinks.

Ab initio trajectory surface-hopping dynamics studies of excited-state proton-coupled electron transfer reactions in trianisoleheptazine-phenol complexes

The excited-state proton-coupled electron-transfer (PCET) reaction in hydrogen-bonded complexes of trianisoleheptazine (TAHz), a chromophore related to polymeric carbon nitrides widely used in hydrogen-evolution photocatalysis, with several phenol derivatives were recently studied by Schlenker and coworkers with time-resolved photoluminescence quenching and pump-probe experiments. A pronounced dependence of the PCET reactivity on the electron-donating/electron-withdrawing character of the substituents on phenol was found, with indications of a barrierless or nearly barrierless PCET reaction for the most strongly electron-donating substituent, methoxy. In the present work, the excited-state PCET dynamics was explored with first-principles nonadiabatic dynamics simulations using the TDDFT/ωB97X-D electronic-structure model for two selected complexes, TAHz-phenol and TAHz-methoxyphenol. The qualitative reliability of the TDDFT/ωB97X-D electronic-structure model was assessed by extensive benchmarking of excitation energies and potential-energy profiles against a wave-function-based ab initio method, the algebraic-diagrammatic construction of second order (ADC(2)). The nonadiabatic dynamics simulations provide temporally and structurally resolved insights into paradigmatic PCET reactions in TAHz-phenol complexes. The radiationless relaxation of the photoexcited bright ¹ππ* state to the long-lived dark S₁ state of TAHz occurs in less than 100 fs. The ensuing PCET reaction on the adiabatic S₁ surface is faster in TAHz-methoxyphenol complexes than in TAHz-phenol complexes due to a lower H-atom-transfer barrier, as observed in the experiments. The relaxation of the complexes to the electronic ground state is found to occur exclusively via PCET within the 250 fs time window covered by the present simulations, confirming the essential role of the hydrogen bond for the fluorescence quenching process. The absolute values of the computed PCET time constants are significantly shorter than those extracted from time-resolved photoluminescence measurements for mixtures of TAHz with phenolic substrates in toluene. The possible origins of this discrepancy are discussed.

Theoretical insights into photo-induced isomerization mechanisms of phenylsulfinyl radical PhSO˙

Sulfinyl radicals (R-SO˙) play important roles in lots of reactions, while the isomer oxathiyl radicals (R-OS˙) and the isomerization between them are rarely observed due to the poor stability of R-OS˙. In this work, the complete active space self-consistent field (CASSCF) and its multi-state second order perturbation (MS-CASPT2) methods were employed to study the photo-induced reaction mechanisms of phenylsulfinyl radical PhSO˙ 1 and its isomer phenoxathiyl radical PhOS˙ 2. Our results show that 1 and 2 have similar singly occupied molecular orbitals in the ground state but different properties in the excited state, which determine their diverse behaviors after irradiation. Radical 1 can generate 2 by light irradiation, but 2 produces isomerization product 3 (2-hydroxyphenylthiyl radical) and ring-opening product 4 (acyclic thioketoketene radical) in two paths via S atom migration intermediate Int1 (2-carbonylcyclohexadienthiyl radical). The former path involves consequent hydrogen shift reactions with a strongly exothermic process while the latter path involves both ring-expansion and ring-opening processes with a high barrier, resulting in a structural and energetic preference for the former path. Moreover, we revealed several conical intersections that participate in the reactions and facilitate the photochemical processes. Our calculations not only remain consistent with and clarify the experimental observations (X. Zeng, et al., J. Am. Chem. Soc., 2018, 140(31), 9972-9978) but also enrich the knowledge of sulfinyl radicals and isomer oxathiyl radicals.

Ab Initio Quasiclassical Simulation of Femtosecond Time-Resolved Two-Dimensional Electronic Spectra of Pyrazine

Two-dimensional (2D) electronic spectroscopy is a powerful nonlinear technique which provides spectroscopic information on two frequency axes as well as dynamical information as a function of the so-called waiting time. Herein, an ab initio theoretical framework for the simulation of electronic 2D spectra has been developed. The method is based on the classical approximation to the doorway-window representation of three-pulse photon-echo signals and the description of nuclear motion by classical trajectories. Nonadiabatic effects are taken into account by a trajectory surface-hopping algorithm. 2D electronic spectra were simulated with ab initio on-the-fly trajectory calculations using the ADC(2) electronic-structure method for the pyrazine molecule, which is a benchmark system for ultrafast radiationless decay through conical intersections. It is demonstrated that 2D spectroscopy with subfemtosecond UV pulses can provide unprecedented detailed information on the ultrafast photodynamics of polyatomic molecules.

Nonradiative relaxation mechanisms of the elusive silole molecule

Silole derivatives have been extensively employed for developing organic optoelectronics, but few studies focused on the photophysical properties of the silole molecule. In this work, we investigate these properties by computing the absorption spectra and performing nonadiabatic molecular dynamics of silole employing the algebraic diagrammatic construction [ADC(2)] and extended multi-state XMS-CASPT2 ab initio electronic structure methods. For vertical excitations and excited state optimizations, the equation of motion coupled-cluster singles and doubles (EOM-CCSD) was also used. The nuclear ensemble and the fewest-switches surface hopping molecular dynamics methods coupled with the first two high-level electronic structure methods were applied to probe the relaxation mechanisms of silole. We could reproduce the experimental first absorption maximum value and found an ultrafast relaxation process occurring exclusively through ring-puckering distortions without breaking ring bonds or hydrogen elimination. Minimum energy conical intersection optimizations were carried out and potential energy curves, including triplet states, were calculated to further elucidate the relaxation process of silole.

Simulation and Analysis of the Transient Absorption Spectrum of 4-(N,N-Dimethylamino)benzonitrile (DMABN) in Acetonitrile

4-(N,N-Dimethylamino)benzonitrile (DMABN) is a well-known model compound for dual fluorescence-in sufficiently polar solvents, it exhibits two distinct fluorescence emission bands. The interpretation of its transient absorption (TA) spectrum in the visible range is the subject of a long-standing controversy. In the present study, we resolve this issue by calculating the TA spectrum on the basis of nonadiabatic molecular dynamics simulations. An unambiguous assignment of spectral signals to specific excited-state structures is achieved by breaking down the calculated spectrum into contributions from twisted and nontwisted molecular geometries. In particular, the much-discussed excited-state absorption band near 1.7 eV (ca. 700 nm) is attributed to the near-planar locally excited (LE) minimum on the S1 state. On the technical side, our study demonstrates that the second-order approximate coupled cluster singles and doubles (CC2) method can be used successfully to calculate the TA spectra of moderately large organic molecules, provided that the system in question does not approach a crossing between the lowest excited state and the singlet ground state within the time frame of the simulation.

Resolving the ultrafast dynamics of the anionic green fluorescent protein chromophore in water

The chromophore of the green fluorescent protein (GFP) is critical for probing environmental influences on fluorescent protein behavior. Using the aqueous system as a bridge between the unconfined vacuum system and a constricting protein scaffold, we investigate the steric and electronic effects of the environment on the photodynamical behavior of the chromophore. Specifically, we apply ab initio multiple spawning to simulate five picoseconds of nonadiabatic dynamics after photoexcitation, resolving the excited-state pathways responsible for internal conversion in the aqueous chromophore. We identify an ultrafast pathway that proceeds through a short-lived (sub-picosecond) imidazolinone-twisted (I-twisted) species and a slower (several picoseconds) channel that proceeds through a long-lived phenolate-twisted (P-twisted) intermediate. The molecule navigates the non-equilibrium energy landscape via an aborted hula-twist-like motion toward the one-bond-flip dominated conical intersection seams, as opposed to following the pure one-bond-flip paths proposed by the excited-state equilibrium picture. We interpret our simulations in the context of time-resolved fluorescence experiments, which use short- and long-time components to describe the fluorescence decay of the aqueous GFP chromophore. Our results suggest that the longer time component is caused by an energetically uphill approach to the P-twisted intersection seam rather than an excited-state barrier to reach the twisted intramolecular charge-transfer species. Irrespective of the location of the nonadiabatic population events, the twisted intersection seams are inefficient at facilitating isomerization in aqueous solution. The disordered and homogeneous nature of the aqueous solvent environment facilitates non-selective stabilization with respect to I- and P-twisted species, offering an important foundation for understanding the consequences of selective stabilization in heterogeneous and rigid protein environments.

Modern Theoretical Approaches to Modeling the Excited-State Intramolecular Proton Transfer: An Overview

The excited-state intramolecular proton transfer (ESIPT) phenomenon is nowadays widely acknowledged to play a crucial role in many photobiological and photochemical processes. It is an extremely fast transformation, often taking place at sub-100 fs timescales. While its experimental characterization can be highly challenging, a rich manifold of theoretical approaches at different levels is nowadays available to support and guide experimental investigations. In this perspective, we summarize the state-of-the-art quantum-chemical methods, as well as molecular- and quantum-dynamics tools successfully applied in ESIPT process studies, focusing on a critical comparison of their specific properties.

Electronic Structure Methods for the Description of Nonadiabatic Effects and Conical Intersections

Nonadiabatic effects are ubiquitous in photophysics and photochemistry, and therefore, many theoretical developments have been made to properly describe them. Conical intersections are central in nonadiabatic processes, as they promote efficient and ultrafast nonadiabatic transitions between electronic states. A proper theoretical description requires developments in electronic structure and specifically in methods that describe conical intersections between states and nonadiabatic coupling terms. This review focuses on the electronic structure aspects of nonadiabatic processes. We discuss the requirements of electronic structure methods to describe conical intersections and nonadiabatic couplings, how the most common excited state methods perform in describing these effects, and what the recent developments are in expanding the methodology and implementing nonadiabatic couplings.

Caveat when using ADC(2) for studying the photochemistry of carbonyl-containing molecules

Several electronic-structure methods are available to study the photochemistry and photophysics of organic molecules. Among them, ADC(2) stands as a sweet spot between computational efficiency and accuracy. As a result, ADC(2) has recently seen its number of applications booming, in particular to unravel the deactivation pathways and photodynamics of organic molecules. Despite this growing success, we demonstrate here that care has to be taken when studying the nonradiative pathways of carbonyl-containing molecules, as ADC(2) appears to suffer from a systematic flaw.

Elucidating an Atmospheric Brown Carbon Species-Toward Supplanting Chemical Intuition with Exhaustive Enumeration and Machine Learning

Brown carbon (BrC) is involved in atmospheric light absorption and climate forcing and can cause adverse health effects. Understanding the formation mechanisms and molecular structure of BrC is of key importance in developing strategies to control its environment and health impact. Structure determination of BrC is challenging, due to the lack of experiments providing molecular fingerprints and the sheer number of molecular candidates with identical mass. Suggestions based on chemical intuition are prone to errors due to the inherent bias. We present an unbiased algorithm, using graph-based molecule generation and machine learning, which can identify all molecular structures of compounds involved in biomass burning and the composition of BrC. We apply this algorithm to C₁₂H₁₂O₇, a light-absorbing "test case" molecule identified in chamber experiments on the aqueous photo-oxidation of syringol, a prevalent marker in wood smoke. Of the 260 million molecular graphs, the algorithm leaves only 36,518 (0.01%) as viable candidates matching the spectrum. Although no unique molecular structure is obtained from only a chemical formula and a UV/vis absorption spectrum, we discuss further reduction strategies and their efficacy. With additional data, the method can potentially more rapidly identify isomers extracted from lab and field aerosol particles without introducing human bias.

2,6-diaminopurine promotes repair of DNA lesions under prebiotic conditions

High-yielding and selective prebiotic syntheses of RNA and DNA nucleotides involve UV irradiation to promote the key reaction steps and eradicate biologically irrelevant isomers. While these syntheses were likely enabled by UV-rich prebiotic environment, UV-induced formation of photodamages in polymeric nucleic acids, such as cyclobutane pyrimidine dimers (CPDs), remains the key unresolved issue for the origins of RNA and DNA on Earth. Here, we demonstrate that substitution of adenine with 2,6-diaminopurine enables repair of CPDs with yields reaching 92%. This substantial self-repairing activity originates from excellent electron donating properties of 2,6-diaminopurine in nucleic acid strands. We also show that the deoxyribonucleosides of 2,6-diaminopurine and adenine can be formed under the same prebiotic conditions. Considering that 2,6-diaminopurine was previously shown to increase the rate of nonenzymatic RNA replication, this nucleobase could have played critical roles in the formation of functional and photostable RNA/DNA oligomers in UV-rich prebiotic environments.

UV-induced photodamage that likely occurred during the prebiotic synthesis of DNA and RNA is still an untackled issue for their origin on early Earth. Here, the authors show that substitution of 2,6-diaminopurine for adenine enables repair of cyclobutane pyrimidine dimers with high yields, and demonstrate that both 2,6-diaminopurine and adenine nucleosides can be formed under the same prebiotic conditions.

Fast and Accurate Computation of Nonadiabatic Coupling Matrix Elements Using the Truncated Leibniz Formula and Mixed-Reference Spin-Flip Time-Dependent Density Functional Theory

We present a fast and accurate numerical algorithm for computing the first-order nonadiabatic coupling matrix element (NACME). The algorithm employs the truncated Leibniz formula (TLF) approximation within the finite-difference method, which makes it easily applicable in connection with any wave function-based methodology. In this work, we used the algorithm in connection with the recently developed mixed-reference spin-flip time-dependent density functional theory (MRSF-TDDFT, MRSF for brevity). The accuracy is assessed for NACME between the singlet electronic states of a dissociating hydrogen molecule. It is demonstrated that an intermediate approximation, TLF(1), affords a negligible numeric error on the order of ∼10^-10 a.u. while enabling a fast computation of NACME. As the MRSF method yields the correct description of the dissociation curves of H₂ for all the electronic states involved, the numeric TLF(1)/MRSF NACME values are in excellent agreement with the reference analytical values obtained by the full configuration interaction. For polyatomic molecules, the MRSF NAC vectors agree very closely with the MRCISD NAC vectors. Hence, the proposed protocol is a promising tool for the evaluation of NACMEs.

Impact of the Dynamic Electron Correlation on the Unusually Long Excited-State Lifetime of Thymine

Non-radiative relaxation of the photoexcited thymine in the gas phase shows an unusually long excited-state lifetime, and, over the years, a number of models, i.e., S₁-trapping, S₂-trapping, and S₁&S₂-trapping, have been put forward to explain its mechanism. Here, we investigate this mechanism using non-adiabatic molecular dynamics (NAMD) simulations in connection with the recently developed mixed-reference spin-flip time-dependent density functional theory (MRSF-TDDFT) method. We show that the previously predicted S₂-trapping model was due to an artifact caused by an insufficient account of the dynamic electron correlation. The current work supports the S₁-trapping mechanism with two lifetimes, τ₁ = 30 ± 1 fs and τ₂ = 6.1 ± 0.035 ps, quantitatively consistent with the recent time-resolved experiments. Upon excitation to the S₂ (ππ*) state, thymine undergoes an ultrafast (ca. 30 fs) S₂→S₁ internal conversion and resides around the minimum on the S₁ (n_Oπ*) surface, slowly decaying to the ground state (ca. 6.1 ps). While the S₂→S₁ internal conversion is mediated by fast bond length alternation distortion, the subsequent S₁→S₀ occurs through several conical intersections, involving a slow puckering motion.

Role of Triplet States in the Photodynamics of Aniline

The dynamics of excited heteroaromatic molecules is a key to understanding the photoprotective properties of many biologically relevant chromophores that dissipate their excitation energy nonreactively and thereby prevent the detrimental effects of ultraviolet radiation. Despite their structural variability, most substituted aromatic compounds share a common feature of a repulsive ¹πσ* potential energy surface. This surface can lead to photoproducts, and it can also facilitate the population transfer back to the ground electronic state by means of a ¹πσ*/S₀ conical intersection. Here, we explore a hidden relaxation route involving the triplet electronic state of aniline, which has recently been discovered by means of time-selected photofragment translational spectroscopy [J. Chem. Phys. 2019, 151, 141101]. By using the recently available analytical gradients for multiconfiguration pair-density functional theory, it is now possible to locate the minimum-energy crossing points between states of different spin and therefore compute the intersystem crossing rates with a multireference method, rather than with the less reliable single-reference methods. Using such calculations, we demonstrate that the population loss of aniline in the T₁(³ππ*) state is dominated by C₆H₅NH₂ → C₆H₅NH· + H· dissociation, and we explain the long nonradiative lifetimes of the T₁(³ππ*) state at the excitation wavelengths of 294-264 nm.

Role of Conical Intersections on the Efficiency of Fluorescent Organic Molecular Crystals

Organic molecular crystals are attractive materials for luminescent applications because of their promised tunability. However, the link between the chemical structure and emissive behavior is poorly understood because of the numerous interconnected factors which are at play in determining radiative and nonradiative behaviors at the solid-state level. In particular, the decay through conical intersection dominates the nonadiabatic regions of the potential energy surface, and thus, their accessibility is a telling indicator of the luminosity of the material. In this study, we investigate the radiative mechanism for five organic molecular crystals which display a solid-state emission, with a focus on the role of conical intersections in their photomechanisms. The objective is to situate the importance of the accessibility of conical intersections with regards to emissive behavior, taking into account other nonradiative decay channels, namely, vibrational decay, and exciton hopping. We begin by giving a brief overview of the structural patterns of the five systems within a larger pool of 13 crystals for a richer comparison. We observe that because of the prevalence of sheet like and herringbone packing in organic molecular crystals, the conformational diversity of crystal dimers is limited. Additionally, similarly spaced dimers have exciton coupling values of a similar order within a 50 meV interval. Next, we focus on three exemplary cases, where we disentangle the role of nonradiative decay mechanisms and show how rotational minimum energy conical intersections in vacuum lead to puckered ones in the crystal, increasing their instability upon crystallization in typical packing motifs. In contrast, molecules with puckered conical intersections in vacuum tend to conserve this trait upon crystallization, and therefore, their quantum yield of fluorescence is determined predominantly by other nonradiative decay mechanisms.

Machine Learning-Assisted Excited State Molecular Dynamics with the State-Interaction State-Averaged Spin-Restricted Ensemble-Referenced Kohn-Sham Approach

We present a machine learning-assisted excited state molecular dynamics (ML-ESMD) based on the ensemble density functional theory framework. Since we represent a diabatic Hamiltonian in terms of generalized valence bond ansatz within the state-interaction state-averaged spin-restricted ensemble-referenced Kohn-Sham (SI-SA-REKS) method, we can avoid singularities near conical intersections, which are crucial in excited state molecular dynamics simulations. We train the diabatic Hamiltonian elements and their analytical gradients with the SchNet architecture to construct machine learning models, while the phase freedom of off-diagonal elements of the Hamiltonian is cured by introducing the phase-less loss function. Our machine learning models show reasonable accuracy with mean absolute errors of ∼0.1 kcal/mol and ∼0.5 kcal/mol/Å for the diabatic Hamiltonian elements and their gradients, respectively, for penta-2,4-dieniminium cation. Moreover, by exploiting the diabatic representation, our models can predict correct conical intersection structures and their topologies. In addition, our ML-ESMD simulations give almost identical result with a direct dynamics at the same level of theory.

Theoretical study on the photochemistry of furoylazides: Curtius rearrangement and subsequent reactions

Organic azides are an efficient source of nitrenes, which serve as vigorous intermediates in many useful organic reactions. In this work, the complete active space self-consistent field (CASSCF) and its second-order perturbation (CASPT2) methods were employed to study the photochemistry of 2-furoylazide 1 and 3-furoylazide 5, including the Curtius rearrangement to two furylisocyanates (3 and 7) and subsequent reactions to the final product cyanoacrolein 9. Our calculations show that the photoinduced Curtius rearrangement of the two furoylazides takes place through similar stepwise mechanisms via two bistable furoylnitrenes 2 and 6. However, the decarbonylation and ring-opening process of 7 to 9 prefers a stepwise mechanism involving the 3-furoylnitrene intermediate 8, while 3 to 9 goes in a concerted asynchronous way without the corresponding 2-furoylnitrene intermediate 4. Importantly, we revealed that several conical intersections play key roles in the photochemistry of furoylazides. Our results are not only consistent and also make clear the experimental observations (X. Zeng, et al., J. Am. Chem. Soc., 2018, 140, 10-13), but additionally provide important information on the chemistry of furoylazides and nitrenes.

Stabilization of the Triplet Biradical Intermediate of 5-Methylcytosine Enhances Cyclobutane Pyrimidine Dimer (CPD) Formation in DNA

Cyclobutane pyrimidine dimer (CPD) is a photoproduct formed by two stacked pyrimidine bases through a cycloaddition reaction upon irradiation. Owing to its close association with skin cancer, the mechanism of CPD formation has been studied thoroughly. Among many aspects of CPD, its formation involving 5-methylcytosine (5mC) has been of special interest because the CPD yield is known to increase with C5-methylation of cytosine. In this work, high-level quantum mechanics/molecular mechanics (QM/MM) calculations are used to examine a previously experimentally detected pathway for CPD formation in hetero (thymine-cytosine and thymine-5mC) dipyrimidines, which is facilitated through intersystem crossing in thymine and formation of a triplet biradical intermediate. A DNA duplex model system containing a core sequence TmCG or TCG is used. The stabilization of a radical center in the biradical intermediate by the methyl group of 5mC can lead to increased CPD yield in TmCG compared with its non-methylated counterpart, TCG, thereby suggesting the existence of a new pathway of CPD formation enhanced by 5mC.

Quantum and Quantum-Classical Studies of the Photoisomerization of a Retinal Chromophore Model

We report an in-depth analysis of the photo-induced isomerization of the 2-cis-penta-2,4-dieniminium cation: a minimal model of the 11-cis retinal protonated Schiff base chromophore of the dim-light photoreceptor rhodopsin. Based on recently developed three-dimensional potentials parametrized on ab initio multi-state multi-configurational second-order perturbation theory data, we perform quantum-dynamical studies. In addition, simulations based on various quantum-classical methods, among which Tully surface hopping and the coupled-trajectory approach derived from the exact factorization, allow us to validate their performance against vibronic wavepacket propagation and, therefore, a purely quantum treatment. Quantum-dynamics results uncover qualitative differences with respect to the two-dimensional Hahn-Stock potentials, widely used as model potentials for the isomerization of the same chromophore, due to the increased dimensionality and three-mode correlation. Quantum-classical simulations show, instead, that three-dimensional model potentials are capable of capturing a number of features revealed by atomistic simulations and experimental observations. In particular, a recently reported vibrational phase relationship between double-bond torsion and hydrogen-out-of-plane modes critical for rhodopsin isomerization efficiency is correctly reproduced.

Molecular Photochemistry: Recent Developments in Theory

Photochemistry is a fascinating branch of chemistry that is concerned with molecules and light. However, the importance of simulating light-induced processes is reflected also in fields as diverse as biology, material science, and medicine. This Minireview highlights recent progress achieved in theoretical chemistry to calculate electronically excited states of molecules and simulate their photoinduced dynamics, with the aim of reaching experimental accuracy. We focus on emergent methods and give selected examples that illustrate the progress in recent years towards predicting complex electronic structures with strong correlation, calculations on large molecules, describing multichromophoric systems, and simulating non-adiabatic molecular dynamics over long time scales, for molecules in the gas phase or in complex biological environments.

First-Principles Quantum and Quantum-Classical Simulations of Exciton Diffusion in Semiconducting Polymer Chains at Finite Temperature

We report on first-principles quantum-dynamical and quantum-classical simulations of photoinduced exciton dynamics in oligothiophene chain segments, representative of intrachain exciton migration in the poly(3-hexylthiophene) (P3HT) polymer. Following up on our recent study (Binder R.; Burghardt, I. Faraday Discuss. 2020, 221, 406), multilayer multiconfiguration time-dependent Hartree calculations for a short oligothiophene segment comprising 20 monomer units (OT-20) are carried out to obtain full quantum-dynamical simulations at finite temperature. These are employed to benchmark mean-field Ehrenfest calculations, which are shown to give qualitatively correct results for the present system. Periodic boundary conditions turn out to significantly improve earlier estimates of diffusion coefficients. Using the Ehrenfest approach, a series of calculations are subsequently carried out for larger lattices (OT-40 to OT-80), leading to estimates for temperature-dependent mean-squared displacements, which are found to exhibit a near-linear dependence as a function of time. The resulting diffusion coefficient estimates are an increasing function of temperature, whose detailed functional form depends on the degree of static disorder. With a realistic static disorder parameter (σ_s ≃ 0.06 eV), the diffusion coefficients decrease from D ∼ 1 × 10^-2 cm² s^-1 to D ∼ 1 × 10^-3 cm² s^-1, in qualitative agreement with experimental data for P3HT. The dynamical scenario obtained from our simulations shows that exciton migration in P3HT-type chains is a largely adiabatic process throughout the temperature regime we investigated (i.e., T = 50-300 K). The resulting picture of exciton migration is a coherent, but not bandlike, motion of an exciton-polaron driven by fluctuations induced by low-frequency modes. This process acquires partial hopping character if static disorder becomes prominent and Anderson localization sets in.

Reduced scaling extended multi-state CASPT2 (XMS-CASPT2) using supporting subspaces and tensor hyper-contraction

We present a reduced scaling formulation of the extended multi-state CASPT2 (XMS-CASPT2) method, which is based on our recently developed state-specific CASPT2 (SS-CASPT2) formulation using supporting subspaces and tensor hyper-contraction. By using these two techniques, the off-diagonal elements of the effective Hamiltonian can be computed with only O(N³) operations and O(N²) memory, where N is the number of basis functions. This limits the overall computational scaling to O(N⁴) operations and O(N²) memory. Thus, excited states can now be obtained at the same reduced (relative to previous algorithms) scaling we achieved for SS-CASPT2. In addition, we also investigate how the energy denominators can be factorized with the Laplace quadrature when some of the denominators are negative, which is critical for excited state calculations. An efficient implementation of the method has been developed using graphical processing units while also exploiting spatial sparsity in tensor operations. We benchmark the accuracy of the new method by comparison to non-THC formulated XMS-CASPT2 for the excited states of various molecules. In our tests, the THC approximation introduces negligible errors (≈0.01 eV) compared to the non-THC reference method. Scaling behavior and computational timings are presented to demonstrate performance. The new method is also interfaced with quantum mechanics/molecular mechanics (QM/MM). In an example study of green fluorescent protein, we show how the XMS-CASPT2 potential energy surfaces and excitation energies are affected by increasing the size of the QM region up to 278 QM atoms with more than 2300 basis functions.

Variational Principles in Quantum Monte Carlo: The Troubled Story of Variance Minimization

We investigate the use of different variational principles in quantum Monte Carlo, namely, energy and variance minimization, prompted by the interest in the robust and accurate estimation of electronic excited states. For two prototypical, challenging molecules, we readily reach the accuracy of the best available reference excitation energies using energy minimization in a state-specific or state-average fashion for states of different or equal symmetry, respectively. On the other hand, in variance minimization, where the use of suitable functionals is expected to target specific states regardless of the symmetry, we encounter severe problems for a variety of wave functions: as the variance converges, the energy drifts away from that of the selected state. This unexpected behavior is sometimes observed even when the target is the ground state and generally prevents the robust estimation of total and excitation energies. We analyze this problem using a very simple wave function and infer that the optimization finds little or no barrier to escape from a local minimum or local plateau, eventually converging to a lower-variance state instead of the target state. For the increasingly complex systems becoming in reach of quantum Monte Carlo simulations, variance minimization with current functionals appears to be an impractical route.

Probing competing relaxation pathways in malonaldehyde with transient X-ray absorption spectroscopy

Excited-state intramolecular hydrogen transfer (ESIHT) is a fundamental reaction relevant to chemistry and biology. Malonaldehyde is the simplest example of ESIHT, yet only little is known experimentally about its excited-state dynamics. Several competing relaxation pathways have been proposed, including internal conversion mediated by ESIHT and C[double bond, length as m-dash]C torsional motion as well as intersystem crossing. We perform an in silico transient X-ray absorption spectroscopy (TRXAS) experiment at the oxygen K-edge to investigate its potential to monitor the proposed ultrafast decay pathways in malonaldehyde upon photoexcitation to its bright S₂(ππ*) state. We employ both restricted active space perturbation theory and algebraic-diagrammatic construction for the polarization propagator along interpolated reaction coordinates as well as representative trajectories from ab initio multiple spawning simulations to compute the TRXAS signals from the lowest valence states. Our study suggests that oxygen K-edge TRXAS can distinctly fingerprint the passage through the H-transfer intersection and the concomitant population transfer to the S₁(nπ*) state. Potential intersystem crossing to T₁(ππ*) is detectable from reappearance of the double pre-edge signature and reversed intensities. Moreover, the torsional deactivation pathway induces transient charge redistribution from the enol side towards the central C-atom and manifests itself as substantial shifts of the pre-edge features. Given the continuous advances in X-ray light sources, our study proposes an experimental route to disentangle ultrafast excited-state decay channels in this prototypical ESIHT system and provides a pathway-specific mapping of the TRXAS signal to facilitate the interpretation of future experiments.

fromage: A library for the study of molecular crystal excited states at the aggregate scale

The study of photoexcitations in molecular aggregates faces the twofold problem of the increased computational cost associated with excited states and the complexity of the interactions among the constituent monomers. A mechanistic investigation of these processes requires the analysis of the intermolecular interactions, the effect of the environment, and 3D arrangements or crystal packing on the excited states. A considerable number of techniques have been tailored to navigate these obstacles; however, they are usually restricted to in‐house codes and thus require a disproportionate effort to adopt by researchers approaching the field. Herein, we present the FRamewOrk for Molecular AGgregate Excitations (fromage), which implements a collection of such techniques in a Python library complemented with ready‐to‐use scripts. The program structure is presented and the principal features available to the user are described: geometrical analysis, exciton characterization, and a variety of ONIOM schemes. Each is illustrated by examples of diverse organic molecules in condensed phase settings. The program is available at https://github.com/Crespo-Otero-group/fromage.

Four resonance structures elucidate double-bond isomerisation of a biological chromophore

Four resonance structures determining the electronic structure of the chromophore’s ground and first excited states. Changing the relative energies of the structures by hydrogen-bonding interactions tunes all chromophore’s photochemical properties.

Conical Intersections in Organic Molecules: Benchmarking Mixed-Reference Spin-Flip Time-Dependent DFT (MRSF-TD-DFT) vs Spin-Flip TD-DFT

The mixed-reference spin-flip time-dependent density functional theory (MRSF-TD-DFT) method eliminates the erroneous spin contamination of the SF-TD-DFT methodology, while retaining the conceptual and practical simplicity of the latter. The availability of the analytic gradient of the energy of the MRSF-TD-DFT response states enables automatic geometry optimization of the targeted states. Here, we apply the new method to optimize the geometry of several S₁/S₀ conical intersections occurring in typical organic molecules. We demonstrate that MRSF-TD-DFT is capable of producing the correct double-cone topology of the intersections and describing the geometry of the lowest-energy conical intersections and their relative energies with accuracy matching that of the best multireference wavefunction ab initio methods. In this regard, MRSF-TD-DFT differs from many popular single-reference methods, such as, e.g., the linear response TD-DFT method, which fail to produce the correct topology of the intersections. As the new methodology completely eliminates the ambiguity with the identification of the response states as proper singlets or triplets, which is plaguing the SF-TD-DFT calculations, it can be used for automatic geometry optimization and molecular dynamic simulations not requiring constant human intervention.