Ab Initio Molecular Dynamics Studies of the Liquid−Vapor Interface of an HCl Solution

Hydronium Ion Transport across the Liquid/Liquid Interface Assisted by a Phase-Transfer Catalyst: Structure and Thermodynamics Using Molecular Dynamics Simulation

Molecular dynamics simulations are used to examine the thermodynamic and structural aspects of the transfer of the classical hydronium ion (H3O+) across a water/1,2-dichloroethane (DCE) interface assisted by the phase-transfer catalyst (PTC) tetrakis(pentafluorophenyl) borate anion (TPFB-). The free energy of transfer from water to DCE of the H3O+-TPFB- ion pair is calculated to be 6 ± 1 kcal/mol, significantly less than that of the free hydronium ion (17 ± 1 kcal/mol). The ion pair is relatively stable at the interface and in the organic phase when it is accompanied by three water molecules with a small barrier to dissociation that supports its utility as a PTC. An examination of the hydration structure that accompanies the transfer of the ion pair shows that the ion pair, like the free hydronium ion, is transferred with the assistance of a finger-like water structure.

Double-Layer Distribution of Hydronium and Hydroxide Ions in the Air-Water Interface

The acid-base nature of the aqueous interface has long been controversial. Most macroscopic experiments suggest that the air-water interface is basic based on the detection of negative charges at the interface that indicates the enrichment of hydroxides (OH-), whereas microscopic studies mostly support the acidic air-water interface with the observation of hydronium (H3O+) accumulation in the top layer of the interface. It is crucial to clarify the interfacial preference of OH- and H3O+ ions for rationalizing the debate. In this work, we perform deep potential molecular dynamics simulations to investigate the preferential distribution of OH- and H3O+ ions at the aqueous interfaces. The neural network potential energy surface is trained based on density functional theory calculations with the SCAN functional, which can accurately describe the diffusion of these two ions both in the interface and in the bulk water. In contrast to the previously reported single ion enrichment, we show that both OH- and H3O+ surprisingly prefer to accumulate in interfaces but at different interfacial depths, rendering a double-layer ionic distribution within ∼1 nm near the Gibbs dividing surface. The H3O+ preferentially resides in the topmost layer of the interface, but the OH-, which is enriched in the deeper interfacial layer, has a higher equilibrium concentration due to the more negative free energy of interfacial stabilization [-0.90 (OH-) vs -0.56 (H3O+) kcal/mol]. The present finding of the ionic double-layer distribution may qualitatively offer a self-consistent explanation for the long-term controversy about the acid-base nature of the air-water interface.

Identifying the acidic or basic behavior of surface water: a QM/MM-MD study

Controversies on the water surface were theoretically addressed with the help of large scale quantum mechanical molecular dynamics (QMMD) simulations on water surface model systems with and without excess hydroniums and hydroxides. It was revealed that the thermodynamic surface structures of these ions strongly depend on their location and dipole orientation. Fast hydronium diffusion by proton transfer establishes a wider kinetic depth distribution (∼6 Å) than that predicted by its thermodynamic affinity for the water surface, while slow hydroxide is shallowly trapped below the outermost molecular layer (3-4 Å). In addition, the anisotropic orientation of surface water dipole can generate a substantial magnitude of surface potential, which extends to a depth of a few molecular layers. With these distinctively different surface properties of two ions and water molecules, the seemingly contradictory observations of acidic and negatively charged water surfaces may be successfully explained. That is, the negative surface charge of neutral water mostly stems from intrinsic water properties such as water dipole orientation and electron density spillage at the surface, rather than surface OH- ions. The enhanced acidity of the water surface can be attributed in large part to the kinetic depth profile of ion density in addition to static thermodynamic origin. Furthermore, the different depth profiles of the two ions may differently affect the surface-sensitive spectroscopic observations.

On the Propensity of Excess Hydroxide Ions at the Alcohol Monolayer-Water Interface

