Ondrej Marsalek

Unraveling the Complex Nature of the Hydrated Electron.

The structure of the hydrated electron, which is a key species in radiative processes in water, has remained elusive. The traditional cavity model has been questioned recently, but the newly suggested picture of an electron delocalized over a region of enhanced water density is controversial. Here, we present results from ab initio molecular dynamics simulations, where not only the excess electron but also the valence electrons of the surrounding water molecules are described quantum mechanically. Unlike in previous one-electron pseudopotential calculations, many-electron interactions are explicitly accounted for. The present approach allows for partitioning of the electron solvated in liquid water into contributions from an inner cavity, neighboring water molecules, and a diffuse tail. We demonstrate that all three of these contributions are sizable and, consequently, important, which underlines the complex nature of the hydrated electron and warns against oversimplified interpretations based on pseudopotential models.

Chasing charge localization and chemical reactivity following photoionization in liquid water

The ultrafast dynamics of the cationic hole formed in bulk liquid water following ionization is investigated by ab initio molecular dynamics simulations and an experimentally accessible signature is suggested that might be tracked by femtosecond pump-probe spectroscopy. This is one of the fastest fundamental processes occurring in radiation-induced chemistry in aqueous systems and biological tissue. However, unlike the excess electron formed in the same process, the nature and time evolution of the cationic hole has been hitherto little studied. Simulations show that an initially partially delocalized cationic hole localizes within ~30 fs after which proton transfer to a neighboring water molecule proceeds practically immediately, leading to the formation of the OH radical and the hydronium cation in a reaction which can be formally written as H(2)O(+) + H(2)O → OH + H(3)O(+). The exact amount of initial spin delocalization is, however, somewhat method dependent, being realistically described by approximate density functional theory methods corrected for the self-interaction error. Localization, and then the evolving separation of spin and charge, changes the electronic structure of the radical center. This is manifested in the spectrum of electronic excitations which is calculated for the ensemble of ab initio molecular dynamics trajectories using a quantum mechanics/molecular mechanics (QM∕MM) formalism applying the equation of motion coupled-clusters method to the radical core. A clear spectroscopic signature is predicted by the theoretical model: as the hole transforms into a hydroxyl radical, a transient electronic absorption in the visible shifts to the blue, growing toward the near ultraviolet. Experimental evidence for this primary radiation-induced process is sought using femtosecond photoionization of liquid water excited with two photons at 11 eV. Transient absorption measurements carried out with ~40 fs time resolution and broadband spectral probing across the near-UV and visible are presented and direct comparisons with the theoretical simulations are made. Within the sensitivity and time resolution of the current measurement, a matching spectral signature is not detected. This result is used to place an upper limit on the absorption strength and/or lifetime of the localized H(2)O(+) ((aq)) species.

Solvation and ion-pairing properties of the aqueous sulfate anion: explicit versus effective electronic polarization

We assessed the relative merits of two approaches for including polarization effects in classical force fields for the sulfate anion. One of the approaches is the explicit shell model for atomic polarization and the other is an implicit dielectric continuum representation of the electronic polarization, wherein the polarizability density is spatially uniform. Both the solvation and ion association properties of sulfate were considered. We carried out an ab initio molecular dynamics simulation for a single sulfate anion in aqueous solution to obtain a benchmark for the solvation structure. For the ion-pairing properties, the models were compared to experimental thermodynamic data through Kirkwood-Buff theory, which relates the integrals of the pair correlation functions to measurable properties. While deficiencies were found for both of the approaches, the continuum polarization model was not systematically worse than the shell model. The shell model was found to give a more structured solution than the continuum polarization model, both with respect to solvation and ion pairing.

Ab initio molecular dynamics with nuclear quantum effects at classical cost: Ring polymer contraction for density functional theory.

Path integral molecular dynamics simulations, combined with an ab initio evaluation of interactions using electronic structure theory, incorporate the quantum mechanical nature of both the electrons and nuclei, which are essential to accurately describe systems containing light nuclei. However, path integral simulations have traditionally required a computational cost around two orders of magnitude greater than treating the nuclei classically, making them prohibitively costly for most applications. Here we show that the cost of path integral simulations can be dramatically reduced by extending our ring polymer contraction approach to ab initio molecular dynamics simulations. By using density functional tight binding as a reference system, we show that our ring polymer contraction scheme gives rapid and systematic convergence to the full path integral density functional theory result. We demonstrate the efficiency of this approach in ab initio simulations of liquid water and the reactive protonated and deprotonated water dimer systems. We find that the vast majority of the nuclear quantum effects are accurately captured using contraction to just the ring polymer centroid, which requires the same number of density functional theory calculations as a classical simulation. Combined with a multiple time step scheme using the same reference system, which allows the time step to be increased, this approach is as fast as a typical classical ab initio molecular dynamics simulation and 35× faster than a full path integral calculation, while still exactly including the quantum sampling of nuclei. This development thus offers a route to routinely include nuclear quantum effects in ab initio molecular dynamics simulations at negligible computational cost.

Electron at the Surface of Water: Dehydrated or Not?

The hydrated electron is a crucial species in radiative processes, and it has been speculated that its behavior at the water surface could lead to specific interfacial chemical properties. Here, we address fundamental questions concerning the structure and energetics of an electron at the surface of water. We use the method of ab initio molecular dynamics, which was shown to provide a faithful description of solvated electrons in large water clusters and in bulk water. The present results clearly demonstrate that the surface electron is mostly buried in the interfacial water layer, with only about 10 % of its density protruding into the vapor phase. Consequently, it has a structure that is very similar to that of an electron solvated in the aqueous bulk. This points to a general feature of charges at the surface of water, namely, that they do not behave as half-dehydrated but rather as almost fully hydrated species.

