Pith. sign in

REVIEW 3 major objections 6 minor 87 references

Nature of Hydrated Electron in Varied Solvation Environments

T0 review · 3 major / 6 minor · reviewed 2026-08-07 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read An excess electron in water localizes only when free OH bonds focus on one spot and a stable second shell surrounds it.

desk verdict Useful comparative hybrid-AIMD study of the hydrated electron in 3D, 2D, and 1D water, but the 'necessary condition' claims outrun the single-trajectory sampling. read the letter →

arxiv 2506.07157 v1 pith:GU4OWWXG submitted 2025-06-08 physics.comp-ph

Ritama Kar , Nisanth N. Nair This is my paper
classification physics.comp-ph
keywords hydratedelectronexcesslocalizationabinitiomoleculardynamicshybriddensityfunctionaltheorywaterhydrogen-bondnetworklow-dimensionalicevacancy
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper asks what makes an excess electron settle into a stable, localized state in water rather than spread out. By simulating liquid water, hexagonal ice with and without a vacancy, a 2D water monolayer, and a 1D water chain at the hybrid density-functional level, the authors isolate two necessary conditions: dangling O–H groups that all point toward one small region, and a second solvation shell that stays intact. Bulk water satisfies both and yields the familiar $\approx 2.7$ Å hydrated-electron cavity; defective ice satisfies them at the vacancy; perfect ice, the monolayer, and the chain do not, and the electron stays delocalized. The point is that under-coordination alone does not explain solvated-electron stability; the orientation of dangling OH groups and a supporting second shell are both required.

What carries the argument

The central observable is the excess electron's spin density, tracked by the radius of gyration $r_g$ computed from its gyration tensor; localization is read as a stably small $r_g$ and a cavity in the electron-oxygen radial distribution. The argument is carried by comparing hydrogen-bonding topology across environments: the count and orientation of dangling OH bonds, first- and second-shell H-bond counts, and distributions of OH-vector and dipole angles relative to the electron. The computational enabler is a multiple-time-step integration scheme for hybrid density-functional forces, called resonance-free adaptively compressed exchange (RF-MTACE), which makes trajectories of tens of picoseconds with 40% exact exchange feasible.

What would settle it

Run the same set of environments with path-integral dynamics that include quantum nuclear motion: if a long-lived localized electron appears in the 2D monolayer or 1D chain, or disappears in defective ice, the claimed necessary conditions are wrong.

Watch

Extended reading notes

Core claim

The central claim is that electron localization in water requires, at the same time, a set of dangling OH bonds whose directions converge on a small spatial region and a dynamically stable second solvation shell around the emerging cavity. In liquid water the excess electron localizes within about 2 ps into a near-spherical cavity with $r_g\approx 2.7$ Å, disrupting the first shell to roughly three H-bonds per molecule while the second shell stays four-coordinated. In hexagonal ice with a molecular vacancy, the dangling OH groups at the defect point inward and the electron locks into the vacancy with $r_g\approx 3.4$ Å. In perfect ice, a 2D monolayer, and a 1D chain, the electron never forms such a state even where OH groups are under-coordinated, with mean $r_g$ values of about 6.35, 5.63, and 6.72 Å respectively.

Load-bearing premise

The conclusions assume that treating the hydrogen atoms as classical point particles leaves the dangling-OH dynamics essentially unchanged, so quantum proton motion could alter both localization and the role of the second shell.

Editorial extensions

If this is right

  • Bulk water with an excess electron should always show a two-shell solvation structure, with the radius of gyration staying near 2.7 Å at room temperature rather than growing over time.
  • A defect-free ice lattice is not a host for a stable hydrated electron, and introducing a single molecular vacancy is enough to create a stable trap with $r_g\approx 3.4$ Å.
  • In a 2D water monolayer between hydrophobic surfaces, excess electrons remain spread across the plane and only transiently compact, so interfacial hydrated-electron chemistry there is dominated by delocalized states.
  • In a 1D water chain, the electron coats the chain uniformly and no compact cage forms, so water in narrow hydrophobic tubes will not solvate electrons the way bulk water does.
  • The two necessary conditions give a structural test: any aqueous environment lacking either focused dangling OH bonds or a stable second shell cannot dynamically stabilize a hydrated electron.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Editorial inference: the necessary-condition claim implies a sharp crossover under confinement, so squeezing water to two or one dimensions should change the solvated electron from a compact quantum droplet into a delocalized band, a signature that could be searched for in absorption spectra of electrons in nanoporous or nanotube-confined water.
  • Editorial inference: the second-shell condition suggests a finite-size threshold, namely that water clusters too small to complete a second solvation shell should not host long-lived interior hydrated electrons, which is checkable in size-resolved cluster anion experiments.
  • Editorial inference: if proton quantum delocalization were included, the strict necessity of the second shell might soften, since quantum nuclei change OH orientational flexibility; this could be tested with path-integral hybrid density-functional simulations of the monolayer and chain.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 6 minor

Summary. The manuscript presents Born–Oppenheimer molecular dynamics simulations of an excess electron in four water environments — bulk liquid water, hexagonal ice (perfect and with a molecular vacancy), a two-dimensional monolayer confined between hydrophobic sheets, and a one-dimensional chain inside a nanotube — using the PBE0 hybrid functional with 40% exact exchange, an 80 Ry plane-wave cutoff, and the RF-MTACE multiple-time-step method. The electron is characterized through the radius of gyration of the spin density, radial distribution functions around the electron center of mass, hydrogen-bond populations, orientational distributions, and asphericity. Based on the observed localization in bulk water and at an ice vacancy, and the lack of stable localization in the perfect ice, monolayer, and chain, the authors conclude that dangling OH bonds all pointing toward a small spatial region are a necessary requirement for electron localization, and that the presence and stability of a second solvation shell is a necessary condition for the dynamic stabilization of the hydrated electron.

