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

Deducing Reaction and Diffusion Depths of Near-Interfacial Solvated Electrons from pH-Dependent Product Evaporation

Near-interfacial electrons in water can be produced by bombarding an aqueous microjet in vacuum with gas-phase sodium atoms. These Na atoms immediately ionize into Na+ and es-, which can then react with surface-active molecules that preferentially populate the surface. We carried out these experiments by reacting es- with the surfactant benzyltrimethylammonium (BTMA+) in a 6.7 M LiBr/H2O microjet at 242 K as a function of pH between 1 and 5. The reaction products, trimethylamine (TMA) and benzyl radical, evaporate into the gas phase where they are detected by a mass spectrometer. We find that TMA evaporation sharply diminishes with increasing H+ concentration and is barely visible at pH = 1, while benzyl evaporation varies much less. These results indicate that TMA protonation overwhelms TMA evaporation at 0.1 M H+. Diffusion-reaction modeling matches the observed trends and predicts that es- reacts with BTMA+ within the top 20 Å at all pH values. However, TMA molecules that evaporate and escape protonation diffuse on average only over 20 Å at pH = 1 but over 1000 Å at pH = 5. These observations emphasize that the near-interfacial region provides a controllable reaction environment that is also an escape route for volatile intermediates, a route that is unavailable deep in the bulk. The competition between evaporation and reaction depends on the solubility of the intermediate, the location of its creation, and the propensity for secondary reactions.

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

Real-time structural dynamics of the ultrafast solvation process around photo-excited aqueous halides

This work investigates and describes the structural dynamics taking place following charge-transfer-to-solvent photo-abstraction of electrons from I- and Br- ions in aqueous solution following single- and 2-photon excitation at 202 nm and 400 nm, respectively. A Time-Resolved X-ray Solution Scattering (TR-XSS) approach with direct sensitivity to the structure of the surrounding solvent as the water molecules adopt a new equilibrium configuration following the electron-abstraction process is utilized to investigate the structural dynamics of solvent shell expansion and restructuring in real-time. The structural sensitivity of the scattering data enables a quantitative evaluation of competing models for the interaction between the nascent neutral species and surrounding water molecules. Taking the I0-O distance as the reaction coordinate, we find that the structural reorganization is delayed by 0.1 ps with respect to the photoexcitation and completes on a time scale of 0.5-1 ps. On longer time scales we determine from the evolution of the TR-XSS difference signal that I0: e- recombination takes place on two distinct time scales of ∼20 ps and 100 s of picoseconds. These dynamics are well captured by a simple model of diffusive evolution of the initial photo-abstracted electron population where the charge-transfer-to-solvent process gives rise to a broad distribution of electron ejection distances, a significant fraction of which are in the close vicinity of the nascent halogen atoms and recombine on short time scales.

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.

Key Role of Cations in Stabilizing Hydrogen Radicals for CO2-to-CO Conversion via a Reverse Water-Gas Shift Reaction

Electrochemically converting CO2 into valuable chemicals and fuels in acidic media is argued as a promising energy- and carbon-efficient route. Although several key roles of alkali cations have been unveiled, the alkali cation trends for CO2 reduction remain largely elusive. With decreasing cation size from Cs+ to Li+, here we show that the apparent proton diffusion coefficient in 3.0 M Li+ is tens-fold lower than in 3.0 M K+ and 3.0 M Cs+ acidic electrolytes. Although Li+ has the strongest inhibition ability for proton transport, it acts the worst for both the CO2-to-CO conversion and partial current density on Au catalysts. Unexpectedly, K+ with a higher proton transport performs the best for CO2-to-CO conversion. We thus revisit the roles of alkali cations and find that hydrated K+ can stabilize hydrogen radicals benefiting CO2 conversion at the electrode interface while for Li+ this is not the case. This study proposes that cation-stabilized atomic hydrogen assists in activating CO2 via a reverse water-gas shift route under electrochemical conditions.

Observation of Nanoconfinement Effect on the Kinetics of Hydrated Electron in the Nanoscale Water Pools of Water-AOT-Cyclohexane Microemulsions by Picosecond Pulse Radiolysis

The decay kinetics of the hydrated electron (eaq-) in aerosol OT (AOT)-based ternary microemulsions with pool sizes ranging from 0.34 to 4.85 nm were studied using picosecond pulse radiolysis coupled with transient absorption UV-vis spectroscopy. Electron transfer from oil to water and the subsequent solvation occurred within a time resolution of 7 ps. The decay kinetics of eaq- were accurately modeled using a double-exponential decay model, revealing the occurrence of two types of reactions, i.e., the recombination reaction at the water-oil interface and the radical-radical reactions in the water pools. The apparent lifetimes of both types of decays decreased significantly as the size of water pools decreased, indicating the influence of nanoconfinement effects. Moreover, the importance of the water-oil interface increased with decreasing water content, regardless of the presence or absence of NO3- as an electron scavenger in the water pools. Our findings provide a comprehensive understanding on the kinetics of the radiation reaction in AOT-based microemulsions.