Atomistic molecular dynamics simulations have been employed to study the self-ion (H⁺ and OH^-) distribution at the interface between long-chain C₁₆-OH alcohol (cetyl alcohol) monolayer and water. It is well known that the free air-water interface is acidic due to accumulation of the hydronium (H₃O⁺) ions at the interface. In the present study, we have observed that contrary to the air-water interface, at the long-chain alcohol monolayer-water interface, it is the hydroxide (OH^-) ion, not the hydronium ion (H₃O⁺) that gets accumulated. By calculating the potential of mean forces, it is confirmed that there is extra stabilization for the OH^- ions at the interface relative to the bulk, but no such stabilization is observed for the H₃O⁺ ions. By analyzing the interaction of the self-ions with other constituents in the medium, it is clearly shown that the favorable interaction of the OH^- ions with the alcoholic -OH groups stabilizes this ion at the interface. By calculating coordination numbers of the self-ions it is observed that around 50% water neighbors are substituted by alcoholic -OH in case of the hydroxide ion at the interface, whereas in the case of hydronium ions, only 15% water neighbors are substituted by the alcoholic -OH. The most interesting observation about the local structure and H-bonding pattern is that the hydroxide ion acts solely as the H-bond acceptor, but the hydronium ion acts only as the H-bond donor.

Molecular origins of asymmetric proton conduction in the influenza M2 channel

The M2 proton channel of influenza A is embedded into the viral envelope and allows acidification of the virion when the external pH is lowered. In contrast, no outward proton conductance is observed when the internal pH is lowered, although outward current is observed at positive voltage. Residues Trp41 and Asp44 are known to play a role in preventing pH-driven outward conductance, but the mechanism for this is unclear. We investigate this issue using classical molecular dynamics simulations with periodic proton hops. When all key His37 residues are neutral, inward proton movement is much more facile than outward movement if the His are allowed to shuttle the proton. The preference for inward movement increases further as the charge on the His37 increases. Analysis of the trajectories reveals three factors accounting for this asymmetry. First, in the outward direction, Asp44 traps the hydronium by strong electrostatic interactions. Secondly, Asp44 and Trp41 orient the hydronium with the protons pointing inward, hampering outward Grotthus hopping. As a result, the effective barrier is lower in the inward direction. Trp41 adds to the barrier by weakly H-bonding to potential H+ acceptors. Finally, for charged His, the H3O+ in the inner vestibule tends to get trapped at lipid-lined fenestrations of the cone-shaped channel. Simulations qualitatively reproduce the experimentally observed higher outward conductance of mutants. The ability of positive voltage, unlike proton gradient, to induce an outward current appears to arise from its ability to bias H3O+ and the waters around it toward more H-outward orientations.

A critique of the capacitor-based "Transmembrane Electrostatically Localized Proton" hypothesis

In his Transmembrane Electrostatically Localized Proton hypothesis (TELP), James W. Lee has modeled the bioenergetic membrane as a simple capacitor. According to this model, the surface concentration of protons is completely independent of proton concentration in the bulk phase, and is linearly proportional to the transmembrane potential. Such a proportionality runs counter to the results of experimental measurements, molecular dynamics simulations, and electrostatics calculations. We show that the TELP model dramatically overestimates the surface concentration of protons, and we discuss the electrostatic reasons why a simple capacitor is not an appropriate model for the bioenergetic membrane.

The Proton in Biochemistry: Impacts on Bioenergetics, Biophysical Chemistry, and Bioorganic Chemistry

The proton is the smallest atomic particle, and in aqueous solution it is the smallest hydrated ion, having only two waters in its first hydration shell. In this article we survey key aspects of the proton in chemistry and biochemistry, starting with the definitions of pH and pK _a and their application inside biological cells. This includes an exploration of pH in nanoscale spaces, distinguishing between bulk and interfacial phases. We survey the Eigen and Zundel models of the structure of the hydrated proton, and how these can be used to explain: a) the behavior of protons at the water-hydrophobic interface, and b) the extraordinarily high mobility of protons in bulk water via Grotthuss hopping, and inside proteins via proton wires. Lastly, we survey key aspects of the effect of proton concentration and proton transfer on biochemical reactions including ligand binding and enzyme catalysis, as well as pH effects on biochemical thermodynamics, including the Chemiosmotic Theory. We find, for example, that the spontaneity of ATP hydrolysis at pH ≥ 7 is not due to any inherent property of ATP (or ADP or phosphate), but rather to the low concentration of H⁺. Additionally, we show that acidification due to fermentation does not derive from the organic acid waste products, but rather from the proton produced by ATP hydrolysis.

