Dynamics of the charge transfer to solvent process in aqueous iodide

Charge-transfer-to-solvent states in aqueous halides are ideal systems for studying the electron-transfer dynamics to the solvent involving a complex interplay between electronic excitation and solvent polarization. Despite extensive experimental investigations, a full picture of the charge-transfer-to-solvent dynamics has remained elusive. Here, we visualise the intricate interplay between the dynamics of the electron and the solvent polarization occurring in this process. Through the combined use of ab initio molecular dynamics and machine learning methods, we investigate the structure, dynamics and free energy as the excited electron evolves through the charge-transfer-to-solvent process, which we characterize as a sequence of states denoted charge-transfer-to-solvent, contact-pair, solvent-separated, and hydrated electron states, depending on the distance between the iodine and the excited electron. Our assignment of the charge-transfer-to-solvent states is supported by the good agreement between calculated and measured vertical binding energies. Our results reveal the charge transfer process in terms of the underlying atomic processes and mechanisms.

Partial Molar Solvation Volume of the Hydrated Electron Simulated Via DFT

Different simulation models of the hydrated electron produce different solvation structures, but it has been challenging to determine which simulated solvation structure, if any, is the most comparable to experiment. In a recent work, Neupane et al. [J. Phys. Chem. B 2023, 127, 5941-5947] showed using Kirkwood-Buff theory that the partial molar volume of the hydrated electron, which is known experimentally, can be readily computed from an integral over the simulated electron-water radial distribution function. This provides a sensitive way to directly compare the hydration structure of different simulation models of the hydrated electron with experiment. Here, we compute the partial molar volume of an ab-initio-simulated hydrated electron model based on density-functional theory (DFT) with a hybrid functional at different simulated system sizes. We find that the partial molar volume of the DFT-simulated hydrated electron is not converged with respect to the system size for simulations with up to 128 waters. We show that even at the largest simulation sizes, the partial molar volume of DFT-simulated hydrated electrons is underestimated by a factor of 2 with respect to experiment, and at the standard 64-water size commonly used in the literature, DFT-based simulations underestimate the experimental solvation volume by a factor of ∼3.5. An extrapolation to larger box sizes does predict the experimental partial molar volume correctly; however, larger system sizes than those explored here are currently intractable without the use of machine-learned potentials. These results bring into question what aspects of the predicted hydrated electron radial distribution function, as calculated by DFT-based simulations with the PBEh-D3 functional, deviate from the true solvation structure.

Tracking Cavity Formation in Electron Solvation: Insights from X-ray Spectroscopy and Theory

We present time-resolved X-ray absorption spectra of ionized liquid water and demonstrate that OH radicals, H3O+ ions, and solvated electrons all leave distinct X-ray-spectroscopic signatures. Particularly, this allows us to characterize the electron solvation process through a tool that focuses on the electronic response of oxygen atoms in the immediate vicinity of a solvated electron. Our experimental results, supported by ab initio calculations, confirm the formation of a cavity in which the solvated electron is trapped. We show that the solvation dynamics are governed by the magnitude of the random structural fluctuations present in water. As a consequence, the solvation time is highly sensitive to temperature and to the specific way the electron is injected into water.

Spectroscopy and dynamics of the hydrated electron at the water/air interface

The hydrated electron, e-(aq), has attracted much attention as a central species in radiation chemistry. However, much less is known about e-(aq) at the water/air surface, despite its fundamental role in electron transfer processes at interfaces. Using time-resolved electronic sum-frequency generation spectroscopy, the electronic spectrum of e-(aq) at the water/air interface and its dynamics are measured here, following photo-oxidation of the phenoxide anion. The spectral maximum agrees with that for bulk e-(aq) and shows that the orbital density resides predominantly within the aqueous phase, in agreement with supporting calculations. In contrast, the chemistry of the interfacial hydrated electron differs from that in bulk water, with e-(aq) diffusing into the bulk and leaving the phenoxyl radical at the surface. Our work resolves long-standing questions about e-(aq) at the water/air interface and highlights its potential role in chemistry at the ubiquitous aqueous interface.

Mechanisms and Opportunities for Rational In Silico Design of Enzymes to Degrade Per- and Polyfluoroalkyl Substances (PFAS)

