Abstract: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.

Molecular Dynamics Simulation Of Fe2+ Hydration Process Using The L-J(12/6) Interaction Potential Model

The Affection to Water Clusters on Temperatures and Electric Fields by Molecular Dynamics Simulation

Molecular dynamics study on ionic hydration

Molecular Dynamics Simulation of Zinc Ion in Water with an Ab Initio Based Neural Network Potential

Local Hydrated Structure of an Fe 2+ /fe 3+ Aqueous Solution: an Investigation Using a Combination of Molecular Dynamics and X-ray Absorption Fine Structure Methods

Importance of Van Der Waals Effects on the Hydration of Metal Ions from the Hofmeister Series

Molecular Dynamics Simulation of Heavy Metal Ions in Aqueous Solution Using Lennard-Jones 12-6 Potential

Revisiting the structure and dynamics of hydrated Cd2+ in aqueous solutions: Insights from the RI‐SCS‐MP2/MM molecular dynamics simulation

Molecular Dynamics Simulations of Electroosmosis in Perfluorosulfonic Acid Polymer

Parameterization of Highly Charged Metal Ions Using the 12-6-4 LJ-Type Nonbonded Model in Explicit Water

Analysis of factors influencing hydration site prediction based on molecular dynamics simulations

Molecular dynamics simulations study for the effect of cations hydration on the surface tension of the electrolyte solutions

Hydration Structure of Sodium and Potassium Ions with DFT-MD

Molecular Dynamics Simulation for Dmf Aqueous Solution

Molecular dynamics simulation of ionic concentration effects on ionic hydration in aqueous KCl solutions

Structural and Dynamic Properties of Solvated Hydroxide and Hydronium Ions in Water from Ab Initio Modeling.

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

Bridging atomistic simulations and thermodynamic hydration models of aqueous electrolyte solutions

Molecular Dynamics Simulation, Ab Initio Calculation, and Size-Selected Anion Photoelectron Spectroscopy Study of Initial Hydration Processes of Calcium Chloride

Quantifying the Structure of Water and Hydrated Monovalent Ions by Density Functional Theory-Based Molecular Dynamics.