Electronic Structure of Adenine and Thymine Base Pairs Studied by Femtosecond Electron−Ion Coincidence Spectroscopy

Toward Consistent Predictions of Core/Valence Ionization Potentials and Valence Excitation Energies by MRSF-TDDFT

Optimizing exchange-correlation functionals for both core/valence ionization potentials (cIPs/vIPs) and valence excitation energies (VEEs) at the same time in the framework of MRSF-TDDFT is self-contradictory. To overcome the challenge, within the previous "adaptive exact exchange" or double-tuning strategy on Coulomb-attenuating XC functionals (CAM), a new XC functional specifically for cIPs and vIPs was first developed by enhancing exact exchange to both short- and long-range regions. The resulting DTCAM-XI functional achieved remarkably high accuracy in its predictions with errors of less than half eV. An additional concept of "valence attenuation", where the amount of exact exchange for the frontier orbital regions is selectively suppressed, was introduced to consistently predict both VEEs and IPs at the same time. The second functional, DTCAM-XIV, exhibits consistent overall prediction accuracy at ∼0.64 eV. By preferentially optimizing VEEs within the same "valence attenuation" concept, a third functional, DTCAM-VAEE, was obtained, which exhibits improved performance as compared to that of the previous DTCAM-VEE and DTCAM-AEE in the prediction of VEEs, making it an attractive alternative to BH&HLYP. As the combination of "adaptive exchange" and "valence attenuation" is operative, it would be exciting to explore its potential with a more tunable framework in the future.

Radiation-Induced Transfer of Charge, Atoms, and Energy within Isolated Biomolecular Systems

In biological tissues, ionizing radiation interacts with a variety of molecules and the consequences include cell killing and the modification of mechanical properties. Applications of biological radiation action are for instance radiotherapy, sterilization, or the tailoring of biomaterial properties. During the first femtoseconds to milliseconds after the initial radiation action, biomolecular systems typically respond by transfer of charge, atoms, or energy. In the condensed phase, it is usually very difficult to distinguish direct effects from indirect effects. A straightforward solution for this problem is the use of gas-phase techniques, for instance from the field of mass spectrometry. In this review, we survey mainly experimental but also theoretical work, focusing on radiation-induced intra- and inter-molecular transfer of charge, atoms, and energy within biomolecular systems in the gas phase. Building blocks of DNA, proteins, and saccharides, but also antibiotics are considered. The emergence of general processes as well as their timescales and mechanisms are highlighted.

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

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

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

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

Photoionization Dynamics and Proton Transfer within the Adenine-Thymine Nucleobase Pair

Studying the stability of hydrogen-bonded nucleobase pairs, at the heart of the genetic code, is of utmost importance for an in-depth understanding of basic mechanisms of life and biomolecular evolution. We present here a VUV single photon ionization dynamic study of the nucleobase pair adenine-thymine (AT), revealing its ionization and dissociative ionization thresholds via double imaging electron/ion coincidence spectroscopy. The experimental data, consisting of cluster mass-resolved threshold photoelectron spectra and photon energy-dependent ion kinetic energy release distributions, allow the unambiguous distinction of the dissociation of AT into protonated adenine AH⁺ and a dehydrogenated thymine radical T(-H) from dissociative ionization processes of other nucleobase clusters. Comparison to high-level ab initio calculations indicates that our experimental observations can be explained by a single hydrogen-bonded conformer present in our molecular beam and allows the estimation of an upper limit of the barrier of the proton transfer in the ionized AT pair.

Photochemistry of Thymine in Solution and DNA Revealed by an Electrostatic Embedding QM/MM Combined with Mixed-Reference Spin-Flip TDDFT

The photochemistry of nucleobases, important for their role as building blocks of DNA, is largely affected by the electrostatic environment in which they are soaked. For example, despite the numerous studies of thymine in solution and DNA, there is still a debate on the photochemical deactivation pathways after UV absorption. Many theoretical models are oversimplified due to the lack of computationally accurate and efficient electronic structure methodologies that capture excited state electron correlation effects when nucleobases are embedded in large electrostatic media. Here, we combine mixed-reference spin-flip time-dependent density functional theory (MRSF-TDDFT) with electrostatic embedding QM/MM using electrostatic potential fittingfitted (ESPF) atomic charges, as a strategy to accurately and efficiently describe the electronic structure of chromophores polarized by an electrostatic medium. In particular, we develop analytic expressions for the energy and gradient of MRSF/MM based on the ESPF coupling using atom-centered grids and total charge conservation. We apply this methodology to the study of solvation effects on thymine photochemistry in water and thymine dimers in DNA. In the former, the combination of trajectory surface hopping (TSH) nonadiabatic molecular dynamics (NAMD) with MRSF/MM remarkably revealed accelerated deactivation decay pathways, which is consistent with the experimental decay time of ∼400 fs. The enhanced hopping rate can be explained by the preferential stabilization of corresponding conical interactions due to their increased dipole moments. Structurally, it is a consequence of characteristic methyl puckered geometries near the conical intersection region. For the thymine dimer in B-DNA, we found new photochemical pathways through conical intersections that could explain the formation of cyclobutadiene dimers and 6-4 photoproducts.