Per and polyfluoroalkyl substances (PFAS) present a unique challenge to remediation techniques because their strong carbon-fluorine bonds make them difficult to degrade. This review explores the use of in silico enzymatic design as a potential PFAS degradation technique. The scope of the enzymes included is based on currently known PFAS degradation techniques, including chemical redox systems that have been studied for perfluorooctanesulfonic acid (PFOS) and perfluorooctanoic acid (PFOA) defluorination, such as those that incorporate hydrated electrons, sulfate, peroxide, and metal catalysts. Bioremediation techniques are also discussed, namely the laccase and horseradish peroxidase systems. The redox potential of known reactants and enzymatic radicals/metal-complexes are then considered and compared to potential enzymes for degrading PFAS. The molecular structure and reaction cycle of prospective enzymes are explored. Current knowledge and techniques of enzyme design, particularly radical-generating enzymes, and application are also discussed. Finally, potential routes for bioengineering enzymes to enable or enhance PFAS remediation are considered as well as the future outlook for computational exploration of enzymatic in situ bioremediation routes for these highly persistent and globally distributed contaminants.

Multifunctional Roles of Clathrate Hydrate Nanoreactors for CO2 Reduction

In this study, we explore a possible platform for the CO2 reduction (CO2 R) in one of water's solid phases, namely clathrate hydrates (CHs), by ab initio molecular dynamics and well-tempered metadynamics simulations with periodic boundary conditions. We found that the stacked H2 O nanocages in CHs help to initialize CO2 R by increasing the electron-binding ability of CO2 . The substantial CO2 R processes are further influenced by the hydrogen bond networks in CHs. The first intermediate CO2 - in this process can be stabilized through cage structure reorganization into the H-bonded [CO2 - ⋅⋅⋅H-OHcage ] complex. Further cooperative structural dynamics enables the complex to convert into a vital transient [CO2 2- ⋅⋅⋅H-OHcage ] intermediate in a low-barrier disproportionation-like process. Such a highly reactive intermediate spontaneously triggers subsequent double proton transfer along its tethering H-bonds, finally converting it into HCOOH. These hydrogen-bonded nanoreactors feature multiple functions in facilitating CO2 R such as confining, tethering, H-bond catalyzing and proton pumping. Our findings have a general interest and extend the knowledge of CO2 R into porous aqueous systems.

Empirically Optimized One-Electron Pseudopotential for the Hydrated Electron: A Proof-of-Concept Study

Mixed quantum-classical molecular dynamics simulations have been important tools for studying the hydrated electron. They generally use a one-electron pseudopotential to describe the interactions of an electron with the water molecules. This approximation shows both the strength and weakness of the approach. On the one hand, it enables extensive statistical sampling and large system sizes that are not possible with more accurate ab initio molecular dynamics methods. On the other hand, there has (justifiably) been much debate about the ability of pseudopotentials to accurately and quantitatively describe the hydrated electron properties. These pseudopotentials have largely been derived by fitting them to ab initio calculations of an electron interacting with a single water molecule. In this paper, we present a proof-of-concept demonstration of an alternative approach in which the pseudopotential parameters are determined by optimizing them to reproduce key experimental properties. Specifically, we develop a new pseudopotential, using the existing TBOpt model as a starting point, which correctly describes the hydrated electron vertical detachment energy and radius of gyration. In addition to these properties, this empirically optimized model displays a significantly modified solvation structure, which improves, for example, the prediction of the partial molar volume.

Influence of H-bonds on the photoionization of aromatic chromophores in water: The aniline molecule

We have conducted time-resolved experiments (pump-probe and pump-repump-probe) on a model aromatic chromophore, aniline, after excitation in water at 267 nm. In the initial spectra recorded, in addition to the absorption corresponding to the bright ππ* excitation, the fingerprint of a transient state with the electron located on the solvent molecule is identified. We postulate that the latter corresponds to the πσ* state along the N-H bond, whose complete relaxation with a ∼500 ps lifetime results in the formation of the fully solvated electron and cation. This ionization process occurs in parallel with the ππ* photophysical channel that yields the characteristic ∼1 ns fluorescence lifetime. The observed branched pathway is rationalized in terms of the different H-bonds that the water establishes with the amino group. The proposed mechanism could be common for aromatics in water containing N-H or O-H bonds and would allow the formation of separated charges after excitation at the threshold of their electronic absorptions.

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.

