Proton transfer through the water gossamer

Vibrational dynamics of aqueous hydroxide solutions probed using broadband 2DIR spectroscopy

We employed ultrafast transient absorption and broadband 2DIR spectroscopy to study the vibrational dynamics of aqueous hydroxide solutions by exciting the O-H stretch vibrations of the strongly hydrogen-bonded hydroxide solvation shell water and probing the continuum absorption of the solvated ion between 1500 and 3800 cm(-1). We observe rapid vibrational relaxation processes on 150-250 fs time scales across the entire probed spectral region as well as slower vibrational dynamics on 1-2 ps time scales. Furthermore, the O-H stretch excitation loses its frequency memory in 180 fs, and vibrational energy exchange between bulk-like water vibrations and hydroxide-associated water vibrations occurs in ∼200 fs. The fast dynamics in this system originate in strong nonlinear coupling between intra- and intermolecular vibrations and are explained in terms of non-adiabatic vibrational relaxation. These measurements indicate that the vibrational dynamics of the aqueous hydroxide complex are faster than the time scales reported for long-range transport of protons in aqueous hydroxide solutions.

Why Is MP2-Water "Cooler" and "Denser" than DFT-Water?

Density functional theory (DFT) with a dispersionless generalized gradient approximation (GGA) needs much higher temperature and pressure than the ambient conditions to maintain water in the liquid phase at the correct (1 g/cm(3)) density during first-principles simulations. Conversely, ab initio second-order many-body perturbation (MP2) calculations of liquid water require lower temperature and pressure than DFT/GGA to keep water liquid. Here we present a unifying explanation of these trends derived from classical water simulations using a polarizable force field with different sets of parameters. We show that the different temperatures and pressures between DFT/GGA and MP2 at which the simulated water displays the experimentally observed liquid structure under the ambient conditions can be largely explained by their differences in polarizability and dispersion interaction, respectively. In DFT/GGA, the polarizability and thus the induced dipole moments and the hydrogen-bond strength are all overestimated. This hinders the rotational motion of molecules and requires a higher temperature for DFT-water to be liquid. MP2 gives a stronger dispersion interaction and thus shorter intermolecular distances than dispersionless DFT/GGA, which is why MP2-water is denser than DFT-water under the same external pressure.

Deprotonation of formic acid in collisions with a liquid water surface studied by molecular dynamics and metadynamics simulations

Formic acid has a lower barrier to deprotonation at the air–water interface than in bulk liquid water.

Ab initio molecular dynamics study of an aqueous NaCl solution under an electric field

Ab initio molecular dynamics simulations of salty water under an electric field reveal two regimes of the relative mobilities of chlorine and sodium ions. When water dissociation and proton transfer are actived at strong field intensities, the presence of the ions hinders the efficiency of the proton transfer mechanism.

Modeling pH variation in reverse osmosis

The transport of hydronium and hydroxide ions through reverse osmosis membranes constitutes a unique case of ionic species characterized by uncommonly high permeabilities. Combined with electromigration, this leads to complex behavior of permeate pH, e.g., negative rejection, as often observed for monovalent ions in nanofiltration of salt mixtures. In this work we employed a rigorous phenomenological approach combined with chemical equilibrium to describe the trans-membrane transport of hydronium and hydroxide ions along with salt transport and calculate the resulting permeate pH. Starting from the Nernst-Planck equation, a full non-linear transport equation was derived, for which an approximate solution was proposed based on the analytical solution previously developed for trace ions in a dominant salt. Using the developed approximate equation, transport coefficients were deduced from experimental results obtained using a spiral wound reverse osmosis module operated under varying permeate flux (2-11 μm/s), NaCl feed concentrations (0.04-0.18 M) and feed pH values (5.5-9.0). The approximate equation agreed well with the experimental results, corroborating the finding that diffusion and electromigration, rather than a priori neglected convection, were the major contributors to the transport of hydronium and hydroxide. The approach presented here has the potential to improve the predictive capacity of reverse osmosis transport models for acid-base species, thereby improving process design/control.

Role of Presolvation and Anharmonicity in Aqueous Phase Hydrated Proton Solvation and Transport

Results from condensed phase ab initio molecular dynamics (AIMD) simulations suggest a proton transfer reaction is facilitated by "presolvation" in which the hydronium is transiently solvated by four water molecules, similar to the typical solvation structure of water, by accepting a weak hydrogen bond from the fourth water molecule. A new version 3.2 multistate empirical valence bond (MS-EVB 3.2) model for the hydrated excess proton incorporating this presolvation behavior is therefore developed. The classical MS-EVB simulations show similar structural properties as those of the previous model but with significantly improved diffusive behavior. The inclusion of nuclear quantum effects in the MS-EVB also provides an even better description of the proton diffusion rate. To quantify the influence of anharmonicity, a second model (aMS-EVB 3.2) is developed using the anharmonic aSPC/Fw water model, which provides similar structural properties but improved spectroscopic responses at high frequencies.

