Toward Accurate Quantum Chemical Methods for Molecules of Increasing Dimension: The New Family of Pisa Composite Schemes

The new versions of the Pisa composite scheme introduced in the present paper are based on the careful selection of different quantum chemical models for energies, geometries, and vibrational frequencies, with the aim of maximizing the accuracy of the overall description while retaining a reasonable cost for all the steps. In particular, the computation of accurate electronic energies has been further improved introducing more reliable complete basis set extrapolations and estimation of core-valence correlation, together with improved basis sets for third-row atoms. Furthermore, the reduced-cost frozen natural orbital (FNO) model has been introduced and validated for large molecules. Accurate molecular structures can be obtained avoiding complete basis set extrapolation and evaluating core-valence correlation at the MP2 level. Unfortunately, analytical gradients are not available for the FNO version of the model. Therefore, for large molecules, an accurate reduced-cost alternative is offered by evaluation of valence contributions with a double-hybrid functional in conjunction with the same MP2 contribution for core-valence correlation or by means of a one-parameter approximation. The same double-hybrid functional and basis set are employed to evaluate zero-point energies and partition functions. After the validation of the new models for small systems, a panel of molecular bricks of life has been used to analyze their performances for problems of current fundamental or technological interest. The fully black-box implementation of the computational workflow paves the way toward the accurate yet not prohibitively expensive study of medium- to large-sized molecules also by experimentally oriented researchers.

Quantum chemistry meets high-resolution spectroscopy for characterizing the molecular bricks of life in the gas-phase

Computation of accurate geometrical structures and spectroscopic properties of large flexible molecules in the gas-phase is tackled at an affordable cost using a general exploration/exploitation strategy. The most distinctive feature of the approach is the careful selection of different quantum chemical models for energies, geometries and vibrational frequencies with the aim of maximizing the accuracy of the overall description while retaining a reasonable cost for all the steps. In particular, a composite wave-function method is used for energies, whereas a double-hybrid functional (with the addition of core-valence correlation) is employed for geometries and harmonic frequencies and a cheaper hybrid functional for anharmonic contributions. A thorough benchmark based on a wide range of prototypical molecular bricks of life shows that the proposed strategy is close to the accuracy of state-of-the-art composite wave-function methods, and is applicable to much larger systems. A freely available web-utility post-processes the geometries optimized by standard electronic structure codes paving the way toward the accurate yet not prohibitively expensive study of medium- to large-sized molecules by experimentally-oriented researchers.

On Prototropy and Bond Length Alternation in Neutral and Ionized Pyrimidine Bases and Their Model Azines in Vacuo

In this review, the complete tautomeric equilibria are derived for disubstituted pyrimidine nucleic acid bases starting from phenol, aniline, and their model compounds-monosubstituted aromatic azines. The differences in tautomeric preferences for isolated (gaseous) neutral pyrimidine bases and their model compounds are discussed in light of different functional groups, their positions within the six-membered ring, electronic effects, and intramolecular interactions. For the discussion of tautomeric preferences and for the analysis of internal effects, recent quantum-chemical results are taken into account and compared to some experimental ones. For each possible tautomer-rotamer of the title compounds, the bond length alternation, measured by means of the harmonic oscillator model of electron delocalization (HOMED) index, is examined. Significant HOMED similarities exist for mono- and disubstituted derivatives. The lack of parallelism between the geometric (HOMED) and energetic (ΔG) parameters for all possible isomers clearly shows that aromaticity is not the main factor that dictates tautomeric preferences for pyrimidine bases, particularly for uracil and thymine. The effects of one-electron loss (positive ionization) and one-electron gain (negative ionization) on prototropy and bond length alternation are also reviewed for pyrimidine bases and their models.

Accurate Thermochemical and Kinetic Parameters at Affordable Cost by Means of the Pisa Composite Scheme (PCS)

