Electronic Structure and Spectroscopy of Nucleic Acid Bases: Ionization Energies, Ionization-Induced Structural Changes, and Photoelectron Spectra

Benchmarking Ionization Potentials from pCCD Tailored Coupled Cluster Models

The ionization potential (IP) is an important parameter providing essential insights into the reactivity of chemical systems. IPs are also crucial for designing, optimizing, and understanding the functionality of modern technological devices. We recently showed that limiting the CC ansatz to the seniority-zero sector proves insufficient in predicting reliable and accurate ionization potentials within an IP equation-of-motion coupled-cluster formalism. Specifically, the absence of dynamical correlation in the seniority-zero pair coupled cluster doubles (pCCD) model led to unacceptably significant errors of approximately 1.5 eV. In this work, we aim to explore the impact of dynamical correlation and the choice of the molecular orbital basis (canonical vs localized) in CC-type methods targeting 230 ionized states in 70 molecules, comprising small organic molecules, medium-sized organic acceptors, and nucleobases. We focus on pCCD-based approaches as well as the conventional IP-EOM-CCD and IP-EOM-CCSD. Their performance is compared to the CCSD(T) or CCSDT equivalent and experimental reference data. Our statistical analysis reveals that all investigated frozen-pair coupled cluster methods exhibit similar performance, with differences in errors typically within chemical accuracy (1 kcal/mol or 0.05 eV). Notably, the effect of the molecular orbital basis, such as canonical Hartree-Fock or natural pCCD-optimized orbitals, on the IPs is marginal if dynamical correlation is accounted for. Our study suggests that triple excitations are crucial in achieving chemical accuracy in IPs when modeling electron detachment processes with pCCD-based methods.

Selective Tautomer Production and Cryogenic Ion Spectroscopy of Radical Cations: The Uracil and Thymine Cases

The vibrational and electronic spectroscopy of the radical cations of two nucleobases (NB) (uracil and thymine) was studied by cryogenic ion photodissociation spectroscopy. The radical cations have been generated from the photodissociation of NB-Ag+ complexes. A charge transfer process from the NB to Ag+ governs the deactivation mechanism, leading to the formation of the radical cation without further tautomerization. Single- and double-resonance spectroscopy allows for structural assignments of both the silver complexes and the radical cations by comparison with DFT-based calculations. Interestingly, a tautomer-dependent fragmentation is observed in the thymine enol form that involves the loss of NCO, a fragment which was never reported before for this NB. This selective photodissociation of silver complexes containing aromatic chromophore greatly expands the current technique to produce isomer-selected radical cations in the gas phase providing benchmark experimental data to assess calculations of open-shell species.

Temozolomide-induced guanine mutations create exploitable vulnerabilities of guanine-rich DNA and RNA regions in drug-resistant gliomas

Temozolomide (TMZ) is a chemotherapeutic agent that has been the first-line standard of care for the aggressive brain cancer glioblastoma (GBM) since 2005. Although initially beneficial, TMZ resistance is universal and second-line interventions are an unmet clinical need. Here, we took advantage of the known mechanism of action of TMZ to target guanines (G) and investigated G-rich G-quadruplex (G4) and splice site changes that occur upon TMZ resistance. We report that TMZ-resistant GBM has guanine mutations that disrupt the G-rich DNA G4s and splice sites that lead to deregulated alternative splicing. These alterations create vulnerabilities, which are selectively targeted by either the G4-stabilizing drug TMPyP4 or a novel splicing kinase inhibitor of cdc2-like kinase. Last, we show that the G4 and RNA binding protein EWSR1 aggregates in the cytoplasm in TMZ-resistant GBM cells and patient samples. Together, our findings provide insight into targetable vulnerabilities of TMZ-resistant GBM and present cytoplasmic EWSR1 as a putative biomarker.

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.

Unimolecular isomerisation of 1,5-hexadiyne observed by threshold photoelectron photoion coincidence spectroscopy

The unimolecular isomerisation of the prompt propargyl + propargyl "head-to-head" adduct, 1,5- hexadiyne, to fulvene and benzene by the 3,4-dimethylenecyclobut-1-ene (DMCB) intermediate (all C6H6) was studied in the high-pressure limit...

Differential Charging in Photoemission from Mercurated DNA Monolayers on Ferromagnetic Films