Efficient Calculation of Free Energy Differences Associated with Isotopic Substitution Using Path-Integral Molecular Dynamics.

The problem of computing free energy differences due to isotopic substitution in chemical systems is discussed. The shift in the equilibrium properties of a system upon isotopic substitution is a purely quantum mechanical effect that can be quantified using the Feynman path integral approach. In this paper, we explore two developments that lead to a highly efficient path integral scheme. First, we employ a mass switching function inspired by the work of Ceriotti and Markland [ J. Chem. Phys. 2013, 138, 014112] that is based on the inverse square root of the mass and which leads to a perfectly constant free energy derivative with respect to the switching parameter in the harmonic limit. We show that even for anharmonic systems, this scheme allows a single-point thermodynamic integration approach to be used in the construction of free energy differences. In order to improve the efficiency of the calculations even further, however, we derive a set of free energy derivative estimators based on the fourth-order scheme of Takahashi and Imada [ J. Phys. Soc. Jpn. 1984, 53, 3765]. The Takahashi-Imada procedure generates a primitive fourth-order estimator that allows the number of imaginary time slices in the path-integral approach to be reduced substantially. However, as with all primitive estimators, its convergence is plagued by numerical noise. In order to alleviate this problem, we derive a fourth-order virial estimator based on a transferring of the difference between second- and fourth-order primitive estimators, which remains relatively constant as a function of the number of configuration samples, to the second-order virial estimator. We show that this new estimator converges as smoothly as the second-order virial estimator but requires significantly fewer imaginary time points.

Hydrogen forms in water by proton transfer to a distorted electron.

Solvated electrons are ubiquitous intermediates in radiation-induced processes, with their lifetime being determined by quenching processes, such as the direct reaction with protons under acidic conditions. Ab initio molecular dynamics simulations allow us to unravel with molecular resolution the ultrafast reaction mechanism by which the electron and proton react in water. The path to a successful reaction involves a distortion and contraction of the hydrated electron and a rapid proton motion along a chain of hydrogen bonds, terminating on the water molecule most protruding into the electron cloud. This fundamental reaction is thus decidedly shown to be of a proton-transfer rather than electron-transfer character. Due to the desolvation penalty connected with breaking of the hydration shells of these charged particles, the reaction is, however, not diffusion-limited, in agreement with the interpretation of kinetics measurements.

Quantum Dynamics and Spectroscopy of Ab Initio Liquid Water: The Interplay of Nuclear and Electronic Quantum Effects

Understanding the reactivity and spectroscopy of aqueous solutions at the atomistic level is crucial for the elucidation and design of chemical processes. However, the simulation of these systems requires addressing the formidable challenges of treating the quantum nature of both the electrons and nuclei. Exploiting our recently developed methods that provide acceleration by up to 2 orders of magnitude, we combine path integral simulations with on-the-fly evaluation of the electronic structure at the hybrid density functional theory level to capture the interplay between nuclear quantum effects and the electronic surface. Here we show that this combination provides accurate structure and dynamics, including the full infrared and Raman spectra of liquid water. This allows us to demonstrate and explain the failings of lower-level density functionals for dynamics and vibrational spectroscopy when the nuclei are treated quantum mechanically. These insights thus provide a foundation for the reliable investigation of spectroscopy and reactivity in aqueous environments.

Arrhenius Analysis Of Anisotropic Surface Self-diffusion On The Prismatic Facet Of Ice

We present an Arrhenius analysis of self-diffusion on the prismatic surface of ice calculated from molecular dynamics simulations. The six-site water model of Nada and van der Eerden was used in combination with a structure-based criterion for determining the number of liquid-like molecules in the quasi-liquid layer. Simulated temperatures range from 230 K-287 K, the latter being just below the melting temperature of the model, 289 K. Calculated surface diffusion coefficients agree with available experimental data to within quoted precision. Our results indicate a positive Arrhenius curvature, implying a change in the mechanism of self-diffusion from low to high temperature, with a concomitant increase in energy of activation from 29.1 kJ mol(-1) at low temperature to 53.8 kJ mol(-1) close to the melting point. In addition, we find that the surface self-diffusion is anisotropic at lower temperatures, transitioning to isotropic in the temperature range of 240-250 K. We also present a framework for self-diffusion in the quasi-liquid layer on ice that aims to explain these observations.

Peptide salt bridge stability:From gas phase via microhydration to bulk water simulations

The salt bridge formation and stability in the terminated lysine-glutamate dipeptide is investigated in water clusters of increasing size up to the limit of bulk water. Proton transfer dynamics between the acidic and basic side chains is described by DFT-based Born-Oppenheimer molecular dynamics simulations. While the desolvated peptide prefers to be in its neutral state, already the addition of a single water molecule can trigger proton transfer from the glutamate side chain to the lysine side chain, leading to a zwitterionic salt bridge state. Upon adding more water molecules we find that stabilization of the zwitterionic state critically depends on the number of hydrogen bonds between side chain termini, the water molecules, and the peptidic backbone. Employing classical molecular dynamics simulations for larger clusters, we observed that the salt bridge is weakened upon additional hydration. Consequently, long-lived solvent shared ion pairs are observed for about 30 water molecules while solvent separated ion pairs are found when at least 40 or more water molecules hydrate the dipeptide. These results have implications for the formation and stability of salt bridges at partially dehydrated surfaces of aqueous proteins.