A new strategy for the computation at an affordable cost of geometrical structures, thermochemical parameters, and rate constants for medium-sized molecules in the gas phase is proposed. The most distinctive features of the new model are the systematic use of cc-pVnZ-F12 basis sets, the addition of MP2 core-valence correlation in geometry optimizations by a double-hybrid functional, the separate extrapolation of MP2 and post-MP2 contributions, and the inclusion of anharmonic contributions in zero-point energies and thermodynamic functions. A thorough benchmark based on a wide range of prototypical systems shows that the new scheme outperforms the most well-known model chemistries without the need for any empirical parameter. Additional tests show that the computed zero-point energies and thermal contributions can be confidently used for obtaining accurate thermochemical and kinetic parameters. Since the whole computational workflow is translated in a black-box procedure, which can be followed with standard electronic structure codes, the way is paved for the accurate yet not prohibitively expensive study of medium- to large-sized molecules also by nonspecialists.

Experimental and theoretical study of the cytosine tautomerism through excited states

CONTEXT

The irradiation of water solution of cytosine with UV light (λmax = 254 nm) shows oxo-hydroxy tautomerism with a rate constant of 6.297 × 10-3 min-1. The order of the reaction implies a tautomeric conversion. After removing the UV light source, we observed a dark reaction with a rate constant of 1.473 × 10-3 min-1 which leads to a restoration of the initial tautomer as before the irradiation. The mechanism of oxo-hydroxy tautomerism of cytosine in water solution was studied in the excited state. It was found that the transformations occur along the 1πσ* excited-state reaction paths which link the Franck-Condon geometries of the tautomers and the conical intersections S0/S1 connected with the H-detachment processes of the corresponding bonds. Furthermore, we established that the conical intersections S0/S1 are also mutually accessible along the 1πσ* excited-state reaction paths.

METHODS

The ground-state equilibrium geometries were optimized at the B3LYP/aug-cc-pVDZ level of theory in water environment according to PCM as well as at the CC2/aug-cc-pVDZ level in the gas phase. The TD B3LYP and CC2 methods were applied for the study of the excited states. The tautomerization mechanisms were studied with the use of the linear interpolation in internal coordinates approach using the optimized geometries of tautomers minima and conical intersections S0/S1 at the CASSCF(6,6)/6-31G* level. All calculations were performed with the GAUSSIAN 16 commercial software.

DFT Meets Wave-Function Composite Methods for Characterizing Cytosine Tautomers in the Gas Phase

A general strategy for the accurate computation of structural and spectroscopic properties of biomolecule building blocks in the gas phase has been further improved and validated with a special reference to tautomeric equilibria. The main improvements concern the use of the cc-pVTZ-F12 basis set in both DFT and CCSD(T)-F12 computations, the inclusion of core-valence correlation in geometry optimizations by double hybrid functionals, and the use of the cc-pVQZ-F12 basis set for complete basis set extrapolation at the MP2-F12 level. The resulting model chemistry is applied to the challenging problem of cytosine tautomers in the gas phase. The results are in remarkable agreement with experiment concerning both rotational and vibrational spectroscopic parameters and permit their unbiased interpretation in terms of structural and thermochemical features. Together with the intrinsic interest of the studied molecule, the accuracy of the results obtained at reasonable cost without any empirical parameter suggests that the proposed composite method can be profitably employed for accurate investigations of other molecular bricks of life.

Characterization of tautomeric forms of anti-cancer drug gemcitabine and their interconversion upon mechano-chemical treatment, using ATR-FTIR spectroscopy and complementary methods