Significance. The study is valuable because it applies a consistent, state-of-the-art hybrid-functional protocol to a wide range of confining environments, and the bulk-water result (rg ≈ 2.7 Å) matches experiment and previous simulations, providing an external benchmark. The RF-MTACE method is a credible acceleration technique, and the paper is clearly written with well-structured analyses. If the central necessary-condition claims hold, the work would give a unified structural/dynamical rationale for electron localization in water across dimensionalities, with implications for radiation chemistry and confined water. However, the claims are stronger than the current evidence: each environment is represented by a single trajectory of 10–20 ps, and the key negative cases lack statistical replication and time-resolved analysis. The paper also does not address nuclear quantum effects, which the authors themselves cite as important in bulk water (Ref. 61).

major comments (3)
  1. [§3.3, Figs. 9 and 10] The conclusion that the 2D monolayer does not sustain a long-lived localized state rests on a single 20 ps trajectory and on a time-averaged projection of the spin density. The compact state at t = 11.77 ps (rg = 3.3 Å, Fig. 9c) shows that focused dangling OH bonds can produce a localized configuration, but the authors characterize this event as 'short-lived and rare' without quantifying its lifetime, its frequency of occurrence, or the uncertainty from a single initial condition. The time-averaged projection in Fig. 10a can appear spread even if the electron repeatedly localizes at different lateral sites, because spatial and temporal averaging cancels localized features. A dwell-time analysis of the electron center of mass and multiple independent trajectories are needed before the 'necessary requirement' claim in Section 4 can be considered established for this system.
  2. [§2 and §4] All simulations treat hydrogen nuclei classically. The authors cite Ref. 61 (Lan et al., Nat. Commun. 2021) as showing a twin-cavity structure in bulk water arising from nuclear quantum effects. Since the paper's central mechanistic arguments concern the dynamics of dangling OH bonds and the stability of the second solvation shell — both properties that involve proton positions — the absence of NQE could change the reported outcomes. The manuscript does not discuss this limitation or provide a sensitivity test (e.g., a path-integral or centroid-molecular-dynamics run for at least one representative system). Please add such a discussion or an explicit test.
  3. [§4] Section 4 states that focused dangling OH bonds are a 'necessary requirement' for electron localization and that a stable second solvation shell is a 'necessary condition' for dynamic stabilization. These universal claims are inferred from four environments, each with one trajectory. The data are consistent with the claims, but the logical form of a necessity claim requires ruling out localization mechanisms in systems where the proposed features are absent. The current set of systems provides only three negative cases, and the negative cases are under-sampled (see the first comment). The authors should either soften the claims to be explicitly about the systems studied, or provide additional evidence, such as a free-energy or mechanistic argument that excludes alternative localization pathways.
minor comments (6)
  1. [§2] The multiple-time-step parameters read 'δt = 0.48 fs and ∆t = 48 ps; i.e., nMTS = 100.' With nMTS = 100 and δt = 0.48 fs, Δt = 48 fs, not 48 ps. Please correct the unit.
  2. [§3.1, Refs. 11–12] The text attributes an experimental rg of 2.48 ± 0.1 Å from moment analysis of the absorption spectrum to Ref. 11 (Kevan, electron spin echo) and a transient terahertz value of 2.7 Å to Ref. 12 (Coe et al., photoelectron spectra). These citations appear mismatched; Ref. 7 (Novelli et al.) is the terahertz study. Please check and correct.
  3. [Fig. 12] The notation '5.6 ± 0.69 Å' and '6.7 ± 0.18 Å' reports the standard deviation over the single trajectory, not a statistical error of the mean. Please clarify this in the caption or text and add an estimate of the standard error from block averaging.
  4. [SI, S4] The description of the 1D chain asphericity — 'spherical charge droplet on every dangling OH group across the chain' — appears inconsistent with the main-text description of a cylindrically symmetric, core-depleted spin density (Section 3.3). Please clarify or correct the wording.
  5. [Introduction] The introduction states that Hart and Boag observed the hydrated electron, but the cited reference (Ref. 1, Weyl, 1864) is about metal-ammonia solutions. Please add the proper Hart and Boag reference.
  6. [Fig. 4] The H-bond distributions are compared across systems with different sizes and geometries, but the H-bond definition is not given in the main text. Please specify the geometric/energetic criterion used.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity found: the paper's causal claims are inductive generalizations from comparative hybrid-DFT simulations, benchmarked against external experiment and independent prior work; the only self-citations are for the RF-MTACE computational method and are not load-bearing.

full rationale