Spin-dependent and enantioselective electron-molecule scattering occurs in photoelectron transmission through chiral molecular films. This spin selectivity leads to electron spin filtering by molecular helices, with increasing magnitude concomitant with increasing numbers of helical turns. Using ultraviolet photoelectron spectroscopy, we measured spin-selective surface charging accompanying photoemission from ferromagnetic substrates functionalized with monolayers of mercurated DNA hairpins that constitute only one helical turn. Mercury ions bind specifically at thymine-thymine mismatches within self-hybridized single-stranded DNA, enabling precise control over the number and position of Hg²⁺ along the helical axis. Differential charging of the organic layers, manifested as substrate-magnetization-dependent photoionization energies, was observed for DNA hairpins containing Hg²⁺; no differences were measured for hairpin monolayers in the absence of Hg²⁺. Inversion of the DNA helical secondary structure at increased metal loading led to complementary inversion in spin selectivity. We attribute these results to increased scattering probabilities from relativistic enhancement of spin-orbit interactions in mercurated DNA.

Electron and hole interactions with P, Z, and P:Z and the formation of mutagenic products by proton transfer reactions

Z would act as an electron acceptor and P would capture a hole in the unnatural DNA. The latter process would produce mutagenic products via a proton transfer reaction.

Valence Photoionization of Thymine: Ionization Energies, Vibrational Structure, and Fragmentation Pathways from the Slow to the Ultrafast

The photoionization of thymine has been studied by using vacuum ultraviolet radiation and imaging photoelectron photoion coincidence spectroscopy after aerosol flash vaporization and bulk evaporation. The two evaporation techniques have been evaluated by comparison of the photoelectron spectra and breakdown diagrams. The adiabatic ionization energies for the first four electronic states were determined to be 8.922±0.008, 9.851±0.008, 10.30±0.02, and 10.82±0.01 eV. Vibrational features have been assigned for the first three electronic states with the help of Franck-Condon factor calculations based on density functional theory and wave function theory vibrational analysis within the harmonic approximation. The breakdown diagram of thymine, as supported by composite method ab initio calculations, suggests that the main fragment ions are formed in sequential HNCO-, CO-, and H-loss dissociation steps from the thymine parent ion, with the first step corresponding to a retro-Diels-Alder reaction. The dissociation rate constants were extracted from the photoion time-of-flight distributions and used together with the breakdown curves to construct a statistical model to determine 0 K appearance energies of 11.15±0.16 and 11.95±0.09 eV for the m/z 83 and 55 fragment ions, respectively. These results have allowed us to revise previously proposed fragmentation mechanisms and to propose a model for the final, nonstatistical H-loss step in the breakdown diagram, yielding the m/z 54 fragment ion at an appearance energy of 13.24 eV.

Approximate Green's Function Coupled Cluster Method Employing Effective Dimension Reduction

The Green's function coupled cluster (GFCC) method, originally proposed in the early 1990s, is a powerful many-body tool for computing and analyzing the electronic structure of molecular and periodic systems, especially when electrons of the system are strongly correlated. However, in order for the GFCC to become a method that may be routinely used in the electronic structure calculations, robust numerical techniques and approximations must be employed to reduce its extremely high computational overhead. In our recent studies, it has been demonstrated that the GFCC equations can be solved directly in the frequency domain using iterative linear solvers, which can be easily distributed in a massively parallel environment. In the present work, we demonstrate a successful application of model-order-reduction (MOR) techniques in the GFCC framework. Briefly speaking, for a frequency regime of interest that requires high-resolution descriptions of spectral function, instead of solving the GFCC linear equation of full dimension for every single frequency point of interest, an efficiently solvable linear system model of a reduced dimension may be built upon projecting the original GFCC linear system onto a subspace. From this reduced order model is obtained a reasonable approximation to the full dimensional GFCC linear equations in both interpolative and extrapolative spectral regions. Here, we show that the subspace can be properly constructed in an iterative manner from the auxiliary vectors of the GFCC linear equations at some selected frequencies within the spectral region of interest. During the iterations, the quality of the subspace, as well as the linear system model, can be systematically improved. The method is tested in this work in terms of the efficiency and accuracy of computing spectral functions for some typical molecular systems such as carbon monoxide, 1,3-butadiene, benzene, and adenine. To reach the same level of accuracy as that of the original GFCC method, the application of MOR in the GFCC method is able to significantly lower the original computational cost for the aforementioned molecules in designated frequency regimes. As a byproduct, the reduced order model obtained by this method is found to provide a high-quality initial guess, which improves the convergence rate for the existing iterative linear solver.