Distortion Correction of Low-Energy Photoelectron Spectra of Liquids Using Spectroscopic Data for Solvated Electrons

Time-resolved photoelectron spectroscopy (TRPES) enables real-time observation of ultrafast electronic dynamics in solutions. When extreme ultraviolet (EUV) probe pulses are employed, they can ionize solutes from all electronic states involved in the dynamics. However, EUV pulses also produce a strong ionization signal from a solvent that is typically 6 orders of magnitude greater than the pump-probe photoelectron signal of solutes. Alternatively, UV probe pulses enable highly sensitive and selective observation of photoexcited solutes because typical solvents such as water are transparent to UV radiation. An obstacle in such UV-TRPES measurements is spectral distortion caused by electron scattering and a yet to be identified mechanism in liquids. We have previously proposed the spectral retrieval (SR) method as an a posteriori approach to removing the distortion and overcoming this difficulty in UV-TRPES; however, its accuracy has not yet been verified by comparison with EUV-TRPES results. In the present study, we perform EUV-TRPES for charge transfer reactions in water, methanol, and ethanol, and verify SR analysis of UV-TRPES. We also estimate a previously undetermined energy-dependent intensity factor and expand the basis sets for SR analysis. The refined SR method is employed for reanalyzing the UV-TRPES data for the formation and relaxation dynamics of solvated electrons in various systems. The electron binding energy distributions for solvated electrons in liquid water, methanol, and ethanol are confirmed to be Gaussian centered at 3.78, 3.39, and 3.25 eV, respectively, in agreement with Nishitani et al. [ Sci. Adv. 2019, 5(8), eaaw6896]. An effective energy gap between the conduction band and the vacuum level at the gas-liquid interface is estimated to be 0.2 eV for liquid water and 0.1 eV for methanol and ethanol.

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.

Birth of the Hydrated Electron via Charge-Transfer-to-Solvent Excitation of Aqueous Iodide

A primary means to generate hydrated electrons in laboratory experiments is excitation to the charge-transfer-to-solvent (CTTS) state of a solute such as I^-(aq), but this initial step in the genesis of e^-(aq) has never been simulated directly using ab initio molecular dynamics. We report the first such simulations, combining ground- and excited-state simulations of I^-(aq) with a detailed analysis of fluctuations in the Coulomb potential experienced by the nascent solvated electron. What emerges is a two-step picture of the evolution of e^-(aq) starting from the CTTS state: I^-(aq) + hν → I^-*(aq) → I^•(aq) + e^-(aq). Notably, the equilibrated ground state of e^-(aq) evolves from I^-*(aq) without any nonadiabatic transitions, simply as a result of solvent reorganization. The methodology used here should be applicable to other photochemical electron transfer processes in solution, an important class of problems directly relevant to photocatalysis and energy transfer.

Observation of a transient intermediate in the ultrafast relaxation dynamics of the excess electron in strong-field-ionized liquid water

A unified picture of the electronic relaxation dynamics of ionized liquid water has remained elusive despite decades of study. Here, we employ sub-two-cycle visible to short-wave infrared pump-probe spectroscopy and ab initio nonadiabatic molecular dynamics simulations to reveal that the excess electron injected into the conduction band (CB) of ionized liquid water undergoes sequential relaxation to the hydrated electron s ground state via an intermediate state, identified as the elusive p excited state. The measured CB and p-electron lifetimes are 0.26 ± 0.02 ps and 62 ± 10 fs, respectively. Ab initio quantum dynamics yield similar lifetimes and furthermore reveal vibrational modes that participate in the different stages of electronic relaxation, with initial relaxation within the dense CB manifold coupled to hindered translational motions whereas subsequent p-to-s relaxation facilitated by librational and even intramolecular bending modes of water. Finally, energetic considerations suggest that a hitherto unobserved trap state resides ~0.3-eV below the CB edge of liquid water. Our results provide a detailed atomistic picture of the electronic relaxation dynamics of ionized liquid water with unprecedented time resolution.

