Quantum Mechanical Studies on the Photophysics and the Photochemistry of Nucleic Acids and Nucleobases

Solvent-Dependent Stabilization of a Charge Transfer State is the Key to Ultrafast Triplet State Formation in an Epigenetic DNA Nucleoside

2'-Deoxy-5-formylcytidine (5fdCyd), a naturally occurring nucleoside found in mammalian DNA and mitochondrial RNA, exhibits important epigenetic functionality in biological processes. Because it efficiently generates triplet excited states, it is an endogenous photosensitizer capable of damaging DNA, but the intersystem crossing (ISC) mechanism responsible for ultrafast triplet state generation is poorly understood. In this study, time-resolved mid-IR spectroscopy and quantum mechanical calculations reveal the distinct ultrafast ISC mechanisms of 5fdCyd in water versus acetonitrile. Our experiment indicates that in water, ISC to triplet states occurs within 1 ps after 285 nm excitation. PCM-TD-DFT computations suggest that this ultrafast ISC is mediated by a singlet state with significant cytosine-to-formyl charge-transfer (CT) character. In contrast, ISC in acetonitrile proceeds via a dark ¹ nπ* state with a lifetime of ∼3 ps. CT-induced ISC is not favored in acetonitrile because reaching the minimum of the gateway CT state is hampered by intramolecular hydrogen bonding, which enforces planarity between the aldehyde group and the aromatic group. Our study provides a comprehensive picture of the non-radiative decay of 5fdCyd in solution and new insights into the factors governing ISC in biomolecules. We propose that the intramolecular CT state observed here is a key to the excited-state dynamics of epigenetic nucleosides with modified exocyclic functional groups, paving the way to study their effects in DNA strands.

Unraveling the Mechanism for Tuning the Fluorescence of Fluorescein Derivatives: The Role of the Conical Intersection and nπ* State

Although a large number of fluorescein derivatives have been developed and applied in many different fields, the general mechanisms for tuning the fluorescence of fluorescein derivatives still remain uncovered. Herein, we found that the fluorescence quenching of neutral form of fluorescein derivatives in acidic medium resulted from a dark nπ* state, whereas the fluorescence of the anionic form of fluorescein derivatives in the gas phase and alkaline solutions was tuned by minimal energy conical intersection (MECI). The formation of MECI involved significant rotation of benzene ring and flip-flop motion of xanthene moiety, which would be restricted by intermolecular hydrogen bonding and lowering temperature. The energy barrier for reaching MECI depended on the substituents in the benzene moiety in accordance with experimentally observed substituent effects. These unprecedented mechanisms would lead to a recognition of fluorescein derivatives and could provide a correct and instructive design strategy for further developing new fluorescein derivatives.

Computational Design of Radical Recognition Assay with the Possible Application of Cyclopropyl Vinyl Sulfides as Tunable Sensors

The processes involving the capture of free radicals were explored by performing DFT molecular dynamics simulations and modeling of reaction energy profiles. We describe the idea of a radical recognition assay, where not only the presence of a radical but also the nature/reactivity of a radical may be assessed. The idea is to utilize a set of radical-sensitive molecules as tunable sensors, followed by insight into the studied radical species based on the observed reactivity/selectivity. We utilize this approach for selective recognition of common radicals-alkyl, phenyl, and iodine. By matching quantum chemical calculations with experimental data, we show that components of a system react differently with the studied radicals. Possible radical generation processes were studied involving model reactions under UV light and metal-catalyzed conditions.

Time Resolved Photoelectron Spectroscopy as a Test of Electronic Structure and Nonadiabatic Dynamics

We compare different levels of theory for simulating excited state molecular dynamics and use time-resolved photoelectron spectroscopy measurements to benchmark the theory. We perform trajectory surface hopping simulations for uracil excited to the first bright state (ππ*) using three different levels of theory (CASSCF, MRCIS, and XMS-CASPT2) in order to understand the role of dynamical correlation in determining the excited state dynamics, with a focus on the coupling between different electronic states and internal conversion back to the ground state. These dynamics calculations are used to simulate the time-resolved photoelectron spectra. The comparison of the calculated and measured spectra allows us to draw conclusions regarding the relative insights and quantitative accuracy of the calculations at the three different levels of theory, demonstrating that detailed quantitative comparisons of time-resolved photoelectron spectra can be used to benchmark methodology.

Spiers Memorial Lecture: Introduction to ultrafast spectroscopy and imaging of photochemical reactions

A brief overview is presented on ultrafast spectroscopy and imaging of photochemical reactions by highlighting several experimental studies reported in the last five years. A particular focus is placed on new experiments performed using high-order harmonic generation, X-ray free electron lasers, and relativistic electron beams. Exploration of fundamental chemical reaction dynamics using these advanced experimental methodologies is in an early stage, and exciting new research opportunities await in this rapidly expanding and advancing research field. At the same time, there is no experimental methodology that provides all aspects of the electronic and structural dynamics in a single experiment, and investigations using different methodologies with various perspectives need to be considered in a comprehensive manner.

Theoretical Study on the Photo-Oxidation and Photoreduction of an Azetidine Derivative as a Model of DNA Repair

Photocycloreversion plays a central role in the study of the repair of DNA lesions, reverting them into the original pyrimidine nucleobases. Particularly, among the proposed mechanisms for the repair of DNA (6-4) photoproducts by photolyases, it has been suggested that it takes place through an intermediate characterized by a four-membered heterocyclic oxetane or azetidine ring, whose opening requires the reduction of the fused nucleobases. The specific role of this electron transfer step and its impact on the ring opening energetics remain to be understood. These processes are studied herein by means of quantum-chemical calculations on the two azetidine stereoisomers obtained from photocycloaddition between 6-azauracil and cyclohexene. First, we analyze the efficiency of the electron-transfer processes by computing the redox properties of the azetidine isomers as well as those of a series of aromatic photosensitizers acting as photoreductants and photo-oxidants. We find certain stereodifferentiation favoring oxidation of the cis-isomer, in agreement with previous experimental data. Second, we determine the reaction profiles of the ring-opening mechanism of the cationic, neutral, and anionic systems and assess their feasibility based on their energy barrier heights and the stability of the reactants and products. Results show that oxidation largely decreases the ring-opening energy barrier for both stereoisomers, even though the process is forecast as too slow to be competitive. Conversely, one-electron reduction dramatically facilitates the ring opening of the azetidine heterocycle. Considering the overall quantum-chemistry findings, N,N-dimethylaniline is proposed as an efficient photosensitizer to trigger the photoinduced cycloreversion of the DNA lesion model.

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.

