Proton transfer and the diffusion of H+ and OH− ions along water wires

Classical Models of Hydroxide for Proton Hopping Simulations

Hydronium (H3O+) and hydroxide (OH-) ions perform structural diffusion in water via sequential proton transfers ("Grotthuss hopping"). This phenomenon can be accounted for by interspersing stochastic proton transfer events in classical molecular dynamics simulations. The implementation of OH--mediated proton hopping is particularly challenging because classical force fields are known to produce overcoordinated solvation structures around the OH- ion. Here, we first explore the ability of two-particle point-charge models to reproduce both the solvation free energy and coordination number in TIP3P water. We find that this is possible only with unphysical changes in the nonbonded parameters which create problems in proton hopping simulations. We then construct a classical OH- model with the charge of oxygen distributed among three auxiliary particles. This model favors a lower coordination number by accepting three hydrogen bonds and weakly donating one. The model was implemented in the MOBHY module of the CHARMM program and was fit to reproduce the experimental aqueous diffusion coefficient of OH-. This parameterization gave reasonable electrophoretic mobilities and the expected accelerated transport under nanoconfinement.

Fluids and Electrolytes under Confinement in Single-Digit Nanopores

Confined fluids and electrolyte solutions in nanopores exhibit rich and surprising physics and chemistry that impact the mass transport and energy efficiency in many important natural systems and industrial applications. Existing theories often fail to predict the exotic effects observed in the narrowest of such pores, called single-digit nanopores (SDNs), which have diameters or conduit widths of less than 10 nm, and have only recently become accessible for experimental measurements. What SDNs reveal has been surprising, including a rapidly increasing number of examples such as extraordinarily fast water transport, distorted fluid-phase boundaries, strong ion-correlation and quantum effects, and dielectric anomalies that are not observed in larger pores. Exploiting these effects presents myriad opportunities in both basic and applied research that stand to impact a host of new technologies at the water-energy nexus, from new membranes for precise separations and water purification to new gas permeable materials for water electrolyzers and energy-storage devices. SDNs also present unique opportunities to achieve ultrasensitive and selective chemical sensing at the single-ion and single-molecule limit. In this review article, we summarize the progress on nanofluidics of SDNs, with a focus on the confinement effects that arise in these extremely narrow nanopores. The recent development of precision model systems, transformative experimental tools, and multiscale theories that have played enabling roles in advancing this frontier are reviewed. We also identify new knowledge gaps in our understanding of nanofluidic transport and provide an outlook for the future challenges and opportunities at this rapidly advancing frontier.

Hydroxide Diffusion in Functionalized Cylindrical Nanopores as Idealized Models of Anion Exchange Membrane Environments: An Ab Initio Molecular Dynamics Study

Anion exchange membranes (AEMs) have attracted significant interest for their applications in fuel cells and other electrochemical devices in recent years. Understanding water distributions and hydroxide transport mechanisms within AEMs is critical to improving their performance as concerns hydroxide conductivity. Recently, nanoconfined environments have been used to mimic AEM environments. Following this approach, we construct nanoconfined cylindrical pore structures using graphane nanotubes (GNs) functionalized with trimethylammonium cations as models of local AEM morphology. These structures were then used to investigate hydroxide transport using ab initio molecular dynamics (AIMD). The simulations showed that hydroxide transport is suppressed in these confined environments relative to the bulk solution although the mechanism is dominated by structural diffusion. One factor causing the suppressed hydroxide transport is the reduced proton transfer (PT) rates due to changes in hydroxide and water solvation patterns under confinement compared to bulk solution as well as strong interactions between hydroxide ions and the tethered cation groups.

Hydroxide ion-conducting viologen-bakelite organic frameworks for flexible solid-state zinc-air battery applications