Predicting Degradation Mechanisms in Lithium Bistriflimide "Water-In-Salt" Electrolytes For Aqueous Batteries

Aqueous solutions are crucial to most domains in biology and chemistry, including in energy fields such as catalysis and batteries. Water-in-salt electrolytes (WISEs), which extend the stability of aqueous electrolytes in rechargeable batteries, are one example. While the hype for WISEs is huge, commercial WISE-based rechargeable batteries are still far from reality, and there remain several fundamental knowledge gaps such as those related to their long-term reactivity and stability. Here, we propose a comprehensive approach to accelerating the study of WISE reactivity by using radiolysis to exacerbate the degradation mechanisms of concentrated LiTFSI-based aqueous solutions. We find that the nature of the degradation species depends strongly on the molality of the electrolye, with degradation routes driven by the water or the anion at low or high molalities, respectively. The main aging products are consistent with those observed by electrochemical cycling, yet radiolysis also reveals minor degradation species, providing a unique glimpse of the long-term (un)stability of these electrolytes.

The primary photolysis of aqueous carbonate di-anions

We study the primary photolysis dynamics of aqueous carbonate, CO₃^2-(aq), and hydrogen carbonate, HCO₃^-(aq), when they are excited at λ = 200 nm. The photolysis is recorded with sub-picosecond time resolution using UV pump-Vis probe and UV pump-IR probe transient absorption spectroscopy and interpreted with the aid of density functional theory calculations. When CO₃^2- is excited via single photon absorption at λ = 200 nm, Φ(t = 20 ps) = 82 ± 5% of the excited di-anions either detach an electron or dissociate. The electron detachment takes place from the excited state in t < 1 ps and forms ground state CO₃˙^- and e_aq^-. Dissociation occurs from both the electronic ground and excited states of CO₃^2-. Dissociation from the CO₃^2- excited state is assisted by water molecules and forms CO₂˙^-, OH˙ and OH^-. The dissociation occurs both directly from the Franck-Condon region in t < 1 ps and indirectly with a time constant of τ = 13.9 ± 0.5 ps as the excited state relaxes. Dissociation of vibrationally excited CO₃^2- molecules in the electronic ground state is also assisted by water molecules and forms CO₂ and two OH^- anions. The dissociation and subsequent vibrational relaxation of CO₂ occur with a time constant of τ = 10.2 ± 0.5 ps. The residual 1 - Φ(t = 20 ps) = 18 ± 5% of the excited CO₃^2- di-anions return by internal conversion to the equilibrated CO₃^2- ground state with a time constant of τ = 4.0 ± 0.4 ps. The extinction coefficient of aqueous hydrogen carbonate HCO₃^-(aq) at λ = 200 nm is an order of magnitude smaller than that of carbonate, so even though the hydrogen carbonate anions dominate the carbonate di-anions in the hydrogen carbonate solution, the primary photolysis of hydrogen carbonate is obscured by the photo-products of carbonate. Hence, we are unable to assess the primary photolysis of hydrogen carbonate. However, the weak one-photon absorption facilitates two-photon ionization of water, which forms hydronium, H₃O⁺, cations. The sudden increase in the acidity induced by two-photon ionization protonates the ground state hydrogen carbonate molecules, thus offering a rare spectroscopic glimpse of aqueous carbonic acid.

How Ions Break Local Symmetry: Simulations of Polarized Transient Hole Burning for Different Models of the Hydrated Electron in Contact Pairs with Na. J Phys Chem Lett 2023;14:3014-3022. [PMID: 36943261 DOI: 10.1021/acs.jpclett.3c00220]

The hydrated electron (e_aq^-) is known via polarized transient hole-burning (pTHB) experiments to have a homogeneously broadened absorption spectrum. Here, we explore via quantum simulation how the pTHB spectroscopy of different e_aq^- models changes in the presence of electrolytes. The idea is that cation-e_aq^- pairing can break the local symmetry and, thus, induce persistent inhomogeneity. We find that a "hard" cavity model shows a modest increase in the pTHB recovery time in the presence of salt, while a "soft" cavity model remains homogeneously broadened independent of the salt concentration. We also explore the orientational anisotropy of a fully ab initio density functional theory-based model of the e_aq^-, which is strongly inhomogeneously broadened without salt and which becomes significantly more inhomogeneously broadened in the presence of salt. The results provide a direct prediction for experiments that can distinguish between different models and, thus, help pin down the hydration structure and dynamics of the e_aq^-.

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.

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.