Evidence for Water Antibonding Orbital Mixing in the Hydrated Electron from Its Oxygen 1s X-ray Absorption Spectrum

The X-ray absorption spectrum (XAS) of the hydrated electron (e_(aq)^-) has been simulated using time-dependent density functional theory with a quantum mechanics/molecular mechanics description. A unique XAS peak at 533 eV is observed with an energy and intensity in quantitative agreement with recent time-resolved experiments, allowing its assignment as arising from water O1s transitions to the singly occupied molecular orbital (SOMO) in which the excess electron resides. The transitions acquire oscillator strength due to the SOMO comprising an admixture of a cavity-localized orbital and water 4a₁ and 2b₂ antibonding orbitals. The mixing of antibonding orbitals has implications for the strength of couplings between e_(aq)^- and intramolecular modes of water.

Investigation of the Failure of Marcus Theory for Hydrated Electron Reactions

Reactions of the hydrated electron with a wide variety of substrates have been found to exhibit unusually similar activation energies in a manner incompatible with Marcus electron transfer theory. Given the fundamental linear response assumption of Marcus theory, one possible explanation for this apparent failure is that the underlying free energy surfaces governing the reactions are not harmonic; i.e., hydrated electron structural fluctuations exhibit non-Gaussian behavior. In this work, we test this hypothesis by using simulations to calculate the hydrated electron vertical detachment energy distribution. We consider both cavity and noncavity models for the hydrated electron, between which the actual hydrated electron behavior is expected to lie. Our results identify a possible origin for non-Gaussian behavior of the hydrated electron but show that it is not of sufficient magnitude to explain the failure of Marcus theory to describe its reactions. Thus, other explanations must be sought.

Temperature Dependent Properties of the Aqueous Electron

The temperature‐dependent properties of the aqueous electron have been extensively studied using mixed quantum‐classical simulations in a wide range of thermodynamic conditions based on one‐electron pseudopotentials. While the cavity model appears to explain most of the physical properties of the aqueous electron, only a non‐cavity model has so far been successful in accounting for the temperature dependence of the absorption spectrum. Here, we present an accurate and efficient description of the aqueous electron under various thermodynamic conditions by combining hybrid functional‐based molecular dynamics, machine learning techniques, and multiple time‐step methods. Our advanced simulations accurately describe the temperature dependence of the absorption maximum in the presence of cavity formation. Specifically, our work reveals that the red shift of the absorption maximum results from an increasing gyration radius with temperature, rather than from global density variations as previously suggested.

Quantum simulations of neutral water clusters and singly-charged water cluster anions

We report a computational study of the structural and energetic properties of water clusters and singly-charged water cluster anions containing from 20 to 573 water molecules. We have used both a classical and a quantum description of the molecular degrees of freedom. Water intra and inter-molecular interactions have been modelled through the SPC/F model, while the water-excess electron interaction has been described via the well-known Turi-Borgis potential. We find that in general the quantum effects of the water degrees of freedom are small, but they do influence the cluster-size at which the excess electron stabilises inside the cluster, which occurs at smaller cluster sizes when quantum effects are taken into consideration.

Coherent polaron dynamics of electrons solvated in polar liquids

An electron solvated in a polar liquid is an elementary quantum system with properties governed by electric interactions with a fluctuating molecular environment. In the prevailing single particle picture, the quantum ground and excited states are determined by a self-consistent potential, as defined by the particular local configuration of the solvation shell. This description neglects collective many-body excitations, which arise from the coupling of electronic degrees of freedom and nuclear motions of the environment. While recent experiments have demonstrated collective nonequilbrium electronic-nuclear motion, i.e. polaron excitations in liquid water, their relevance in the broader context of polar liquids has remained unexplored. Here, we study the nonequilibrium dielectric response of the, compared to water, less polar alcohols isopropanol, and ethylene glycol, that also display a different hydrogen bond pattern. We demonstrate that ultrafast relaxation of photogenerated electrons impulsively induces coherent charge oscillations, which persist for some 10 ps. They emit electric waves in a frequency range from 0.1 to 2 THz, depending on electron concentration. Oscillation frequencies and line shapes are reproduced by a unified polaron picture for alcohols and water, which is based on a Clausius-Mossotti local field approach for the THz dielectric function. The analysis suggests a longitudinal character of many-body polaron excitations and a weak coupling to transverse excitations, supported by the underdamped character of charge oscillations. Polaron dynamics are governed by the long-range Coulomb interaction between an excess electron and several thousands of polar solvent molecules, while local electron solvation geometries play a minor role.