The paper's derivation chain is self-contained. The localization measure rg (Eqs. 1-2 of Section 2) is a direct descriptor of the computed spin density, and the dangling-OH orientations, H-bond counts, and second-shell structure are independent geometric observables computed from the same trajectories; the Section 4 conclusions (dangling OH bonds all pointing to a small region are necessary for localization; a stable second solvation shell is necessary for dynamic stabilization) are inductive generalizations across the five simulated systems rather than identities between inputs and outputs. No parameter is fitted to the results being explained: the PBE0 functional with 40% Fock exchange is adopted from refs 56, 58, 59 (Pasquarello and Schwartz groups), the SIN(R) thermostat parameters from refs 75-77, the ice structure from ref 82, and the Lennard-Jones parameters for the hydrophobic confinements from ref 85; all are external to this paper. The bulk comparison rg = 2.7 Å versus experimental 2.48 ± 0.1 Å (ref 11) and 2.7 Å (ref 12) is an external benchmark used for validation, not an input to any fit. The self-citations that do exist (refs 68-71, all including the present authors) support only the computational acceleration claim for RF-MTACE; the physical conclusions are cross-checked against independent hybrid-functional AIMD studies by other groups (refs 51, 56-59, 61-62), so the central theses do not reduce to the authors' own method. The legitimate weaknesses here are robustness issues, not circularity: the 2D monolayer 'no stable localization' verdict rests on a single 20 ps trajectory and time-averaged spin-density projections (Figure 10a) that can average over transient events, as the paper itself reports a compact rg = 3.3 Å state at t = 11.77 ps (Figure 9), and the second-solvation-shell 'necessary condition' is tested only in geometries that preclude a second shell by construction, leaving confounds such as dimensionality and confinement uncontrolled. These are sampling and inference-strength limitations, which are correctness risks rather than circular steps.

Assumptions & free parameters 4 free parameters · 5 assumptions · 0 invented entities

The paper does not fit any parameter to its own results; the free parameters listed are imported modeling choices. The central claim rests on the transferability of PBE0-40 to non-bulk environments, on classical nuclei, on the representativeness of single short trajectories, and on the fidelity of the fixed Lennard-Jones confining walls. No new physical entities are postulated. The largest hidden cost is the functional and thermostat parameter choice, which is inherited from bulk-water benchmarks.

free parameters (4)
  • PBE0 exact exchange fraction = 0.40
    Chosen on the basis of prior benchmark studies of the hydrated electron (refs 56, 57); not fitted to the data reported here, but the central conclusions depend on this functional choice.
  • Graphene/CNT-dummy Lennard-Jones parameters = From Deshmukh et al. 2013 (ref 85)
    Define the hydrophobic confinement used for the 2D monolayer and 1D chain; not fitted to electron data, but they set the environment whose dynamics drive the delocalization result.
  • RF-MTACE outer timestep = delta t=0.48 fs, nMTS=100; Delta t=48 fs (typo in text: 48 ps)
    Chosen for efficiency and stability in RF-MTACE; the reported 48 ps appears to be a typo for 48 fs. The timestep influences the sampled dynamics.
  • SIN(R) thermostat parameters = L=4, tau=9.7 fs, gamma=0.01 fs^-1
    Stochastic resonance-free thermostat settings; standard but hand-chosen, and they affect dynamical properties such as localization stability.
assumptions (5)
  • domain assumption PBE0 with 40% exact exchange accurately describes the hydrated electron in all environments studied
    Adopted from prior benchmark studies in bulk water (refs 56, 57); transferability to ice, monolayer, and chain is assumed, not validated.
  • domain assumption Classical treatment of nuclei is adequate for electron localization dynamics
    Born-Oppenheimer MD with classical thermostat (Section 2); nuclear quantum effects are known to alter bulk hydrated-electron structure (ref 61).
  • domain assumption Spin-density-based radius of gyration and center of mass are robust descriptors even for delocalized periodic systems
    Used in Eqs. 1-2; for a delocalized electron spanning the simulation cell, periodic images can bias rg. The chain rg = 6.72 Å is close to half the chain length, so this assumption is load-bearing.
  • domain assumption One 20 ps trajectory per system is sufficient to characterize localization and stability
    Single trajectory, no replicate initial conditions (Section 2); slow reorganization could be missed in rigid ice.
  • domain assumption Fixed dummy atoms interacting via Lennard-Jones potentials mimic hydrophobic confinement
    Used for monolayer and chain (Section 2); neglects electronic polarization and flexibility of the confining walls.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Nature of Hydrated Electron in Varied Solvation Environments." pith.science (2026). https://pith.science/paper/GU4OWWXG

@misc{pith2026250607157,
  author       = {Pith},
  title        = {Pith review of: Nature of Hydrated Electron in Varied Solvation Environments},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/GU4OWWXG}},
  note         = {Machine review of arXiv:2506.07157}
}
read the original abstract

Understanding the nature of solvated electrons is important in studying a range of chemical and biological phenomena. This study investigates the structural and dynamical behavior of an excess electron in water, examining different solvation environments, including liquid water, ice, monolayer, and chain. To accurately model these systems, we carry out molecular dynamics (MD) simulations using hybrid density functionals, employing the computationally efficient resonance-free multiple time-stepping based adaptively compressed exchange operator method. Through these simulations, we create a comprehensive and detailed picture of how excess electrons are solvated across different aqueous environments. We report the factors influence the localization and dynamic stability of the hydrated electron. The determinants include the presence and reorganization flexibility of the dangling OH groups and the spatial arrangement of the surrounding water molecules.

Figures

Figures reproduced from arXiv: 2506.07157 by the authors.

