Atomistic simulation of water percolation and proton hopping in Nafion fuel cell membrane.

We have performed a detailed analysis of water clustering and percolation in hydrated Nafion configurations generated by classical molecular dynamics simulations. Our results show that at low hydration levels H(2)O molecules are isolated and a continuous hydrogen-bonded network forms as the hydration level is increased. Our quantitative analysis has established a hydration level (λ) between 5 and 6 H(2)O/SO(3)(-) as the percolation threshold of Nafion. We have also examined the effect of such a network on proton transport by studying the structural diffusion of protons using the quantum hopping molecular dynamics method. The mean residence time of the proton on a water molecule decreases by 2 orders of magnitude when the λ value is increased from 5 to 15. The proton diffusion coefficient in Nafion at a λ value of 15 is about 1.1 × 10(-5) cm(2)/s in agreement with experiment. The results provide quantitative atomic-level evidence of water network percolation in Nafion and its effect on proton conductivity.

Dynamics of electron localization in warm versus cold water clusters.

The process of electron localization on a cluster of 32 water molecules at 20, 50, and 300 K is unraveled using ab initio molecular dynamics simulations. In warm, liquid clusters, the excess electron relaxes from an initial diffuse and weakly bound structure to an equilibrated, strongly bound species within 1.5 ps. In contrast, in cold, glassy clusters the relaxation processes is not completed and the electron becomes trapped in a metastable surface state with an intermediate binding energy. These results question the validity of extrapolations of the properties of solvated electrons from cold clusters of increasing size to the liquid bulk.

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.

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.

We present a computational study of the structure and dynamics of an excess electron in a medium-sized water cluster aimed at addressing the question of interior vs exterior solvation. Ab initio Born-Oppenheimer molecular dynamics simulations were performed within the DFT framework, employing a hybrid Gaussian and plane-wave formalism together with the PBE exchange-correlation functional and norm-conserving pseudopotentials. Analysis of a 15-ps trajectory allowed us to reach the following conclusions: (i) the excess electron is predominantly located at the cluster surface (even if it is initially placed in the interior), (ii) the computed electron binding energies correlate with the electron localization rather than with its bulk vs surface location, and (iii) a dynamical interconversion between two different H-bond patterns around the electron occurs. The computed electron binding energies and the most relevant features of the IR spectrum are in a very good agreement with results of previous experimental studies.

Energy barriers of proton transfer reactions between amino acid side chain analogs and water from ab initio calculations.

Proton transfer reactions were studied in all titratable pairs of amino acid side chains where, under physiologically reasonable conditions, one amino acid may function as a donor and the other one as an acceptor. Energy barriers for shifting the proton from donor to acceptor atom were calculated by electronic structure methods at the MP2/6-31++G(d,p) level, and the well-known double-well potentials were characterized. The energy difference between both minima can be expressed by a parabola using as argument the donor-acceptor distance R(DA). In this work, the fit parameters of the quadratic expression are determined for each donor-acceptor pair. Moreover, it was found previously that the energy barriers of the reactions can be expressed by an analytical expression depending on the distance between donor and acceptor and the energy difference between donor and acceptor bound states. The validity of this approach is supported by the extensive new data set. This new parameterization of proton transfer barriers between titratable amino acid side chains allows us to very efficiently estimate proton transfer probabilities in molecular modelling studies or during classical molecular dynamics simulation of biomolecular systems.

Valence- and dipole-bound anions of the thymine-water complex: ab initio characterization of the potential energy surfaces.

The potential energy surfaces of the neutral and anionic thymine-water complexes are investigated using high-level ab initio calculations. Both dipole-bound (DB) and valence-bound (VB) anionic forms are considered. Four minima and three first-order stationary points are located, and binding energies are computed. All minima, for both anions, are found to be vertically and adiabatically stable. The binding energies are much higher for valence-bound than for dipole-bound anions. Adiabatic electron affinities are in the 66-287 meV range for VB anions and the 4-60 meV range for DB anions, and vertical detachment energies are in the 698-977 meV and 10-70 meV range for VB and DB anions, respectively. For cases where literature data are available, the computed values are in good agreement with previous experimental and theoretical studies. It is observed that electron attachment modifies the shape of the potential energy surfaces of the systems, especially for the valence-bound anions. Moreover, for both anions the size of the energy barrier between the two lowest energy minima is strongly reduced, rendering the coexistence of different structures more probable.

Unified molecular picture of the surfaces of aqueous acid, base, and salt solutions.

The molecular structure of the interfacial regions of aqueous electrolytes is poorly understood, despite its crucial importance in many biological, technological, and atmospheric processes. A long-term controversy pertains between the standard picture of an ion-free surface layer and the strongly ion specific behavior indicating in many cases significant propensities of simple inorganic ions for the interface. Here, we present a unified and consistent view of the structure of the air/solution interface of aqueous electrolytes containing monovalent inorganic ions. Molecular dynamics calculations show that in salt solutions and bases the positively charged ions, such as alkali cations, are repelled from the interface, whereas the anions, such as halides or hydroxide, exhibit a varying surface propensity, correlated primarily with the ion polarizability and size. The behavior of acids is different due to a significant propensity of hydronium cations for the air/solution interface. Therefore, both cations and anions exhibit enhanced concentrations at the surface and, consequently, these acids (unlike bases and salts) reduce the surface tension of water. The results of the simulations are supported by surface selective nonlinear vibrational spectroscopy, which reveals among other things that the hydronium cations are present at the air/solution interface. The ion specific propensities for the air/solution interface have important implications for a whole range of heterogeneous physical and chemical processes, including atmospheric chemistry of aerosols, corrosion processes, and bubble coalescence.

Ab initio electronic structure of thymine anions.

Dipole-bound and valence-bound anions of thymine were characterised by means of state-of-the-art ab initio calculations. The dipole-bound anion of thymine is both vertically and adiabatically stable with adiabatic electron affinity of 51 meV. The valence-bound anion is also adiabatically stable, albeit only marginally. Its vertical stability of 501 meV is, however, higher than that of the dipole-bound anion. All these values agree well with experimental findings.

The mechanism of proton exclusion in the aquaporin-1 water channel.

Aquaporins are efficient, yet strictly selective water channels. Remarkably, proton permeation is fully blocked, in contrast to most other water-filled pores which are known to conduct protons well. Blocking of protons by aquaporins is essential to maintain the electrochemical gradient across cellular and subcellular membranes. We studied the mechanism of proton exclusion in aquaporin-1 by multiple non-equilibrium molecular dynamics simulations that also allow proton transfer reactions. From the simulations, an effective free energy profile for the proton motion along the channel was determined with a maximum-likelihood approach. The results indicate that the main barrier is not, as had previously been speculated, caused by the interruption of the hydrogen-bonded water chain, but rather by an electrostatic field centered around the fingerprint Asn-Pro-Ala (NPA) motif. Hydrogen bond interruption only forms a secondary barrier located at the ar/R constriction region. The calculated main barrier height of 25-30 kJ mol(-1) matches the barrier height for the passage of protons across pure lipid bilayers and, therefore, suffices to prevent major leakage of protons through aquaporins. Conventional molecular dynamics simulations additionally showed that negatively charged hydroxide ions are prevented from being trapped within the NPA region by two adjacent electrostatic barriers of opposite polarity.