Gemcitabine is a widely used anti-cancer drug of pyrimidine structure, which can exist as a free base molecular form in crystals. Tautomers are structural isomers of molecules, which interconvert via proton transfer. Mechano-chemistry studies reactions of solids under mechanical impact. We investigated gemcitabine free base for the presence of specific molecular tautomers, using ATR-FTIR spectroscopic analysis, powder XRD, optical microscopy and HPLC. The amino-keto tautomer has the characteristic infrared (IR) peak of the amino group at 3390 cm-1. For the first time, the imino-keto tautomer of gemcitabine free base was detected. The imino-keto tautomer has the characteristic IR peak of the =N-H group, and its peak due to the CO group in pyrimidine ring is shifted vs. that of the amino-keto tautomer. This serves as the unique spectroscopic "fingerprints" of these tautomers. The ATR-FTIR spectroscopic analysis shows that gemcitabine free base can be enriched with the amino-keto or the imino-keto tautomer. Further, we studied the transformation of gemcitabine free base in crystals between its tautomers under conditions of liquid-assisted grinding (LAG). The imino-keto tautomer undergoes tautomerization to the amino-keto tautomer, while the amino-keto tautomer remains stable. No destruction of molecules of gemcitabine free base, when present as either tautomer, occurs during LAG as was verified by the HPLC-UV analysis. LAG is a new, straightforward, facile and fast method to interconvert tautomers in crystals, and ATR-FTIR spectroscopy is a method of choice to study tautomerization reactions of pharmaceuticals. The presented approach is promising for analysis of crystals of drugs containing one or more than one tautomer, and the knowledge-driven design of composite materials, which contain specific tautomeric molecular forms of pyrimidines, purines and other biologically active heterocyclic compounds.

The Resonance Raman Spectrum of Cytosine in Water: Analysis of the Effect of Specific Solute-Solvent Interactions and Non-Adiabatic Couplings

In this contribution, we report a computational study of the vibrational Resonance Raman (vRR) spectra of cytosine in water, on the grounds of potential energy surfaces (PES) computed by time-dependent density functional theory (TD-DFT) and CAM-B3LYP and PBE0 functionals. Cytosine is interesting because it is characterized by several close-lying and coupled electronic states, challenging the approach commonly used to compute the vRR for systems where the excitation frequency is in quasi-resonance with a single state. We adopt two recently developed time-dependent approaches, based either on quantum dynamical numerical propagations of vibronic wavepackets on coupled PES or on analytical correlation functions for cases in which inter-state couplings were neglected. In this way, we compute the vRR spectra, considering the quasi-resonance with the eight lowest-energy excited states, disentangling the role of their inter-state couplings from the mere interference of their different contributions to the transition polarizability. We show that these effects are only moderate in the excitation energy range explored by experiments, where the spectral patterns can be rationalized from the simple analysis of displacements of the equilibrium positions along the different states. Conversely, at higher energies, interference and inter-state couplings play a major role, and the adoption of a fully non-adiabatic approach is strongly recommended. We also investigate the effect of specific solute-solvent interactions on the vRR spectra, by considering a cluster of cytosine, hydrogen-bonded by six water molecules, and embedded in a polarizable continuum. We show that their inclusion remarkably improves the agreement with the experiments, mainly altering the composition of the normal modes, in terms of internal valence coordinates. We also document cases, mostly for low-frequency modes, in which a cluster model is not sufficient, and more elaborate mixed quantum classical approaches, in explicit solvent models, need to be applied.

Prototropically Controlled Dynamics of Cytosine Photodecay

The effect of the existence of several prototropic tautomers of cytosine on its UV/vis spectra and the excited state decay dynamics is studied by spectral and nonadiabatic molecular dynamics (NAMD) simulations in connection with the mixed-reference spin-flip time-dependent density functional theory (MRSF-TDDFT) method. Simulated UV/vis spectra provide a strong indication that the H3N keto-amino cytosine tautomer (the least anticipated species) may be present under experimental conditions. The NAMD simulations yield a wide range of excited state decay constants for various tautomers of cytosine, ranging from ∼1.3 ps for the biologically relevant H1N keto-amino tautomer to ∼0.1 ps for the keto-imino tautomer. The slowness of the H1N decay dynamics follows from the presence of a barrier on the excited state energy surface separating the Franck-Condon structure from the major decay funnel, the conical intersection seam. It is suggested that the experimentally observed photodecay dynamics may result from a combination of the decay processes of various tautomers (H3N in particular) present simultaneously under the experimental conditions.

Non-radiative deactivation of excited cytosine: probing of different DFT functionals and basis sets in solvents with different polarity