Field-Induced Electron Generation in Water: Solvation Dynamics and Many-Body Interactions

The solvated electron represents an elementary quantum system in a liquid environment. Electrons solvated in water have raised strong interest because of their prototypical properties, their role in radiation chemistry, and their relevance for charge separation and transport. Nonequilibrium dynamics of photogenerated electrons in water occur on ultrafast time scales and include charge transfer, localization, and energy dissipation processes. We present new insight into the role of fluctuating electric fields of the liquid for generating electrons in the presence of an external terahertz field and address polaronic many-body properties of solvated electrons. This Perspective combines a review of recent results from experiment and theory with a discussion of basic electric interactions of electrons in water.

Multicomponent Quantum Mechanics/Molecular Mechanics Study of Hydrated Positronium

We propose a model for solvated positronium (Ps) atoms in water, based on the sequential quantum mechanics/molecular mechanics (s-QM/MM) protocol. We developed a Lennard-Jones force field to account for Ps-water interactions in the MM step. The repulsive term was obtained from a previously reported model for the solvated electron, while the dispersion constant was derived from the Slater-Kirkwood formula. The force field was employed in classical Monte Carlo (MC) simulations to generate Ps-solvent configurations in the NpT ensemble, while the quantum properties were computed with the any-particle molecular orbital method in the subsequent QM step. Our approach is general, as it can be applied to other liquids and materials. One basically needs to describe the solvated electron in the environment of interest to obtain the Ps solvation model. The thermodynamical properties computed from the MC simulations point out similarities between the solvation of Ps and noble gas atoms, hydrophobic solutes that form clathrate structures. We performed convergence tests for the QM step, with particular attention to the choice of basis set and expansion centers for the positronic and electronic subsystems. Our largest model was composed of the Ps atom and 22 water molecules in the QM region, corresponding to the first solvation shell, surrounded by 128 molecules described as point charges. The mean electronic and positronic vertical detachment energies were (4.73 ± 0.04) eV and (5.33 ± 0.04) eV, respectively. The latter estimates were computed with Koopmans' theorem corrected by second-order self-energies, for a set of statistically uncorrelated MC configurations. While the Hartree-Fock wave functions do not properly account for the annihilation rates, they were useful for numerical tests, pointing out that annihilation is more sensitive to the choice of basis sets and expansion centers than the detachment energies. We further explored a model with reduced solute cavity size by changing the Ps-solvent force field. Although the pick-off annihilation lifetimes were affected by the cavity size, essentially the same conclusions were drawn from both models.

Nonlinear electronic excitation in water under proton irradiation: a first principles study

Nonlinear dynamics of electronic excitation bridge physical and physicochemical stages of water radiolysis under proton irradiation, a multi-scale physicochemical process that is fundamental to a wide range of technological and medical applications of high-energy protons. We study the spatial and temporal changes to the excited holes generated in this ionization event using first-principles theory simulation. A significant majority of holes are formed in the immediate vicinity of the irradiating proton paths, and these holes decay rapidly, while secondary excitations are simultaneously induced in regions farther away. While the hole population remains constant, the observed spatially spreading hole distribution cannot be described as concentration-dependent diffusion current. Compared to the primary excitation induced by the irradiating protons, the secondary excitation farther away is somewhat less energetic. The first-principles theory simulation here provides a detailed description of how the primary excitation by proton radiation precedes the formation of cationic holes, which undergo ultrafast chemical processes in water radiolysis.

Shallow and deep trap states of solvated electrons in methanol and their formation, electronic excitation, and relaxation dynamics