Figure 1
Figure 1. Time series plot of the rg of the excess e − in liquid water. Red dots indicate the specific frames for which the corresponding spin density is shown. The histogram is shown separately for the first 2 ps (left panel; red) and after localization (right panel; green). The black dotted line indicates the average rg computed after 2 ps. Spin density plots for selected time frames are shown with an isovalue of 0.001 au. … view at source ↗
Figure 2
Figure 2. RDF of oxygen (green) and hydrogen (red) with respect to COM of solvated [PITH_FULL_IMAGE:figures/full_fig_p009_2.png] view at source ↗
Figure 3
Figure 3. Hydrogen atom coordination around COM of solvated electron. Spin densities for [PITH_FULL_IMAGE:figures/full_fig_p011_3.png] view at source ↗
Figures from the paper (9 more)
Figure 4
Figure 4. Figure 4: Average H-bond distribution per water molecule for different systems. [PITH_FULL_IMAGE:figures/full_fig_p012_4.png]
Figure 5
Figure 5. Figure 5: Average H-bond distribution per water molecule along the distance from the COM [PITH_FULL_IMAGE:figures/full_fig_p013_5.png]
Figure 6
Figure 6. Figure 6: (a) Average angular distribution of COM of the electron to O-H vectors of the first [PITH_FULL_IMAGE:figures/full_fig_p014_6.png]
Figure 7
Figure 7. Figure 7: Time series plot of the rg of the excess e − in ice, with and without defect. the rigidity of the lattice, which restricts reorganization of hydrogen bond network, thereby inhibiting the cavity formation. Projected spin density analysis in Figure 8a-b showed bright reg…
Figure 8
Figure 8. Figure 8: (a) Snapshot from MD simulation of perfect hexagonal ice with an excess electron, [PITH_FULL_IMAGE:figures/full_fig_p016_8.png]
Figure 9
Figure 9. Figure 9: Time evolution of rg of spin density of excess e − for the 2D monolayer and 1D chain water molecules. Spin densities are shown for both cases as mesh with isovalue 0.0001 au for selected frames. (a)-(c) For 2D monolayer: dangling O-H bonds interacting with e − are visu…
Figure 10
Figure 10. Figure 10: Two-dimensional projection of normalized spin density on [PITH_FULL_IMAGE:figures/full_fig_p020_10.png]
Figure 11
Figure 11. Figure 11: Average distribution of OH vector and dipole angle for (a) 2D system with [PITH_FULL_IMAGE:figures/full_fig_p021_11.png]
Figure 12
Figure 12. Figure 12: Distribution of rg for 3D (liquid, ice; with and without defect), 2D monolayer and 1D water chain systems. The most probable values of rg are mentioned in the plot. Systems containing monolayer and chain of confined water molecules have under-coordinated 23 [PITH_FUL…

Discussion (0). Sign in to comment.

Reference graph

Works this paper leans on

87 extracted references · 79 canonical work pages

  1. [1]

    Ueber Metallammonium-Verbindungen

    Weyl, W. Ueber Metallammonium-Verbindungen. Annalen der Physik 1864, 197, 601--612

  2. [2]

    Analytical investigations of an electron–water molecule pseudopotential

    Turi, L.; Borgis, D. Analytical investigations of an electron–water molecule pseudopotential. II. Development of a new pair potential and molecular dynamics simulations . J. Chem. Phys. 2002, 117, 6186--6195

    2002

  3. [3]

    Turi, L.; Sheu, W.-S.; Rossky, P. J. Characterization of excess electrons in water-cluster anions by quantum simulations. Science 2005, 309, 914--917

    2005

  4. [4]

    E.; Glover, W

    Larsen, R. E.; Glover, W. J.; Schwartz, B. J. Does the hydrated electron occupy a cavity? Science 2010, 329, 65--69

    2010

  5. [5]

    W.; Karashima, S.; Suzuki, T

    Nishitani, J.; Yamamoto, Y.-i.; West, C. W.; Karashima, S.; Suzuki, T. Binding energy of solvated electrons and retrieval of true UV photoelectron spectra of liquids. Sci. Adv. 2019, 5, eaaw6896

    2019

  6. [6]

    Herbert, J. M. Structure of the aqueous electron. Phys. Chem. Chem. Phys. 2019, 21, 20538--20565

    2019

  7. [7]

    The birth and evolution of solvated electrons in the water

    Novelli, F.; Chen, K.; Buchmann, A.; Ockelmann, T.; Hoberg, C.; Head-Gordon, T.; Havenith, M. The birth and evolution of solvated electrons in the water. Proc. Natl. Acad. Sci. 2023, 120, e2216480120

    2023

  8. [8]

    Schnitker, J.; Rossky, P. J. Quantum simulation study of the hydrated electron . J. Chem. Phys. 1987, 86, 3471--3485

    work page 1987

Show all 87 references
  1. [9]

    C.; Stienkemeier, F.; Slav \' c ek, P.; W \"o rner, H

    Svoboda, V.; Michiels, R.; LaForge, A. C.; Stienkemeier, F.; Slav \' c ek, P.; W \"o rner, H. J. Real-time observation of water radiolysis and hydrated electron formation induced by extreme-ultraviolet pulses. Sci. Adv. 2020, 6, eaaz0385

    2020

  2. [10]

    E.; Ban, L.; Yoder, B

    Gartmann, T. E.; Ban, L.; Yoder, B. L.; Hartweg, S.; Chasovskikh, E.; Signorell, R. Relaxation dynamics and genuine properties of the solvated electron in neutral water clusters. J. Phys. Chem. Lett. 2019, 10, 4777--4782

    2019

  3. [11]

    Electron spin echo studies of solvation structure

    Kevan, L. Electron spin echo studies of solvation structure. J. Phys. Chem. 1981, 85, 1628--1636

    1981

  4. [12]

    V.; Williams, S

    Coe, J. V.; Williams, S. M.; and, K. H. B. Photoelectron spectra of hydrated electron clusters vs. cluster size: connecting to bulk. Int. Rev. Phys. Chem. 2008, 27, 27--51

    2008

  5. [13]

    Hart, E. J. Research potentials of the hydrated electron. Acc. Chem. Res. 1969, 2, 161--167

    1969

  6. [14]

    Rybkin, V. V. Mechanism of aqueous carbon dioxide reduction by the solvated electron. J. Phys. Chem. B 2020, 124, 10435--10441

    2020

  7. [15]

    C.; Dixon, D

    Garrett, B. C.; Dixon, D. A.; Camaioni, D. M.; Chipman, D. M.; Johnson, M. A.; Jonah, C. D.; Kimmel, G. A.; Miller, J. H.; Rescigno, T. N.; Rossky, P. J., et al. Role of water in electron-initiated processes and radical chemistry: Issues and scientific advances. Chem. Rev. 200...

    2004

  8. [16]

    Kumar, A.; Becker, D.; Adhikary, A.; Sevilla, M. D. Reaction of electrons with DNA : radiation damage to radiosensitization. International Journal of Molecular Sciences 2019, 20, 3998

    2019

  9. [17]

    Schnitker, J.; Rossky, P. J. An electron–water pseudopotential for condensed phase simulation. The Journal of Chemical Physics 1987, 86, 3462--3470

    1987

  10. [18]

    D.; Williams, C

    Jacobson, L. D.; Williams, C. F.; Herbert, J. M. The static-exchange electron-water pseudopotential, in conjunction with a polarizable water model: A new Hamiltonian for hydrated-electron simulations. J. Chem. Phys. 2009, 130, 124115

    2009

  11. [19]

    Does the Hydrated Electron Occupy a Cavity?

    Turi, L.; Ádám Madarász, Comment on “Does the Hydrated Electron Occupy a Cavity?”. Science 2011, 331, 1387--1387

    2011

  12. [20]

    Does the Hydrated Electron Occupy a Cavity?

    Jacobson, L. D.; Herbert, J. M. Comment on “Does the Hydrated Electron Occupy a Cavity?”. Science 2011, 331, 1387--1387

    2011

  13. [21]

    Does the Hydrated Electron Occupy a Cavity?

    Larsen, R. E.; Glover, W. J.; Schwartz, B. J. Response to Comments on “Does the Hydrated Electron Occupy a Cavity?”. Science 2011, 331, 1387--1387

    2011

  14. [22]

    R.; Kahros, A.; Schwartz, B

    Casey, J. R.; Kahros, A.; Schwartz, B. J. To be or not to be in a cavity: the hydrated electron dilemma. J. Phys. Chem. B 2013, 117, 14173--14182

    2013

  15. [23]

    D.; Romero, C.; Rahman, A

    Jonah, C. D.; Romero, C.; Rahman, A. Hydrated electron revisited via the Feynman path integral route. Chem. Phys. letters 1986, 123, 209--214

    1986

  16. [24]

    Schnitker, J.; Rossky, P. J. Quantum simulation study of the hydrated electron. J. Chem. Phys. 1987, 86, 3471--3485

    1987

  17. [25]

    N.; Landman, U.; Cleveland, C

    Barnett, R. N.; Landman, U.; Cleveland, C. L.; Jortner, J. Surface states of excess electrons on water clusters. Phys. Rev. Lett. 1987, 59, 811--814

    1987

  18. [26]

    On the calculation of dynamical properties of solvated electrons by maximum entropy analytic continuation of path integral Monte Carlo data

    Gallicchio, E.; Berne, B. On the calculation of dynamical properties of solvated electrons by maximum entropy analytic continuation of path integral Monte Carlo data. J. Chem. Phys. 1996, 105, 7064--7078

    1996

  19. [27]

    Sprik, M.; Klein, M. L. Optimization of a distributed Gaussian basis set using simulated annealing: application to the adiabatic dynamics of the solvated electron. J. Chem. Phys. 1988, 89, 1592--1607

    1988

  20. [28]

    J.; Schnitker, J

    Rossky, P. J.; Schnitker, J. The hydrated electron: quantum simulation of structure, spectroscopy, and dynamics. J. Phys. Chem. 1988, 92, 4277--4285

    1988

  21. [29]

    D.; Herbert, J

    Jacobson, L. D.; Herbert, J. M. A one-electron model for the aqueous electron that includes many-body electron-water polarization: Bulk equilibrium structure, vertical electron binding energy, and optical absorption spectrum. J. Chem. Phys. 2010, 133

    2010

  22. [30]

    J.; Schwartz, B

    Glover, W. J.; Schwartz, B. J. Short-range electron correlation stabilizes noncavity solvation of the hydrated electron. J. Chem. Theory Comput. 2016, 12, 5117--5131

    2016

  23. [31]

    B.; Kim, K

    Kim, J.; Suh, S. B.; Kim, K. S. Water dimer to pentamer with an excess electron: Ab initio study. The Journal of Chemical Physics 1999, 111, 10077--10087

    1999

  24. [32]

    I.; Shin, J.-W.; Headrick, J

    Hammer, N. I.; Shin, J.-W.; Headrick, J. M.; Diken, E. G.; Roscioli, J. R.; Weddle, G. H.; Johnson, M. A. How Do Small Water Clusters Bind an Excess Electron? Science 2004, 306, 675--679

    2004

  25. [33]

    M.; Suh, S

    Lee, H. M.; Suh, S. B.; Tarakeshwar, P.; Kim, K. S. Origin of the magic numbers of water clusters with an excess electron. J. Chem. Phys. 2005, 122, 044309

    2005

  26. [34]

    J.; Turi, L

    Madar \'a sz, \'A .; Rossky, P. J.; Turi, L. Interior- and surface-bound excess electron states in large water cluster anions. J. Chem. Phys. 2009, 130, 124319

    2009

  27. [35]

    Impact of Electron Solvation on Ice Structures at the Molecular Scale

    Bertram, C.; Auburger, P.; Bockstedte, M.; Stähler, J.; Bovensiepen, U.; Morgenstern, K. Impact of Electron Solvation on Ice Structures at the Molecular Scale. J. Phys. Chem. Lett. 2020, 11, 1310--1316

    2020

  28. [36]

    K.; Inam, F.; Scandolo, S

    Bhattacharya, S. K.; Inam, F.; Scandolo, S. Excess electrons in ice: a density functional theory study. Phys. Chem. Chem. Phys. 2014, 16, 3103--3107

    2014

  29. [37]

    de Koning, M.; Fazzio, A.; da Silva, A. J. R.; Antonelli, A. On the nature of the solvated electron in ice I h. Phys. Chem. Chem. Phys. 2016, 18, 4652--4658

    2016

  30. [38]

    Surface Trapped Excess Electrons on Ice

    Baletto, F.; Cavazzoni, C.; Scandolo, S. Surface Trapped Excess Electrons on Ice. Phys. Rev. Lett. 2005, 95, 176801

    2005

  31. [39]

    Crystal imperfections in ice I h

    de Koning, M. Crystal imperfections in ice I h. J. Chem. Phys. 2020, 153, 110902

    2020

  32. [40]

    Ab Initio Molecular Dynamics: Basic Theory and Advanced Methods; Cambridge University Press: Cambridge, 2009

    Marx, D.; Hutter, J. Ab Initio Molecular Dynamics: Basic Theory and Advanced Methods; Cambridge University Press: Cambridge, 2009

    2009

  33. [41]

    Adiabatic connection approach to density functional theory of electronic systems

    Savin, A.; Colonna, F.; Pollet, R. Adiabatic connection approach to density functional theory of electronic systems. Int. J. Quantum Chem. 2003, 93, 166--190

    2003

  34. [42]

    Boero, M.; Parrinello, M.; Terakura, K.; Ikeshoji, T.; Liew, C. C. First-Principles Molecular-Dynamics Simulations of a Hydrated Electron in Normal and Supercritical Water. Phys. Rev. Lett. 2003, 90, 226403

    2003

  35. [43]

    Excess Electron in Water at Different Thermodynamic Conditions

    Boero, M. Excess Electron in Water at Different Thermodynamic Conditions. J. Phys. Chem. A 2007, 111, 12248--12256

    2007

  36. [44]

    Lee, C.; Yang, W.; Parr, R. G. Development of the C olle- S alvetti correlation-energy formula into a functional of the electron density. Phys. Rev. B 1988, 37, 785--789

    1988

  37. [45]

    M.; Head-Gordon, M

    Herbert, J. M.; Head-Gordon, M. First-principles, quantum-mechanical simulations of electron solvation by a water cluster. Proc. Natl. Acad. Sci. 2006, 103, 14282--14287

    2006

  38. [46]

    Ab Initio Molecular Dynamics Simulation of a Medium-Sized Water Cluster Anion: From an Interior to a Surface-Located Excess Electron via a Delocalized State

    Frigato, T.; VandeVondele, J.; Schmidt, B.; Schütte, C.; Jungwirth, P. Ab Initio Molecular Dynamics Simulation of a Medium-Sized Water Cluster Anion: From an Interior to a Surface-Located Excess Electron via a Delocalized State. J. Phys. Chem. A 2008, 112, 6125--6133

    2008

  39. [47]

    Dynamics of Electron Localization in Warm versus Cold Water Clusters

    Marsalek, O.; Uhlig, F.; Frigato, T.; Schmidt, B.; Jungwirth, P. Dynamics of Electron Localization in Warm versus Cold Water Clusters. Phys. Rev. Lett. 2010, 105, 043002

    2010

  40. [48]

    E.; Schmidt, B.; Schütte, C.; Jungwirth, P

    Marsalek, O.; Frigato, T.; VandeVondele, J.; Bradforth, S. E.; Schmidt, B.; Schütte, C.; Jungwirth, P. Hydrogen Forms in Water by Proton Transfer to a Distorted Electron. J. Phys. Chem. B 2010, 114, 915--920

    2010

  41. [49]

    Electrons in Cold Water Clusters: An ab Initio Molecular Dynamics Study of Localization and Metastable States

    Marsalek, O.; Uhlig, F.; Jungwirth, P. Electrons in Cold Water Clusters: An ab Initio Molecular Dynamics Study of Localization and Metastable States. J. Phys. Chem. C 2010, 114, 20489--20495

    2010

  42. [50]

    N.; Giniger, R.; Cheshnovsky, O.; Landman, U

    Barnett, R. N.; Giniger, R.; Cheshnovsky, O.; Landman, U. Dielectron Attachment and Hydrogen Evolution Reaction in Water Clusters. J. Phys. Chem. A 2011, 115, 7378--7391

    2011

  43. [51]

    Unraveling the Complex Nature of the Hydrated Electron

    Uhlig, F.; Marsalek, O.; Jungwirth, P. Unraveling the Complex Nature of the Hydrated Electron. J. Phys. Chem. Lett. 2012, 3, 3071--3075

    2012

  44. [52]

    C.; Rana, B.; Herbert, J

    Holden, Z. C.; Rana, B.; Herbert, J. M. Analytic gradient for the QM/MM -Ewald method using charges derived from the electrostatic potential: Theory, implementation, and application to ab initio molecular dynamics simulation of the aqueous electron . J. Chem. Phys. 2019, 150, 144115

    2019

  45. [53]

    Structure, Dynamics, and Reactivity of Hydrated Electrons by Ab Initio Molecular Dynamics

    Marsalek, O.; Uhlig, F.; VandeVondele, J.; Jungwirth, P. Structure, Dynamics, and Reactivity of Hydrated Electrons by Ab Initio Molecular Dynamics. Acc. Chem. Res. 2012, 45, 23--32

    2012

  46. [54]

    Direct observation of the collapse of the delocalized excess electron in water

    Savolainen, J.; Uhlig, F.; Ahmed, S.; Hamm, P.; Jungwirth, P. Direct observation of the collapse of the delocalized excess electron in water. Nat. Chem. 2014, 6, 697--701

    2014

  47. [55]

    J.; Mori-S \'a nchez, P.; Yang, W

    Cohen, A. J.; Mori-S \'a nchez, P.; Yang, W. Insights into Current Limitations of Density Functional Theory. Science 2008, 321, 792--794

    2008

  48. [56]

    Electronic Levels of Excess Electrons in Liquid Water

    Ambrosio, F.; Miceli, G.; Pasquarello, A. Electronic Levels of Excess Electrons in Liquid Water. J. Phys. Chem. Lett. 2017, 8, 2055--2059

    2017

  49. [57]

    R.; Liu, X.; Schwartz, B

    Borrelli, W. R.; Liu, X.; Schwartz, B. J. How the choice of exchange–correlation functional affects DFT-based simulations of the hydrated electron. J. Chem. Phys. 2025, 162, 110901

    2025

  50. [58]

    V.; Pasquarello, A

    Lan, J.; Rybkin, V. V.; Pasquarello, A. Temperature Dependent Properties of the Aqueous Electron. Angew. Chem. 2022, 61, e202209398

    2022

  51. [59]

    J.; Schwartz, B

    Park, S. J.; Schwartz, B. J. Understanding the Temperature Dependence and Finite Size Effects in Ab Initio MD Simulations of the Hydrated Electron. J. Chem. Theory Comput. 2022, 18, 4973--4982

    2022

  52. [60]

    Wilhelm, J.; VandeVondele, J.; Rybkin, V. V. Dynamics of the Bulk Hydrated Electron from Many-Body Wave-Function Theory. Angew. Chem. 2019, 58, 3890--3893

    2019

  53. [61]

    Lan, J.; Kapil, V.; Gasparotto, P.; Ceriotti, M.; Iannuzzi, M.; Rybkin, V. V. Simulating the ghost: quantum dynamics of the solvated electron. Nat. Commun. 2021, 12, 766

    2021

  54. [62]

    Enhanced deep potential model for fast and accurate molecular dynamics: application to the hydrated electron

    Gao, R.; Li, Y.; Car, R. Enhanced deep potential model for fast and accurate molecular dynamics: application to the hydrated electron. Phys. Chem. Chem. Phys. 2024, 26, 23080--23088

    2024

  55. [63]

    P.; You, Z.-Q.; Herbert, J

    Coons, M. P.; You, Z.-Q.; Herbert, J. M. The Hydrated Electron at the Surface of Neat Liquid Water Appears To Be Indistinguishable from the Bulk Species. J. Am. Chem. Soc. 2016, 138, 10879--10886

    2016

  56. [64]

    Tuckerman, M. E. Statistical Mechanics: Theory and Molecular Simulation; Oxford University Press: Oxford, 2010

    2010

  57. [65]

    Adaptively Compressed Exchange Operator

    Lin, L. Adaptively Compressed Exchange Operator. J. Chem. Theory Comput. 2016, 12, 2242--2249

    2016

  58. [66]

    S.; Vecharynski, E.; Yang, C

    Hu, W.; Lin, L.; Banerjee, A. S.; Vecharynski, E.; Yang, C. Adaptively Compressed Exchange Operator for Large-Scale Hybrid Density Functional Calculations with Applications to the Adsorption of Water on Silicene. J. Chem. Theory Comput. 2017, 13, 1188--1198

    2017

  59. [67]

    Low-rank approximations for accelerating plane-wave hybrid functional calculations in unrestricted and noncollinear spin density functional theory

    Chen, S.; Wu, K.; Hu, W.; Yang, J. Low-rank approximations for accelerating plane-wave hybrid functional calculations in unrestricted and noncollinear spin density functional theory. J. Chem. Phys. 2023, 158, 134106

    2023

  60. [68]

    Mandal, S.; Nair, N. N. Speeding-up ab initio molecular dynamics with hybrid functionals using adaptively compressed exchange operator based multiple timestepping. J. Chem. Phys. 2019, 151, 151102

    2019

  61. [69]

    Mandal, S.; Kar, R.; Klöffel, T.; Meyer, B.; Nair, N. N. Improving the scaling and performance of multiple time stepping-based molecular dynamics with hybrid density functionals. J. Comput. Chem. 2022, 43, 588--597

    2022

  62. [70]

    Mandal, S.; Thakkur, V.; Nair, N. N. Achieving an Order of Magnitude Speedup in Hybrid-Functional- and Plane-Wave-Based Ab Initio Molecular Dynamics: Applications to Proton-Transfer Reactions in Enzymes and in Solution. J. Chem. Theory Comput. 2021, 17, 2244--2255

    2021

  63. [71]

    Kar, R.; Mandal, S.; Thakkur, V.; Meyer, B.; Nair, N. N. Speeding-up Hybrid Functional-Based Ab Initio Molecular Dynamics Using Multiple Time-stepping and Resonance-Free Thermostat. J. Chem. Theory Comput. 2023, 19, 8351--8364

    2023

  64. [72]

    Copyright 2000-2021 jointly by IBM Corp

    CPMD , version 4.3. Copyright 2000-2021 jointly by IBM Corp. and by Max Planck Institute, Stuttgart., 2021; http://www.cpmd.org (accessed on Oct 18, 2021)

    2000

  65. [73]

    Kl \"o ffel, T.; Mathias, G.; Meyer, B. Integrating state of the art compute, communication, and autotuning strategies to multiply the performance of ab initio molecular dynamics on massively parallel multi-core supercomputers. Comput. Phys. Commun. 2021, 260, 107745

    2021

  66. [74]

    E.; Martyna, G

    Tuckerman, M. E.; Martyna, G. J.; Berne, B. J. Molecular dynamics algorithm for condensed systems with multiple time scales. J. Chem. Phys. 1990, 93, 1287--1291

    1990

  67. [75]

    T.; Tuckerman, M

    Leimkuhler, B.; Margul, D. T.; Tuckerman, M. E. Stochastic, resonance-free multiple time-step algorithm for molecular dynamics with very large time steps. Mol. Phys. 2013, 111, 3579--3594

    2013

  68. [76]

    R.; Tuckerman, M

    Abreu, C. R.; Tuckerman, M. E. Molecular Dynamics with Very Large Time Steps for the Calculation of Solvation Free Energies. J. Chem. Theory Comput. 2020, 16, 7314--7327

    2020

  69. [77]

    Abreu, C. R. A.; Tuckerman, M. E. Hamiltonian based resonance-free approach for enabling very large time steps in multiple time-scale molecular dynamics. Mol. Phys. 2021, 119, e1923848

    2021

  70. [78]

    Troullier, N.; Martins, J. L. Efficient pseudopotentials for plane-wave calculations. Phys. Rev. B 1991, 43, 1993--2006

    1991

  71. [79]

    Convergence acceleration of iterative sequences

    Pulay, P. Convergence acceleration of iterative sequences. the case of scf iteration. Chem. Phys. Lett. 1980, 73, 393 -- 398

    1980

  72. [80]

    P.; Parrinello, M

    Hutter, J.; Lüthi, H. P.; Parrinello, M. Electronic structure optimization in plane-wave-based density functional calculations by direct inversion in the iterative subspace. Comput. Mater. Sci. 1994, 2, 244 -- 248

    1994

  73. [81]

    Time-reversible always stable predictor–corrector method for molecular dynamics of polarizable molecules

    Kolafa, J. Time-reversible always stable predictor–corrector method for molecular dynamics of polarizable molecules. J. Comput. Chem. 2004, 25, 335--342

    2004

  74. [82]

    Band gaps of liquid water and hexagonal ice through advanced electronic-structure calculations

    Bischoff, T.; Reshetnyak, I.; Pasquarello, A. Band gaps of liquid water and hexagonal ice through advanced electronic-structure calculations. Phys. Rev. Res. 2021, 3, 023182

    2021

  75. [83]

    D.; Fowler, R

    Bernal, J. D.; Fowler, R. H. A Theory of Water and Ionic Solution, with Particular Reference to Hydrogen and Hydroxyl Ions . J. Chem. Phys. 1933, 1, 515--548

    1933

  76. [84]

    C.; Schwegler, E.; Gygi, F.; Galli, G

    Cicero, G.; Grossman, J. C.; Schwegler, E.; Gygi, F.; Galli, G. Water Confined in Nanotubes and between Graphene Sheets: A First Principle Study. J. Am. Chem. Soc. 2008, 130, 1871--1878

    2008

  77. [85]

    A.; Kamath, G.; Baker, G

    Deshmukh, S. A.; Kamath, G.; Baker, G. A.; Sumant, A. V.; Sankaranarayanan, S. K. The interfacial dynamics of water sandwiched between graphene sheets are governed by the slit width. Surface Science 2013, 609, 129--139

    2013

  78. [86]

    Picture of the wet electron: a localized transient state in liquid water

    Pizzochero, M.; Ambrosio, F.; Pasquarello, A. Picture of the wet electron: a localized transient state in liquid water. Chem. Sci. 2019, 10, 7442--7448

    2019

  79. [87]

    Eu v,! inS <v1ӢBk˜oՖ # [3

    Tsurusawa, T.; Iwata, S. Dipole-bound and interior electrons in water dimer and trimer anions: ab initio MO studies. Chem. Phys. Lett. 1998, 287, 553--562 mcitethebibliography achemso-demo.tex0000664000000000000000000023126315021314064012644 0ustar rootroot [journal=jctcce,man...

    1998

Pith tools

Reviewed August 7, 2026 · model on record in the stance chip above.