The vertical excitation and emission energies of the cytosine oxo-amino form were calculated at the TD DFT level of theory with several functionals, basis sets, and solvents with different polarity (PCM). They were compared with the experimental UV absorption and fluorescence maxima, revealing that the minimal deviation of the vertical excitation energies from the of UV absorption maxima can be achieved when the hybrid functional B3LYP is applied within the calculations. Regular correlations like curves with saturation between the vertical excitation/emission energies and the dielectric constants were registered. The relaxation of the ¹ππ^* excited state through an ethylene-like conical intersection S₀/S₁ should occur through decrease of the decay rate (commented qualitatively here) with the rise of the solvent polarity.

Excited state absorption of DNA bases in the gas phase and in chloroform solution: a comparative quantum mechanical study

We study the excited state absorption (ESA) properties of the four DNA bases (thymine, cytosine, adenine, and guanine) by different single reference quantum mechanical methods, namely, equation of motion coupled cluster singles and doubles (EOM-CCSD), singles, doubles and perturbative triples (EOM-CC3), and time-dependent density functional theory (TD-DFT), with the long-range corrected CAM-B3LYP functional. Preliminary results at the Tamm-Dancoff (TDA) CAM-B3LYP level using the maximum overlap method (MOM) are reported for thymine. In the gas phase, the three methods predict similar One Photon Absorption (OPA) spectra, which are consistent with the experimental results and with the most accurate computational studies available in the literature. The ESA spectra are then computed for the ππ* states (one for pyrimidine, two for purines) associated with the lowest-energy absorption band, and for the close-lying nπ* state. The EOM-CC3, EOM-CCSD and CAM-B3LYP methods provide similar ESA spectral patterns, which are also in qualitative agreement with literature RASPT2 results. Once validated in the gas phase, TD-CAM-B3LYP has been used to compute the ESA in chloroform, including solvent effects by the polarizable continuum model (PCM). The predicted OPA and ESA spectra in chloroform are very similar to those in the gas phase, most of the bands shifting by less than 0.1 eV, with a small increase of the intensities and a moderate destabilization of the nπ* state. Finally, ESA spectra have been computed from the minima of the lowest energy ππ* state, and found in line with the available experimental transient absorption spectra of the nucleosides in solution, providing further validation of our computational approach.

do Casal M, Pinheiro M, Toldo JM, Barbatti M. Modeling the heating and cooling of a chromophore after photoexcitation

The heating of a chromophore due to internal conversion and its cooling down due to energy dissipation to the solvent are crucial phenomena to characterize molecular photoprocesses. In this work, we simulated the ab initio nonadiabatic dynamics of cytosine, a prototypical chromophore undergoing ultrafast internal conversion, in three solvents—argon matrix, benzene, and water—spanning an extensive range of interactions. We implemented an analytical energy-transfer model to analyze these data and extract heating and cooling times. The model accounts for nonadiabatic effects, and excited- and ground-state energy transfer, and can analyze data from any dataset containing kinetic energy as a function of time. Cytosine heats up in the subpicosecond scale and cools down within 25, 4, and 1.3 ps in argon, benzene, and water, respectively. The time constants reveal that a significant fraction of the benzene and water heating occurs while cytosine is still electronically excited.

An analytical energy-transfer model is implemented to obtain a chromophore's heating and cooling times in a given solvent by using quantities available in nonadiabatic dynamics simulations.

Cryogenic IR and UV spectroscopy of isomer-selected cytosine radical cation

The UV photodissociation of cryogenic-cooled isomer-selected cytosine–silver complex leads to the production of cytosine radical cation without isomerization.

Nucleobases thin films deposited on nanostructured transparent conductive electrodes for optoelectronic applications

Environmentally-friendly bio-organic materials have become the centre of recent developments in organic electronics, while a suitable interfacial modification is a prerequisite for future applications. In the context of researches on low cost and biodegradable resource for optoelectronics applications, the influence of a 2D nanostructured transparent conductive electrode on the morphological, structural, optical and electrical properties of nucleobases (adenine, guanine, cytosine, thymine and uracil) thin films obtained by thermal evaporation was analysed. The 2D array of nanostructures has been developed in a polymeric layer on glass substrate using a high throughput and low cost technique, UV-Nanoimprint Lithography. The indium tin oxide electrode was grown on both nanostructured and flat substrate and the properties of the heterostructures built on these two types of electrodes were analysed by comparison. We report that the organic-electrode interface modification by nano-patterning affects both the optical (transmission and emission) properties by multiple reflections on the walls of nanostructures and the electrical properties by the effect on the organic/electrode contact area and charge carrier pathway through electrodes. These results encourage the potential application of the nucleobases thin films deposited on nanostructured conductive electrode in green optoelectronic devices.