Adaptable polymer-based solid-state electrolytes can be a game-changer toward safe, lightweight flexible batteries. We present a robust Bakelite-type organic polymer covalently decked with viologen, triazine, and phenolic moieties. Its flexible structure with cationic viologen centers incorporates counter-balancing free hydroxide ions into the polymeric framework. By design, the aromatic groups and heteroatoms in the framework can be activated under an applied potential to prompt a push-pull drive, setting off the towing of hydroxide ions via weak electrostatic, van der Waals, and hydrogen-bond interactions. The frontier orbitals from a DFT-modeled structure certify this. The hydroxyl-polymer requires minimal KOH wetting to maintain a humid environment for Grotthuss-type transport. The hydroxide ion conductivity reaches a value of 1.4 × 10^-2 S cm^-1 at 80 °C and 95% RH, which is retained for over 15 h. We enhanced its practical utility by coating it as a thin solid-state separator-cum-electrolyte on readily available filter paper. The composite exhibits a conductivity of 4.5 × 10^-3 S cm^-1 at 80 °C and 95% RH. A zinc-air battery (ZAB) constructed using this polymer-coated paper as electrolyte yields a maximum power density of 115 mW cm^-2 and high specific capacitance of 435 mA h g^-1. The power density recorded for our ZAB is among the best reported for polymer electrolyte-based batteries. Subsequently, the flexible battery fabricated with IISERP-POF11_OH@FilterPaper exhibits an OCV of 1.44 V, and three batteries in series power a demo traffic signal. To underscore the efficiency of hydroxide ion transport through the complex multifunctional backbone of the polymer, we calculated the diffusion coefficient for OH^- (Exp: 2.9 × 10^-5 cm² s^-1; Comp. 5.2 × 10^-6 cm² s^-1) using electrochemical methods and MD simulations. Climbing-edge NEB calculations reveal a large energy barrier of 2.11 eV for Zn²⁺ to penetrate the polymer and identify hydroxide ions within the polymer, suggesting no undesirable Zn²⁺ crossover. Our findings assert the readily accessible C-C-linked cationic polymer's capacity as a solid-state electrolyte for ZABs and any anion-conducting membrane.

Putting together the puzzle of ion transfer in single-digit carbon nanotubes: mean-field meets ab initio

Nature employs channel proteins to selectively pass water across cell membranes, which inspires the search for bio-mimetic analogues. Carbon nanotube porins (CNTPs) are intriguing mimics of water channels, yet ion transport in CNTPs still poses questions. As an alternative to continuum models, here we present a molecular mean-field model that transparently describes ion coupling, yet unlike continuum models, computes ab initio all required thermodynamic quantities for the KCl salt and H⁺ and OH^- ions present in water. Starting from water transfer, the model considers the transfer of free ions, along with ion-pair formation as a proxy of non-mean-field ion-ion interactions. High affinity to hydroxide, suggested by experiments, making it a dominant charge carrier in CNTPs, is revealed as an exceptionally favorable transfer of KOH pairs. Nevertheless, free ions, coexisting with less mobile ion-pairs, apparently control ion transport. The model well explains the observed effects of salt concentration and pH on conductivity, transport numbers, anion permeation and its activation energies, and current rectification. The proposed approach is extendable to other sub-nanochannels and helps design novel osmotic materials and devices.

Systematic Evaluation of Ion Diffusion and Water Exchange

As a fundamental property of all fluids, diffusion plays myriad roles in both science and our daily lives. Diffusive properties of many liquids including water have been extensively studied both experimentally and theoretically, while for transition metal ions, there exist significant experimental data that have not been extensively studied theoretically. Hence, high-confidence predictions for challenging systems like radioactive ions that are biohazardous cannot be reliably generated. In this work, a workflow named ISAIAH (Ion Simulation using AMBER for dIffusion Action when Hydrated) was designed to accurately simulate the diffusion coefficients of 15 monoatomic ions with charges varying from -1 to +3 in four water models. As the results indicate, good agreement with experimental values was achieved, leading us to select ²³⁹Pu⁴⁺ (for which no experimental data are available) as a candidate ion to make a theoretical prediction of its diffusion coefficient in water. Among all the force field parameter sets, the ones parametrized using an augmented 12-6-4 Lennard-Jones (LJ) potential showed lower average unsigned errors (AUE) for ions of various radii and electron configurations relative to some 12-6 LJ parameters. This observation agrees well with the fact that diffusion is affected by both the hydration free energy (HFE) and the ion-oxygen distance (IOD) between solute and solvent molecules, both of which are handled well by the 12-6-4 model.