We present condensed-phase first-principles molecular dynamics simulations to elucidate the presence of different electron trapping sites in liquid methanol and their roles in the formation, electronic transitions, and relaxation of solvated electrons (e_met⁻) in methanol. Excess electrons injected into liquid methanol are most likely trapped by methyl groups, but rapidly diffuse to more stable trapping sites with dangling OH bonds. After localization at the sites with one free OH bond (1OH trapping sites), reorientation of other methanol molecules increases the OH coordination number and the trap depth, and ultimately four OH bonds become coordinated with the excess electrons under thermal conditions. The simulation identified four distinct trapping states with different OH coordination numbers. The simulation results also revealed that electronic transitions of e_met⁻ are primarily due to charge transfer between electron trapping sites (cavities) formed by OH and methyl groups, and that these transitions differ from hydrogenic electronic transitions involving aqueous solvated electrons (e_aq⁻). Such charge transfer also explains the alkyl-chain-length dependence of the photoabsorption peak wavelength and the excited-state lifetime of solvated electrons in primary alcohols.

Condensed-phase first-principles molecular dynamics simulations elucidate the presence of different electron trapping sites in liquid methanol and their roles in the formation, electronic transitions, and relaxation of solvated electrons.

Flexible boundary layer using exchange for embedding theories. I. Theory and implementation

Embedding theory is a powerful computational chemistry approach to exploring the electronic structure and dynamics of complex systems, with Quantum Mechanical/Molecular Mechanics (QM/MM) being the prime example. A challenge arises when trying to apply embedding methodology to systems with diffusible particles, e.g., solvents, if some of them must be included in the QM region, for example, in the description of solvent-supported electronic states or reactions involving proton transfer or charge-transfer-to-solvent: without a special treatment, inter-diffusion of QM and MM particles will eventually lead to a loss of QM/MM separation. We have developed a new method called Flexible Boundary Layer using Exchange (FlexiBLE) that solves the problem by adding a biasing potential to the system that closely maintains QM/MM separation. The method rigorously preserves ensemble averages by leveraging their invariance to an exchange of identical particles. With a careful choice of the biasing potential and the use of a tree algorithm to include only important QM and MM exchanges, we find that the method has an MM-forcefield-like computational cost and thus adds negligible overhead to a QM/MM simulation. Furthermore, we show that molecular dynamics with the FlexiBLE bias conserves total energy, and remarkably, sub-diffusional dynamical quantities in the inner QM region are unaffected by the applied bias. FlexiBLE thus widens the range of chemistry that can be studied with embedding theory.

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.

Software for the frontiers of quantum chemistry: An overview of developments in the Q-Chem 5 package

This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange-correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear-electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an "open teamware" model and an increasingly modular design.

Probing Interfacial Effects on Ionization Energies: The Surprising Banality of Anion-Water Hydrogen Bonding at the Air/Water Interface

Liquid microjet photoelectron spectroscopy is an increasingly common technique to measure vertical ionization energies (VIEs) of aqueous solutes, but the interpretation of these experiments is subject to questions regarding sensitivity to bulk versus interfacial solvation environments. We have computed aqueous-phase VIEs for a set of inorganic anions, using a combination of molecular dynamics simulations and electronic structure calculations, with results that are in excellent agreement with experiment regardless of whether the simulation data are restricted to ions at the air/water interface or to those in bulk aqueous solution. Although the computed VIEs are sensitive to ion-water hydrogen bonding, we find that the short-range solvation structure is sufficiently similar in both environments that it proves impossible to discriminate between the two on the basis of the VIE, a conclusion that has important implications for the interpretation of liquid-phase photoelectron spectroscopy. More generally, analysis of the simulation data suggests that the surface activity of soft anions is largely a second or third solvation shell effect, arising from disruption of water-water hydrogen bonds and not from significant changes in first-shell anion-water hydrogen bonding.

Electron Attachment to Isolated Molecules as a Probe to Understand Mitochondrial Reductive Processes

This chapter describes the complementary experimental techniques Electron Transmission Spectroscopy and Dissociative Electron Attachment Spectroscopy, two of the most suitable means for investigating interactions between electrons and gas-phase molecules, resonance formation of temporary molecular negative ions, and their possible decay through the dissociative electron attachment (DEA) mechanism. The latter can be seen as the gas-phase counterpart of the transfer of a solvated electron in solution, accompanied by dissociation of the molecular anion, referred to as dissociative electron transfer (DET). DET takes place in vivo under reductive conditions, for instance, in the intermembrane space of mitochondria under interaction of xenobiotic molecules possessing high electron affinity with electrons "leaked" from the mitochondrial respiratory chain. A likely mechanism of the toxic activity of dichlorodiphenyltrichloroethane based on its DEA properties is briefly outlined, and compared with the well-established harmful effects of the model toxicant carbon tetrachloride ascribed to reductive dechlorination in a cellular ambient. A possible mechanism of the antioxidant activity of polyphenolic compounds present near the main site of superoxide anion production in mitochondria is also briefly discussed.