The thermochemical, structural, and spectroscopic analyses of the tautomers of sulfur and selenium modified emissive nucleobases

In this study, density functional theory (DFT) and time dependent density functional theory (TD-DFT) are used to investigate the stabilities and spectral properties [IR, UV–vis, and two-photon absorption (2PA)] of two sets of modified RNA nucleobase tautomers. The modifications introduce either a sulfur or selenium atom to form an isothiazolo[4,3-d]pyrimidine or isoselenazolo[4,3-d]pyrimidine heterocyclic core, respectively. The relative stabilities of both sets of modified tautomers determined with B3LYP/6-31++G(d,p) reveal that in water (with a polarizable continuum model), the 6-keto-2-amino tautomer of guanine and the rare 4-imino-2-keto tautomer of cytosine may be present at significant populations, whereas the 6-enol-2-amino tautomer of guanine is more common in the gas phase. The identification of these modified tautomers is possible due to the natural differences in their vibrational modes and hence IR spectra. Furthermore, the photophysical properties of both these sets of modified tautomers indicate that excitation and emission energies are shifted relative to their more abundant form in both one photon absorption and emission and 2PA spectra, as determined at the B3LYP/6-31++G(d,p) and CAM-B3LYP/aug-cc-pVDZ levels of theory, respectively. Even though the 2PA cross sections in water for all of the species are small (0.3–2.3 GM), the modified cytosine tautomer shows promise, as its cross section is larger than the more dominant form. The spectral separation between the dominant form and the tautomers of isoselenazole and isothiazole modified cytosine and guanine are relatively similar, suggesting both modifications could be useful in elucidating the tautomers from their more abundant counterparts.

First-principles characterization of the singlet excited state manifold in DNA/RNA nucleobases

TD-DFT characterization of the high-energy singlet excited state manifold of the canonical DNA/RNA nucleobasesin vacuumis assessed against RASPT2 reference computations for reliable simulations of linear and non-linear electronic spectra.

Ultrafast Excited-State Dynamics of Hydrogen-Bonded Cytosine Microsolvated Clusters with Protic and Aprotic Polar Solvents

Microsolvation effects on the ultrafast excited-state deactivation dynamics of cytosine (Cy) were studied in hydrogen-bonded Cy clusters with protic and aprotic solvents using mass-resolved femtosecond pump-probe ionization spectroscopy. Two protic solvents, water (H₂O) and methanol (MeOH), and one aprotic solvent, tetrahydrofuran (THF), were investigated, and transients of Cy·(H₂O)_1-6, Cy·(MeOH)_1-3, and Cy·THF microsolvated clusters produced in supersonic expansions were measured. With the aid of electronic structure calculations, we assigned the observed dynamics to the low-energy isomers of various Cy clusters and discussed the microsolvation effect on the excited-state deactivation dynamics. With the protic solvents only the microsolvated clusters of Cy keto tautomer were observed. The observed decay time constants of Cy·(H₂O) _n are 0.5 ps for n = 1 and ∼0.2-0.25 ps for n = 2-6. For Cy·(MeOH) _n clusters, the decay time constant for n = 1 cluster is similar to that of the Cy monohydrate, but for n = 2 and 3 the decays are about a factor of 2 slower than the corresponding microhydrates. With the aprotic solvent, THF, hydrogen-bonded complexes of both keto and enol tautomers are present in the beam. The keto-Cy·THF shows a decay similar to that of the keto-Cy monomer, whereas the enol-Cy·THF exhibits a 2-fold slower decay than the enol-Cy monomer, suggesting an increase in the barrier to excited-state deactivation upon binding of one THF molecule to the enol form of Cy.