Velocity Adjustment in Surface Hopping: Ethylene as a Case Study of the Maximum Error Caused by Direction Choice

The most common surface hopping dynamics algorithms require velocity adjustment after hopping to ensure total-energy conservation. Based on the semiclassical analysis, this adjustment must be made parallel to the nonadiabatic coupling vector's direction. Nevertheless, this direction is not always known, and the common practice has been to adjust the velocity in either the linear momentum or velocity directions. This paper benchmarks surface hopping dynamics of photoexcited ethylene with velocity adjustment in several directions, including those of the nonadiabatic coupling vector, the momentum, and the energy gradient difference. It is shown that differences in time constants and structural evolution fall within the statistical uncertainty of the method considering up to 500 trajectories in each dynamics set, rendering the three approaches statistically equivalent. For larger ensembles beyond 1000 trajectories, significant differences between the results arise, limiting the validity of adjustment in alternative directions. Other possible adjustment directions (velocity, single-state gradients, angular momentum) are evaluated as well. Given the small size of ethylene, the results reported in this paper should be considered an upper limit for the error caused by the choice of the velocity-adjustment direction on surface hopping dynamics.

The effect of flanking bases on direct and triplet sensitized cyclobutane pyrimidine dimer formation in DNA depends on the dipyrimidine, wavelength and the photosensitizer

Cyclobutane pyrimidine dimers (CPDs) are the major products of DNA produced by direct absorption of UV light, and result in C to T mutations linked to human skin cancers. Most recently a new pathway to CPDs in melanocytes has been discovered that has been proposed to arise from a chemisensitized pathway involving a triplet sensitizer that increases mutagenesis by increasing the percentage of C-containing CPDs. To investigate how triplet sensitization may differ from direct UV irradiation, CPD formation was quantified in a 129-mer DNA designed to contain all 64 possible NYYN sequences. CPD formation with UVB light varied about 2-fold between dipyrimidines and 12-fold with flanking sequence and was most frequent at YYYR and least frequent for GYYN sites in accord with a charge transfer quenching mechanism. In contrast, photosensitized CPD formation greatly favored TT over C-containing sites, more so for norfloxacin (NFX) than acetone, in accord with their differing triplet energies. While the sequence dependence for photosensitized TT CPD formation was similar to UVB light, there were significant differences, especially between NFX and acetone that could be largely explained by the ability of NFX to intercalate into DNA.

Ultrafast Dynamics of the Isolated Adenosine-5'-triphosphate Dianion Probed by Time-Resolved Photoelectron Imaging

The excited state dynamics of the doubly deprotonated dianion of adenosine-5'-triphosphate, [ATP-H₂]^2-, has been spectroscopically explored by time-resolved photoelectron spectroscopy following excitation at 4.66 eV. Time-resolved photoelectron spectra show that two competing processes occur for the initially populated ¹ππ* state. The first is rapid electron emission by tunneling through a repulsive Coulomb barrier as the ¹ππ* state is a resonance. The second is nuclear motion on the ¹ππ* state surface leading to an intermediate that no longer tunnels and subsequently decays by internal conversion to the ground electronic state. The spectral signatures of the features are similar to those observed for other adenine-derivatives, suggesting that this nucleobase is quite insensitive to the nearby negative charges localized on the phosphates, except of course for the appearance of the additional electron tunneling channel, which is open in the dianion.

Simulation of UV absorption spectra and relaxation dynamics of uracil and uracil-water clusters

Despite many studies, the mechanisms of nonradiative relaxation of uracil in the gas phase and in aqueous solution are still not fully resolved. Here we combine theoretical UV absorption spectroscopy with nonadiabatic dynamics simulations to identify the photophysical mechanisms that can give rise to experimentally observed decay time constants. We first compute and theoretically assign the electronic spectra of uracil using the second-order algebraic-diagrammatic-construction (ADC(2)) method. The obtained electronic states, their energy differences and state-specific solvation effects are the prerequisites for understanding the photodynamics. We then use nonadiabatic trajectory-surface-hopping dynamics simulations to investigate the photoinduced dynamics of uracil and uracil-water clusters. In contrast to previous studies, we found that a single mechanism - the ethylenic twist around the C[double bond, length as m-dash]C bond - is responsible for the ultrafast component of the nonradiative decay, both in the gas phase and in solution. Very good agreement with the experimentally determined ultrashort decay time constants is obtained.

Quantum and semiclassical dynamical studies of nonadiabatic processes in solution: achievements and perspectives

We concisely review the main methodological approaches to model nonadiabatic dynamics in isotropic solutions and their applications. Three general classes of models are identified as the most used to include solvent effects in the simulations. The first model describes the solvent as a set of harmonic collective modes coupled to the solute degrees of freedom, and the second as a continuum, while the third explicitly includes solvent molecules in the calculations. The issues related to the use of these models in semiclassical and quantum dynamical simulations are discussed, as well as the main limitations and perspectives of each approach.

Theoretical insights on the excited-state-deactivation mechanisms of protonated thymine and cytosine

Ab initio and surface-hopping nonadiabatic dynamics simulation methods were employed to investigate relaxation mechanisms in protonated thymine (TH⁺) and cytosine (CH⁺). A few conical intersections were located between ¹ππ* and S₀ states for each system with the CASSCF (8,8) theoretical model and relevant contributions to the deactivation mechanism of titled systems were addressed by the determination of potential energy profiles at the CASPT2 (12,10) theoretical level. It was revealed that the relaxation of the ¹ππ* state of the most stable conformer of both systems to the ground state is mostly governed by the accessible S₁/S₀ conical intersection resulting from the barrier-free out-of-plane deformation. Interestingly, it was exhibited that the ring puckering coordinate driven from the C₆ position of the heterocycle ring in TH⁺ and CH⁺ plays the most prominent role in the deactivation mechanism of considered systems. Our ab initio results are also supported by excited-state nonadiabatic dynamics simulations based on ADC(2), describing the ultrashort S₁ lifetime of TH⁺/CH⁺ by analyzing trajectories leading excited systems to the ground. It was confirmed that the excited-state population mostly relaxes to the ground via the ring puckering coordinate from the C₆ moiety. Overall, the theoretical results of this study shed light on the deactivation mechanism of protonated DNA bases.

Nonadiabatic Absorption Spectra and Ultrafast Dynamics of DNA and RNA Photoexcited Nucleobases