Dynamics and Surface Propensity of H+ and OH- within Rigid Interfacial Water: Implications for Electrocatalysis

Facile solvent reorganization promoting ion transfer across the solid-liquid interface is considered a prerequisite for efficient electrocatalysis. We provide first-principles insight into this notion by examining water self-ion dynamics at a highly rigid NaCl(100)-water interface. Through extensive density functional theory molecular dynamics simulations, we demonstrate for both acidic and alkaline solutions that Grotthuss dynamics is not impeded by a rigid water structure. Conversely, decreased proton transfer barriers and a striking propensity of H₃O⁺ and OH^- for stationary interfacial water are found. Differences in the ideal hydration structure of the ions, however, distinguish their behavior at the water contact layer. While hydronium can maintain its optimal solvation, the preferentially hypercoordinated hydroxide is repelled from the immediate vicinity of the surface due to interfacial coordination reduction. This has implications for alkaline hydrogen electrosorption in which the formation of undercoordinated OH^- at the surface is proposed to contribute to the observed sluggish kinetics.

Molecular Dynamics Studies on the Effect of Surface Roughness and Surface Tension on the Thermodynamics and Dynamics of Hydronium Ion Transfer Across the Liquid/Liquid Interface

Molecular dynamics simulations are used to examine the effect of surface roughness and surface tension on the transfer of the classical hydronium ion (H₃O⁺) across the water/1,2-dichloroethane interface. Free energy of transfer, hydration structure, and dynamics as a function of the ion location along the interface normal are calculated with six different values of a control parameter whose variation modifies the surface tension without impacting the bulk properties of the two solvents. Transfer of the classical hydronium ion across the water/1,2-dichloroethan interface involves the cotransfer of three hydration shell water molecules, independent of the surface tension. However, as the interaction between the two liquids weakens, a rise in interfacial tension and decrease in intrinsic water fingering and capillary fluctuations result in fewer water molecules cotransported with the ion in the second shell and a reduction in the length of the finger that the ion is attached to, consistent with the reduced size of the second hydration shell. First shell water residence time and lateral ion diffusion constants vary with the surface tension in a way that is consistent with the abovementioned structural insight.

Revisiting the Volmer-Heyrovský mechanism of hydrogen evolution on a nitrogen doped carbon nanotube: constrained molecular dynamics versus the nudged elastic band method

Density functional theory (DFT) based computational electrochemistry has the potential to serve as a tool with predictive power in the rational development and screening of electrocatalysts for renewable energy technologies. It is, however, of paramount importance that simulations are conducted rigorously at a level of theory that is sufficiently accurate in order to obtain physicochemically sensible results. Herein, we present a comparative study of the performance of the static climbing image nudged elastic band method (CI-NEB) vs. DFT based constrained molecular dynamics simulations with thermodynamic integration in estimating activation and reaction (free) energies of the Volmer-Heyrovský mechanism on a nitrogen doped carbon nanotube. Due to cancellation of errors within the CI-NEB calculations, static and dynamic activation barriers are observed to be surprisingly similar, while a substantial decrease in reaction energies is seen upon incorporation of solvent dynamics. This finding is attributed to two competing effects; (1) solvent reorganization that stabilizes the transition and, in particular, the product states with respect to the reactant state and (2) destabilizing entropic contributions due to solvent fluctuations. Our results highlight the importance of explicitly sampling the interfacial solvent dynamics when studying hydrogen evolution at solid-liquid interfaces.

Nature of Excess Hydrated Proton at the Water-Air Interface

Understanding the interfacial molecular structure of acidic aqueous solutions is important in the context of, e.g., atmospheric chemistry, biophysics, and electrochemistry. The hydration of the interfacial proton is necessarily different from that in the bulk, given the lower effective density of water at the interface, but has not yet been elucidated. Here, using surface-specific vibrational spectroscopy, we probe the response of interfacial protons at the water-air interface and reveal the interfacial proton continuum. Combined with spectral calculations based on ab initio molecular dynamics simulations, the proton at the water-air interface is shown to be well-hydrated, despite the limited availability of hydration water, with both Eigen and Zundel structures coexisting at the interface. Notwithstanding the interfacial hydrated proton exhibiting bulk-like structures, a substantial interfacial stabilization by -1.3 ± 0.2 kcal/mol is observed experimentally, in good agreement with our free energy calculations. The surface propensity of the proton can be attributed to the interaction between the hydrated proton and its counterion.