Photoexcitation of iodide ion-pyrimidine clusters above the electron detachment threshold: Intracluster electron transfer versus nucleobase-centred excitations

Laser photodissociation spectroscopy of the I^-·thymine (I^-·T) and I^-·cytosine (I^-·C) nucleobase clusters has been conducted for the first time across the regions above the electron detachment thresholds to explore the excited states and photodissociation channels. Although photodepletion is strong, only weak ionic photofragment signals are observed, indicating that the clusters decay predominantly by electron detachment. The photodepletion spectra of the I^-·T and I^-·C clusters display a prominent dipole-bound excited state (I) in the vicinity of the vertical detachment energy (∼4.0 eV). Like the previously studied I^-·uracil (I^-·U) cluster [W. L. Li et al., J. Chem. Phys. 145, 044319 (2016)], the I^-·T cluster also displays a second excited state (II) centred at 4.8 eV, which we similarly assign to a π-π* nucleobase-localized transition. However, no distinct higher-energy absorption bands are evident in the spectra of the I^-·C. Time-dependent density functional theory (TDDFT) calculations are presented, showing that while each of the I^-·T and I^-·U clusters displays a single dominant π-π* nucleobase-localized transition, the corresponding π-π* nucleobase transitions for I^-·C are split across three separate weaker electronic excitations. I^- and deprotonated nucleobase anion photofragments are observed upon photoexcitation of both I^-·U and I^-·T, with the action spectra showing bands (at 4.0 and 4.8 eV) for both the I^- and deprotonated nucleobase anion production. The photofragmentation behaviour of the I^-·C cluster is distinctive as its I^- photofragment displays a relatively flat profile above the expected vertical detachment energy. We discuss the observed photofragmentation profiles of the I^-·pyrimidine clusters, in the context of the previous time-resolved measurements, and conclude that the observed photoexcitations are primarily consistent with intracluster electron transfer dominating in the near-threshold region, while nucleobase-centred excitations dominate close to 4.8 eV. TDDFT calculations suggest that charge-transfer transitions [Iodide n (5p⁶) → Uracil σ*] may contribute to the cluster absorption profile across the scanned spectral region, and the possible role of these states is also discussed.

A 'bottom up', ab initio computational approach to understanding fundamental photophysical processes in nitrogen containing heterocycles, DNA bases and base pairs

The availability of non-radiative decay mechanisms by which photoexcited molecules can revert to their ground electronic state, without experiencing potentially deleterious chemical transformation, is fundamental to molecular photostability. This Perspective Article combines results of new ab initio electronic structure calculations and prior experimental data in an effort to systematise trends in the non-radiative decay following UV excitation of selected families of heterocyclic molecules. We start with the prototypical uni- and bicyclic molecules phenol and indole, and explore the structural and photophysical consequences of incorporating progressively more nitrogen atoms within the respective ring structures en route to the DNA bases thymine, cytosine, adenine and guanine. For each of the latter, we identify low energy non-radiative decay pathways via conical intersections with the ground state potential energy surface accessed by out-of-plane ring deformations. This is followed by summary descriptions and illustrations of selected rival (electron driven H atom transfer) non-radiative excited state decay processes that demand consideration once the nucleobases are merely components in larger biomolecular systems like nucleosides, and both individual and stacked base-pairs.

Cooperativity effect involving drug-DNA/RNA intermolecular interaction: A B3LYP-D3 and MP2 theoretical investigation on ketoprofen⋯cytosine⋯H2O system

