Multiple Decay Mechanisms and 2D-UV Spectroscopic Fingerprints of Singlet Excited Solvated Adenine-Uracil Monophosphate

Galván I, Rivalta I, Aquilante F, Garavelli M, Lindh R, Segarra-Martí J. Characterizing Conical Intersections in DNA/RNA Nucleobases with Multiconfigurational Wave Functions of Varying Active Space Size

We characterize the photochemically relevant conical intersections between the lowest-lying accessible electronic excited states of the different DNA/RNA nucleobases using Cholesky decomposition-based complete active space self-consistent field (CASSCF) algorithms. We benchmark two different basis set contractions and several active spaces for each nucleobase and conical intersection type, measuring for the first time how active space size affects conical intersection topographies in these systems and the potential implications these may have toward their description of photoinduced phenomena. Our results show that conical intersection topographies are highly sensitive to the electron correlation included in the model: by changing the amount (and type) of correlated orbitals, conical intersection topographies vastly change, and the changes observed do not follow any converging pattern toward the topographies obtained with the largest and most correlated active spaces. Comparison across systems shows analogous topographies for almost all intersections mediating population transfer to the dark 1nO/Nπ* states, while no similarities are observed for the "ethylene-like" conical intersection ascribed to mediate the ultrafast decay component to the ground state in all DNA/RNA nucleobases. Basis set size seems to have a minor effect, appearing to be relevant only for purine-based derivatives. We rule out structural changes as a key factor in classifying the different conical intersections, which display almost identical geometries across active space and basis set change, and we highlight instead the importance of correctly describing the electronic states involved at these crossing points. Our work shows that careful active space selection is essential to accurately describe conical intersection topographies and therefore to adequately account for their active role in molecular photochemistry.

Coherent Two-Dimensional and Broadband Electronic Spectroscopies

Over the past few decades, coherent broadband spectroscopy has been widely used to improve our understanding of ultrafast processes (e.g., photoinduced electron transfer, proton transfer, and proton-coupled electron transfer reactions) at femtosecond resolution. The advances in femtosecond laser technology along with the development of nonlinear multidimensional spectroscopy enabled further insights into ultrafast energy transfer and carrier relaxation processes in complex biological and material systems. New discoveries and interpretations have led to improved design principles for optimizing the photophysical properties of various artificial systems. In this review, we first provide a detailed theoretical framework of both coherent broadband and two-dimensional electronic spectroscopy (2DES). We then discuss a selection of experimental approaches and considerations of 2DES along with best practices for data processing and analysis. Finally, we review several examples where coherent broadband and 2DES were employed to reveal mechanisms of photoinitiated ultrafast processes in molecular, biological, and material systems. We end the review with a brief perspective on the future of the experimental techniques themselves and their potential to answer an even greater range of scientific questions.

Modelling Photoionisation in Isocytosine: Potential Formation of Longer-Lived Excited State Cations in its Keto Form

Studying the effects of UV and VUV radiation on non-canonical DNA/RNA nucleobases allows us to compare how they release excess energy following absorption with respect to their canonical counterparts. This has attracted much research attention in recent years because of its likely influence on the origin of our genetic lexicon in prebiotic times. Here we present a CASSCF and XMS-CASPT2 theoretical study of the photoionisation of non-canonical pyrimidine nucleobase isocytosine in both its keto and enol tautomeric forms. We analyse their lowest energy cationic excited states including 2 π + , 2 n O + and 2 n N + and compare these to the corresponding electronic states in cytosine. Investigating lower-energy decay pathways we find - unexpectedly - that keto-isocytosine+ presents a sizeable energy barrier potentially inhibiting decay to its cationic ground state, whereas enol-isocytosine+ features a barrierless and consequently ultrafast pathway analogous to the one previously found for the canonical (keto) form of cytosine+ . Dynamic electron correlation reduces the energy barrier in the keto form substantially (by ∼1 eV) but it is nevertheless still present. We additionally compute the UV/Vis absorption signals of the structures encountered along these decay channels to provide spectroscopic fingerprints to assist future experiments in monitoring these intricate photo-processes.

Ultrafast Spectroscopy of Photoactive Molecular Systems from First Principles: Where We Stand Today and Where We Are Going