Thermochemistry of Guanine Tautomers Re-Examined by Means of High-Level CCSD(T) Composite Ab Initio Methods

We obtained accurate gas-phase tautomerization energies for a set of 14 guanine tautomers by means of high-level thermochemical procedures approximating the CCSD(T) energy at the complete basis set (CBS) limit. For the five low-lying tautomers, we use the computationally demanding W1-F12 composite method for obtaining the tautomerization energies. The relative W1-F12 tautomerization enthalpies at 298K are: 0.00 (1), 2.37 (2), 2.63 (3), 4.03 (3′), and 14.31 (4) kJmol−1. Thus, as many as four tautomers are found within a small energy window of less than 1.0kcalmol−1 (1kcalmol−1=4.184kJmol−1). We use these highly accurate W1-F12 tautomerization energies to evaluate the performance of a wide range of lower-level composite ab initio procedures. The Gn composite procedures (G4, G4(MP2), G4(MP2)-6X, G3, G3B3, G3(MP2), and G3(MP2)B3) predict that the enol tautomer (3) is more stable than the keto tautomer (2) by amounts ranging from 0.36 (G4) to 1.28 (G3(MP2)) kJmol−1. We also find that an approximated CCSD(T)/CBS energy calculated as HF/jul-cc-pV{D,T}Z+CCSD/jul-cc-pVTZ+(T)/jul-cc-pVDZ results in a root-mean-square deviation (RMSD) of merely 0.11kJmol−1 relative to the W1-F12 reference values. We use this approximated CCSD(T)/CBS method to obtain the tautomerization energies of 14 guanine tautomers. The relative tautomerization enthalpies at 298K are: 0.00 (1), 2.20 (2), 2.51 (3), 4.06 (3′), 14.30 (4), 25.65 (5), 43.78 (4′), 53.50 (6′), 61.58 (6), 77.37 (7), 82.52 (8′), 86.02 (9), 100.70 (10), and 121.01 (8) kJmol−1. Using these tautomerization enthalpies, we evaluate the performance of standard and composite methods for the entire set of 14 guanine tautomers. The best-performing procedures emerge as (RMSDs are given in parentheses): G4(MP2)-6X (0.51), CCSD(T)+ΔMP2/CBS (0.52), and G4(MP2) (0.64kJmol−1). The worst performers are CCSD(T)/AVDZ (1.05), CBS-QB3 (1.24), and CBS-APNO (1.38kJmol−1).

Green's function coupled cluster formulations utilizing extended inner excitations

In this paper, we analyze new approximations of the Green's function coupled cluster (GFCC) method where locations of poles are improved by extending the excitation level of inner auxiliary operators. These new GFCC approximations can be categorized as the GFCC-i(n, m) method, where the excitation level of the inner auxiliary operators (m) used to describe the ionization potential and electron affinity effects in the N - 1 and N + 1 particle spaces is higher than the excitation level (n) used to correlate the ground-state coupled cluster wave function for the N-electron system. Furthermore, we reveal the so-called "n + 1" rule in this category [or the GFCC-i(n, n + 1) method], which states that in order to maintain size-extensivity of the Green's function matrix elements, the excitation level of inner auxiliary operators X _p (ω) and Y _q (ω) cannot exceed n + 1. We also discuss the role of the moments of coupled cluster equations that in a natural way assures these properties. Our implementation in the present study is focused on the first approximation in this GFCC category, i.e., the GFCC-i(2,3) method. As our first practice, we use the GFCC-i(2,3) method to compute the spectral functions for the N₂ and CO molecules in the inner and outer valence regimes. In comparison with the Green's function coupled cluster singles, doubles results, the computed spectral functions from the GFCC-i(2,3) method exhibit better agreement with the experimental results and other theoretical results, particularly in terms of providing higher resolution of satellite peaks and more accurate relative positions of these satellite peaks with respect to the main peak positions.

Systematic High-Accuracy Prediction of Electron Affinities for Biological Quinones