We have recently proposed a protocol for Quantum Dynamics (QD) calculations, which is based on a parameterisation of Linear Vibronic Coupling (LVC) Hamiltonians with Time Dependent (TD) Density Functional Theory (TD-DFT), and exploits the latest developments in multiconfigurational TD-Hartree methods for an effective wave packet propagation. In this contribution we explore the potentialities of this approach to compute nonadiabatic vibronic spectra and ultrafast dynamics, by applying it to the five nucleobases present in DNA and RNA. For all of them we computed the absorption spectra and the dynamics of ultrafast internal conversion (100 fs timescale), fully coupling the first 2-3 bright states and all the close by dark states, for a total of 6-9 states, and including all the normal coordinates. We adopted two different functionals, CAM-B3LYP and PBE0, and tested the effect of the basis set. Computed spectra are in good agreement with the available experimental data, remarkably improving over pure electronic computations, but also with respect to vibronic spectra obtained neglecting inter-state couplings. Our QD simulations indicate an effective population transfer from the lowest energy bright excited states to the close-lying dark excited states for uracil, thymine and adenine. Dynamics from higher-energy states show an ultrafast depopulation toward the more stable ones. The proposed protocol is sufficiently general and automatic to promise to become useful for widespread applications.

Long living excited state of protonated adenosine unveiled by ultrafast fluorescence spectroscopy and density functional theoretical study

Adenosine (Ado) possesses ultrafast nonradiative dynamics accounting for its remarkably high photostability. The deactivation dynamics of Ado after protonation in an aqueous solution remains an elusive issue. Herein we report an investigation of the excited state dynamics of protonated Ado (AdoH⁺) performed using ultrafast time-resolved fluorescence spectroscopy combined with density functional theoretical calculation. The result obtained from comparison of conformers with protonation at different sites revealed that the syn-conformer with protonation occurring at the N3 position (syn-N3) is the predominant form of AdoH⁺ in the ground state, similar to that of Ado. In contrast, the fluorescence of AdoH⁺ with maximum intensity at 385 nm, significantly red-shifted from that of Ado, displaying decay dynamics composed of an ultrafast component with the lifetime of ∼0.5 ps and a slower one of ∼2.9 ns. The former is because of the decay of the syn-N3 conformer, similar to that reported for AdoH⁺ under the gas phase condition. The latter is due to the syn-N1 conformer formed via ultrafast proton transfer of the syn-N3. The excited state of syn-N1 has a peculiar nonplanar conformation over the purine molecule, which is responsible for the substantial Stokes shift showed in the fluorescence spectrum and correlates with a large energy barrier for nonradiative decay likely involving a reversed proton transfer. This study demonstrates the importance of protonation and solvent environment in altering dramatically the excited states of Ado, providing insight for better understanding nonradiative dynamics of both the monomeric bases and the oligomeric or polymeric DNAs.

Fundamentals of the Intrinsic DNA Fluorescence

The intrinsic fluorescence of nucleic acids is extremely weak compared to that of the fluorescent labels used to probe their structural and functional behavior. Thus, for technical reasons, the investigation of the intrinsic DNA fluorescence was limited for a long time. But with the improvement in spectroscopic techniques, the situation started to change around the turn of the century. During the past two decades, various factors modulating the static and dynamic properties of the DNA fluorescence have been determined; it was shown that, under certain conditions, quantum yields may be up 100 times higher than what was known so far. The ensemble of these studies opened up new paths for the development of label-free DNA fluorescence for biochemical applications. In parallel, these studies have shed new light on the primary processes leading to photoreactions that damage DNA when it absorbs UV radiation.We have been studying a variety of DNA systems, ranging from the monomeric nucleobases to double-stranded and four-stranded structures using fluorescence spectroscopy. The specificity of our work resides in the quantitative association of the steady-state fluorescence spectra with time-resolved data recorded from the femtosecond to the nanosecond timescales, made possible by the development of specific methodologies.Among others, our fluorescence studies provide information on the energy and the polarization of electronic transitions. These are valuable indicators for the evolution of electronic excitations in complex systems, where the electronic coupling between chromophores plays a key role. Highlighting collective effects that originate from electronic interactions in DNA multimers is the objective of the present Account.In contrast to the monomeric chromophores, whose fluorescence decays within a few picoseconds, that of DNA multimers persists on the nanosecond timescale. Even if long-lived states represent only a small fraction of electronic excitations, they may be crucial to the DNA photoreactivity because the probability to reach reactive conformations increases over time, owing to the incessant structural dynamics of nucleic acids.Our femtosecond studies have revealed that an ultrafast excitation energy transfer takes place among the nucleobases within duplexes and G-quadruplexes. Such an ultrafast process is possible when collective states are populated directly upon photon absorption. At much longer times, we discovered an unexpected long-lived high-energy emission stemming from what was coined "HELM excitons". These collective states, whose emission increases with the duplex size, could be responsible for the delayed fluorescence of ππ* states observed for genomic DNA.Most studies dealing with excited-state relaxation in DNA were carried out with excitation in the absorption band peaking at around 260 nm. We went beyond this and also performed the first time-resolved study with excitation in the UVA spectral range, where a very weak absorption tail is present. The resulting fluorescence decays are much slower and the fluorescence quantum yields are much higher than for UVC excitation. We showed that the base pairing of DNA strands enhances the UVA fluorescence and, in parallel, increases the photoreactivity because it modifies the nature of the involved collective excited states.

Excited-State Properties and Relaxation Pathways of Selenium-Substituted Guanine Nucleobase in Aqueous Solution and DNA Duplex

The excited-state properties and relaxation mechanisms after light irradiation of 6-selenoguanine (6SeG) in water and in DNA have been investigated using a quantum mechanics/molecular mechanics (QM/MM) approach with the multistate complete active space second-order perturbation theory (MS-CASPT2) method. In both environments, the S₁ ¹(n_Seπ₅^*) and S₂ ¹(π_Seπ₅^*) states are predicted to be the spectroscopically dark and bright states, respectively. Two triplet states, T₁ ³(π_Seπ₅^*) and T₂ ³(n_Seπ₅^*), are found energetically below the S₂ state. Extending the QM region to include the 6SeG-Cyt base pair slightly stabilizes the S₂ state and destabilizes the S₁, due to hydrogen-bonding interactions, but it does not affect the order of the states. The optimized minima, conical intersections, and singlet-triplet crossings are very similar in water and in DNA, so that the same general mechanism is found. Additionally, for each excited state geometry optimization in DNA, three kind of structures ("up", "down", and "central") are optimized which differ from each other by the orientation of the C═Se group with respect to the surrounding guanine and thymine nucleobases. After irradiation to the S₂ state, 6SeG evolves to the S₂ minimum, near to a S₂/S₁ conical intersection that allows for internal conversion to the S₁ state. Linear interpolation in internal coordinates indicate that the "central" orientation is less favorable since extra energy is needed to surmount the high barrier in order to reach the S₂/S₁ conical intersection. From the S₁ state, 6SeG can further decay to the T₁ ³(π_Seπ₅^*) state via intersystem crossing, where it will be trapped due to the existence of a sizable energy barrier between the T₁ minimum and the T₁/S₀ crossing point. Although this general S₂ → T₁ mechanism takes place in both media, the presence of DNA induces a steeper S₂ potential energy surface, that it is expected to accelerate the S₂ → S₁ internal conversion.

