A Simple ab Initio Model for the Hydrated Electron That Matches Experiment

Ionization Patterns and Chemical Reactivity of Cytosine-Guanine Watson-Crick Pairs

Gas-phase and aqueous vertical ionization potentials, vIPgas and vIPaq respectively and measurements of the molecular electrostatic and local ionization maps calculated at the DFT/B3LYP-D3/ 6-311+G** level of theory and the C-PCM reaction field model for single- and double-stranded CpG and 5MeCpG pairs show that the vIPaq for single- and double-stranded pairs of C-G and 5MeC-G are practically the same, in the range of 5.79 to 5.81 eV. The aqueous adiabatic ionization potentials for single-stranded CpG and 5MeCpG are 5.52 eV and 5.51 eV respectively and they reflect the nuclear reorganization that takes place after the abstraction of the electron. The aqueous adiabatic ionization energy values that correspond to the CpG+. radical cation and the hydrated electron, e-,, being at infinite distance, adIPaq+Vo, are 3.92 eV and 3.91 eV respectively with (Vo=-1.6 eV) Analysis of data suggest that the HOMO-LUMO energy gap in the hard/soft-acid/base (HSAB) concept cannot be used a priori to determine the effect of cytosine methylation on the guanine enhanced oxidative damage in DNA.

Photo-excited charge transfer from adamantane to electronic bound states in water

Aqueous nanodiamonds illuminated by UV light produce free solvated electrons, which may drive high-energy reduction reactions in water. However, the influence of water conformations on the excited-state electron-transfer mechanism are still under debate. In this work, we offer a theoretical study of charge-transfer states in adamantane-water structures obtained by linear-response time-dependent density-functional theory. Small water clusters with broken hydrogen bonds are found to efficiently bind the electron from adamantane. A distinction is made with respect to the nature of the water clusters: some bind the electron in a water cavity, others along a strong permanent total dipole. These two types of bound states are more strongly binding, the higher their electron affinity and their positive electrostatic potential, the latter being dominated by the energy of the lowest unoccupied molecular orbital of the isolated water clusters. Structural sampling in a thermal equilibrium at room temperature via molecular dynamics snapshots confirms under which conditions the underlying waters clusters can occur and verifies that broken hydrogen bonds in the water network close to adamantane can create traps for the solvated electron.

Immobilization of Bi2WO6 on Polymer Membranes for Photocatalytic Removal of Micropollutants from Water - A Stable and Visible Light Active Alternative

In this work, bismuth tungstate Bi2WO6 is immobilized on polymer membranes to photocatalytically remove micropollutants from water as an alternative to titanium dioxide TiO2. A synthesis method for Bi2WO6 preparation and its immobilization on a polymer membrane is developed. Bi2WO6 is characterized using X-ray diffraction and UV-vis reflectance spectroscopy, while the membrane undergoes analysis through scanning electron microscopy, X-ray photoelectron spectroscopy, and degradation experiments. The density of states calculations for TiO2 and Bi2WO6, along with PVDF reactions with potential reactive species, are investigated by density functional theory. The generation of hydroxyl radicals OH• is investigated via the reaction of coumarin to umbelliferone via fluorescence probe detection and electron paramagnetic resonance. Increasing reactant concentration enhances Bi2WO6 crystallinity. Under UV light at pH 7 and 11, the Bi2WO6 membrane completely degrades propranolol in 3 and 1 h, respectively, remaining stable and reusable for over 10 cycles (30 h). Active under visible light with a bandgap of 2.91 eV, the Bi2WO6 membrane demonstrates superior stability compared to a TiO2 membrane during a 7-day exposure to UV light as Bi2WO6 does not generate OH• radicals. The Bi2WO6 membrane is an alternative for water pollutant degradation due to its visible light activity and long-term stability.

A fresh perspective on metal ammonia molecular complexes and expanded metals: opportunities in catalysis and quantum information

Recent advances in our comprehension of the electronic structure of metal ammonia complexes have opened avenues for novel materials with diffuse electrons. These complexes in their ground state can host peripheral "Rydberg" electrons which populate a hydrogenic-type shell model imitating atoms. Aggregates of such complexes form the so-called expanded or liquid metals. Expanded metals composed of d- and f-block metal ammonia complexes offer properties, such as magnetic moments and larger numbers of diffuse electrons, not present for alkali and alkaline earth (s-block) metals. In addition, tethering metal ammonia complexes via hydrocarbon chains (replacement of ammonia ligands with diamines) yields materials that can be used for redox catalysis and quantum computing, sensing, and optics. This perspective summarizes the recent findings for gas-phase isolated metal ammonia complexes and projects the obtained knowledge to the condensed phase regime. Possible applications for the newly introduced expanded metals and linked solvated electrons precursors are discussed and future directions are proposed.

Theoretical investigation of Aryl/Alkyl halide reduction with hydrated electrons from energy and AIMD aspects

CONTEXT

