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

Simulating the Competitive Ion Pairing of Hydrated Electrons with Chaotropic Cations

Experiments show that the absorption spectrum of the hydrated electron (ehyd-) blue-shifts in electrolyte solutions compared with what is seen in pure water. This shift has been assigned to the ehyd-'s competitive ion-pairing interactions with the salt cation relative to the salt anion based on the ions' positions on the Hofmeister series. Remarkably, little work has been done investigating the ehyd-'s behavior when the salts have chaotropic cations, which should greatly change the ion-pairing interactions given that the ehyd- is a champion chaotrope. In this work, we remedy this by using mixed quantum/classical simulations to analyze the behavior of two different models of the ehyd- in aqueous RbF and RbI electrolyte solutions as a function of salt concentration. We find that the magnitude of the salt-induced spectral blue-shift is determined by a combination of the number of chaotropic Rb+ cations near the ehyd- and the number of salt anions near those cations so that the spectrum of the ehyd- directly reflects its local environment. We also find that the use of a soft-cavity ehyd- model predicts stronger competitive interactions with Rb+ relative to I- than a more traditional hard cavity model, leading to different predicted spectral shifts that should provide a way to distinguish between the two models experimentally. Our simulations predict that at the same concentration, salts with chaotropic cations should produce larger spectral blue-shifts than salts with kosmotropic cations. We also found that at high salt concentrations with chaotropic cations, the predicted blue-shift is greater when the salt anion is kosmotropic instead of chaotropic. Our goal is for this work to inspire experimentalists to make such measurements, which will help provide a spectroscopic means to distinguish between simulations models that predict different hydration structures for the ehyd-.

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.

Low thermal expansion of layered electrides predicted by density-functional theory

Layered electrides are a unique class of materials with anionic electrons bound in interstitial regions between thin, positively charged atomic layers. While density-functional theory is the tool of choice for computational study of electrides, there has to date been no systematic comparison of density functionals or dispersion corrections for their accurate simulation. There has also been no research into the thermomechanical properties of layered electrides, with computational predictions considering only static lattices. In this work, we investigate the thermomechanical properties of five layered electrides using density-functional theory to evaluate the magnitude of thermal effects on their lattice constants and cell volumes. We also assess the accuracy of five popular dispersion corrections with both planewave and numerical atomic orbital calculations.

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.

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.

Competitive Ion Pairing and the Role of Anions in the Behavior of Hydrated Electrons in Electrolytes

Experiments have shown that in the presence of electrolytes, the hydrated electron's absorption spectrum experiences a blue shift whose magnitude depends on both the salt concentration and chemical identity. Previous computer simulations have suggested that the spectral blue shift results from the formation of (cation, electron) contact pairs and that the concentration dependence arises because the number of cations simultaneously paired with the electron increases with increasing concentration. In this work, we perform new simulations to build an atomistic picture that explains the effect of salt identity on the observed hydrated electron spectral shifts. We simulate hydrated electrons in the presence of both monovalent (Na⁺) and divalent (Ca²⁺) cations paired with both Cl^- and a spherical species representing ClO₄^- anions. Our simulations reproduce the experimental observations that divalent ions produce larger blue shifts of the hydrated electron's spectrum than monovalent ions with the same anion and that perchlorate salts show enhanced blue shifts compared to chloride salts with the same cation. We find that these observations can be explained by competitive ion pairing. With small kosmotropic cations such as Na⁺ and Ca²⁺, aqueous chloride salts tend to form (cation, anion) contact pairs, whereas there is little ion pairing between these cations and chaotropic perchlorate anions. Hydrated electrons also strongly interact with these cations, but if the cations are also paired with anions, this affects the free energy of the electron-cation interaction. With chloride salts, hydrated electrons end up in complexes containing multiple cations plus a few anions as well as the electron. Repulsive interactions between the electron and the nearby Cl^- anions reduce the cation-induced spectral blue shift of the hydrated electron. With perchlorate salts, hydrated electrons pair with multiple cations without any associated anions, leading to the largest possible cation-induced spectral blue shift. We also see that the reason multivalent cations produce larger spectral blue shifts than monovalent cations is because hydrated electrons are able to simultaneously pair with a larger number of multivalent cations due to a larger free energy of interaction. Overall, the interaction of hydrated electrons with electrolytes fits well with the Hofmeister series, where the electron behaves as an anion that is slightly more able to break water's H-bond structure than chloride.

Trap-Seeking or Trap-Digging? Photoinjection of Hydrated Electrons into Aqueous NaCl Solutions

It is well-known that when excess electrons are injected into an aqueous solution, they localize and solvate in ∼1 ps. Still debated is whether localization occurs via "trap-digging", in which the electron carves out a suitable localization site, or by "trap-seeking", where the electron prefers to localize at pre-existing low-energy trap sites in solution. To distinguish between these two possible mechanisms, we study the localization dynamics of excess electrons in aqueous NaCl solutions using both ultrafast spectroscopy and mixed quantum-classical molecular dynamics simulations. By introducing pre-existing traps in the form of Na⁺ ions, we can use the cation-induced blue-shift of the hydrated electron's absorption spectrum to directly monitor the site of electron localization. Our experimental and computational results show that the electron prefers to localize directly at the sites of Na⁺ traps; the presence of concentrated electrolytes otherwise has little impact on the way trap-seeking hydrated electrons relax following injection.

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.

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.

Hydrated Electrons in High-Concentration Electrolytes Interact with Multiple Cations: A Simulation Study

Experimental studies have demonstrated that the hydrated electron's absorption spectrum undergoes a concentration-dependent blue-shift in the presence of electrolytes such as NaCl. The blue-shift increases roughly linearly at low salt concentration but saturates as the solubility limit of the salt is approached. Previous attempts to understand the origin of the concentration-dependent spectral shift using molecular simulation have only examined the interaction between the hydrated electron and a single sodium cation, and these simulations predicted a spectral blue-shift that was an order of magnitude larger than that seen experimentally. Thus, in this paper, we first explore the reasons for the exaggerated spectral blue-shift when a simulated hydrated electron interacts with a single Na⁺. We find that the issue arises from nonpairwise additivity of the Na⁺-e^- and H₂O-e^- pseudopotentials used in the simulation. This effect arises because the solvating water molecules donate charge into the empty orbitals of Na⁺, lowering the effective charge of the cation and thus reducing the excess electron-cation interaction. Careful analysis shows, however, that although this nonpairwise additivity changes the energetics of the electron-Na⁺ interaction, the forces between the electron, Na⁺, and water are unaffected, so that mixed quantum/classical (MQC) simulations produce the correct structure and dynamics. With this in hand, we then use MQC simulations to explore the behavior of the hydrated electron as an explicit function of NaCl salt concentration. We find that the simulations correctly reproduce the observed experimental spectral shifting behavior. The reason for the spectral shift is that as the electrolyte concentration increases, the average number of cations simultaneously interacting in contact pairs with the hydrated electron increases from 1.0 at low concentrations to ∼2.5 near the saturation limit. As the number of cations that interact with the electron increases, the cation/electron interactions becomes slightly weaker, so that the corresponding Na⁺-e^- distance increases with increasing salt concentration. We also find that the dielectric constant of the solution plays little role in the observed spectroscopy, so that the electrolyte-dependent spectral shifts of the hydrated electron are directly related to the concentration-dependent number of closely interacting cations.

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.