Modeling solvation effects on absorption and fluorescence spectra of indole in aqueous solution

Modeling the optical spectra of molecules in solution presents a challenge, so it is important to understand which of the solvation effects (i.e., electrostatics, mutual polarization, and hydrogen bonding interactions between solute and solvent molecules) are crucial in reproducing the various features of the absorption and fluorescence spectra and to identify a sufficient theoretical model that accurately captures these effects with minimal computational cost. In this study, we use various implicit and explicit solvation models, such as molecular dynamics coupled with non-polarizable and polarizable force fields, as well as Car-Parrinello molecular dynamics, to model the absorption and fluorescence spectra of indole in aqueous solution. The excited states are computed using the equation of motion coupled cluster with single and double excitations combined with the effective fragment potential to represent water molecules, which we found to be a computationally efficient approach for modeling large solute-solvent clusters at a high level of quantum theory. We find that modeling mutual polarization, compared to other solvation effects, is a dominating factor for accurately reproducing the position of the peaks and spectral line shape of the absorption spectrum of indole in solution. We present an in-depth analysis of the influence that different solvation models have on the electronic excited states responsible for the features of the absorption spectra. Modeling fluorescence is more challenging since it is hard to reproduce even the correct emitting state, and force field parameters need to be re-evaluated.

Excited-State Absorption of Uracil in the Gas Phase: Mapping the Main Decay Paths by Different Electronic Structure Methods

We present a computational study of the one-photon and excited-state absorption (ESA) from the two lowest energy excited states of uracil in the gas phase: an nπ* dark state (1n) and the lowest energy bright ππ* state (1π). The predictions of six different linear response electronic structure methods, namely, TD-CAM-B3LYP, EOM-CCSD, EOM-CC3, ADC(2), ADC(2)-x, and ADC(3) are critically compared. In general, the spectral shapes predicted by TD-CAM-B3LYP, EOM-CCSD, EOM-CC3, and ADC(3) are fairly similar, though the quality of TD-CAM-B3LYP slightly deteriorates in the high-energy region. By computing the spectra at some key structures on different potential energy surfaces (PES), that is, the Franck-Condon point, the 1n minimum, and structures representative of different regions of the 1π PES, we obtain important insights into the shift of the ESA spectra, following the motion of the wavepacket on the excited-state PES. Though 1π has larger ESA than 1n, some spectral regions are dominated by these latter signals. Aside from its methodological interest, we thus obtain interesting indications to interpret transient absorption spectra to disentangle the photoactivated dynamics of nucleobases.

Controlled Fluorescence Enhancement of DNA-Binding Dye Through Chain Length Match between Oligoguanine and TOTO

Fluorescent DNA-binding dyes are extensively employed as probe and biosensing in biological detection and imaging. Experiments and theoretical calculations of thiazole orange homodimeric (TOTO) dye binding to a single-strand DNA (ssDNA), poly(dG)_n (n = 2, 4, 6, 8), reveal that the n = 6 complex shows about 300-fold stronger fluorescence than n = 2, 4 and a slightly stronger one than n = 8 complexes, which is benefited from the length match between TOTO and poly(dG)₆. The machine learning, based on molecular dynamics trajectories, indicates that TOTO is featured by the dihedral angle along its backbone and its end-to-end distance, in which the latter one defines the stretch and hairpin structures of TOTO, respectively. The time-dependent density functional theory calculations on the low-lying excited states show that the stretched TOTO with π-π end-stacking binding mode can bring about strong fluorescence with localized π-π^* transitions. For the n = 2, 4, and 8 complexes, the linear scaling quantum mechanics calculations indicate that the dominant hairpin TOTO with intercalative binding modes have relatively larger binding energies, leading to fluorescence quenching by intramolecular charge transfer. Our results may provide an insight for modulating the DNA-dye binding modes to tune the degree of charge transfer and designing fluorescent probes for the recognition of specific DNA sequences.

Molecular Rotors as Guest Fluorophores Probing the Local Environment inside Host G4 Supramolecular Hydrogels

Fluorescent molecular rotors with a high binding affinity toward the guanosine quartet (G4) were incorporated as guest fluorophores into host supramolecular hydrogels based on the self-assembly of G4 units, to probe the local environment. Torsional dynamics of the rotors were severely inhibited inside the hydrogels in comparison with aqueous solutions, although the hydrogels were composed of >95% water. Moreover, even though all the gels were rigid bodies with no spontaneous deformation or flow property at room temperature, torsional dynamics in G4 borate gels was found to be consistently several orders of magnitude slower than those in the other G4 gels, irrespective of the identity of the molecular rotor probe. This clear difference in the molecular mobilities of the guest fluorophore could be attributed to systematic differences in the internal structure between the two categories of host G4 hydrogels. In specific terms, the borate groups in G4 borate hydrogels serve as bridging units between separate G4 quadruplex strands, generating additional cross-links that reinforce the network structure of the gel. The results demonstrate that molecular rotors act as efficient fluorescent probes for the quantitative assessment of the molecular-level environment and dynamics inside the hydrogels, an aspect that is missed out by most other analytical methods that are routinely employed for studying them.

Excitonic Model for Strongly Coupled Multichromophoric Systems: The Electronic Circular Dichroism Spectra of Guanine Quadruplexes as Test Cases