In this study, the reactions of hydrated electron (e^-(aq)) with alkyl and aryl halides were simulated with an ab initial molecular dynamics (AIMD) method to reveal the underlying mechanism. An original protocol was developed for preparing the proper initial wavefunction guess of AIMD, in which a single electron was curled in a tetrahedral cavity of four water molecules. Our results show that the stability of e^-(aq) increases with the hydrogen bond grid integrity. The organic halides prefer to react with e^-(aq) in neutral or alkaline environment, while they are more likely to react with hydrogen radical (the product of e^-(aq) and proton) under acidic conditions. The reaction between fluorobenzene/fluoromethane and hydrogen radical is considered as the least favorable reaction due to the highest reaction barriers. The bond dissociation energy (BDE) suggested that the cleavage of the carbon-halogen bond of their anion radical might be a thermodynamically favorable reaction. AIMD results indicated that the LUMO or higher orbitals were the e^-(aq) migration destination. The transplanted electron enhanced carbon-halogen bond vibration intensively, leading to bond cleavage. The solvation process of the departing halogen anions was observed in both fluorobenzene and fluoromethane AIMD simulation, indicating that it might have a significant effect on enthalpy. Side reactions and byproducts obtained during the AIMD simulation suggested the complexity of the e^-(aq) reactions and further investigation was needed to fully understand the reaction mechanisms. This study provided theoretical insight into the pollutant environmental fate and constructed a methodological foundation for AIMD simulation of analogous free radical reactions.

METHODS

The theoretical calculation was conducted on the combination of Gaussian16 and ORCA5.0.3 software packages. The initial geometries, as well as the wavefunction initial guesses, were obtained at PBE0/ma-def2-TZVP/IEFPCM-water level in Gaussian16 unless otherwise stated. AIMD simulations were performed at the same level in ORCA. Wavefunction analysis was carried out with Multiwfn. The details methods were described in the section "Computational details" section.

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.

Hydrated electrons as nodes in porous clathrate hydrates

We investigate the structures of hydrated electrons (e^- _aq) in one of water's solid phases, namely, clathrate hydrates (CHs). Using density functional theory (DFT) calculations, DFT-based ab initio molecular dynamics (AIMD), and path-integral AIMD simulations with periodic boundary conditions, we find that the structure of the e^- _aq@node model is in good agreement with the experiment, suggesting that an e^- _aq could form a node in CHs. The node is a H₂O defect in CHs that is supposed to be composed of four unsaturated hydrogen bonds. Since CHs are porous crystals that possess cavities that can accommodate small guest molecules, we expect that these guest molecules can be used to tailor the electronic structure of the e^- _aq@node, and it leads to experimentally observed optical absorption spectra of CHs. Our findings have a general interest and extend the knowledge of e^- _aq into porous aqueous systems.

Photoionization of 3-Methylindole Embedded in Sodium Dodecyl Sulfate and Dodecyltrimethylammonium Chloride Micelles: Migration of Electrons Generated in Micelle Cores and Their Solvation in Outside Water

Electrons were generated in the core of micelles formed by negatively charged sodium dodecyl sulfate (SDS) or positively charged dodecyltrimethylammonium chloride (DTAC) by photoionization of 3-methylindole embedded in the core. The electrons were hydrated after they moved out of the core to the outer aqueous phase. These processes were monitored with femtosecond time-resolved absorption spectroscopy. The migration of electrons from the micelle core to the outer aqueous phase was faster than the instrumental response time of 200 fs. Hot electrons in the aqueous phase were produced in ≤320 fs. There was no significant difference observed for the micellar solutions of negatively charged SDS and positively charged DTAC, or for water. The geminate recombination between the electrons and the radical cations was hindered to a large extent once the electrons hydrated at the outer aqueous phase were separated from the radical cations remaining in the micelle core.

A General Framework of Scanning Tunneling Microscopy Based on Bardeen's Approximation for Isolated Molecules

Scanning tunneling microscopy (STM) is one of the most popular techniques for precise characterization. Yet, its current theoretical implementation is often based on the periodic boundary condition with the Tersoff-Hamann approximation, which is inefficient to explore the tip states other than the s-wave and to treat properly the charged molecules that are ubiquitous in chemistry. In this work, we establish a general theoretical framework for STM image simulations, which is based on the Bardeen's approximation and utilizes the boundary condition of the cluster model. We develop an analytic algorithm for the precise evaluation of the transfer Hamiltonian matrix, addressing correctly the asymptotic behaviors of the tip states. Numerical results demonstrate that the molecular images under different STM tip states and mapping modes can be quantitatively simulated in the present framework, which paves the avenue for the conclusive investigation of the ground state electronic structures for either neutral or charged molecules.

Ab Initio Studies of Hydrated Electron/Cation Contact Pairs: Hydrated Electrons Simulated with Density Functional Theory Are Too Kosmotropic

We have performed the first DFT-based ab initio MD simulations of a hydrated electron (e_aq^-) in the presence of Na⁺, a system chosen because ion-pairing behavior in water depends sensitively on the local hydration structure. Experiments show that e_aq^-'s interact weakly with Na⁺; the e_aq^-'s spectrum blue shifts by only a few tens of meV upon ion pairing without changing shape. We find that the spectrum of the DFT-simulated e_aq^- red shifts and changes shape upon interaction with Na⁺, in contrast with experiment. We show that this is because the hydration structure of the DFT-simulated e_aq^- is too ordered or kosmotropic. Conversely, simulations that produce e_aq^-'s with a less ordered or chaotropic hydration structure form weaker ion pairs with Na⁺, yielding predicted spectral blue shifts in better agreement with experiment. Thus, ab initio simulations based on hybrid GGA DFT functionals fail to produce the correct solvation structure for the hydrated electron.

Understanding the Temperature Dependence and Finite Size Effects in Ab Initio MD Simulations of the Hydrated Electron