Quinones play vital roles as electron carriers in fundamental biological processes; therefore, the ability to accurately predict their electron affinities is crucial for understanding their properties and function. The increasing availability of cost-effective implementations of correlated wave function methods for both closed-shell and open-shell systems offers an alternative to density functional theory approaches that have traditionally dominated the field despite their shortcomings. Here, we define a benchmark set of quinones with experimentally available electron affinities and evaluate a range of electronic structure methods, setting a target accuracy of 0.1 eV. Among wave function methods, we test various implementations of coupled cluster (CC) theory, including local pair natural orbital (LPNO) approaches to canonical and parameterized CCSD, the domain-based DLPNO approximation, and the equations-of-motion approach for electron affinities, EA-EOM-CCSD. In addition, several variants of canonical, spin-component-scaled, orbital-optimized, and explicitly correlated (F12) Møller-Plesset perturbation theory are benchmarked. Achieving systematically the target level of accuracy is challenging and a composite scheme that combines canonical CCSD(T) with large basis set LPNO-based extrapolation of correlation energy proves to be the most accurate approach. Methods that offer comparable performance are the parameterized LPNO-pCCSD, the DLPNO-CCSD(T₀ ), and the orbital optimized OO-SCS-MP2. Among DFT methods, viable practical alternatives are only the M06 and the double hybrids, but the latter should be employed with caution because of significant basis set sensitivity. A highly accurate yet cost-effective DLPNO-based coupled cluster approach is used to investigate the methoxy conformation effect on the electron affinities of ubiquinones found in photosynthetic bacterial reaction centers. © 2018 Wiley Periodicals, Inc.

Unveiling the complex vibronic structure of the canonical adenine cation

Adenine, a DNA base, exists as several tautomers and isomers that are closely lying in energy and that may form a mixture upon vaporization of solid adenine. Indeed, it is challenging to bring adenine into the gas phase, especially as a unique tautomer. The experimental conditions were tuned to prepare a jet-cooled canonical adenine (9H-adenine). This isolated DNA base was ionized by single VUV photons from a synchrotron beamline and the corresponding slow photoelectron spectrum was compared to ab initio computations of the neutral and ionic species. We report the vibronic structure of the X+ 2A'' (D0), A+ 2A' (D1) and B+ 2A'' (D2) electronic states of the 9H adenine cation, from the adiabatic ionization energy (AIE) up to AIE + 1.8 eV. Accurate AIEs are derived for the 9H-adenine (X[combining tilde] 1A') + hν → 9H-adenine+ (X+ 2A'', A+ 2A', B+ 2A'') + e- transitions. Close to the AIE, we fully assign the rich vibronic structure solely to the 9H-adenine (X 1A') + hν → 9H-adenine+ (X+ 2A'') transition. Importantly, we show that the lowest cationic electronic states of canonical adenine are coupled vibronically. The present findings are important for understanding the effects of ionizing radiation and the charge distribution on this elementary building block of life, at ultrafast, short, and long timescales.

Green's Function Coupled-Cluster Approach: Simulating Photoelectron Spectra for Realistic Molecular Systems

In this paper, we present an efficient implementation for the analytical energy-dependent Green's function coupled-cluster with singles and doubles (GFCCSD) approach with our first practice being computing spectral functions of realistic molecular systems. Because of its algebraic structure, the presented method is highly scalable and is capable of computing spectral function for a given molecular system in any energy region. Several typical examples have been given to demonstrate its capability of computing spectral functions not only in the valence band but also in the core-level energy region. Satellite peaks have been observed in the inner valence band and core-level energy region where a many-body effect becomes significant and the single particle picture of ionization often breaks down. The accuracy test has been carried out by extensively comparing the computed spectral functions by our GFCCSD method with experimental photoelectron spectra as well as the theoretical ionization potentials obtained from other methods. It turns out the GFCCSD method is able to provide a qualitative or semiquantitative level of description of ionization processes in both the core and valence regimes. To significantly improve the GFCCSD results for the main ionic states, a larger basis set can usually be employed, whereas the improvement of the GFCCSD results for the satellite states needs higher-order many-body terms to be included in the GFCC implementation.

Effect of solvation on the ionization of guanine nucleotide: A hybrid QM/EFP study