Confinement-Controlled Aqueous Chemistry within Nanometric Slit Pores

In this Focus Review, we put the spotlight on very recent insights into the fascinating world of wet chemistry in the realm offered by nanoconfinement of water in mechanically rather rigid and chemically inert planar slit pores wherein only monolayer and bilayer water lamellae can be hosted. We review the effect of confinement on different aspects such as hydrogen bonding, ion diffusion, and charge defect migration of H⁺(aq) and OH^-(aq) in nanoconfined water depending on slit pore width. A particular focus is put on the strongly modulated local dielectric properties as quantified in terms of anisotropic polarization fluctuations across such extremely confined water films and their putative effects on chemical reactions therein. The stunning findings disclosed only recently extend wet chemistry in particular and solvation science in general toward extreme molecular confinement conditions.

Minimizing the electrosorption of water from humid ionic liquids on electrodes

In supercapacitors based on ionic liquid electrolytes, small amounts of absorbed water could potentially reduce the electrochemical window of electrolytes and cause performance degradation. The same would take place if ionic liquids are used as solvents for electrocatalysis involving the dissolved molecular species. In this work, we carry out molecular dynamics simulations, with gold and carbon electrodes in typical ionic liquids, hydrophobic and hydrophilic, to study electrosorption of water. We investigate the effects of hydrophobicity/hydrophilicity of ionic liquids and electrodes on interfacial distribution of ions and electrosorbed water. Results reveal that using hydrophilic ionic liquids would help to keep water molecules away from the negatively charged electrodes, even at large electrode polarizations. This conclusion is supported by electrochemical cyclic voltammetry measurements on gold and carbon electrodes in contact with humid ionic liquids. Thereby, our findings suggest potential mechanisms for protection of electrodes from water electrosorption.

Ionic liquid electrolytes can impart increased operational voltage and energy density in supercapacitors, but water may diminish performance. Here the authors show that the hydrophilicity/hydrophobicity of ionic liquids can influence electrosorption of water and ultimately the supercapacitor performance.

Role of Extracellular Carbonic Anhydrase in Dissolved Inorganic Carbon Uptake in Alkaliphilic Phototrophic Biofilm

Alkaline Soda Lakes are extremely productive ecosystems, due to their high dissolved inorganic carbon (DIC) concentrations. Here, we studied the dynamics of the carbonate system, in particular, the role of extracellular carbonic anhydrase (eCA) of an alkaliphilic phototrophic biofilm composed of bacteria enriched from soda lake benthic mats. By using measurements with microsensors and membrane inlet mass spectrometry, combined with mathematical modeling, we show how eCA controls DIC uptake. In our experiments, the activity of eCA varied four-fold, and was controlled by the bicarbonate concentration during growth: a higher bicarbonate concentration led to lower eCA activity. Inhibition of eCA decreased both the net and the gross photosynthetic productivities of the investigated biofilms. After eCA inhibition, the efflux of carbon dioxide (CO2) from the biofilms increased two- to four-fold. This could be explained by the conversion of CO2, leaking from cyanobacterial cells, by eCA, to bicarbonate. Bicarbonate is then taken up again by the cyanobacteria. In suspensions, eCA reduced the CO2 leakage to the bulk medium from 90 to 50%. In biofilms cultivated at low bicarbonate concentration (~0.13 mM), the oxygen production was reduced by a similar ratio upon eCA inhibition. The role of eCA in intact biofilms was much less significant compared to biomass suspensions, as CO2 loss to the medium is reduced due to mass transfer resistance.

One-Dimensional Confinement Inhibits Water Dissociation in Carbon Nanotubes