Ultrafast 2D IR spectroscopy of the excess proton in liquid water

Despite decades of study, the structures adopted to accommodate an excess proton in water and the mechanism by which they interconvert remain elusive. We used ultrafast two-dimensional infrared (2D IR) spectroscopy to investigate protons in aqueous hydrochloric acid solutions. By exciting O–H stretching vibrations and detecting the spectral response throughout the mid-IR region, we observed the interaction between the stretching and bending vibrations characteristic of the flanking waters of the Zundel complex, [H(H2O)2]+, at 3200 and 1760 cm−1, respectively. From time-dependent shifts of the stretch-bend cross peak, we determined a lower limit on the lifetime of this complex of 480 femtoseconds. These results suggest a key role for the Zundel complex in aqueous proton transfer.

How amide hydrogens exchange in native proteins

Amide hydrogen exchange (HX) is widely used in protein biophysics even though our ignorance about the HX mechanism makes data interpretation imprecise. Notably, the open exchange-competent conformational state has not been identified. Based on analysis of an ultralong molecular dynamics trajectory of the protein BPTI, we propose that the open (O) states for amides that exchange by subglobal fluctuations are locally distorted conformations with two water molecules directly coordinated to the N-H group. The HX protection factors computed from the relative O-state populations agree well with experiment. The O states of different amides show little or no temporal correlation, even if adjacent residues unfold cooperatively. The mean residence time of the O state is ∼100 ps for all examined amides, so the large variation in measured HX rate must be attributed to the opening frequency. A few amides gain solvent access via tunnels or pores penetrated by water chains including native internal water molecules, but most amides access solvent by more local structural distortions. In either case, we argue that an overcoordinated N-H group is necessary for efficient proton transfer by Grotthuss-type structural diffusion.

Proton transfer in a short hydrogen bond caused by solvation shell fluctuations: an ab initio MD and NMR/UV study of an (OHO)(-) bonded system

We present a joint experimental and quantum chemical study on the influence of solvent dynamics on the protonation equilibrium in a strongly hydrogen bonded phenol-acetate complex in CD2Cl2. Particular attention is given to the correlation of the proton position distribution with the internal conformation of the complex itself and with fluctuations of the aprotic solvent. Specifically, we have focused on a complex formed by 4-nitrophenol and tetraalkylammonium-acetate in CD2Cl2. Experimentally we have used combined low-temperature (1)H and (13)C NMR and UV-vis spectroscopy and showed that a very strong OHO hydrogen bond is formed with proton tautomerism (PhOH···(-)OAc and PhO(-)···HOAc forms, both strongly hydrogen bonded). Computationally, we have employed ab initio molecular dynamics (70 and 71 solvent molecules, with and without the presence of a counter-cation, respectively). We demonstrate that the relative motion of the counter-cation and the "free" carbonyl group of the acid plays the major role in the OHO bond geometry and causes proton "jumps", i.e. interconversion of PhOH···(-)OAc and PhO(-)···HOAc tautomers. Weak H-bonds between CH(CD) groups of the solvent and the oxygen atom of carbonyl stabilize the PhOH···(-)OAc type of structures. Breaking of CH···O bonds shifts the equilibrium towards PhO(-)···HOAc form.

Ab initio molecular dynamics simulations of aqueous triflic acid confined in carbon nanotubes

Ab initio molecular dynamics simulations were performed to investigate the effects of nanoscale confinement on the structural and dynamical properties of aqueous triflic acid (CF3SO3H). Single-walled carbon nanotubes (CNTs) with diameters ranging from ∼11 to 14 Å were used as confinement vessels, and the inner surface of the CNT were either left bare or fluorinated to probe the influence of the confined environment on structural and dynamical properties of the water and triflic acidic. The systems were simulated at hydration levels of n = 1-3 H2O/CF3SO3H. Proton dissociation expectedly increased with increasing hydration. Along with the level of hydration, hydrogen bond connectivity between the triflic acid molecules, both directly and via a single water molecule, played a role on proton dissociation. Direct hydrogen bonding between the CF3SO3H molecules, most commonly found in the larger bare CNT, also promoted interactions between water molecules allowing for greater separation of the dissociated protons from the CF3SO3(-) as the hydration level was increased. However, this also resulted in a decrease in the overall proportion of dissociated protons. The confinement dimensions altered both the hydrogen bond network and the distribution of water molecules where the H2O in the fluorinated CNTs tended to form small clusters with less proton dissociation at n = 1 and 2 but the highest at n = 3. In the absence of nearby hydrogen bond accepting sites from H2O or triflic acid SO3H groups, the water molecules formed weak hydrogen bonds with the fluorine atoms. In the bare CNT systems, these involved the CF3 groups of triflic acid and were more frequently observed when direct hydrogen bonding between CF3SO3H hindered potential hydrogen bonding sites. In the fluorinated tubes, interactions with the covalently bound fluorine atoms of the CNT wall dominated which appear to stabilize the hydrogen bond network. Increasing the hydration level increased the frequency of the OH···F (CNT) hydrogen bonding which was highly pronounced in the smaller fluorinated CNT indicating an influence on the confinement dimensions on these interactions.

