Reliability and performances of real-time time-dependent auxiliary density functional theory

Use of Gaussian-Type Functions for Describing Fast Ion-Matter Irradiation with Time-Dependent Density Functional Theory

The electronic stopping power is an observable property that quantifies the ability of swift ions to penetrate matter to transfer energy to the electron cloud. The recent literature has proven the value of Real-Time Time-Dependent Density Functional Theory to accurately evaluate this property from first-principles, but questions remain regarding the capability of computer codes relying on atom-centered basis functions to capture the physics at play. In this Perspective, we draw attention to the fact that irradiation by swift ions triggers electron emission into the continuum, especially at the Bragg peak. We investigate the ability of Gaussian atomic orbitals (AOC), which were fitted to mimic continuum wave functions, to improve electronic stopping power predictions. AOC are added to standard correlation-consistent basis sets or STO minimal basis sets. Our benchmarks for water irradiation by fast protons clearly advocate for the use of AOC, especially near the Bragg peak. We show that AOC only need to be placed on the molecules struck by the ion. The number of AOC that are added to the usual basis set is relatively small compared to the total number of atomic orbitals, making the use of such a basis set an excellent choice from a computational cost point of view. The optimum basis set combination is applied for the calculation of the stopping power of a proton in water with encouraging agreement with experimental data.

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.

Following the density evolution using real time density functional theory and density based indexes: Application to model push-pull molecules

Considering as test case a family of organic rod like push-pull molecules, we derived and applied density based index enabling the description and diagnostic of the electronic density evolution in real time-time dependent density functional theory (RT-TDDFT) simulations. In particular, both the charge transfer (CT) distance and a diagnostic index, the D_CT and M_AC ^RT respectively, were computed on the fly from the density distribution obtained at a given time and the reference ground state density and their mean values were compared with what obtained at Linear Response-TDDFT level. Besides giving a way of analyzing the density redistribution occurring in time, these tools allowed to show how RT-TDDFT, which is definitely a powerful method to model the evolution of the density in CT or charge separation processes, can be affected by the same artifacts known for LR-TDDFT approaches and, particularly, to those related to the use of approximate exchange correlation functionals. The analysis here performed allowed to identify and discard on fly the electronic configurations corresponding to spurious situations.

Ab Initio Prediction of a Negative Barkas Coefficient for Slow Protons and Antiprotons in LiF

We report the ab initio prediction of a negative Barkas coefficient in lithium fluoride (LiF) insulator at low velocity (v<0.25  a.u., E_{kin}∼2  keV). The electronic stopping power of protons in LiF has been extensively studied both experimentally and theoretically because of a controversial threshold effect. While our time-dependent density-functional theory simulations confirm the presence of a velocity threshold below which the proton stopping power vanishes, our calculations demonstrate that the antiprotons do not experience such a threshold. The combination of those two contrasting behaviors gives rise to an unprecedented negative Barkas effect: the stopping power of antiprotons is larger than that of protons. We identify that the slow antiproton at close encounter destabilizes a p orbital of the F^{-} anion pointing toward the antiproton. This particular orbital becomes highly polarizable and hence contributes much to the stopping power.

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

The early mechanisms by which ionizing rays damage biological structures by so-called direct effects are largely elusive. In a recent picosecond pulse radiolysis study of concentrated uridine monophosphate solutions [J. Ma, S. A. Denisov, J.-L. Marignier, P. Pernot, A. Adhikary, S. Seki and M. Mostafavi, J. Phys. Chem. Lett., 2018, 9, 5105], unexpected results were found regarding the oxidation of the nucleobase. The signature of the oxidized nucleobase could not be detected 5 ps after the electron pulse, but only the oxidized phosphate, raising intriguing questions about the identity of charge-transfer mechanisms that could explain the absence of U⁺. We address here this question by means of advanced first-principles atomistic simulations of solvated uridine monophosphate, combining Density Functional Theory (DFT) with polarizable embedding schemes. We contrast three very distinct mechanisms of charge transfer covering the atto-, femto- and pico-second timescales. We first investigate the ionization mechanism and subsequent hole/charge migrations on a timescale of attoseconds to a few femtoseconds under the frozen nuclei approximation. We then consider a nuclear-driven phosphate-to-oxidized-nucleobase electron transfer, showing that it is an uncompetitive reaction channel on the sub-picosecond timescale, despite its high exothermicity and significant electronic coupling. Finally, we show that non-adiabatic charge transfer is enabled by femtosecond nuclear relaxation after ionization. We show that electronic decoherence and the electronic coupling strength are the key parameters that determine the hopping probabilities. Our results provide important insight into the interplay between electronics and nuclear motions in the early stages of the multiscale responses of biological matter subjected to ionizing radiation.