Natural Charge-Transfer Analysis: Eliminating Spurious Charge-Transfer States in Time-Dependent Density Functional Theory via Diabatization, with Application to Projection-Based Embedding

For many types of vertical excitation energies, linear-response time-dependent density functional theory (LR-TDDFT) offers a useful degree of accuracy combined with unrivaled computational efficiency, although charge-transfer excitation energies are often systematically and dramatically underestimated, especially for large systems and those that contain explicit solvent. As a result, low-energy electronic spectra of solution-phase chromophores often contain tens to hundreds of spurious charge-transfer states, making LR-TDDFT needlessly expensive in bulk solution. Intensity borrowing by these spurious states can affect intensities of the valence excitations, altering electronic bandshapes. At higher excitation energies, it is difficult to distinguish spurious charge-transfer states from genuine charge-transfer-to-solvent (CTTS) excitations. In this work, we introduce an automated diabatization that enables fast and effective screening of the CTTS acceptor space in bulk solution. Our procedure introduces "natural charge-transfer orbitals" that provide a means to isolate orbitals that are most likely to participate in a CTTS excitation. Projection of these orbitals onto solvent-centered virtual orbitals provides a criterion for defining the most important solvent molecules in a given excitation and be used as an automated subspace selection algorithm for projection-based embedding of a high-level description of the CTTS state in a lower-level description of its environment. We apply this method to an ab initio molecular dynamics trajectory of I^-(aq) and report the lowest-energy CTTS band in the absorption spectrum. Our results are in excellent agreement with the experiment, and only one-third of the water molecules in the I^-(H₂O)₉₆ simulation cell need to be described with LR-TDDFT to obtain excitation energies that are converged to <0.1 eV. The tools introduced herein will improve the accuracy, efficiency, and usability of LR-TDDFT in solution-phase environments.

Two-dimensional terahertz spectroscopy of condensed-phase molecular systems

Nonlinear terahertz (THz) spectroscopy relies on the interaction of matter with few-cycle THz pulses of electric field amplitudes up to megavolts/centimeter (MV/cm). In condensed-phase molecular systems, both resonant interactions with elementary excitations at low frequencies such as intra- and intermolecular vibrations and nonresonant field-driven processes are relevant. Two-dimensional THz (2D-THz) spectroscopy is a key method for following nonequilibrium processes and dynamics of excitations to decipher the underlying interactions and molecular couplings. This article addresses the state of the art in 2D-THz spectroscopy by discussing the main concepts and illustrating them with recent results. The latter include the response of vibrational excitations in molecular crystals up to the nonperturbative regime of light-matter interaction and field-driven ionization processes and electron transport in liquid water.

Photoelectron Spectroscopy of Liquid Water with Tunable Extreme-Ultraviolet Radiation: Effects of Electron Scattering

We report the first systematic photoelectron measurements of the three outer-valence bands of liquid water as a function of the ionizing photon energy in the near-threshold region. We use extreme-ultraviolet (XUV) radiation tunable between ∼17.1 and 35.6 eV, obtained through monochromatization of a high-harmonic source. We show that the absolute values of the apparent vertical ionization energies and their respective peak widths show a decreasing trend of their magnitudes with increasing photon energy close to the ionization threshold. We find that the observed effects do not only depend on the electron kinetic energy but are also different for the various outer-valence bands. These observations are consistent with, but not fully explained by, the effects of inelastic electron scattering.

Terahertz Polaron Oscillations of Electrons Solvated in Liquid Water