We here propose a general and flexible approach, based on fragment diabatization, which incorporates charge transfer states and significantly increases the reliability of excitonic Hamiltonians for systems where the chromophores are very close. This model (FrDEx) is used to compute the electronic circular dichroism and absorption spectra of two prototype guanine-rich DNA sequences folded in quadruple helices (GQs), i.e., a fragment of the human telomeric sequence (Tel21, antiparallel), and (TGGGGT)₄ (TG4T, parallel). Calculations on different subsets of Tel21 and TG4T, from dimers to tetramers, show that FrDEx provides spectra close to the reference full quantum mechanical (QM) ones (obtained with time-dependent density functional theory), with significant improvements with respect to "standard" excitonic Hamiltonians. Furthermore, these tests enable the most cost-effective procedure for the whole GQ to be determined. FrDEx spectra of Tel21 and TG4T are also in good agreement with the QM and experimental ones and give access to interesting insights into the chemical-physical effects modulating the spectral signals. FrDEx could be profitably used to investigate many other biological and nanotechnological materials, from DNA to (opto)electronic polymers.

Accurate Molecular Geometries in Complex Excited-State Potential Energy Surfaces from Time-Dependent Density Functional Theory

The interplay of electronic excitations and structural changes in molecules impacts nonradiative decay and charge transfer in the excited state, thus influencing excited-state lifetimes and photocatalytic reaction rates in optoelectronic and energy devices. To capture such effects requires computational methods providing an accurate description of excited-state potential energy surfaces and geometries. We suggest time-dependent density functional theory using optimally tuned range-separated hybrid (OT-RSH) functionals as an accurate approach to obtain excited-state molecular geometries. We show that OT-RSH provides accurate molecular geometries in excited-state potential energy surfaces that are complex and involve an interplay of local and charge-transfer excitations, for which conventional semilocal and hybrid functionals fail. At the same time, the nonempirical OT-RSH approach maintains the high accuracy of parametrized functionals (e.g., B3LYP) for predicting excited-state geometries of small organic molecules showing valence excited states.

Nonadiabatic molecular dynamics simulation for the ultrafast photoisomerization of dMe-OMe-NAIP based on TDDFT on-the-fly potential energy surfaces

Global switching on-the-fly trajectory surface hopping molecular dynamics simulation was performed on the accurate TD-B3LYP/6-31G* potential energy surfaces for E-to-Z and Z-to-E photoisomerization of dMe-OMe-NAIP up to S1(ππ*) excitation. The present TD(DFT) simulation provides accurate calculation for conical intersections between the first-excited and ground states. Thus, simulated quantum yield and lifetime of 0.23 and 620 fs (0.15 and 600 fs) for E-to-Z (Z-to-E) isomerization are in good (relatively good) agreement with experimental observation of 0.25 and 480 fs (0.24 and 430 fs), respectively. Simulated results reveal that photoisomerization pathways are initially uphill to conical intersection zones on the S1 potential energy surface and then downhill to product zones. Three types of representative conical intersections are found for determining photoisomerization mechanisms: one is the rotation type responsible for reactive isomerization and the other two are close to E and Z configurations, respectively, only for nonreactive isomerization. The present conclusions can be held in general for similar large NAIP systems of photoinduced isomerization based on E and Z configurations.

Effect of dynamic correlation on the ultrafast relaxation of uracil in the gas phase

High level multi-reference non-adiabatic dynamics simulations reveal that uracil’s photoexcited S₂ state decays very quickly without any significant trapping.

Benchmarking Quantum Mechanical Methods for the Description of Charge-Transfer States in π-Stacked Nucleobases

Charge-transfer (CT) states are of special interest in photochemical research because they can facilitate chemical reactions through the rearrangement of electrons and subsequently chemical bonds in a molecular system. Of particular importance to this research is the transfer of electrons between π-stacked nucleobases in DNA because they play an important role in its photophysics and photochemistry. Computational methods are paramount for the study of CT states because of the current inability of experimental methods to easily detect such states. However, many ab-initio wavefunction-based and density functional theory (DFT) methods fail to accurately describe these CT states. Here, we benchmark how 40 different quantum mechanical methods describe the excited states of a guanine-thymine π-stacked nucleobase dimer system, both in 5'-TG-3' and 5'-GT-3' conformations. We find that the distance between the nucleobases plays a major role in the energy of the CT state and in the difference of the dipole moments between the CT and ground state. There is a larger range of values (and errors) for the energies of CT states compared to those of states localized on one nucleobase. Wavefunction-based methods have similar errors for the CT and localized valence states, while DFT methods are very sensitive to the amount of Hartree-Fock exchange. Long-range-corrected functionals with a careful balance of the Hartree-Fock exchange included can predict very accurate CT states and a balanced description with the localized states.

Carbonyl Stretch as a Franck-Condon Active Mode and Driving Force for Excited-State Decay of 8-Methoxy-4-methyl-2H-benzo[g]chromen-2-one from nπ* State

The fluorescence of most organic chromophore is emitted from the ππ* state, whereas the nπ* state, as a dark state, plays an important role in quenching the fluorescence when its energy is close to the ππ* state. Herein, we report a theoretical study on the fluorescence quenching of 8-methoxy-4-methyl-2H-benzo[g]chromen-2-one by the nπ* state and propose a new mechanism for describing the vibronic coupling between the ππ* and nπ* states. By applying extended multistate complete-active-space second-order perturbation theory (XMS-CASPT2) to optimize the geometries, the geometry distortion of the ππ* state along the out-of-plane mode is observed. This geometry distortion causes the stretching vibration of the carbonyl group to be coupled with the C-C bonds of the pyran ring, which become a Franck-Condon active mode upon photoexcitation and provides a driving force for nonradiative decay from the nπ* state, even if it is energetically unfavorable. This mechanism is significantly different from the previously proposed "proximity effect" and cannot be captured by the popularly used time-dependent density functional theory (TDDFT) and complete-active-space self-consistent field (CASSCF) methods.

Electronic couplings for photo-induced processes from subsystem time-dependent density-functional theory: The role of the diabatization

Subsystem time-dependent density-functional theory (sTDDFT) making use of approximate non-additive kinetic energy (NAKE) functionals is known to be capable of describing excitation energy transfer processes in a variety of applications. Here, we show that sTDDFT, especially when combined with projection-based embedding (PbE), can be employed for the entire range of photo-induced electronic couplings essential for modeling photophysical properties of complex chemical and biological systems and therefore represents a complete toolbox for this class of problems. This means that it is capable of capturing the interaction/coupling associated with local- and charge-transfer (CT) excitons. However, this requires the choice of a reasonable diabatic basis. We therefore propose different diabatization strategies of the virtual orbital space in PbE-sTDDFT and show how CT excitations can be included in sTDDFT using NAKE functionals via a phenomenological approach. Finally, these electronic couplings are compared to couplings from a multistate fragment excitation difference (FED)-fragment charge difference (FCD) diabatization procedure. We show that both procedures, multistate FED-FCD and sTDDFT (with the right diabatization procedure chosen), lead to an overall good agreement for the electronic couplings, despite differences in their general diabatization strategy. We conclude that the entire range of photo-induced electronic couplings can be obtained using sTDDFT (with the right diabatization procedure chosen) in a black-box manner.

