The mystery of sub-picosecond charge transfer following irradiation of hydrated uridine monophosphate

Addressing electronic and dynamical evolution of molecules and molecular clusters: DFTB simulations of energy relaxation in polycyclic aromatic hydrocarbons

We present a review of the capabilities of the density functional based Tight Binding (DFTB) scheme to address the electronic relaxation and dynamical evolution of molecules and molecular clusters following energy deposition via either collision or photoabsorption. The basics and extensions of DFTB for addressing these systems and in particular their electronic states and their dynamical evolution are reviewed. Applications to PAH molecules and clusters, carbonaceous systems of major interest in astrochemical/astrophysical context, are reported. A variety of processes are examined and discussed such as collisional hydrogenation, fast collisional processes and induced electronic and charge dynamics, collision-induced fragmentation, photo-induced fragmentation, relaxation in high electronic states, electronic-to-vibrational energy conversion and statistical versus non-statistical fragmentation. This review illustrates how simulations may help to unravel different relaxation mechanisms depending on various factors such as the system size, specific electronic structure or excitation conditions, in close connection with experiments.

Irradiation of Plutonium Tributyl Phosphate Complexes by Ionizing Alpha Particles: A Computational Study

The PUREX solvent extraction process, widely used for recovering uranium and plutonium from spent nuclear fuel, utilizes an organic solvent composed of tributyl phosphate (TBP). The emission of ionizing particles such as alpha particles, resulting from the decay of plutonium, makes the organic solvent vulnerable to degradation. Here, we study the ultrashort time alpha irradiation of tributylphosphate (TBP) and Pu(NO3)4(TBP)2 complex formed in the PUREX process. Electron dynamics is propagated by Real-Time-Dependent Auxiliary Density Functional Theory (RT-TD-ADFT). We investigate the use of previously proposed absorption boundary conditions (ABC) in the molecular orbital space to treat secondary electron emission. Basis set and exchange correlation functional effects with ABC are reported as well as a detailed analysis of the ABC parametrization. Preliminary results on the water molecule and then on TBP show that the phenomenological nature of the ABC parameters necessitates selecting appropriate values for each system under study. Irradiation of free and complexed TBP shows an influence of the ligands on the variation of atomic charges on the femtosecond time scale. An accumulation of atomic charges in the alkyl chains of TBP is observed in the case where the nitrate groups are predominantly irradiated. In addition, we find that the Pu atom regains its electric charge very rapidly after being hit by the projectile, with the coordination sphere serving as an electron reservoir to preserve its formal redox state. This study paves the road toward a full understanding of the degradation of organic extracants employed in the nuclear industry.

Secondary Electron Attachment-Induced Radiation Damage to Genetic Materials

Reactions of radiation-produced secondary electrons (SEs) with biomacromolecules (e.g., DNA) are considered one of the primary causes of radiation-induced cell death. In this Review, we summarize the latest developments in the modeling of SE attachment-induced radiation damage. The initial attachment of electrons to genetic materials has traditionally been attributed to the temporary bound or resonance states. Recent studies have, however, indicated an alternative possibility with two steps. First, the dipole-bound states act as a doorway for electron capture. Subsequently, the electron gets transferred to the valence-bound state, in which the electron is localized on the nucleobase. The transfer from the dipole-bound to valence-bound state happens through a mixing of electronic and nuclear degrees of freedom. In the presence of aqueous media, the water-bound states act as the doorway state, which is similar to that of the presolvated electron. Electron transfer from the initial doorway state to the nucleobase-bound state in the presence of bulk aqueous media happens on an ultrafast time scale, and it can account for the decrease in DNA strand breaks in aqueous environments. Analyses of the theoretically obtained results along with experimental data have also been discussed.

Multiscale molecular modelling: from electronic structure to dynamics of nanosystems and beyond

Important contemporary biological and materials problems often depend on interactions that span orders of magnitude differences in spatial and temporal dimensions. This Tutorial Review attempts to provide an introduction to such fascinating problems through a series of case studies, aimed at beginning researchers, graduate students, postdocs and more senior colleagues who are changing direction to focus on multiscale aspects of their research. The choice of specific examples is highly personal, with examples either chosen from our own work or outstanding multiscale efforts from the literature. I start with various embedding schemes, as exemplified by polarizable continuum models, 3-D RISM, molecular DFT and frozen-density embedding. Next, QM/MM (quantum mechanical/molecular mechanical) techniques are the workhorse of pm-to-nm/ps-to-ns simulations; examples are drawn from enzymes and from nanocatalysis for oil-sands upgrading. Using polarizable force-fields in the QM/MM framework represents a burgeoning subfield; with examples from ion channels and electron dynamics in molecules subject to strong external fields, probing the atto-second dynamics of the electrons with RT-TDDFT (real-time - time-dependent density functional theory) eventually coupled with nuclear motion through the Ehrenfest approximation. This is followed by a section on coarse graining, bridging dimensions from atoms to cells. The penultimate chapter gives a quick overview of multiscale approaches that extend into the meso- and macro-scales, building on atomistic and coarse-grained techniques to enter the world of materials engineering, on the one hand, and cell biology, on the other. A final chapter gives just a glimpse of the burgeoning impact of machine learning on the structure-dynamics front. I aim to capture the excitement of contemporary leading-edge breakthroughs in the description of physico-chemical systems and processes in complex environments, with only enough historical content to provide context and aid the next generation of methodological development. While I aim also for a clear description of the essence of methodological breakthroughs, equations are kept to a minimum and detailed formalism and implementation details are left to the references. My approach is very selective (case studies) rather than exhaustive. I think that these case studies should provide fodder to build as complete a reference tree on multiscale modelling as the reader may wish, through forward and backward citation analysis. I hope that my choices of cases will excite interest in newcomers and help to fuel the growth of multiscale modelling in general.

Modulation of the Directionality of Hole Transfer between the Base and the Sugar-Phosphate Backbone in DNA with the Number of Sulfur Atoms in the Phosphate Group

This work shows that S atom substitution in phosphate controls the directionality of hole transfer processes between the base and sugar-phosphate backbone in DNA systems. The investigation combines synthesis, electron spin resonance (ESR) studies in supercooled homogeneous solution, pulse radiolysis in aqueous solution at ambient temperature, and density functional theory (DFT) calculations of in-house synthesized model compound dimethylphosphorothioate (DMTP(O^-)═S) and nucleotide (5'-O-methoxyphosphorothioyl-2'-deoxyguanosine (G-P(O^-)═S)). ESR investigations show that DMTP(O^-)═S reacts with Cl₂^•- to form the σ²σ*¹ adduct radical -P-S[Formula: see text]Cl, which subsequently reacts with DMTP(O^-)═S to produce [-P-S[Formula: see text]S-P-]^-. -P-S[Formula: see text]Cl in G-P(O^-)═S undergoes hole transfer to Gua, forming the cation radical (G^•+) via thermally activated hopping. However, pulse radiolysis measurements show that DMTP(O^-)═S forms the thiyl radical (-P-S^•) by one-electron oxidation, which did not produce [-P-S[Formula: see text]S-P-]^-. Gua in G-P(O^-)═S is oxidized unimolecularly by the -P-S^• intermediate in the sub-picosecond range. DFT thermochemical calculations explain the differences in ESR and pulse radiolysis results obtained at different temperatures.