The effect of nanoconfinement on the self-dissociation of water constitutes an open problem whose elucidation poses a serious challenge to experiments and simulations alike. In slit pores of width ≈1 nm, recent first-principles calculations have predicted that the dissociation constant of H₂O increases by almost 2 orders of magnitude [ Muñoz-Santiburcio and Marx, Phys. Rev. Lett. 2017 , 119 , 056002 ]. In the present study, quantum mechanics-molecular mechanics simulations are employed to compute the dissociation free-energy profile of water in a (6,6) carbon nanotube. According to our results, the equilibrium constant K_w drops by 3 orders of magnitude with respect to the bulk phase value, at variance with the trend predicted for confinement in two dimensions. The higher barrier to dissociation can be ascribed to the undercoordination of the hydroxide and hydronium ions in the nanotube and underscores that chemical reactivity does not exhibit a monotonic behavior with respect to pore size but may vary substantially with the characteristic length scale and dimensionality of the confining media.

Classical Molecular Dynamics with Mobile Protons

An important limitation of standard classical molecular dynamics simulations is the inability to make or break chemical bonds. This restricts severely our ability to study processes that involve even the simplest of chemical reactions, the transfer of a proton. Existing approaches for allowing proton transfer in the context of classical mechanics are rather cumbersome and have not achieved widespread use and routine status. Here we reconsider the combination of molecular dynamics with periodic stochastic proton hops. To ensure computational efficiency, we propose a non-Boltzmann acceptance criterion that is heuristically adjusted to maintain the correct or desirable thermodynamic equilibria between different protonation states and proton transfer rates. Parameters are proposed for hydronium, Asp, Glu, and His. The algorithm is implemented in the program CHARMM and tested on proton diffusion in bulk water and carbon nanotubes and on proton conductance in the gramicidin A channel. Using hopping parameters determined from proton diffusion in bulk water, the model reproduces the enhanced proton diffusivity in carbon nanotubes and gives a reasonable estimate of the proton conductance in gramicidin A.

Nanoconfinement in Slit Pores Enhances Water Self-Dissociation

We investigate the self-dissociation of water that is nanoconfined between the sheets of a realistic layered mineral, FeS mackinawite, as well as between Lennard-Jones walls via ab initio simulations. By comparing it with the same reaction in bulk water under various thermodynamic conditions, we show that such strong two-dimensional confinement between hard surfaces greatly enhances the self-dissociation process of water-thus increasing its ionic product K_{w} due to nanoconfinement. In addition to providing free energies, we analyze in detail the underlying dielectric properties in terms of dipole moment distributions, and thus the polarity of the liquid, as well as local polarization fluctuations as quantified by dielectric tensor profiles perpendicular to the lamella.

Nonequilibrium Self-Assembly of π-Conjugated Oligopeptides in Solution

Supramolecular assembly is a powerful method that can be used to generate materials with well-defined structures across multiple length scales. Supramolecular assemblies consisting of biopolymer-synthetic polymer subunits are specifically known to exhibit exceptional structural and functional diversity as well as programmable control of noncovalent interactions through hydrogen bonding in biopolymer subunits. Despite recent progress, there is a need to control and quantitatively understand assembly under nonequilibrium conditions. In this work, we study the nonequilibrium self-assembly of π-conjugated synthetic oligopeptides using a combination of experiments and analytical modeling. By isolating an aqueous peptide solution droplet within an immiscible organic layer, the rate of peptide assembly in the aqueous solution can be controlled by tuning the transport rate of acid that is used to trigger assembly. Using this approach, peptides are guided to assemble under reaction-dominated and diffusion-dominated conditions, with results showing a transition from a diffusion-limited reaction front to spatially homogeneous assembly as the transport rate of acid decreases. Interestingly, our results show that the morphology of self-assembled peptide fibers is controlled by the assembly kinetics such that increasingly homogeneous structures of self-assembled synthetic oligopeptides were generally obtained using slower rates of assembly. We further developed an analytical reaction-diffusion model to describe oligopeptide assembly, and experimental results are compared to the reaction-diffusion model across a range of parameters. Overall, this work highlights the importance of molecular self-assembly under nonequilibrium conditions, specifically showing that oligopeptide assembly is governed by a delicate balance between reaction kinetics and transport processes.