Multiscale Conformational Sampling Reveals Excited-State Locality in DNA Self-Repair Mechanism

Ultraviolet (UV) irradiation is known to be responsible for DNA damage. However, experimental studies in DNA oligonucleotides have shown that UV light can also induce sequence-specific self-repair. Following charge transfer from a guanine adenine sequence adjacent to a cyclobutane pyrimidine dimer (CPD), the covalent bond between the two thymines could be cleaved, recovering the intact base sequence. Mechanistic details promoting the self-repair remained unclear, however. In our theoretical study, we investigated whether optical excitation could directly lead to a charge-transfer state, thereby initiating the repair, or whether the initial excited state remains localized on a single nucleobase. We performed conformational sampling of 200 geometries of the damaged DNA double strand solvated in water and used a hybrid quantum and molecular mechanics approach to compute excited states at the complete active space perturbation level of theory. Analysis of the conformational data set clearly revealed that the excited-state properties are uniformly distributed across the fluctuations of the nucleotide in its natural environment. From the electronic wavefunction, we learned that the electronic transitions remained predominantly local on either adenine or guanine, and no direct charge transfer occurred in the experimentally accessed energy range. The investigated base sequence is not only specific to the CPD repair mechanism but ubiquitously occurs in nucleic acids. Our results therefore give a very general insight into the charge locality of UV-excited DNA, a property that is regarded to have determining relevance in the structural consequences following absorption of UV photons.

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.

What Makes Thienoguanosine an Outstanding Fluorescent DNA Probe?

Thienoguanosine (^thG) is an isomorphic guanosine (G) surrogate that almost perfectly mimics G in nucleic acids. To exploit its full potential and lay the foundation for future applications, 20 DNA duplexes, where the bases facing and neighboring ^thG were systematically varied, were thoroughly studied using fluorescence spectroscopy, molecular dynamics simulations, and mixed quantum mechanical/molecular mechanics calculations, yielding a comprehensive understanding of its photophysics in DNA. In matched duplexes, ^thG's hypochromism was larger for flanking G/C residues but its fluorescence quantum yield (QY) and lifetime values were almost independent of the flanking bases. This was attributed to high duplex stability, which maintains a steady orientation and distance between nucleobases, so that a similar charge transfer (CT) mechanism governs the photophysics of ^thG independently of its flanking nucleobases. ^thG can therefore replace any G residue in matched duplexes, while always maintaining similar photophysical features. In contrast, the local destabilization induced by a mismatch or an abasic site restores a strong dependence of ^thG's QY and lifetime values on its environmental context, depending on the CT route efficiency and solvent exposure of ^thG. Due to this exquisite sensitivity, ^thG appears ideal for monitoring local structural changes and single nucleotide polymorphism. Moreover, ^thG's dominant fluorescence lifetime in DNA is unusually long (9-29 ns), facilitating its selective measurement in complex media using a lifetime-based or a time-gated detection scheme. Taken together, our data highlight ^thG as an outstanding emissive substitute for G with good QY, long fluorescence lifetimes, and exquisite sensitivity to local structural changes.

An nπ* gated decay mediates excited-state lifetimes of isolated azaindoles

Aiming to serve as a guide to understand the relaxation mechanisms of more complex aza-aromatic compounds, such as purine bases, we have studied the non-radiative channels of a set of azaindole structural isomers: 4-, 5-, 6- and 7-azaindole (AI). The relaxation of the isolated molecules, after excitation at the low energy portion of their spectra, has been tracked by femtosecond time-resolved ionization, and the decay paths have been obtained with MS-CASPT2//TD-DFT calculations. Although the ultrashort measured lifetimes for 5- and 6-AI are in contrast to the long-living excited state found in 7-AI, the calculations describe a common relaxation pathway. Along it, the initially excited ππ* states decay to the ground state through a conical intersection accessed through an nπ* state that functions as a gate state. The work reveals that the position of the nitrogen atoms in the purine ring determines the barrier to access the gate state and therefore, the rate of the non-radiative relaxation.

The Bethe-Salpeter Equation Formalism: From Physics to Chemistry

The Bethe-Salpeter equation (BSE) formalism is steadily asserting itself as a new efficient and accurate tool in the ensemble of computational methods available to chemists in order to predict optical excitations in molecular systems. In particular, the combination of the so-called GW approximation, giving access to reliable ionization energies and electron affinities, and the BSE formalism, able to model UV/vis spectra, has shown to provide accurate singlet excitation energies with a typical error of 0.1-0.3 eV. With a similar computational cost as time-dependent density-functional theory (TD-DFT), BSE is able to provide an accuracy on par with the most accurate global and range-separated hybrid functionals without the unsettling choice of the exchange-correlation functional, resolving further known issues (e.g., charge-transfer excitations). In this Perspective, we provide a historical overview of BSE, with a particular focus on its condensed-matter roots. We also propose a critical review of its strengths and weaknesses in different chemical situations.

Embedding-theory-based simulations using experimental electron densities for the environment

For the first time, the use of experimentally derived molecular electron densities as ρ_B(r) in calculations based on frozen-density embedding theory (FDET) of environment-induced shifts of electronic excitations for chromophores in clusters is demonstrated. ρ_B(r) was derived from X-ray restrained molecular wavefunctions of glycylglycine to obtain environment densities for simulating electronic excitations in clusters.

The basic idea of frozen-density embedding theory (FDET) is the constrained minimization of the Hohenberg–Kohn density functional E^HK[ρ] performed using the auxiliary functional , where Ψ_A is the embedded N_A-electron wavefunction and ρ_B(r) is a non-negative function in real space integrating to a given number of electrons N_B. This choice of independent variables in the total energy functional makes it possible to treat the corresponding two components of the total density using different methods in multi-level simulations. The application of FDET using ρ_B(r) reconstructed from X-ray diffraction data for a molecular crystal is demonstrated for the first time. For eight hydrogen-bonded clusters involving a chromophore (represented as Ψ_A) and the glycylglycine molecule [represented as ρ_B(r)], FDET is used to derive excitation energies. It is shown that experimental densities are suitable for use as ρ_B(r) in FDET-based simulations.