The Ultrafast Quantum Dynamics of Photoexcited Adenine-Thymine Basepair Investigated with a Fragment-based Diabatization and a Linear Vibronic Coupling Model

In this contribution we present a quantum dynamical study of the photoexcited hydrogen bonded base pair adenine–thymine (AT) in a Watson–Crick arrangement. To that end, we parametrize Linear Vibronic Coupling (LVC) models with Time-Dependent Density Functional Theory (TD-DFT) calculations, exploiting a fragment diabatization scheme (FrD) we have developed to define diabatic states on the basis of individual chromophores in a multichromophoric system. Wavepacket propagations were run with the multilayer extension of the Multiconfiguration Time-Dependent Hartree method. We considered excitations to the three lowest bright states, a ππ* state of thymine and two ππ* states (L_a and L_b) of adenine, and we found that on the 100 fs time scale the main decay pathways involve intramonomer population transfers toward nπ* states of the same nucleobase. In AT this transfer is less effective than in the isolated nucleobases, because hydrogen bonding destabilizes the nπ* states. The population transfer to the A → T charge transfer state is negligible, making the ultrafast (femtosecond) decay through the proton coupled electron transfer mechanism unlikely, in line with experimental results in apolar solvents. The excitation energy transfer is also very small. We carefully compare the predictions of LVC Hamiltonians obtained with different sets of diabatic states, defined so to match either local states of the two separated monomers or the base pair adiabatic states in the Franck–Condon region. To that end we also extend the flexibility of the FrD-LVC approach, introducing a new strategy to define fragments diabatic states that account for the effect of the rest of the multichromohoric system through a Molecular Mechanics potential.

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.

Excited state dynamics of aniline homoclusters

We have investigated the relaxation, following excitation in the 290-235 nm region, of neutral aniline homoclusters (An)n formed in a supersonic expansion by femtosecond time resolved ionization. The applied method permits isolation of the dynamics of the dimer from that originated in bigger species of the generated distribution. Interestingly, and differently from the monomer and (An)n≥3 clusters, the dimer does not present a N-H dissociative 1πσ* channel. This fact can be explained in terms of the symmetric structure adopted, in which each molecule establishes two N-Hπ interactions, destabilizing the H dissociation channel. The observations permit relating the photophysics to the interactions established by the aniline units and confirming previous observations and theoretical predictions on the structure of aniline aggregates.

Internal conversion and intersystem crossing pathways in UV excited, isolated uracils and their implications in prebiotic chemistry

The photodynamic properties of molecules determine their ability to survive in harsh radiation environments. As such, the photostability of heterocyclic aromatic compounds to electromagnetic radiation is expected to have been one of the selection pressures influencing the prebiotic chemistry on early Earth. In the present study, the gas-phase photodynamics of uracil, 5-methyluracil (thymine) and 2-thiouracil-three heterocyclic compounds thought to be present during this era-are assessed in the context of their recently proposed intersystem crossing pathways that compete with internal conversion to the ground state. Specifically, time-resolved photoelectron spectroscopy measurements evidence femtosecond to picosecond timescales for relaxation of the singlet (1)ππ* and (1)nπ* states as well as for intersystem crossing to the triplet manifold. Trapping in the excited triplet state and intersystem crossing back to the ground state are investigated as potential factors contributing to the susceptibility of these molecules to ultraviolet photodamage.

New Insights into the State Trapping of UV-Excited Thymine

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

Femtosecond Excited State Dynamics of Size Selected Neutral Molecular Clusters

The work describes a novel experimental approach to track the relaxation dynamics of an electronically excited distribution of neutral molecular clusters formed in a supersonic expansion, by pump-probe femtosecond ionization. The introduced method overcomes fragmentation issues and makes possible to retrieve the dynamical signature of a particular cluster from each mass channel, by associating it to an IR transition of the targeted structure. We have applied the technique to study the nonadiabatic relaxation of pyrrole homoclusters. The results obtained exciting at 243 nm, near the origin of the bare pyrrole electronic absorption, allow us to identify the dynamical signature of the dimer (Py)2, which exhibits a distinctive lifetime of τ1 ∼ 270 fs, considerably longer than the decays recorded for the monomer and bigger size clusters (Py)n>2. A possible relationship between the measured lifetime and the clusters geometries is tentatively discussed.

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

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

Photoinduced processes in nucleic acids

