We present a first-principles molecular dynamics study of an excess electron in condensed phase models of solvated DNA bases. Calculations on increasingly large microsolvated clusters taken from liquid phase simulations show that adiabatic electron affinities increase systematically upon solvation, as for optimized gas-phase geometries. Dynamical simulations after vertical attachment indicate that the excess electron, which is initially found delocalized, localizes around the nucleobases within a 15 fs time scale. This transition requires small rearrangements in the geometry of the bases.

When biological matter is subjected to ionizing radiation, a wealth of secondary low-energy (<20 eV) electrons are produced. These electrons propagate inelastically, losing energy to the medium until they reach energies low enough to localize in regions of high electron affinity. We have recently shown that in fully solvated DNA fragments, nucleobases are particularly attractive for such excess electrons. The next question is what is their longer-term effect on DNA. It has been advocated that they can lead to strand breaks by cleavage of the phosphodiester C3′–O3′ bond. Here we present a first-principles study of free energy barriers for the cleavage of this bond in fully solvated nucleotides. We have found that except for dAMP, the barriers are on the order of 6 kcal/mol, suggesting that bond cleavage is a regular feature at 300 K. Such low barriers are possible only as a result of solvent and thermal fluctuations. These findings support the notion that low-energy electrons can indeed lead to strand breaks in DNA.

Reactions that can damage DNA have been simulated using a combination of molecular dynamics and density functional theory. In particular, the damage caused by the attachment of a low energy electron to the nucleobase. Simulations of anionic single nucleotides of DNA in an aqueous environment that was modelled explicitly have been performed. This has allowed us to examine the role played by the water molecules that surround the DNA in radiation damage mechanisms. Our simulations show that hydrogen bonding and protonation of the nucleotide by the water can have a significant effect on the barriers to strand breaking reactions. Furthermore, these effects are not the same for all four of the bases.

We present an efficient and accurate method to study electron detachment from negative ions by a few-cycle linearly polarized laser pulse. The adiabatic saddle-point method of Gribakin and Kuchiev [Phys. Rev. A 55, 3760 (1997)] is adapted to calculate the transition amplitude for a short laser pulse. Its application to a pulse with N optical cycles produces 2(N + 1) saddle points in complex time, which form a characteristic "smile." Numerical calculations are performed for H − in a 5-cycle pulse with frequency 0.0043 a.u. and intensities of 10 10 , 5 × 10 10 , and 10 11 W/cm 2 , and for various carrier-envelope phases. We determine the spectrum of the photoelectrons as a function of both energy and emission angle, as well as the angle-integrated energy spectra and total detachment probabilities. Our calculations show that the dominant contribution to the transition amplitude is given by 5-6 central saddle points, which correspond to the strongest part of the pulse. We examine the dependence of the photoelectron angular distributions on the carrier-envelope phase and show that measuring such distributions can provide a way of determining this phase.

Low-energy electron-impact hydrogen loss due to dissociative electron attachment (DEA) to the uracil and thymine molecules in a water cluster environment is investigated theoretically. Only the A(')-resonance contribution, describing the near-threshold behavior of DEA, is incorporated. Calculations are based on the nonlocal complex potential theory and the multiple scattering theory, and are performed for a model target with basic properties of uracil and thymine, surrounded by five water molecules. The DEA cross section is strongly enhanced when the attaching molecule is embedded in a water cluster. This growth is due to two effects: the increase of the resonance lifetime and the negative shift in the resonance position due to interaction of the intermediate negative ion with the surrounding water molecules. A similar effect was earlier found in DEA to chlorofluorocarbons.