Solvation Structures and Deactivation Pathways of Luminescent Isothiazole-Derived Nucleobases: tz, tz, and tz

The photophysical relaxation pathways of ^tzA, ^tzG, and ^tzI luminescent nucleobases were investigated with the MS-CASPT2 quantum-chemical method and double-ζ basis sets (cc-pVDZ) in gas and condensed phases (1,4-dioxane and water) with the sequential Monte Carlo/CASPT2 and free energy gradient (FEG) methods. Solvation shell structures, in the ground and excited states, were examined with the pairwise radial distribution function (G(r)) and solute-solvent hydrogen-bond networks. Site-specific hydrogen bonding analysis evidenced relevant changes between both electronic states. The three luminescent nucleobases share a common photophysical pattern, summarized as the lowest-lying ¹(ππ*) bright state that is populated directly after the absorption of radiation and evolves barrierless to the minimum energy structure, from where the excess of energy is released by fluorescence. From the ¹(ππ*)_min region, the conical intersection with the ground state ((ππ*/GS)_CI) is not accessible due to the presence of high energetic barriers. By combining the present results with those reported earlier by us for the pyrimidine fluorescent nucleobases, we present a comprehensive description of the photophysical properties of this important class of new fluorescent nucleosides.

Ultrafast Dynamics of the Two Lowest Bright Excited States of Cytosine and 1-Methylcytosine: A Quantum Dynamical Study

The nonadiabatic quantum dynamics (QD) of cytosine and 1-methylcytosine in the gas phase is simulated for 250 fs after a photoexcitation to one of the first two bright states. The nuclear wavepacket is propagated on the coupled diabatic potential energy surfaces of the lowest seven excited states, including ππ*, nπ*, and Rydberg states along all the vibrational degrees of freedom. We focus in particular on the interplay between the bright and the dark nπ* states, not considering the decay to the ground electronic state. To run these simulations, we implemented an automatic general procedure to parametrize linear vibronic coupling (LVC) models with time-dependent density functional theory (DFT) computations and interfaced it with Gaussian package. The wavepacket was propagated with the multilayer version of the multiconfigurational time dependent Hartree method. Two different density functionals, PBE0 and CAM-B3LYP, which provide a different description of the relative stability of the lowest energy dark states, were used to parametrize the LVC Hamiltonian. Part of the photoexcited population on lowest HOMO-LUMO transition (π_Hπ_L^*) decays within less than 100 fs to a nπ* state which mainly involves a promotion of an electron from the oxygen lone pair to the LUMO (n_Oπ_L^*). The population of the second ππ* state decays almost completely, in <100 fs, not only to π_Hπ_L^* and to n_Oπ_L^* states but also to another nπ_L^* state involving the nitrogen lone pair. The efficiency of the adopted protocol allowed us to check the accuracy of the predictions by repeating the QD simulations with different LVC Hamiltonians parametrized either at the ground-state minimum or at stationary structures of different relevant excited states.

Observation of Enhanced Dissociative Photochemistry in the Non-Native Nucleobase 2-Thiouracil

We present the first study to measure the dissociative photochemistry of 2-thiouracil (2-TU), an important nucleobase analogue with applications in molecular biology and pharmacology. Laser photodissociation spectroscopy is applied to the deprotonated and protonated forms of 2-TU, which are produced in the gas-phase using electrospray ionization mass spectrometry. Our results show that the deprotonated form of 2-thiouracil ([2-TU-H]⁻) decays predominantly by electron ejection and hence concomitant production of the [2-TU-H]· free-radical species, following photoexcitation across the UVA-UVC region. Thiocyanate (SCN⁻) and a m/z 93 fragment ion are also observed as photodecay products of [2-TU-H]⁻ but at very low intensities. Photoexcitation of protonated 2-thiouracil ([2-TU·H]⁺) across the same UVA-UVC spectral region produces the m/z 96 cationic fragment as the major photofragment. This ion corresponds to ejection of an HS· radical from the precursor ion and is determined to be a product of direct excited state decay. Fragment ions associated with decay of the hot ground state (i.e., the ions we would expect to observe if 2-thiouracil was behaving like UV-dissipating uracil) are observed as much more minor products. This behaviour is consistent with enhanced intersystem crossing to triplet excited states compared to internal conversion back to the ground state. These are the first experiments to probe the effect of protonation/deprotonation on thionucleobase photochemistry, and hence explore the effect of pH at a molecular level on their photophysical properties.

Analytic gradients for state-averaged multiconfiguration pair-density functional theory

Analytic gradients are important for efficient calculations of stationary points on potential energy surfaces, for interpreting spectroscopic observations, and for efficient direct dynamics simulations. For excited electronic states, as are involved in UV-Vis spectroscopy and photochemistry, analytic gradients are readily available and often affordable for calculations using a state-averaged complete active space self-consistent-field (SA-CASSCF) wave function. However, in most cases, a post-SA-CASSCF step is necessary for quantitative accuracy, and such calculations are often too expensive if carried out by perturbation theory or configuration interaction. In this work, we present the analytic gradients for multiconfiguration pair-density functional theory based on SA-CASSCF wave functions, which is a more affordable alternative. A test set of molecules has been studied with this method, and the stationary geometries and energetics are compared to values in the literature as obtained by other methods. Excited-state geometries computed with state-averaged pair-density functional theory have similar accuracy to those from complete active space perturbation theory at the second-order.

Rethinking Causation for Data-intensive Biology: Constraints, Cancellations, and Quantized Organisms: Causality in complex organisms is sculpted by constraints rather than instigators, with outcomes perhaps better described by quantized patterns than rectilinear pathways

Complex organisms thwart the simple rectilinear causality paradigm of "necessary and sufficient," with its experimental strategy of "knock down and overexpress." This Essay organizes the eccentricities of biology into four categories that call for new mathematical approaches; recaps for the biologist the philosopher's recent refinements to the causation concept and the mathematician's computational tools that handle some but not all of the biological eccentricities; and describes overlooked insights that make causal properties of physical hierarchies such as emergence and downward causation straightforward. Reviewing and extrapolating from similar situations in physics, it is suggested that new mathematical tools for causation analysis incorporating feedback, signal cancellation, nonlinear dependencies, physical hierarchies, and fixed constraints rather than instigative changes will reveal unconventional biological behaviors. These include "eigenisms," organisms that are limited to quantized states; trajectories that steer a system such as an evolving species toward optimal states; and medical control via distributed "sheets" rather than single control points.