Collective vibrations of water-solvated hydroxide ions investigated with broadband 2DIR spectroscopy

The infrared spectra of aqueous solutions of NaOH and other strong bases exhibit a broad continuum absorption for frequencies between 800 and 3500 cm(-1), which is attributed to the strong interactions of the OH(-) ion with its solvating water molecules. To provide molecular insight into the origin of the broad continuum absorption feature, we have performed ultrafast transient absorption and 2DIR experiments on aqueous NaOH by exciting the O-H stretch vibrations and probing the response from 1350 to 3800 cm(-1) using a newly developed sub-70 fs broadband mid-infrared source. These experiments, in conjunction with harmonic vibrational analysis of OH(-)(H2O)n (n = 17) clusters, reveal that O-H stretch vibrations of aqueous hydroxides arise from coupled vibrations of multiple water molecules solvating the ion. We classify the vibrations of the hydroxide complex by symmetry defined by the relative phase of vibrations of the O-H bonds hydrogen bonded to the ion. Although broad and overlapping spectral features are observed for 3- and 4-coordinate ion complexes, we find a resolvable splitting between asymmetric and symmetric stretch vibrations, and assign the 2850 cm(-1) peak infrared spectra of aqueous hydroxides to asymmetric stretch vibrations.

Promoting the "water-wire" mechanism of double proton transfer in [2,2'-bipyridyl]-3,3'-diol by porous gold nanoparticles

The effect of nanopores in porous gold nanoparticles (Au NPs) on excited-state double proton transfer (DPT) in [2,2'-bipyridyl]-3,3'-diol (BP(OH)2) in an aqueous environment is the main focus of the present work. DPT in BP(OH)2 is known to take place through two mechanisms. In a bulk environment, an open solvated molecule facilitates the process and emits at 460 nm whereas, in a confined situation, formation of a "water wire" between the prototropic centers leads to the transfer of protons. It has been shown spectroscopically in the present study that in the nanovessels provided by nanoporous Au NPs, the unconventional mechanism of DPT through the formation of a "water wire" is promoted due to the presence of a limited number of water molecules around the probe. Experiments in the presence of solid pure Au, Ag and Au/Ag NPs were performed to support our proposition. Time-resolved fluorescence spectral changes confirm our findings.

Developing a general interaction potential for hydrophobic and hydrophilic interactions

We review direct force measurements on a broad class of hydrophobic and hydrophilic surfaces. These measurements have enabled the development of a general interaction potential per unit area, W(D) = -2γ(i)Hy exp(-D/D(H)) in terms of a nondimensional Hydra parameter, Hy, that applies to both hydrophobic and hydrophilic interactions between extended surfaces. This potential allows one to quantitatively account for additional attractions and repulsions not included in the well-known combination of electrostatic double layer and van der Waals theories, the so-called Derjaguin-Landau-Verwey-Overbeek (DLVO) theory. The interaction energy is exponentially decaying with decay length D(H) ≈ 0.3-2 nm for both hydrophobic and hydrophilic interactions, with the exact value of D(H) depending on the precise system and conditions. The pre-exponential factor depends on the interfacial tension, γ(i), of the interacting surfaces and Hy. For Hy > 0, the interaction potential describes interactions between partially hydrophobic surfaces, with the maximum hydrophobic interaction (i.e., two fully hydrophobic surfaces) corresponding to Hy = 1. Hydrophobic interactions between hydrophobic monolayer surfaces measured with the surface forces apparatus (SFA) are shown to be well described by the proposed interaction potential. The potential becomes repulsive for Hy < 0, corresponding to partially hydrophilic (hydrated) interfaces. Hydrated surfaces such as mica, silica, and lipid bilayers are discussed and reviewed in the context of the values of Hy appropriate for each system.

Unveiling the Janus-Like Properties of OH(-)

Using ab initio simulations, we explore the glassy landscape of the OH(-)(H2O)20 cluster and its infrared spectrum. We show that the OH(-) has an amphiphilic Janus-type behavior like the hydronium ion induced by the ability of its O-H bond to be buried inside of the cluster or exposed at the surface with different coordination numbers. Recent infrared experiments of aqueous NaOH have found two pronounced peaks at 2000 and 2850 cm(-1) [Mandal, A.; J. Chem. Phys. 2014, 140, 1-12]. The microscopic origins of these spectral features remain elusive. Herein, we disentangle the contribution of the spectra between 1700 and 3000 cm(-1) in terms of the microscopic solvation structure of OH(-) and dub this as the amphiphilic band. The delocalized nature of OH(-) results in a red shift to the O-H stretch, which mixes with bend-vibrations, the extent to which is tuned by the local coordination number. These results have important bearing on understanding the spectroscopic signatures of OH(-) in environments like the air-water interface.

A THz/FTIR fingerprint of the solvated proton: evidence for Eigen structure and Zundel dynamics

Zundel (orange), Eigen (red) and hydration water (light blue) contributions to the THz/FIR extinction of the solvated proton.

Unraveling the intramolecular cyclization mechanism of oxidized tryptophan in aqueous solution as a function of pH

pH tunes the mechanism of the intramolecular cyclization of 3a-substituted tryptophan derivatives.

Mechanism for formation of atmospheric Cl atom precursors in the reaction of dinitrogen oxides with HCl/Cl− on aqueous films

Formation of atmospheric chlorine atom precursors ClNO₂ and ClNO in the reaction of HCl with oxides of nitrogen on a water film: left – formation of N–Cl bond as N–O bond breaks; right – concurrent changes in Mulliken charges.

Local and collective reaction coordinates in the transport of the aqueous hydroxide ion

We investigate local and collective reaction coordinates for the structural diffusion of the hydroxide ion in dilute aqueous NaOH solution using a multistate empirical valence bond (MS-EVB) simulation. We characterize a 15 fs time scale associated with shifting of the equally shared proton within a Zundel-like H3O2(-) ion to form a water molecule, a 550 fs relaxation from this transition state largely guided by electrostatic fluctuations of the surrounding environment, and a 9.6 ps time scale that corresponds to the solvation of the water molecule formed by the proton transfer event. When individual proton transfer events are examined, we are unable to identify a unique transition state solely on the basis of a decrease in the hydroxide ion's coordination number. Instead, we find that the collective electric field along the proton transfer direction is better suited to describe the creation and relaxation of Zundel-like transition states that allow structural diffusion of the hydroxide ion.

Philosophy of voltage-gated proton channels

In this review, voltage-gated proton channels are considered from a mainly teleological perspective. Why do proton channels exist? What good are they? Why did they go to such lengths to develop several unique hallmark properties such as extreme selectivity and ΔpH-dependent gating? Why is their current so minuscule? How do they manage to be so selective? What is the basis for our belief that they conduct H(+) and not OH(-)? Why do they exist in many species as dimers when the monomeric form seems to work quite well? It is hoped that pondering these questions will provide an introduction to these channels and a way to logically organize their peculiar properties as well as to understand how they are able to carry out some of their better-established biological functions.

Anomalous water diffusion in salt solutions

The dynamics of water exhibits anomalous behavior in the presence of different electrolytes. Recent experiments [Kim JS, Wu Z, Morrow AR, Yethiraj A, Yethiraj A (2012) J Phys Chem B 116(39):12007-12013] have found that the self-diffusion of water (Dw) can either be enhanced or suppressed around CsI and NaCl, respectively, relative to that of neat water. Here we show that unlike classical empirical potentials, ab initio molecular dynamics simulations successfully reproduce the qualitative trends observed experimentally. These types of phenomena have often been rationalized in terms of the "structure-making" or "structure-breaking" effects of different ions on the solvent, although the microscopic origins of these features have remained elusive. Rather than disrupting the network in a significant manner, the electrolytes studied here cause rather subtle changes in both structural and dynamical properties of water. In particular, we show that water in the ab initio molecular dynamics simulations is characterized by dynamic heterogeneity, which turns out to be critical in reproducing the experimental trends.

Water under inner pressure: a dielectric spectroscopy study

Water is the most studied substance on Earth. However, it is not completely understood why its structural and dynamical properties give rise to some anomalous behaviors. Some of them emerge when experiments at low temperatures and/or high pressures are performed. Here we report dielectric measurements on cold water under macroscopically constrained conditions, i.e., water in a large container at constant volume that cannot freeze below the melting point. The inner pressure in these conditions shifts the α relaxation peak to similar frequencies as seen in ice Ih. At 267 K we observe a peculiar response possibly due to the Grotthuss mechanism. At 251 K (the triple point) ice III forms.

Competition between excited state proton and OH− transport via a short water wire: solvent effects open the gate

We investigate the acid–base proton exchange reaction in a microsolvated bifunctional chromophore by means of quantum chemical calculations.