In order to examine the origin of the drug action and design new DNA/RNA-targeted drugs, the cooperativity effect involving drug-DNA/RNA intermolecular interaction in ketoprofen⋯cytosine⋯H₂O ternary system were investigated by the B3LYP, B3LYP-D3, and MP2 methods with the 6-311++G(2d,p) basis set. The thermodynamic cooperativity was also evaluated at 310.15 K. The N-H⋯O, O-H⋯O, O-H⋯N, C-H⋯N, and C-H⋯O H bonds coexist in ternary complexes. The intermolecular interactions obtained by B3LYP-D3 are close to those calculated by MP2. The steric effects and van der Waals interactions have little influence on the cooperativity effects. The anti-cooperativity effect in ket⋯cyt⋯H₂O is far more notable than the cooperativity effect, and the stability of the cyclic structure with anti-cooperativity effect is higher than that of the linear structure with cooperativity effect, as is confirmed by the AIM (atoms in molecules) and RDG (reduced density gradient) analysis. Thus, it can be inferred that, in the presence of H₂O, the anti-cooperativity effect plays a dominant role in the drug-DNA/RNA interaction, and the nature of the hydration in the binding of drugs to DNA/RNA bases is the H-bonding anti-cooperativity effect. Furthermore, the drug always links simultaneously with DNA/RNA base and H₂O, and only in this way can the biological activity of drugs play a role. In most cases, the enthalpy change is the major factor driving the cooperativity, as is different from most of biomacromolecule complexes.

Solid State Separation and Isolation of Tautomers of Fused-Ring Triazolotriazoles

Fine control of the tautomeric forms of [1,2,4]triazolo[3,2-c][1,2,4]triazole derivatives in acidic conditions has been achieved by acting on the electronic character of the substituent at position 7 of the heterobicycle and on the counterion. Strong electron releasing or electron withdrawing substituents lead almost exclusively to a single tautomeric form, the 1H-3H or the 2H-3H, respectively. In the case of the phenol substituent, both tautomeric forms are present in comparable amount in solution; the two tautomers can also be selectively precipitated in different crystalline salts using suitable counterions.

Revealing Deactivation Pathways Hidden in Time-Resolved Photoelectron Spectra

Time-resolved photoelectron spectroscopy is commonly employed with the intention to monitor electronic excited-state dynamics occurring in a neutral molecule. With the help of theory, we show that when excited-state processes occur on similar time scales the different relaxation pathways are completely obscured in the total photoionization signal recorded in the experiment. Using non-adiabatic molecular dynamics and Dyson norms, we calculate the photoionization signal of cytosine and disentangle the transient contributions originating from the different deactivation pathways of its tautomers. In the simulations, the total signal from the relevant keto and enol tautomers can be decomposed into contributions either from the neutral electronic state populations or from the distinct mechanistic pathways across the multiple potential surfaces. The lifetimes corresponding to these contributions cannot be extracted from the experiment, thereby illustrating that new experimental setups are necessary to unravel the intricate non-adiabatic pathways occurring in polyatomic molecules after irradiation by light.

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

The photophysics and photochemistry of DNA is of great importance due to the potential damage of the genetic code by UV light. Quantum mechanical studies have played a key role in interpretating the results of modern time-resolved pump-probe spectroscopy, and in elucidating the main photoactivated reactive paths. This review provides a concise, complete picture of the computational studies carried out, approximately, in the past decade. We start with an overview of the photophysics of the nucleobases in the gas phase and in solution. We discuss the proposed mechanisms for ultrafast decay to the ground state, that involve conical intersections, consider the role of triplet states, and analyze how the solvent modulates the photophysics. Then we move to larger systems, from dinucleotides to single- and double-stranded oligonucleotides. We focus on the possible role of charge transfer and delocalized or excitonic states in the photophysics of these systems and discuss the main photochemical paths. We finish with an outlook on the current challenges in the field and future directions of research.

Excitation Spectra of Nucleobases with Multiconfigurational Density Functional Theory

Range-separated hybrid methods between wave function theory and density functional theory (DFT) can provide high-accuracy results, while correcting some of the inherent flaws of both the underlying wave function theory and DFT. We here assess the accuracy for excitation energies of the nucleobases thymine, uracil, cytosine, and adenine, using a hybrid between complete active space self-consistent field (CASSCF) and DFT methods. The method is based on range separation, thereby avoiding all double-counting of electron correlation and is denoted long-range CASSCF short-range DFT (CAS-srDFT). Using a linear response extension of CAS-srDFT, we compare the first 7-8 excited states of the nucleobases with perturbative multireference approaches as well as coupled cluster based methods. Our results show that the CAS-srDFT method can provide accurate excitation energies in good correspondence with the computationally more expensive methods.