Formation and migration of H3O+ and OH− ions at the water/silica and water/vapor interfaces under the influence of a static electric field: a molecular dynamics study

Water ‘layers’ 1 and 2 in pink; ‘layer’ 3 in blue and green over portion of glass surface (grey). +90° field causes water migration and clustering.

Protons and Hydroxide Ions in Aqueous Systems

Understanding the structure and dynamics of water's constituent ions, proton and hydroxide, has been a subject of numerous experimental and theoretical studies over the last century. Besides their obvious importance in acid-base chemistry, these ions play an important role in numerous applications ranging from enzyme catalysis to environmental chemistry. Despite a long history of research, many fundamental issues regarding their properties continue to be an active area of research. Here, we provide a review of the experimental and theoretical advances made in the last several decades in understanding the structure, dynamics, and transport of the proton and hydroxide ions in different aqueous environments, ranging from water clusters to the bulk liquid and its interfaces with hydrophobic surfaces. The propensity of these ions to accumulate at hydrophobic surfaces has been a subject of intense debate, and we highlight the open issues and challenges in this area. Biological applications reviewed include proton transport along the hydration layer of various membranes and through channel proteins, problems that are at the core of cellular bioenergetics.

Surface Propensities of the Self-Ions of Water

The surface charge of water, which is important in a wide range of chemical, biological, material, and environmental contexts, has been a subject of lengthy and heated debate. Recently, it has been shown that the highly efficient LEWIS force field, in which semiclassical, independently mobile valence electron pairs capture the amphiproticity, polarizability and H-bonding of water, provides an excellent description of the solvation and dynamics of hydroxide and hydronium in bulk water. Here we turn our attention to slabs, cylinders, and droplets. In extended simulations with 1000 molecules, we find that hydroxide consistently prefers the surface, hydronium consistently avoids the surface, and the two together form an electrical double layer until neutralization occurs. The behavior of hydroxide can largely be accounted for by the observation that hydroxide moving to the surface loses fewer hydrogen bonds than are gained by the water molecule that it displaces from the surface. At the same time, since the orientation of the hydroxide increases the ratio of dangling hydrogens to dangling lone pairs, the proton activity of the exposed surface may be increased, rather than decreased. Hydroxide also moves more rapidly in the surface than in the bulk, likely because the proton donating propensity of neighboring water molecules is focused on the one hydrogen that is not dangling from the surface.

Propensity of Hydrated Excess Protons and Hydroxide Anions for the Air-Water Interface

Significant effort has been undertaken to better understand the molecular details governing the propensity of ions for the air-water interface. Facilitated by computationally efficient reactive molecular dynamics simulations, new and statistically conclusive molecular-scale results on the affinity of the hydrated excess proton and hydroxide anion for the air-water interface are presented. These simulations capture the dynamic bond breaking and formation processes (charge defect delocalization) that are important for correctly describing the solvation and transport of these complex species. The excess proton is found to be attracted to the interface, which is correlated with a favorable enthalpic contribution and consistent with reducing the disruption in the hydrogen bond network caused by the ion complex. However, a recent refinement of the underlying reactive potential energy function for the hydrated excess proton shows the interfacial attraction to be weaker, albeit nonzero, a result that is consistent with the experimental surface tension measurements. The influence of a weak hydrogen bond donated from water to the protonated oxygen, recently found to play an important role in excess hydrated proton transport in bulk water, is seen to also be important for this study. In contrast, the hydroxide ion is found to be repelled from the air-water interface. This repulsion is characterized by a reduction of the energetically favorable ion-water interactions, which creates an enthalpic penalty as the ion approaches the interface. Finally, we find that the fluctuation in the coordination number around water sheds new light on the observed entropic trends for both ions.

Toward a unified picture of the water self-ions at the air-water interface: a density functional theory perspective