Computational spectroscopy is becoming a mandatory tool for the interpretation of the complex, and often congested, spectral maps delivered by modern non-linear multi-pulse techniques. The fields of Electronic Structure Methods, Non-Adiabatic Molecular Dynamics, and Theoretical Spectroscopy represent the three pillars of the virtual ultrafast optical spectrometer, able to deliver transient spectra in silico from first principles. A successful simulation strategy requires a synergistic approach that balances between the three fields, each one having its very own challenges and bottlenecks. The aim of this Perspective is to demonstrate that, despite these challenges, an impressive agreement between theory and experiment is achievable now regarding the modeling of ultrafast photoinduced processes in complex molecular architectures. Beyond that, some key recent developments in the three fields are presented that we believe will have major impacts on spectroscopic simulations in the very near future. Potential directions of development, pending challenges, and rising opportunities are illustrated.

Exploring ultraviolet photoinduced charge-transfer dynamics in a model dinucleotide of guanine and thymine

An understanding of the initial photoexcited states of DNA is essential to unravelling deleterious photoinduced chemical reactions and the intrinsic ultrafast photoprotection of the genetic code for all life. In our combined experimental and theoretical study, we have elucidated the primary non-radiative relaxation dynamics of a model nucleotide of guanine and thymine (2'-deoxyguanosine 3'-monophosphate 5'-thymidine, d(GpT)) in buffered aqueous solution. Experimentally, we unequivocally demonstrate that the Franck-Condon excited states of d(GpT) are significantly delocalised across both nucleobases, and mediate d(G⁺pT^-) exciplex product formation on an ultrafast (<350 fs) timescale. Theoretical studies show that the nature of the vertical excited states is very dependent on the specific geometry of the dinucleotide, and dictate the degree of delocalised, charge-transfer or localised character. Our mechanism for prompt exciplex formation involves a rapid change in electronic structure and includes a diabatic surface crossing very close to the Franck-Condon region mediating fast d(G⁺pT^-) formation. Exciplexes are quickly converted back to neutral ground state molecules on a ∼10 ps timescale with a high quantum yield, ensuring the photostability of the nucleotide sequence.

Ultrafast photoinduced energy and charge transfer

After presenting the basic theoretical models of excitation energy transfer and charge transfer, I describe some of the novel experimental methods used to probe them. Finally, I discuss recent results concerning ultrafast energy and charge transfer in biological systems, in chemical systems and in photovoltaics based on sensitized transition metal oxides.

Ultrafast and radiationless electronic excited state decay of uracil and thymine cations: computing the effects of dynamic electron correlation

In this article we characterise the radiationless decay of the first few electronic excited states of the cations of DNA/RNA nucleobases uracil and thymine, including the effects of dynamic electron correlation on energies and geometries (optimised with XMS-CASPT2). In both systems, we find that one state of ²n and another two of ²π⁺ character can be populated following photoionisation, and their different minima and interstate crossings are located. We find strong similarities between uracil and thymine cations: with accessible conical intersections suggesting that depopulation of their electronic excited states takes place on ultrafast timescales in both systems, suggesting that they are photostable in agreement with previous theoretical (uracil⁺) evidence. We find that dynamic electron correlation separates the energy levels of the "3-state" conical intersection (D₂/D₁/D₀)_CI previously located with CASSCF for uracil⁺, which will therefore have a different geometry and higher energy. Simulating the electronic and vibrational absorptions allows us to characterise spectral fingerprints that could be used to monitor these cation photo-processes experimentally.

Two-dimensional UV spectroscopy: a new insight into the structure and dynamics of biomolecules

Two-dimensional ultraviolet spectroscopy has the potential to deliver rich structural and dynamical information on biomolecules such as DNA and proteins.

Similar chemical structures, dissimilar triplet quantum yields: a CASPT2 model rationalizing the trend of triplet quantum yields in nitroaromatic systems

S₁/S₀ accessibility strongly influences the triplet quantum yields of nitronaphthalenes.

UV-induced long-lived decays in solvated pyrimidine nucleosides resolved at the MS-CASPT2/MM level

The most relevant 'dark' electronic excited states in DNA/RNA pyrimidine nucleosides are mapped in water employing hybrid MS-CASPT2/MM optimisations with explicit solvation and including the sugar. Conical intersections (CIs) between initially accessed bright ¹ππ* and the lowest energy dark ¹nπ* excited states, involving the lone pair localised on the oxygen and/or nitrogen atoms are characterised. They are found in the vicinities of the Franck-Condon (FC) region and are shown to facilitate non-adiabatic population transfer. The excited state population of the ¹n_Oπ* state, localised in the carbonyl moiety on all pyrimidine nucleosides, is predicted to rapidly evolve to its minimum, displaying non-negligible potential energy barriers along its non-radiative decay, and accounting for the ps signal registered in pump-probe experiments as well as for an efficient population of the triplet state. Cytidine displays an additional ¹n_Nπ* state localised in the N3 atom and that leads to its excited state minimum displaying large potential energy barriers in the pathway connecting to the CI with the ground state. Sugar-to-base hydrogen/proton transfer processes are assessed in solution for the first time, displaying a sizable barrier along its decay and thus being competitive with other slow decay channels in the ps and ns timescales. A unified deactivation scheme for the long-lived channels of pyrimidine nucleosides is delivered, where the ¹n_Oπ* state is found to mediate the long-lived decay in the singlet manifold and act as the doorway for triplet population and thus accounting for the recorded phosphorescence and, more generally, for the transient/photoelectron spectral signals registered up to the ns timescale.

The effect of solvent relaxation in the ultrafast time-resolved spectroscopy of solvated benzophenone

Modeling time-resolved spectra to unravel ultra fast solvent reorganization.

Two-dimensional electronic spectroscopy as a tool for tracking molecular conformations in DNA/RNA aggregates

A computational strategy to simulate two-dimensional electronic spectra (2DES) is introduced, which allows characterising ground state conformations of flexible nucleobase aggregates that play a crucial role in nucleic acid photochemistry.

On the Simulation of Two-dimensional Electronic Spectroscopy of Indole-containing Peptides

A benchmark study of low-cost multiconfigurational CASSCF/CASPT2 schemes for computing the electronic structure of indole is presented. This facilitates the simulation of near-ultraviolet (UV) pump visible (VIS) probe (i.e. two-color) two-dimensional electronic spectra (2DES) of homo- and hetero-aggregates as well as for processing of multiple snapshots from molecular dynamics simulations. Fingerprint excited-state absorption signatures of indole are identified in a broad spectral window between 10 and 25 k cm^-1 . The 18-24 k cm^-1 spectral window which has no absorption of the monomer and noninteracting aggregates is ideally suited to embed charge-transfer signatures in stacked aggregates. The small peptide Trp-cage, containing a tryptophan and a tyrosine amino acids, having indole and phenol as side chains, respectively, serves to prove the concept. Clear charge-transfer signatures are found in the proposed spectral window for an interchromophore distance of 5 Å making near-UV pump VIS probe 2DES a suitable technique for resolving closely packed aggregates. We demonstrate that 2DES utilizing ultra-short pulses has the potential to resolve the nature of the spectroscopically resolved electronic states and that the line shapes of the excited-state absorption signals can be correlated to the polarity of the relevant states.

Linear photoabsorption spectra and vertical excitation energies of microsolvated DNA nucleobases in aqueous solution

Employing density functional theory (DFT) and time-dependent density functional theory (TDDFT) calculations in combination with the semiclassical nuclear ensemble method, we have simulated the photoabsorption spectra of the four canonical DNA nucleobases in aqueous solution. In order to model the effects of solvation, for each nucleobase, a number of solvating water molecules were explicitly included in the simulations, and additionally, the bulk solvent was represented by a continuous polarizable medium. We find that the effect of the solvation shell in general is significant, and its inclusion improves the realism of the spectral simulations. The involvement of lone electron pairs in the hydrogen bonding with the solvating water molecules has the effect of systematically increasing the energies of vertical excitation into the [Formula: see text]-type states. Apart from a systematic blue shift of around [Formula: see text][Formula: see text]eV observed in the absorption peaks, the calculated photoabsorption spectra reproduce the measured ones with good accuracy. The photoabsorption spectra are dominated by excited states with [Formula: see text] and partial [Formula: see text] character. No low-energy charge transfer states are observed with the use of the CAM-B3LYP and M06-2X functionals.

Insights into the Thermal Eliminations and Photoeliminations of B,N-Heterocycles: A Theoretical Study

Understanding the photochemistry of organoboron compounds is essential to expand optoelectronic applications. In this work, the complete active space self-consistent field (CASSCF) and its second-order perturbation (CASPT2) methods combining with density functional theory (DFT) have been employed to investigate the elimination mechanisms of compound 6,7-dihydro-5⁴-benzo[d]pyrido[2,1-f][1,2]azaborininr (B4) on the ground state (S₀) and the first excited state (S₁). B4 is one of the 1,2-B,N-heterocycles that undergo competitive thermal elimination and photoelimination depending on the substitution groups on the B atom and the chelate backbone, thus providing a high-selectivity synthesis strategy for luminescent compounds. Since the energy barrier from B4 to BH₃-pyrido[1,2-a]isoindole (D1) and pyrido[1,2-a]isoindole (A1) on the ground state is lower than that from B4 to 5⁴-benzo[d]pyrido[2,1-f][1,2]azaborininr (C4), the retraction ring reaction is expected to proceed with larger probability than the H₂ elimination upon heating. On the contrary, photoelimination of H₂ may take place easily due to the almost barrierless pathway on the S₁ state. Remarkably, we have located an energetically available conical intersection (S₁/S₀)_X-1, which allows for ultrafast S₁ → S₀ decay and subsequently generation of C4. Our results not only throw light on the experimental observations of the selectivity of thermal elimination and photoelimination but also provide detailed information on the excited state as instructional implications for further synthesis and application of B,N-embedded aromatics.

Challenges in Simulating Light-Induced Processes in DNA

In this contribution, we give a perspective on the main challenges in performing theoretical simulations of photoinduced phenomena within DNA and its molecular building blocks. We distinguish the different tasks that should be involved in the simulation of a complete DNA strand subject to UV irradiation: (i) stationary quantum chemical computations; (ii) the explicit description of the initial excitation of DNA with light; (iii) modeling the nonadiabatic excited state dynamics; (iv) simulation of the detected experimental observable; and (v) the subsequent analysis of the respective results. We succinctly describe the methods that are currently employed in each of these steps. While for each of them, there are different approaches with different degrees of accuracy, no feasible method exists to tackle all problems at once. Depending on the technique or combination of several ones, it can be problematic to describe the stacking of nucleobases, bond breaking and formation, quantum interferences and tunneling or even simply to characterize the involved wavefunctions. It is therefore argued that more method development and/or the combination of different techniques are urgently required. It is essential also to exercise these new developments in further studies on DNA and subsystems thereof, ideally comprising simulations of all of the different components that occur in the corresponding experiments.

Two-dimensional electronic spectroscopy in the ultraviolet by a birefringent delay line

We introduce a 2D electronic spectroscopy setup in the UV spectral range in the partially collinear pump-probe geometry. The required interferometrically phase-locked few-optical-cycle UV pulse pair is generated by combining a passive birefringent interferometer in the visible and nonlinear phase transfer. This is achieved by sum-frequency generation between the phase-locked visible pulse pair and narrowband infrared pulses. We demonstrate a pair of 16-fs, 330-nm pulses whose delay is interferometrically stable with an accuracy better than λ/450. 2DUV maps of pyrene solution probed in the UV and visible spectral ranges are demonstrated.

Assessment of the Potential Energy Hypersurfaces in Thymine within Multiconfigurational Theory: CASSCF vs. CASPT2

The present study provides new insights into the topography of the potential energy hypersurfaces (PEHs) of the thymine nucleobase in order to rationalize its main ultrafast photochemical decay paths by employing two methodologies based on the complete active space self-consistent field (CASSCF) and the complete active space second-order perturbation theory (CASPT2) methods: (i) CASSCF optimized structures and energies corrected with the CASPT2 method at the CASSCF geometries and (ii) CASPT2 optimized geometries and energies. A direct comparison between these strategies is drawn, yielding qualitatively similar results within a static framework. A number of analyses are performed to assess the accuracy of these different computational strategies under study based on a variety of numerical thresholds and optimization methods. Several basis sets and active spaces have also been calibrated to understand to what extent they can influence the resulting geometries and subsequent interpretation of the photochemical decay channels. The study shows small discrepancies between CASSCF and CASPT2 PEHs, displaying a shallow planar or twisted ¹(ππ*) minimum, respectively, and thus featuring a qualitatively similar scenario for supporting the ultrafast bi-exponential deactivation registered in thymine upon UV-light exposure. A deeper knowledge of the PEHs at different levels of theory provides useful insight into its correct characterization and subsequent interpretation of the experimental observations. The discrepancies displayed by the different methods studied here are then discussed and framed within their potential consequences in on-the-fly non-adiabatic molecular dynamics simulations, where qualitatively diverse outcomes are expected.