Ionization of nucleobases is affected by their biological environment, which includes both the effect of adjacent nucleotides as well as the presence of water around it. Guanine and its nucleotide have the lowest ionization potentials among the various DNA bases. Therefore, the threshold of ionization is dependent on that of guanine and its characterization is crucial to the prediction of interaction of light with DNA. We investigate the effect of solvation on the vertical ionization energies (VIEs) of guanine and its nucleotide. In this work, we have used hybrid quantum mechanics/molecular mechanics (QM/MM) approach with effective fragment potential as the MM method of choice and equation-of-motion coupled-cluster for ionization potential with singles and doubles (EOM-IP-CCSD) as the QM method. The performance of the hybrid scheme with respect to the full QM method shows an accuracy of ≤ 0.02-0.04 eV. The lowest few ionizations of the nucleotide are found to be from different parts of the moiety, that is, the nucleic acid base, phosphate, or sugar, and these ionization energies are very closely spaced giving rise to a very complicated spectrum. Furthermore, microsolvation has large effects on these ionizations and can lead to red or blue shift depending on the position of the water molecule. Even a single water molecule can change the order of ionized states in the nucleotide. The VIEs of the bulk solvated chromophores are predicted and compared to existing experimental spectra. The predominant role of polarization in the solvatochromic shift is noticed. © 2017 Wiley Periodicals, Inc.

The effect of sequence on the ionization of guanine in DNA

The accurate estimation of the ionization energies and understanding the nature of the ionized states of the nucleic acid bases (NABs) are crucial to the understanding of the DNA damage mechanism. The vertical ionization energy (VIE) of guanine is the lowest among the NABs and the ionization energies are strongly affected by the environment, such as solvation and characteristics of nearby NABs. Therefore, we investigate the sequence dependence of the VIEs of guanine in B-DNA. We use the equation of motion coupled cluster method for the estimation of ionization potential with single and double excitations (EOM-IP-CCSD) and density functional theory with dispersion corrected ωB97x-D for the estimation of VIEs. A significant amount of non-additivity or cooperativity, directly proportional to charge delocalization, is noticed in the change in VIE due to the interaction with the nearby NABs. While the change in VIE due to base pairing originates predominantly from charge-dipole interactions, stacking between base pairs is a more complicated balance of dispersion and charge-dipole interactions as well as stabilization due to the delocalization of the positive charge. The long range interactions are however dominated by 1/r(3) distance dependence which shows the major role played by charge-dipole interactions. The extent of localization of positive holes on guanine is also estimated for various sequences.

The Importance of Short- and Long-Range Exchange on Various Excited State Properties of DNA Monomers, Stacked Complexes, and Watson-Crick Pairs

We present a detailed analysis of several time-dependent DFT (TD-DFT) methods, including conventional hybrid functionals and two types of nonempirically tuned range-separated functionals, for predicting a diverse set of electronic excitations in DNA nucleobase monomers and dimers. This large and extensive set of excitations comprises a total of 50 different transitions (for each tested DFT functional) that includes several n → π and π → π* valence excitations, long-range charge-transfer excitations, and extended Rydberg transitions (complete with benchmark calculations from high-level EOM-CCSD(T) methods). The presence of localized valence excitations as well as extreme long-range charge-transfer excitations in these systems poses a serious challenge for TD-DFT methods that allows us to assess the importance of both short- and long-range exchange contributions for simultaneously predicting all of these various transitions. In particular, we find that functionals that do not have both short- and full long-range exchange components are unable to predict the different types of nucleobase excitations with the same accuracy. Most importantly, the current study highlights the importance of both short-range exchange and a nonempirically tuned contribution of long-range exchange for accurately predicting the diverse excitations in these challenging nucleobase systems.

Electron accommodation dynamics in the DNA base thymine

The dynamics of electron attachment to the DNA base thymine are investigated using femtosecond time-resolved photoelectron imaging of the gas phase iodide-thymine (I(-)T) complex. An ultraviolet pump pulse ejects an electron from the iodide and prepares an iodine-thymine temporary negative ion that is photodetached with a near-IR probe pulse. The resulting photoelectrons are analyzed with velocity-map imaging. At excitation energies ranging from -120 meV to +90 meV with respect to the vertical detachment energy (VDE) of 4.05 eV for I(-)T, both the dipole-bound and valence-bound negative ions of thymine are observed. A slightly longer rise time for the valence-bound state than the dipole-bound state suggests that some of the dipole-bound anions convert to valence-bound species. No evidence is seen for a dipole-bound anion of thymine at higher excitation energies, in the range of 0.6 eV above the I(-)T VDE, which suggests that if the dipole-bound anion acts as a "doorway" to the valence-bound anion, it only does so at excitation energies near the VDE of the complex.

The role of the environment in the ion induced fragmentation of uracil