The hydrated electron is of interest to both theorists and experimentalists as a paradigm solution-phase quantum system. Although the bulk of the theoretical work studying the hydrated electron is based on mixed quantum/classical (MQC) methods, recent advances in computer power have allowed several attempts to study this object using ab initio methods. The difficulty with employing ab initio methods for this system is that even with relatively inexpensive quantum chemistry methods such as density functional theory (DFT), such calculations are still limited to at most a few tens of water molecules and only a few picoseconds duration, leaving open the question as to whether the calculations are converged with respect to either system size or dynamical fluctuations. Moreover, the ab initio simulations of the hydrated electron that have been published to date have provided only limited analysis. Most works calculate the electron's vertical detachment energy, which can be compared to experiment, and occasionally the electronic absorption spectrum is also computed. Structural features, such as pair distribution functions, are rare in the literature, with the majority of the structural analysis being simple statements that the electron resides in a cavity, which are often based only on a small number of simulation snapshots. Importantly, there has been no ab initio work examining the temperature-dependent behavior of the hydrated electron, which has not been satisfactorily explained by MQC simulations. In this work, we attempt to remedy this situation by running DFT-based ab initio simulations of the hydrated electron as a function of both box size and temperature. We show that the calculated properties of the hydrated electron are not converged even with simulation sizes up to 128 water molecules and durations of several tens of picoseconds. The simulations show significant changes in the water coordination and solvation structure with box size. Our temperature-dependent simulations predict a red-shift of the absorption spectrum (computed using TD-DFT with an optimally tuned range-separated hybrid functional) with increasing temperature, but the magnitude of the predicted red-shift is larger than that observed experimentally, and the absolute position of the calculated spectra are off by over half an eV. The spectral red-shift at high temperatures is accompanied by both a partial loss of structure of the electron's central cavity and an increased radius of gyration that pushes electron density onto and beyond the first solvation shell. Overall, although ab initio simulations can provide some insights into the temperature-dependent behavior of the hydrated electron, the simulation sizes and level of quantum chemistry theory that are currently accessible are inadequate for correctly describing the experimental properties of this fascinating object.

Degradation of Per- and Polyfluoroalkyl Substances with Hydrated Electrons: A New Mechanism from First-Principles Calculations

Per- and polyfluoroalkyl substances (PFASs) are synthetic contaminants found in drinking groundwater sources and a wide variety of consumer products. Because of their adverse environmental and human health effects, remediation of these persistent compounds has attracted significant recent attention. To gain mechanistic insight into their remediation, we present the first ab initio study of PFAS degradation via hydrated electrons─a configuration that has not been correctly considered in previous computational studies up to this point. To capture these complex dynamical effects, we harness ab initio molecular dynamics (AIMD) simulations to probe the reactivities of perfluorooctanoic (PFOA) and perfluorooctane sulfonic acid (PFOS) with hydrated electrons in explicit water. We complement our AIMD calculations with advanced metadynamics sampling techniques to compute free energy profiles and detailed statistical analyses of PFOA/PFOS dynamics. Although our calculations show that the activation barrier for C-F bond dissociation in PFOS is three times larger than that in PFOA, all the computed free energy barriers are still relatively low, resulting in a diffusion-limited process. We discuss our results in the context of recent studies on PFAS degradation with hydrated electrons to give insight into the most efficient remediation strategies for these contaminants. Most importantly, we show that the degradation of PFASs with hydrated electrons is markedly different from that with excess electrons/charges, a common (but largely incomplete) approach used in several earlier computational studies.

Nitro rotation tuned dissociative electron attachment upon targeted radiosensitizer 4-substituted Z bases

In this work, a set of new potential radiation sensitizers (4-substituted Z-bases: 4XZ, X = F, Cl, Br, and I) are designed based on the artificial 6-amino-5-nitro-3-(1'-β-D-2'-deoxyribofuranosyl)-2(1H)-pyridone (Z), which can selectively bind to breast cancer cells. The calculated electron affinities in water solution show that the halogenated Z-bases are efficient electron acceptors which possess significant electron-withdrawing characters following the order of 4XZ > Z ≫ U. To ensure the effective electron attachment induced dissociation, we constructed the energy profiles related to the X-C bond cleavage of neutral and anionic bases. The results show that the X-C bond becomes relatively weak after the electron attachment. In particular, the electron induced dehalogenations of (4BrZ)^- and (4IZ)^- are low-barrier and exothermic, which support a high radiosensitivity. Furthermore, we characterized the vibrational excitation effect on the dissociative electron attachment, which demonstrates that the charge distribution can be regulated by the rotation-induced structural distortion accompanied by the electron localization on the nitro group. Also examined is the influence of base pairing on the dehalogenation, which is not only conducive to the electron-driven dissociation but is also beneficial to the stabilization of related products. The current study suggests 4BrZ and 4IZ can be regarded as potential targeted radiosensitizers with possible applications in reducing the side effects in radiotherapy.

Proton-Transfer Reactions in One-Electron-Oxidized G-Quadruplexes: A Density Functional Theory Study