The terahertz (THz) response of solvated electrons in liquid water is studied in nonlinear ultrafast pump-probe experiments. Free electrons with concentrations from c_{e}=4 to 140×10^{-6} moles/liter are generated by high-field THz or near-infrared multiphoton excitation. The time-resolved change of the dielectric function as mapped by broadband THz pulses exhibits pronounced oscillations persisting up to 30 ps. Their frequency increases with electron concentration from 0.2 to 1.5 THz. The oscillatory response is assigned to impulsively excited coherent polarons involving coupled electron and water shell motions with a frequency set by the local electric field.

Simulating the ghost: quantum dynamics of the solvated electron

The nature of the bulk hydrated electron has been a challenge for both experiment and theory due to its short lifetime and high reactivity, and the need for a high-level of electronic structure theory to achieve predictive accuracy. The lack of a classical atomistic structural formula makes it exceedingly difficult to model the solvated electron using conventional empirical force fields, which describe the system in terms of interactions between point particles associated with atomic nuclei. Here we overcome this problem using a machine-learning model, that is sufficiently flexible to describe the effect of the excess electron on the structure of the surrounding water, without including the electron in the model explicitly. The resulting potential is not only able to reproduce the stable cavity structure but also recovers the correct localization dynamics that follow the injection of an electron in neat water. The machine learning model achieves the accuracy of the state-of-the-art correlated wave function method it is trained on. It is sufficiently inexpensive to afford a full quantum statistical and dynamical description and allows us to achieve accurate determination of the structure, diffusion mechanisms, and vibrational spectroscopy of the solvated electron.

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.

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.

Probing the Birth and Ultrafast Dynamics of Hydrated Electrons at the Gold/Liquid Water Interface via an Optoelectronic Approach

The hydrated electron has fundamental and practical significance in radiation and radical chemistry, catalysis, and radiobiology. While its bulk properties have been extensively studied, its behavior at solid/liquid interfaces is still unclear due to the lack of effective tools to characterize this short-lived species in between two condensed matter layers. In this study, we develop a novel optoelectronic technique for the characterization of the birth and structural evolution of solvated electrons at the metal/liquid interface with a femtosecond time resolution. Using this tool, we record for the first time the transient spectra (in a photon energy range from 0.31 to 1.85 eV) in situ with a time resolution of 50 fs revealing several novel aspects of their properties at the interface. Especially the transient species show state-dependent optical transition behaviors from being isotropic in the hot state to perpendicular to the surface in the trapped and solvated states. The technique will enable a better understanding of hot electron driven reactions at electrochemical interfaces.

Field-Induced Tunneling Ionization and Terahertz-Driven Electron Dynamics in Liquid Water

Liquid water at ambient temperature displays ultrafast molecular motions and concomitant fluctuations of very strong electric fields originating from the dipolar H₂O molecules. We show that such random intermolecular fields induce the tunnel ionization of water molecules, which becomes irreversible if an external terahertz (THz) pulse imposes an additional directed electric field on the liquid. Time-resolved nonlinear THz spectroscopy maps charge separation, transport, and localization of the released electrons on a few-picosecond time scale. The highly polarizable localized electrons modify the THz absorption spectrum and refractive index of water, a manifestation of a highly nonlinear response. Our results demonstrate how the interplay of local electric field fluctuations and external electric fields allows for steering charge dynamics and dielectric properties in aqueous systems.

Impact of Electron Solvation on Ice Structures at the Molecular Scale

Electron attachment and solvation at ice structures are well-known phenomena. The energy liberated in such events is commonly understood to cause temporary changes at such ice structures, but it may also trigger permanent modifications to a yet unknown extent. We determine the impact of electron solvation on D₂O structures adsorbed on Cu(111) with low-temperature scanning tunneling microscopy, two-photon photoemission, and ab initio theory. Solvated electrons, generated by ultraviolet photons, lead not only to transient but also to permanent structural changes through the rearrangement of individual molecules. The persistent changes occur near sites with a high density of dangling OD groups that facilitate electron solvation. We conclude that energy dissipation during solvation triggers permanent molecular rearrangement via vibrational excitation.

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.

Role of hemibonding in the structure and ultraviolet spectroscopy of the aqueous hydroxyl radical

The presence of a two-center, three-electron hemibond in the solvation structure of the aqueous hydroxl radical has long been debated, as its appearance can be sensitive to self-interaction error in density functional theory.