The fragmentation of uracil molecules and pure and nano-hydrated uracil clusters by ¹²C⁴⁺ ion impact is investigated.

Modeling photoionization of aqueous DNA and its components

Radiation damage to DNA is usually considered in terms of UVA and UVB radiation. These ultraviolet rays, which are part of the solar spectrum, can indeed cause chemical lesions in DNA, triggered by photoexcitation particularly in the UVB range. Damage can, however, be also caused by higher energy radiation, which can ionize directly the DNA or its immediate surroundings, leading to indirect damage. Thanks to absorption in the atmosphere, the intensity of such ionizing radiation is negligible in the solar spectrum at the surface of Earth. Nevertheless, such an ionizing scenario can become dangerously plausible for astronauts or flight personnel, as well as for persons present at nuclear power plant accidents. On the beneficial side, ionizing radiation is employed as means for destroying the DNA of cancer cells during radiation therapy. Quantitative information about ionization of DNA and its components is important not only for DNA radiation damage, but also for understanding redox properties of DNA in redox sensing or labeling, as well as charge migration along the double helix in nanoelectronics applications. Until recently, the vast majority of experimental and computational data on DNA ionization was pertinent to its components in the gas phase, which is far from its native aqueous environment. The situation has, however, changed for the better due to the advent of photoelectron spectroscopy in liquid microjets and its most recent application to photoionization of aqueous nucleosides, nucleotides, and larger DNA fragments. Here, we present a consistent and efficient computational methodology, which allows to accurately evaluate ionization energies and model photoelectron spectra of aqueous DNA and its individual components. After careful benchmarking, the method based on density functional theory and its time-dependent variant with properly chosen hybrid functionals and polarizable continuum solvent model provides ionization energies with accuracy of 0.2-0.3 eV, allowing for faithful modeling and interpretation of DNA photoionization. The key finding is that the aqueous medium is remarkably efficient in screening the interactions within DNA such that, unlike in the gas phase, ionization of a base, nucleoside, or nucleotide depends only very weakly on the particular DNA context. An exception is the electronic interaction between neighboring bases which can lead to sequence-specific effects, such as a partial delocalization of the cationic hole upon ionization enabled by presence of adjacent bases of the same type.

Vibrationally resolved photoelectron spectroscopy of electronic excited states of DNA bases: application to the ã state of thymine cation

For fully understanding the light-molecule interaction dynamics at short time scales, recent theoretical and experimental studies proved the importance of accurate characterizations not only of the ground (D0) but also of the electronic excited states (e.g., D1) of molecules. While ground state investigations are currently straightforward, those of electronic excited states are not. Here, we characterized the à electronic state of ionic thymine (T(+)) DNA base using explicitly correlated coupled cluster ab initio methods and state-of-the-art synchrotron-based electron/ion coincidence techniques. The experimental spectrum is composed of rich and long vibrational progressions corresponding to the population of the low frequency modes of T(+)(Ã). This work challenges previous numerous works carried out on DNA bases using common synchrotron and VUV-based photoelectron spectroscopies. We provide hence a powerful theoretical and experimental framework to study the electronic structure of ionized DNA bases that could be generalized to other medium-sized biologically relevant systems.

Automated quantum chemistry based molecular dynamics simulations of electron ionization induced fragmentations of the nucleobases Uracil, Thymine, Cytosine, and Guanine

The gas-phase decomposition pathways of electron ionization (EI)-induced radical cations of the nucleobases uracil, thymine, cytosine, and guanine are investigated by means of mixed quantum-classical molecular dynamics. No preconceived fragmentation channels are used in the calculations. The results compare well to a plethora of experimental and theoretical data for these important biomolecules. With our combined stochastic and dynamic approach, one can access in an unbiased way the energetically available decomposition mechanisms. Additionally, we are able to separate the EI mass spectra of different tautomers of cytosine and guanine. Our method (previously termed quantum chemistry electron ionization mass spectra) reproduces free nucleobase experimental mass spectra well and provides detailed mechanistic in-sight into high-energy unimolecular decomposition processes.

Proton irradiation of DNA nucleosides in the gas phase

Charge localization within nucleosides after proton irradiation is strongly influenced by the ionization energy of the base.

Ultrafast non-radiative decay of gas-phase nucleosides

De-excitation of DNA nucleosides on picosecond timescales was measured and found to be twice as fast as the equivalent nucleobases.