Recently, G-quadruplexes (Gq) formed in B-DNA as secondary structures are found to be important therapeutic targets and material for developing nanodevices. Gq are guanine-rich and thus susceptible to oxidative damage by producing short-lived intermediate radicals via proton-transfer reactions. Understanding the mechanisms of radical formation in Gq is of fundamental interest to understand the early stages of DNA damage. Herein, we used density functional theory including aqueous phase (ωB97XD-PCM/6-31++G**) and considered single layer of Gq [G-quartets (G4): association of four guanines in a cyclic Hoogsteen hydrogen-bonded arrangement (Scheme 1)] to unravel the mechanisms of formation of intermediates by calculating the relative Gibbs free energies and spin density distributions of one-electron-oxidized G4 and its various proton-transfer states: G^•+, G(N₁-H)^•, G(N₂-H')^•, G(N₂-H″)^•, G(N₁-H)^•-(H⁺O₆)G, and G(N₂-H)^•-(H⁺N₇)G. The present calculation predicts the formation of G(N₂-H)^•-(H⁺N₇)G, which is only ca. 0.8 kcal/mol higher in energy than the initially formed G^•+. The formation of G(N₂-H)^•-(H⁺N₇)G plays a key role in explaining the formation of 8-OG along with G(N₁-H)^• formation via tautomerization from G(N₂-H)^•, as proposed recently.

Insights into the hydrated electron generation from UV/aniline: Mechanism and quantum efficiency

The insight into the neglected reduction process accompanied by UV-based oxidation process may provide new ideas for the development of advanced oxidation and reduction technologies. In this study, aniline was comprehensively investigated as an unexpected indicator of hydrated electron (e_aq^-) under UV irradiation. Monochloroacetic acid (MCAA) was selected as the probe of e_aq^- and the balance of chloride ions indicated the reduction of MCAA. Further, laser flash photolysis experiments demonstrated the generation of e_aq^- in the UV/aniline process and the half-life period of formed e_aq^- was demonstrated to be 0.13 μs. The photolysis of aniline along with the decay of the excited state of aniline was responsible for the e_aq^- generation. Besides, the hydrogen atom (H) generated from the photolysis can subsequently reacted with OH^- to generate e_aq^-. The photolysis pathways of aniline were proposed by the results of GC-MS. Aniline was abstracted of H in solution to the formation of aniline radical (PhNH) or form aminophenol in three different isomers (orto-, meta- and para-aminophenol). Moreover, UV/aniline showed a higher reducing capacity of MCAA compared with other organic electron donors and sustained a highly reducing ability in a wide pH. And the calculation results of quantum efficiency (Φ) showed that excessive aniline was not conducive to the elevation of Φ. This study introduced a novel pathway of e_aq^- generation during the photolysis of aniline and provided a new perspective for e_aq^--based advanced reduction processes.

Flexible boundary layer using exchange for embedding theories. II. QM/MM dynamics of the hydrated electron

The FlexiBLE embedding method introduced in Paper I [Z. Shen and W. J. Glover, J. Chem. Phys. 155, 224112 (2021)] is applied to explore the structure and dynamics of the aqueous solvated electron at an all-electron density functional theory Quantum Mechanics/Molecular Mechanics level. Compared to a one-electron mixed quantum/classical description, we find the dynamics of the many-electron model of the hydrated electron exhibits enhanced coupling to water OH stretch modes. Natural bond orbital analysis reveals this coupling is due to significant population of water OH σ^* orbitals, reaching 20%. Based on this, we develop a minimal frontier orbital picture of the hydrated electron involving a cavity orbital and important coupling to 4-5 coordinating OH σ^* orbitals. Implications for the interpretation of the spectroscopy of this interesting species are discussed.

Electric fields and potentials in condensed phases

The electric fields and potentials inside and at the interface of matter are relevant to many branches of physics, chemistry, and biology. Accurate quantification of these fields and/or potentials is essential to control and exploit chemical and physical transformations. Before we understand the response of matter to external fields, it is first important to understand the intrinsic interior and interfacial fields and potentials, both classically and quantum mechanically, as well as how they are probed experimentally. Here we compare and contrast, beginning with the hydrogen atom in vacuum and ending with concentrated aqueous NaCl electrolyte, both classical and quantum mechanical electric potentials and fields. We make contact with experimental vibrational Stark, electrochemical, X-ray, and electron spectroscopic probes of these potentials and fields, outline relevant conceptual difficulties, and underscore the advantage of electron holography as a basis to better understand electrostatics in matter.

Electronic Structure of Water from Koopmans-Compliant Functionals

Obtaining a precise theoretical description of the spectral properties of liquid water poses challenges for both molecular dynamics (MD) and electronic structure methods. The lower computational cost of the Koopmans-compliant functionals with respect to Green's function methods allows the simulations of many MD trajectories, with a description close to the state-of-art quasi-particle self-consistent GW plus vertex corrections method (QSGW + f_xc). Thus, we explore water spectral properties when different MD approaches are used, ranging from classical MD to first-principles MD, and including nuclear quantum effects. We have observed that different MD approaches lead to up to 1 eV change in the average band gap; thus, we focused on the band gap dependence with the geometrical properties of a system to explain such spread. We have evaluated the changes in the band gap due to variations in the intramolecular O-H bond distance and HOH angle, as well as the intermolecular hydrogen bond O···O distance and the OHO angles. We have observed that the dominant contribution comes from the O-H bond length; the O···O distance plays a secondary role, and the other geometrical properties do not significantly influence the gap. Furthermore, we analyze the electronic density of states (DOS), where the KIPZ functional shows good agreement with the DOS obtained with state-of-art approaches employing quasi-particle self-consistent GW plus vertex corrections. The O-H bond length also significantly influences the DOS. When nuclear quantum effects are considered, broadening of the peaks driven by the broader distribution of the O-H bond lengths is observed, leading to a closer agreement with the experimental photoemission spectra.

Radiation-Induced Graft Immobilization (RIGI): Covalent Binding of Non-Vinyl Compounds on Polymer Membranes

Radiation-induced graft immobilization (RIGI) is a novel method for the covalent binding of substances on polymeric materials without the use of additional chemicals. In contrast to the well-known radiation-induced graft polymerization (RIGP), RIGI can use non-vinyl compounds such as small and large functional molecules, hydrophilic polymers, or even enzymes. In a one-step electron-beam-based process, immobilization can be performed in a clean, fast, and continuous operation mode, as required for industrial applications. This study proposes a reaction mechanism using polyvinylidene fluoride (PVDF) and two small model molecules, glycine and taurine, in aqueous solution. Covalent coupling of single molecules is achieved by radical recombination and alkene addition reactions, with water radiolysis playing a crucial role in the formation of reactive solute species. Hydroxyl radicals contribute mainly to the immobilization, while solvated electrons and hydrogen radicals play a minor role. Release of fluoride is mainly induced by direct ionization of the polymer and supported by water. Hydrophobic chains attached to cations appear to enhance the covalent attachment of solutes to the polymer surface. Computational work is complemented by experimental studies, including X-ray photoelectron spectroscopy (XPS) and fluoride high-performance ion chromatography (HPIC).

Electron-Induced Repair of 2'-Deoxyribose Sugar Radicals in DNA: A Density Functional Theory (DFT) Study

In this work, we used ωB97XD density functional and 6-31++G** basis set to study the structure, electron affinity, populations via Boltzmann distribution, and one-electron reduction potentials (E°) of 2′-deoxyribose sugar radicals in aqueous phase by considering 2′-deoxyguanosine and 2′-deoxythymidine as a model of DNA. The calculation predicted the relative stability of sugar radicals in the order C4′^• > C1′^• > C5′^• > C3′^• > C2′^•. The Boltzmann distribution populations based on the relative stability of the sugar radicals were not those found for ionizing radiation or OH-radical attack and are good evidence the kinetic mechanisms of the processes drive the products formed. The adiabatic electron affinities of these sugar radicals were in the range 2.6–3.3 eV which is higher than the canonical DNA bases. The sugar radicals reduction potentials (E°) without protonation (−1.8 to −1.2 V) were also significantly higher than the bases. Thus the sugar radicals will be far more readily reduced by solvated electrons than the DNA bases. In the aqueous phase, these one-electron reduced sugar radicals (anions) are protonated from solvent and thus are efficiently repaired via the “electron-induced proton transfer mechanism”. The calculation shows that, in comparison to efficient repair of sugar radicals by the electron-induced proton transfer mechanism, the repair of the cyclopurine lesion, 5′,8-cyclo-2′-dG, would involve a substantial barrier.

A review of simulation codes and approaches for radiation chemistry

Historically, the field of radiation chemistry began shortly after the discovery of radioactivity, and its development has been closely related to discoveries in other related fields such as radiation and nuclear physics. Radiolysis of water and radiation chemistry have been very important in elucidating how radiation affects living matter and how it induces DNA damage. Nowadays, we recognize the importance of chemistry to understanding the effects of radiation on cells; however, it took several decades to obtain this insight, and much is still unknown. The radiolysis of water and aqueous solutions have been the subject of much experimental and theoretical research for many decades. One important concept closely related to radiation chemistry is radiation track structure. Track structure results from early physical and physicochemical events that lead to a highly non-homogenous distribution of radiolytic species. Because ionizing radiation creates unstable species that are distributed non-homogenously, the use of conventional reaction kinetics methods does not describe this chemistry well. In recent years, several methods have been developed for simulating radiation chemistry. In this review, we give a brief history of the field and the development of the simulation codes. We review the current methods used to simulate radiolysis of water and radiation chemistry, and we describe several radiation chemistry codes and their applications.

Mechanism of Aqueous Carbon Dioxide Reduction by the Solvated Electron

Aqueous solvated electron (e_aq^-), a key species in radiation and plasma chemistry, can efficiently reduce CO₂ in a potential green chemistry application. Here, the mechanism of this reaction is unravelled by condensed-phase molecular dynamics based on the correlated wave function and an accurate density functional theory (DFT) approximation. Here, we design and apply the holistic protocol for solvated electron's reactions encompassing all relevant reaction stages starting from diffusion. The carbon dioxide reduction proceeds via a cavity intermediate, which is separated from the product (CO₂^-) by an energy barrier due to the bending of CO₂ and the corresponding solvent reorganization energy. The formation of the intermediate is caused by solvated electron's diffusion, whereas the intermediate transformation to CO₂^- is triggered by hydrogen bond breaking in the second solvation shell of the solvated electron. This picture of an activation-controlled e_aq^- reaction is very different from both rapid barrierless electron transfer and proton-coupled electron transfer, where key transformations are caused by proton migration.

5-Nitro-2,4-Dichloropyrimidine as an Universal Model for Low-Energy Electron Processes Relevant for Radiosensitization

We report experimental results of low-energy electron interactions with.

Evaluating Simple Ab Initio Models of the Hydrated Electron: The Role of Dynamical Fluctuations

Despite its importance in electron transfer reactions and radiation chemistry, there has been disagreement over the fundamental nature of the hydrated electron, such as whether or not it resides in a cavity. Mixed quantum/classical simulations of the hydrated electron give different structures depending on the pseudopotential employed, and ab initio models of computational necessity use small numbers of water molecules and/or provide insufficient statistics to compare to experimental observables. A few years ago, Kumar et al. (J. Phys. Chem. A 2015, 119, 9148) proposed a minimalist ab initio model of the hydrated electron with only a small number of explicitly treated water molecules plus a polarizable continuum model (PCM). They found that the optimized geometry had four waters arranged tetrahedrally around a central cavity, and that the calculated vertical detachment energy and radius of gyration agreed well with experiment, results that were largely independent of the level of theory employed. The model, however, is based on a fixed structure at 0 K and does not explicitly incorporate entropic contributions or the thermal fluctuations that should be associated with the room-temperature hydrated electron. Thus, in this paper, we extend the model of Kumar et al. by running Born-Oppenheimer molecular dynamics (BOMD) of a small number of water molecules with an excess electron plus PCM at room temperature. We find that when thermal fluctuations are introduced, the level of theory chosen becomes critical enough when only four waters are used that one of the waters dissociates from the cluster with certain density functionals. Moreover, even with an optimally tuned range-separated hybrid functional, at room temperature the tetrahedral orientation of the 0 K first-shell waters is entirely lost and the central cavity collapses, a process driven by the fact that the explicit water molecules prefer to make H-bonds with each other more than with the excess electron. The resulting average structure is quite similar to that produced by a noncavity mixed quantum/classical model, so that the minimalist 4-water BOMD models suffer from problems similar to those of noncavity models, such as predicting the wrong sign of the hydrated electron's molar solvation volume. We also performed BOMD with 16 explicit water molecules plus an extra electron and PCM. We find that the inclusion of an entire second solvation shell of explicit water leads to little change in the outcome from when only four waters were used. In fact, the 16-water simulations behave much like those of water cluster anions, in which the electron localizes at the cluster surface, showing that PCM is not acceptable for use in minimalist models to describe the behavior of the bulk hydrated electron. For both the 4- and 16-water models, we investigate how the introduction of thermal motions alters the predicted absorption spectrum, vertical detachment energy, and resonance Raman spectrum of the simulated hydrated electron. We also present a set of structural criteria that can be used to numerically determine how cavity-like (or not) a particular hydrated electron model is. All of the results emphasize that the hydrated electron is a statistical object whose properties are inadequately captured using only a small number of explicit waters, and that a proper treatment of thermal fluctuations is critical to understanding the hydrated electron's chemical and physical behavior.

Gold Nanoparticle Formation via X-ray Radiolysis Investigated with Time-Resolved X-ray Liquidography

We report the generation of gold nanoparticles (AuNPs) from the aqueous solution of chloro(2,2',2″-terpyridine)gold(III) ion ([Au(tpy)Cl]²⁺) through X-ray radiolysis and optical excitation at a synchrotron. The original purpose of the experiment was to investigate the photoinduced structural changes of [Au(tpy)Cl]²⁺ upon 400 nm excitation using time-resolved X-ray liquidography (TRXL). Initially, the TRXL data did not show any signal that would suggest structural changes of the solute molecule, but after an induction time, the TRXL data started to show sharp peaks and valleys. In the early phase, AuNPs with two types of morphology, dendrites, and spheres, were formed by the reducing action of hydrated electrons generated by the X-ray radiolysis of water, thereby allowing the detection of TRXL data due to the laser-induced lattice expansion and relaxation of AuNPs. Along with the lattice expansion, the dendritic and spherical AuNPs were transformed into smaller, raspberry-shaped AuNPs of a relatively uniform size via ablation by the optical femtosecond laser pulse used for the TRXL experiment. Density functional theory calculations confirm that the reduction potential of the metal complex relative to the hydration potential of X-ray-generated electrons determines the facile AuNP formation observed for [Au(tpy)Cl]²⁺.

Fate of Dipole-Bound Anion States when Hydrated

Many strongly polar molecules can form an anion by attaching an electron to either an empty or half-filled valence-bound (VB) orbital or a so-called dipole-bound (DB) orbital. These two families of orbitals can be very different in their radial extent (the former are usually more compact, while the latter are quite diffuse) and in the degree to which they are affected by surrounding solvent molecules. In this study, the effects of hydration (representative of strong solvation) on the DB state of a model polar species are investigated with an eye toward determining whether this state is stabilized or even persists when a few to 100 water molecules surround the polar molecule. It is found that in the presence of up to ca. 10-12 water molecules, the excess electron can remain in a DB orbital. However, once there are enough water molecules to form a complete first hydration shell (or more), the excess electron migrates into an orbital localized on the outer surface of the water solvent cage. These findings have implications on the possible role of DB states as doorways to facilitating electron attachment and subsequent electron transfer to VB states. It is shown that even when the electron is bound to the surface of the surrounding solvent, the dipole potential of the solute molecule can influence where on the surface the electron binds. It is also illustrated that using continuum dielectric methods to describe the hydration of DB states is fraught with danger because much of the outermost electron density in such states penetrates outside the boundary of the cavity used in these methods.

The Fluxional Nature of the Hydrated Electron: Energy and Entropy Contributions to Aqueous Electron Free Energies

There has been a great deal of recent controversy over the structure of the hydrated electron and whether it occupies a cavity or contains a significant number of interior waters (noncavity). The questions we address in this work are, from a free energy perspective, how different are these proposed structures? Do the different structures all lie along a single continuum, or are there significant differences (i.e., free energy barriers) between them? To address these questions, we have performed a series of one-electron calculations using umbrella sampling with quantum biased molecular dynamics along a coordinate that directly reflects the number of water molecules in the hydrated electron's interior. We verify that a standard cavity model of the hydrated electron behaves essentially as a hard sphere: the model is dominated by repulsion at short range such that water is expelled from a local volume around the electron, leading to a water solvation shell like that of a pseudohalide ion. The repulsion is much larger than thermal energies near room temperature, explaining why such models exhibit properties with little temperature dependence. On the other hand, our calculations reveal that a noncavity model is highly fluxional, meaning that thermal motions cause the number of interior waters to fluctuate from effectively zero (i.e., a cavity-type electron) to potentially above the bulk water density. The energetic contributions in the noncavity model are still repulsive in the sense that they favor cavity formation, so the fluctuations in structure are driven largely by entropy: the entropic cost for expelling water from a region of space is large enough that some water is still driven into the electron's interior. As the temperature is lowered and entropy becomes less important, the noncavity electron's structure is predicted to become more cavity-like, consistent with the observed temperature dependence of the hydrated electron's properties. Thus, we argue that although the specific noncavity model we study overestimates the preponderance of fluctuations involving interior water molecules, with appropriate refinements to correctly capture the true average number of interior waters and molar solvation volume, a fluxional model likely makes the most sense for understanding the various experimental properties of the hydrated electron.

Real-time observation of water radiolysis and hydrated electron formation induced by extreme-ultraviolet pulses

The dominant pathway of radiation damage begins with the ionization of water. Thus far, however, the underlying primary processes could not be conclusively elucidated. Here, we directly study the earliest steps of extreme ultraviolet (XUV)-induced water radiolysis through one-photon excitation of large water clusters using time-resolved photoelectron imaging. Results are presented for H2O and D2O clusters using femtosecond pump pulses centered at 133 or 80 nm. In both excitation schemes, hydrogen or proton transfer is observed to yield a prehydrated electron within 30 to 60 fs, followed by its solvation in 0.3 to 1.0 ps and its decay through geminate recombination on a ∼10-ps time scale. These results are interpreted by comparison with detailed multiconfigurational non-adiabatic ab-initio molecular dynamics calculations. Our results provide the first comprehensive picture of the primary steps of radiation chemistry and radiation damage and demonstrate new approaches for their study with unprecedented time resolution.

Ring Formation and Hydration Effects in Electron Attachment to Misonidazole

We study the reactivity of misonidazole with low-energy electrons in a water environment combining experiment and theoretical modelling. The environment is modelled by sequential hydration of misonidazole clusters in vacuum. The well-defined experimental conditions enable computational modeling of the observed reactions. While the NO 2 - dissociative electron attachment channel is suppressed, as also observed previously for other molecules, the OH - channel remains open. Such behavior is enabled by the high hydration energy of OH - and ring formation in the neutral radical co-fragment. These observations help to understand the mechanism of bio-reductive drug action. Electron-induced formation of covalent bonds is then important not only for biological processes but may find applications also in technology.

Reaction of Electrons with DNA: Radiation Damage to Radiosensitization

This review article provides a concise overview of electron involvement in DNA radiation damage. The review begins with the various states of radiation-produced electrons: Secondary electrons (SE), low energy electrons (LEE), electrons at near zero kinetic energy in water (quasi-free electrons, (e-qf)) electrons in the process of solvation in water (presolvated electrons, e-pre), and fully solvated electrons (e-aq). A current summary of the structure of e-aq, and its reactions with DNA-model systems is presented. Theoretical works on reduction potentials of DNA-bases were found to be in agreement with experiments. This review points out the proposed role of LEE-induced frank DNA-strand breaks in ion-beam irradiated DNA. The final section presents radiation-produced electron-mediated site-specific formation of oxidative neutral aminyl radicals from azidonucleosides and the evidence of radiosensitization provided by these aminyl radicals in azidonucleoside-incorporated breast cancer cells.

Partial Molar Volume of the Hydrated Electron

The partial molar volume of the hydrated electron was investigated with pulse radiolysis and transient absorption by measuring the pressure dependence of the equilibrium constant for e^-_aq + NH₄⁺ ⇔ H + NH₃. At 2 kbar pressure, the equilibrium constant decreases relative to 1 bar by only 6%. Using tabulated molar volumes for ammonia and ammonium, we have the result V̅(e^-_aq) - V̅(H) = 11.3 cm³/mol at 25 °C, confirming that V̅(e^-_aq) is positive and even larger than the hydrophobic H atom. Assuming on the basis of recent molecular dynamics simulations that the molar volume of the H atom is somewhat less than that of H₂, we estimate V̅(e^-_aq) = 26 ± 6 cm³/mol. The positive molar volume is consistent with an electron that exists largely in a small solvent void (cavity), ruling out a recent model ( Larsen , R. E. ; Glover , W. J. ; Schwartz , B. J. Science 2010 , 329 , 65 - 69 ) that suggests a noncavity structure with negative molar volume.

Selectivity in Electron Attachment to Water Clusters

Electron attachment onto water clusters to form water cluster anions is studied by varying the point of electron attachment along a molecular beam axis and probing the produced cluster anions using photoelectron spectroscopy. The results show that the point of electron attachment has a clear effect on the final distribution of isomers for a cluster containing 78 water molecules, with isomer I formed preferentially near the start of the expansion and isomer II formed preferentially once the molecular beam has progressed for several millimeters. These changes can be accounted for by the cluster growth rate along the beam. Near the start of the expansion, cluster growth is proceeding rapidly with condensing water molecules solvating the electron, while further along the expansion, the growth has terminated and electrons are attached to large and cold preformed clusters, leading to the isomer associated with a loosely bound surface state.

Quantitative structure-activity relationship models for predicting reaction rate constants of organic contaminants with hydrated electrons and their mechanistic pathways

The hydrated electron (e_aq^-)-based reduction processes are promising for removing organic pollutants in water engineering systems. The reductive kinetics, especially the second order rate constants ( [Formula: see text] ) of e_aq^- with organic compounds, is important for evaluating and modeling the advanced reduction processes. In this study, the [Formula: see text] values for aliphatic compounds and phenyl-based compounds are, for the first time, modeled by the quantitative structure-activity relationship (QSAR) method. The structural features governing the reactivity of two classes of organic compounds toward e_aq^- were revealed, and the energy of the lowest unoccupied molecular orbital (E_LUMO), one-electron reduction potential (E_RED) and polarizability (α) were found to be the important molecular parameters in both two models. The built QSAR models provide robust predictive tools for estimating the removal of emerging pollutants using e_aq^- during wastewater treatment processes. Additionally, quantum chemical calculations were employed to probe into the mechanism and feasibility of the single electron transfer (SET) pathway in the e_aq^--based reduction process. The thermodynamic investigation suggests that the compounds with electron-withdrawing groups tend to possess higher [Formula: see text] and lower Gibbs free energy (ΔG_SET) and Gibbs free energies of activation (∆^‡G_SET^∘) than the ones with electron-donating groups, indicating the SET process occurs more readily. It is also found that the refractory halogenated compounds can achieve dehalogenation via the SET pathway.

Dynamics of the Bulk Hydrated Electron from Many-Body Wave-Function Theory

The structure of the hydrated electron is a matter of debate as it evades direct experimental observation owing to the short life time and low concentrations of the species. Herein, the first molecular dynamics simulation of the bulk hydrated electron based on correlated wave‐function theory provides conclusive evidence in favor of a persistent tetrahedral cavity made up by four water molecules, and against the existence of stable non‐cavity structures. Such a cavity is formed within less than a picosecond after the addition of an excess electron to neat liquid water, with less regular cavities appearing as intermediates. The cavities are bound together by weak H−H bonds, the number of which correlates well with the number of coordinated water molecules, each type of cavity leaving a distinct spectroscopic signature. Simulations predict regions of negative spin density and a gyration radius that are both in agreement with experimental data.

Multistep and multiscale electron transfer and localization dynamics at a model electrolyte/metal interface

The lifetime, coupling, and localization dynamics of electronic states in molecular films near metal electrodes fundamentally determine their propensity to act as precursors or reactants in chemical reactions, crucial for a detailed understanding of charge transport and degradation mechanisms in batteries. In the current study, we investigate the formation dynamics of small polarons and their role as intermediate electronic states in thin films of dimethyl sulfoxide (DMSO) on Cu(111) using time- and angle-resolved two-photon photoemission spectroscopy. Upon photoexcitation, a delocalized DMSO electronic state is initially populated two monolayers from the Cu surface, becoming a small polaron on a 200 fs time scale, consistent with localization due to vibrational dynamics of the DMSO film. The small polaron is a precursor state for an extremely long-lived and weakly coupled multilayer electronic state, with a lifetime of several seconds, thirteen orders of magnitude longer than the small polaron. Although the small polaron in DMSO has a lifetime of 140 fs, its role as a precursor state for long-lived electronic states could make it an important intermediate in multistep battery reactivity.

Structure of the aqueous electron

A cavity or excluded-volume structure best explains the experimental properties of the aqueous or “hydrated” electron.

Electron affinity of liquid water

Understanding redox and photochemical reactions in aqueous environments requires a precise knowledge of the ionization potential and electron affinity of liquid water. The former has been measured, but not the latter. We predict the electron affinity of liquid water and of its surface from first principles, coupling path-integral molecular dynamics with ab initio potentials, and many-body perturbation theory. Our results for the surface (0.8 eV) agree well with recent pump-probe spectroscopy measurements on amorphous ice. Those for the bulk (0.1-0.3 eV) differ from several estimates adopted in the literature, which we critically revisit. We show that the ionization potential of the bulk and surface are almost identical; instead their electron affinities differ substantially, with the conduction band edge of the surface much deeper in energy than that of the bulk. We also discuss the significant impact of nuclear quantum effects on the fundamental gap and band edges of the liquid.

Reactivity and energy level of a localized hole in liquid water

Reaction and redox level of hole capture in liquid water from first principles.

A DFT study on the mechanism of photoselective catalytic reduction of 4-bromobenzaldehyde in different solvents employing an OH-defected TiO2 cluster model

Density functional theory calculations are employed to study the mechanism of photoselective catalytic reduction of 4-bromobenzaldehyde (4-BBA) in acetonitrile and in ethanol solvents. A totally relaxed Ti₃O₉H₆ cluster model is proposed to represent titanium dioxide (TiO₂) surfaces. The reduction selectivity of an adsorbed 4-BBA molecule on Ti₃O₉H₆ has been investigated. Owing to the difference in the proton and H atom donating capabilities between explicit CH₃CN and C₂H₅OH solvent molecules, the photocatalytic reduction of 4-BBA is the debromination process in acetonitrile, whereas in ethanol it is the carbonyl reduction process. Therefore 4-BBA can be selectively reduced to benzaldehyde in acetonitrile and 4-bromobenzyl alcohol in ethanol, respectively. Our computational results have verified the reaction mechanism proposed by experiments and show that the debromination of 4-BBA would be efficient if both 4-BBA and Ti₃O₉H₆ have an extra photoelectron. The Ti₃O₉H₆ cluster, playing a role as a hydrogen source and a bridge to transfer photoelectrons from bulk TiO₂, would have potential to be an ideal molecular model for understanding photocatalytic reactions on the TiO₂ surface.