On the complex structural diffusion of proton holes in nanoconfined alkaline solutions within slit pores

The hydroxide anion OH⁻(aq) in homogeneous bulk water, that is, the solvated proton hole, is known to feature peculiar properties compared with excess protons solvated therein. In this work, it is disclosed that nanoconfinement of such alkaline aqueous solutions strongly affects the key structural and dynamical properties of OH⁻(aq) compared with the bulk limit. The combined effect of the preferred hypercoordinated solvation pattern of OH⁻(aq), its preferred perpendicular orientation relative to the confining surfaces, the pronounced layering of nanoconfined water and the topology of the hydrogen bond network required for proton hole transfer lead to major changes of the charge transport mechanism, in such a way that the proton hole migration mechanism depends exquisitely on the width of the confined space that hosts the water film. Moreover, the anionic Zundel complex, which is of transient nature in homogeneous bulk solutions, can be dynamically trapped as a shallow intermediate species by suitable nanoconfinement conditions.

Confined liquids can exhibit very different properties compared with the bulk. Here, the authors show that the migration mechanism of the hydroxide anion in water is greatly affected by bi-dimensional nano-confinement and strongly depends on the width of the confined space.

Nuclear Quantum Effects in H(+) and OH(-) Diffusion along Confined Water Wires

The diffusion of protons and hydroxide ions along water wires provides an efficient mechanism for charge transport that is exploited by biological membrane channels and shows promise for technological applications such as fuel cells. However, what is lacking for a better control and design of these systems is a thorough theoretical understanding of the diffusion process at the atomic scale. Here we focus on two aspects of this process that are often disregarded because of their high computational cost: the use of first-principles potential energy surfaces and the treatment of the nuclei as quantum particles. We consider proton and hydroxide ions in finite water wires using density functional theory augmented with an apolar cylindrical confining potential. We employ machine learning techniques to identify the charged species, thus obtaining an agnostic definition that takes explicitly into account the delocalization of the charge in the Grotthus-like mechanism. We include nuclear quantum effects (NQEs) through the thermostated ring polymer molecular dynamics method and model finite system size effects by considering Langevin dynamics on the potential of mean force of the charged species, allowing us to extract the same "universal" diffusion coefficient from simulations with different wire sizes. In the classical case, diffusion coefficients depend significantly on the potential energy surface, in particular on how dispersion forces modulate water-water distances. NQEs, however, make the diffusion less sensitive to the underlying potential and geometry of the wire.

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.

Hydrated Excess Protons Can Create Their Own Water Wires

Grotthuss shuttling of an excess proton charge defect through hydrogen bonded water networks has long been the focus of theoretical and experimental studies. In this work we show that there is a related process in which water molecules move ("shuttle") through a hydrated excess proton charge defect in order to wet the path ahead for subsequent proton charge migration. This process is illustrated through reactive molecular dynamics simulations of proton transport through a hydrophobic nanotube, which penetrates through a hydrophobic region. Surprisingly, before the proton enters the nanotube, it starts "shooting" water molecules into the otherwise dry space via Grotthuss shuttling, effectively creating its own water wire where none existed before. As the proton enters the nanotube (by 2-3 Å), it completes the solvation process, transitioning the nanotube to the fully wet state. By contrast, other monatomic cations (e.g., K(+)) have just the opposite effect, by blocking the wetting process and making the nanotube even drier. As the dry nanotube gradually becomes wet when the proton charge defect enters it, the free energy barrier of proton permeation through the tube via Grotthuss shuttling drops significantly. This finding suggests that an important wetting mechanism may influence proton translocation in biological systems, i.e., one in which protons "create" their own water structures (water "wires") in hydrophobic spaces (e.g., protein pores) before migrating through them. An existing water wire, e.g., one seen in an X-ray crystal structure or MD simulations without an explicit excess proton, is therefore not a requirement for protons to transport through hydrophobic spaces.