Photoinduced processes in nucleic acids are phenomena of fundamental interest in diverse fields, from prebiotic studies, through medical research on carcinogenesis, to the development of bioorganic photodevices. In this contribution we survey many aspects of the research across the boundaries. Starting from a historical background, where the main milestones are identified, we review the main findings of the physical-chemical research of photoinduced processes on several types of nucleic-acid fragments, from monomers to duplexes. We also discuss a number of different issues which are still under debate.

A new double imaging velocity focusing coincidence experiment: i2PEPICO

The vacuum ultraviolet (VUV) beamline of the Swiss Light Source has been upgraded after two years of operation. A new, turntable-type monochromator was constructed at the Paul Scherrer Institut, which allows for fast yaw-alignment as well as quick grating change and exchange. In addition to the original imaging photoelectron photoion coincidence endstation (iPEPICO), a second, complementary double imaging setup (i(2)PEPICO) has been built. Volatile samples can be introduced at room temperature or in a molecular beam, a pyrolysis source allows for radical production, and non-volatile solids can be evaporated in a heated cell. Monochromatic VUV radiation ionizes the sample and both photoelectrons and photoions are velocity map imaged onto two fast position sensitive detectors and detected in delayed coincidence. High intensity synchrotron radiation leads to ionization rates above 10(5) s(-1). New data acquisition and processing approaches are discussed for recording coincidence processes at high rates. The setup is capable of resolving pulsed molecular beam profiles and the synchrotron time structure temporally. The latter is shown by photoelectron autocorrelation, which displays both the 1.04 MHz ring clock frequency as well as resolving the micro-pulses with a separation of 2 ns. Kinetic energy release analysis on the dissociative photoionization of CF(4) indicates a dissociation mechanism change in the Franck-Condon allowed energy range of the first ion state.

Photoreaction channels of the guanine-cytosine base pair explored by long-range corrected TDDFT calculations

Photoinduced processes in the Watson-Crick guanine-cytosine base pair are comprehensively studied by means of long-range corrected (LC) TDDFT calculations of potential energy profiles using the LC-BLYP and CAM-B3LYP functionals. The ab initio CC2 method and the conventional TDDFT method with the B3LYP functional are also employed to assess the reliability of the LC-TDDFT method. The present approach allows us to compare the potential energy profiles at the same computational level for excited-state reactions of the base pair, including single and double proton transfer between the bases and nonradiative decay via ring puckering in each base. In particular, long-range correction to the TDDFT method is critical for a qualitatively correct description of the proton transfer reactions. The calculated energy profiles exhibit low barriers for out-of-plane deformation of the guanine moiety in the locally-excited state, which is expected to lead to a conical intersection with the ground state, as well as for single proton transfer from guanine to cytosine with the well-known electron-driven proton transfer mechanism. Thus the present results suggest that both processes can compete in hydrogen-bonded base pairs and play a significant role in the mechanism of photostability.

Specific photodynamics in thymine clusters: the role of hydrogen bonding

A photoionization detected IR study of thymine and 1-methylthymine monohydrates and of their homodimers was carried out to shed some light on the structure of the thymine clusters whose complex photodynamics has recently been the subject of great interest. Under supersonic jet conditions, thymine forms doubly H-bonded cyclic clusters with water or another base preferentially via its N1-H group and the adjacent carbonyl group. This hydrate is of no biological relevance since the N1-H group is the sugar binding site in thymidine. On the other hand, 1-methylthymine forms the donor H-bonds only via the N3-H group. Hence, properties of the N1-H and the N3-H bound clusters of thymine can be studied using thymine and 1-methylthymine molecules, respectively. No biologically relevant conformations of the dimers and hydrates of thymine, contrary to those of 1-methylthymine, are observed under supersonic jet conditions. Thymine homodimer, which extensively fragments upon UV ionization by formation of a protonated monomer, exhibits two N1-H···O═C2 hydrogen bonds. The photodynamics of hydrated thymines is found to be extremely sensitive to the hydration site: ranging from an ultrafast relaxation in less than 100 fs up to formation of a dark state with the lifetime on the microsecond time scale.

Ultrafast relaxation dynamics observed through time-resolved photoelectron angular distributions

Time-resolved photoelectron imaging of the 7,7,8,8-tetracyanoquinodimethane (TCNQ) radical anion is presented. Photoelectron angular distributions (PADs) are qualitatively analyzed in terms of the simple s-p model that is based on symmetry arguments. The internal conversion dynamics from the first excited state (1(2)B(3u)) to the ground state ((2)B(2g)) may be observed through temporal changes in the PADs of the spectrally overlapping photoelectron features arising from photodetachment of the ground state and the excited state. A formulism for extracting the population dynamics from the β(2) anisotropy parameter of overlapping spectroscopic features is presented. This is used to extract the lifetime of the first excited state, which is in good agreement with that observed in the time-resolved photoelectron spectra.