The propensities of the water self-ions, H3O(+) and OH(-), for the air-water interface have implications for interfacial acid-base chemistry. Despite numerous experimental and computational studies, no consensus has been reached on the question of whether or not H3O(+) and/or OH(-) prefer to be at the water surface or in the bulk. Here we report a molecular dynamics simulation study of the bulk vs interfacial behavior of H3O(+) and OH(-) that employs forces derived from density functional theory with a generalized gradient approximation exchange-correlation functional (specifically, BLYP) and empirical dispersion corrections. We computed the potential of mean force (PMF) for H3O(+) as a function of the position of the ion in the vicinity of an air-water interface. The PMF suggests that H3O(+) has equal propensity for the interface and the bulk. We compare the PMF for H3O(+) to our previously computed PMF for OH(-) adsorption, which contains a shallow minimum at the interface, and we explore how differences in solvation of each ion at the interface vs in the bulk are connected with interfacial propensity. We find that the solvation shell of H3O(+) is only slightly dependent on its position in the water slab, while OH(-) partially desolvates as it approaches the interface, and we examine how this difference in solvation behavior is manifested in the electronic structure and chemistry of the two ions.

Negligible water surface charge determined using Kelvin probe and total reflection X-ray fluorescence techniques

The water surface charge has been extensively debated in recent decades. Electrophoretic mobilities of air bubbles in water and disjoining pressures between the surfaces of aqueous films suggest that the surface of water exhibits a significant negative charge. This is commonly attributed to a strong adsorption of hydroxide ions at the interface, though spectroscopic measurements and simulation studies suggest surface depletion of hydroxide ions. Alternatively, the negative surface charge could arise from surface contamination with trace charged surfactants. We have probed the variation in the surface charge of water with pH by measuring surface potentials using the Kelvin probe technique. Independently, the abundance in the interfacial layer of "reporter ions" (Rb(+) and Br(-)), which must be affected by a charged surface, has been monitored using the total reflection X-ray fluorescence (TRXF) technique. Special care was taken to prove the high sensitivity of this technique as well as to avoid surface contaminants. The magnitude of the surface charge was found to be below 1 e per 500 nm(2) (TRXF). No evidence of variations in the surface potential between pH 2-3 and pH 9-12 was detected within the accuracies of the methods (5 mV for Kelvin probe and 2 mV for TRXF). Hence, our findings suggest that the clean water surface exhibits negligible charge in a wide pH range.

Water at hydrophobic interfaces delays proton surface-to-bulk transfer and provides a pathway for lateral proton diffusion

Fast lateral proton migration along membranes is of vital importance for cellular energy homeostasis and various proton-coupled transport processes. It can only occur if attractive forces keep the proton at the interface. How to reconcile this high affinity to the membrane surface with high proton mobility is unclear. Here, we tested whether a minimalistic model interface between an apolar hydrophobic phase (n-decane) and an aqueous phase mimics the biological pathway for lateral proton migration. The observed diffusion span, on the order of tens of micrometers, and the high proton mobility were both similar to the values previously reported for lipid bilayers. Extensive ab initio simulations on the same water/n-decane interface reproduced the experimentally derived free energy barrier for the excess proton. The free energy profile G(H(+)) adopts the shape of a well at the interface, having a width of two water molecules and a depth of 6 ± 2RT. The hydroniums in direct contact with n-decane have a reduced mobility. However, the hydroniums in the second layer of water molecules are mobile. Their in silico diffusion coefficient matches that derived from our in vitro experiments, (5.7 ± 0.7) 10(-5) cm(2) s(-1). Conceivably, these are the protons that allow for fast diffusion along biological membranes.

Analysis of anisotropic local field in sum frequency generation spectroscopy with the charge response kernel water model

A new flexible and polarizable water model based on the charge response kernel (CRK) theory is developed for the analysis of sum frequency generation (SFG) spectroscopy. The CRK model well describes several bulk water properties and SFG spectrum by molecular dynamics (MD) calculations. While the flexible and polarizable MD simulation generally adopts the short-range damping of intermolecular interaction, it is found that the same procedure is not adequate for the calculation of transition dipole in strongly hydrogen bonding environment. Accordingly, the improved calculation of the nonlinear susceptibility of water surface results in the positive imaginary part in the 3000-3200 cm(-1) region, which is consistent with recent phase-sensitive experiments. The mechanism of the positive region is attributed to the anisotropic local field effect induced by the orientational correlation of surface water.