Simultaneous observation of nuclear and electronic dynamics by ultrafast electron diffraction

Simultaneous observation of nuclear and electronic motion is crucial for a complete understanding of molecular dynamics in excited electronic states. It is challenging for a single experiment to independently follow both electronic and nuclear dynamics at the same time. Here we show that ultrafast electron diffraction can be used to simultaneously record both electronic and nuclear dynamics in isolated pyridine molecules, naturally disentangling the two components. Electronic state changes (S₁→S₀ internal conversion) were reflected by a strong transient signal in small-angle inelastic scattering, and nuclear structural changes (ring puckering) were monitored by large-angle elastic diffraction. Supported by ab initio nonadiabatic molecular dynamics and diffraction simulations, our experiment provides a clear view of the interplay between electronic and nuclear dynamics of the photoexcited pyridine molecule.

Selenium substitution effects on excited-state properties and photophysics of uracil: a MS-CASPT2 study

The photophysics of selenium-substituted nucleobases has attracted recent experimental attention because they could serve as potential photosensitizers in photodynamic therapy. Herein, we present a comprehensive MS-CASPT2 study on the spectroscopic and excited-state properties, and photophysics of 2-selenouracil (2SeU), 4-selenouracil (4SeU), and 2,4-selenouracil (24SeU). Relevant minima, conical intersections, crossing points, and excited-state relaxation paths in the lowest five electronic states (i.e., S₀, S₁, S₂, T₂, and T₁) are explored. On the basis of these results, their photophysical mechanisms are proposed. Upon photoirradiation to the bright S₂ state, 2SeU quickly relaxes to its S₂ minimum and then moves in an essentially barrierless way to a nearby S₂/S₁ conical intersection near which the S₁ state is populated. Next, the S₁ system arrives at an S₁/T₂/T₁ intersection where a large S₁/T₁ spin-orbit coupling of 430.8 cm^-1 makes the T₁ state populated. In this state, a barrier of 6.8 kcal mol^-1 will trap 2SeU for a while. In parallel, for 4SeU or 24SeU, the system first relaxes to the S₂ minimum and then overcomes a small barrier to approach an S₂/S₁ conical intersection. Once hopping to the S₁ state, there exists an extended region with very close S₁, T₂, and T₁ energies. Similarly, a large S₁/T₁ spin-orbit coupling of 426.8 cm^-1 drives the S₁→ T₁ intersystem crossing process thereby making the T₁ state populated. Similarly, an energy barrier heavily suppresses electronic transition to the S₀ state. The present work manifests that different selenium substitutions on uracil can lead to a certain extent of different vertical and adiabatic excitation energies, excited-state properties, and relaxation pathways. These insights could help understand the photophysics of selenium-substituted nucleobases.

TURBOMOLE: Modular program suite for ab initio quantum-chemical and condensed-matter simulations

TURBOMOLE is a collaborative, multi-national software development project aiming to provide highly efficient and stable computational tools for quantum chemical simulations of molecules, clusters, periodic systems, and solutions. The TURBOMOLE software suite is optimized for widely available, inexpensive, and resource-efficient hardware such as multi-core workstations and small computer clusters. TURBOMOLE specializes in electronic structure methods with outstanding accuracy-cost ratio, such as density functional theory including local hybrids and the random phase approximation (RPA), GW-Bethe-Salpeter methods, second-order Møller-Plesset theory, and explicitly correlated coupled-cluster methods. TURBOMOLE is based on Gaussian basis sets and has been pivotal for the development of many fast and low-scaling algorithms in the past three decades, such as integral-direct methods, fast multipole methods, the resolution-of-the-identity approximation, imaginary frequency integration, Laplace transform, and pair natural orbital methods. This review focuses on recent additions to TURBOMOLE's functionality, including excited-state methods, RPA and Green's function methods, relativistic approaches, high-order molecular properties, solvation effects, and periodic systems. A variety of illustrative applications along with accuracy and timing data are discussed. Moreover, available interfaces to users as well as other software are summarized. TURBOMOLE's current licensing, distribution, and support model are discussed, and an overview of TURBOMOLE's development workflow is provided. Challenges such as communication and outreach, software infrastructure, and funding are highlighted.

Locating Cytosine Conical Intersections by Laser Experiments and Ab Initio Calculations

The decay mechanism of S₀ → S₁ excited cytosine (Cyt) and the effect of substitution are studied combining jet-cooled spectroscopy (nanosecond resonant two-photon ionization (R2PI) and picosecond lifetime measurements) with CASPT2//CASSCF computations for eight derivatives. For Cyt and five derivatives substituted at N1, C5, and C6, rapid internal conversion sets in at 250-1200 cm^-1 above the 0₀⁰ bands. The break-off in the spectra correlates with the calculated barriers toward the "C5-C6 twist" conical intersection, which unambiguously establishes the decay mechanism at low S₁ state vibrational energies. The barriers increase with substituents that stabilize the charge shifts at C4, C5, and C6 following (¹ππ*) excitation. The R2PI spectra of the clamped derivatives 5,6-trimethyleneCyt (TMCyt) and 1-methyl-TMCyt (1M-TMCyt), which decay along an N3 out-of-plane coordinate, extend up to +3500 and +4500 cm^-1.

O-H vibrational motions promote sub-50 fs nonadiabatic dynamics in 3-hydroxypyran-4-one: interplay between internal conversion and ESIPT

A theoretical study is used to explore the involvement of O-H vibrational motions in the S0 → S2 photoinduced dynamics of 3-hydroxypyran-4-one (3-HOX). Two transitions, S0 → S1 and S0 → S2, are attributed to the experimentally observed electronic absorption spectral features in the range of 3.5-5.5 eV. We compute model potential energy surfaces of vibronically coupled S1 (nπ*) and S2 (ππ*) states with the aid of extensive electronic structure calculations. The S1-S2 conical intersection is characterized in the O-H bend and O-H stretch vibrational coordinate space. Quantum wavepacket dynamics simulations reveal an ultrafast S2 → S1 internal conversion decay, where about 90% of the S2 population disappears within the first 50 fs of the propagation time. The participation of O-H vibrational motions in the early events of nonadiabatic dynamics is analyzed based on the time evolution of nuclear densities on S2. We discuss the implications of these observations to provide fundamental insights into the nonadiabatic excited-state intramolecular proton transfer in 3-HOX and its derivatives.