Ion Transport in Electrically Imperfect Nanopores

Scaling of ionic conductance in a fluctuating single-layer nanoporous membrane

Single-layer membranes have emerged as promising candidates for applications requiring high transport rates due to their low resistance to molecular transport. Owing to their atomically thin structure, these membranes experience significant microscopic fluctuations, emphasizing the need to explore their impact on ion transport processes. In this study, we investigate the effects of membrane fluctuations on the elementary scaling behavior of ion conductance [Formula: see text] as a function of ion concentration [Formula: see text], represented as [Formula: see text], using molecular dynamics simulations. Our findings reveal that membrane fluctuations not only alter the conductance coefficient [Formula: see text] but also the power-law exponent [Formula: see text]. We identify two distinct frequency regimes of membrane fluctuations, GHz-scale and THz-scale fluctuations, and examine their roles in conductance scaling. Furthermore, we demonstrate that the alteration of conductance scaling arises from the non-linearity between ion conductance and membrane shape. This work provides a fundamental understanding of ion transport in fluctuating membranes.

First-Principles Prediction of Amorphous Silica Nanoparticle Surface Charge: Effect of Size, pH, and Ionic Strength

The quantification of surface charge properties of silica nanoparticles is essential for several applications. To determine these properties, many experimental and theoretical methods have been introduced, which are time-consuming and/or challenging to use. In this study, a first-principles approach is developed to determine the surface charge properties of amorphous silica nanoparticles against the nanoparticle size, pH, and ionic strength without relying on experimental data. An amorphous silica nanoparticle of 1.34 nm diameter is simulated by using integrated molecular dynamics and Monte Carlo methods. A detailed analysis of the nanoparticle structure is provided by analyzing the types of silanol groups on the surface. Moreover, a model is developed to estimate the probability distribution of the surface silanol groups based on the nearest neighbor distances and the diameter of the nanoparticle to determine the number of surface silanols on larger nanoparticles. Thereafter, a computational chemistry approach is used to calculate the acid dissociation constants of the corresponding deprotonation reactions. The calculated constants and the point of zero charge value are in excellent agreement with experiments. The surface charge properties of the nanoparticle with various diameters are then estimated by using a mean-field model at different pH and ionic strength values. The results of the developed model are compared to the Poisson-Boltzmann equation as a reference model. The developed model predictions agree well with the reference model for low and mid-electrolyte concentrations (1 and 10 mM) and small nanoparticles (smaller than 100 nm). However, the developed model seems to qualitatively predict the surface charge properties more accurately than the Poisson-Boltzmann model for high electrolyte concentrations.

Ion transport in two-dimensional flexible nanoporous membranes

Ion transport is a fundamental mechanism in living systems that plays a role in cell proliferation, energy conversion, and maintaining homeostasis. This has inspired various nanofluidic applications such as electricity harvesting, molecular sensors, and molecular separation. Two dimensional (2D) nanoporous membranes are particularly promising for these applications due to their ultralow transport barriers. We investigated ion conduction across flexible 2D membranes via extensive molecular dynamics simulations. We found that the microscopic fluctuations of these membranes can significantly increase ion conductance, for example, by 320% in Cu-HAB with 0.5 M KCl. Our analysis of ion dynamics near the flexible membranes revealed that ion hydration is destabilized when the membrane fluctuated within a specific frequency range leading to improved ion conduction. Our results show that the dynamic coupling between the fluctuating membrane and ions can play a crucial role in ion conduction across 2D nanoporous membranes.

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.

Characterization of the surface charge property and porosity of track-etched polymer membranes

As an important property of porous membranes, the surface charge property determines many ionic behaviors of nanopores, such as ionic conductance and selectivity. Based on the dependence of electric double layers on bulk concentrations, ionic conductance through nanopores at high and low concentrations is governed by the bulk conductance and surface charge density, respectively. Here, through the investigation of ionic conductance inside track-etched single polyethylene terephthalate (PET) nanopores under various concentrations, the surface charge density of PET membranes is extracted as ∼-0.021 C/m² at pH 10 over measurements with 40 PET nanopores. Simulations show that surface roughness can cause underestimation in surface charge density due to the inhibited electroosmotic flow. Then, the averaged pore size and porosity of track-etched multipore PET membranes are characterized by the developed ionic conductance method. Through coupled theoretical predictions in ionic conductance under high and low concentrations, the averaged pore size and porosity of porous membranes can be obtained simultaneously. Our method provides a simple and precise way to characterize the pore size and porosity of multipore membranes, especially for those with sub-100 nm pores and low porosities.

Effective Modulation of Ion Mobility through Solid-State Single-Digit Nanopores

Many experimental studies have proved that ion dynamics in a single-digit nanopore with dimensions comparable to the Debye length deviate from the bulk values, but we still have critical knowledge gaps in our understanding of ion transport in nanoconfinement. For many energy devices and sensor designs of nanoporous materials, ion mobility is a key parameter for the performance of nanofluidic equipment. However, investigating ion mobility remains an experimental challenge. This study experimentally investigated the monovalent ion dynamics of single-digit nanopores from the perspective of ionic conductance. In this article, we present a theory that is sufficient for a basic understanding of ion transport through a single-digit nanopore, and we subdivided and separately analyzed the contribution of each conductance component. These conclusions will be useful not only in understanding the behavior of ion migration but also in the design of high-performance nanofluidic devices.

Electrical Conductance of Charged Nanopores

A nanopore's response to an electrical potential drop is characterized by its electrical conductance, . For the last two decades, it has been thought that at low electrolyte concentrations, , the conductance is concentration-independent such that . It has been recently demonstrated that surface charge regulation changes the dependency to , whereby the slope typically takes the values α = 1/3 or 1/2. However, experiments have observed slopes of 2/3 and 1 suggesting that additional mechanisms, such as convection and slip-lengths, appear. Here, we elucidate the interplay between three mechanisms: surface charge regulation, convection, and slip lengths. We show that the inclusion of convection does not change the slope, and when the effects of hydrodynamic slip are included, the slope is doubled. We show that when all effects are accounted for, α can take any value between 0 and 1 where the exact value of the slope depends on the material properties. This result is of utmost importance in designing any electro-kinetically driven nanofluidic system characterized by its conductance.

High-Performance Osmotic Power Generators Based on the 1D/2D Hybrid Nanochannel System

Extracting clean energy by converting the salinity gradient between river and sea into energy is an effective way to reduce the global pollution and carbon emissions. Reverse electrodialysis (RED) is of great importance to realize the energy conversion assisting the ion-selective membrane. However, its higher ion resistance and lower conversion efficiency results in the undesirable power conversion performance. Here, we demonstrate a 1D/2D hybrid nanochannel system to achieve high osmotic energy conversion and output power. This heterogeneous structure is composed of two structures, in which the subnanometer nanochannels in graphene oxide membrane (GOM) can serve as a selective layer and reduce the ion diffusion energy barrier, while the nanochannel in the polymer can introduce asymmetry to enhance ionic rectification and conversion efficiency. This heterogeneous membrane exhibits excellent cation selectivity and enhanced ionic current rectification (ICR) performance. The application of the GOM/PET hybrid nanochannel system in osmotic energy harvesting is evaluated, and the output power can reach up to 118.2 pW with the energy conversion efficiency of 40.3%. Theoretical calculation indicates that the 1D/2D hybrid system can effectively take the advantage of excellent cation selectivity of 2D lamellar nanochannels to improve its RED performance significantly.

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.

Ion Transport in pH-Regulated Double-Barreled Nanopores

Single-molecule detection and characterization with nanopores is a powerful technique that does not require labeling. Multinanopore systems, especially double nanopores, have attracted wide attention and have been applied in many fields. However, theoretical studies of electrokinetic ion transport in nanopores mainly focus on single nanopores. In this paper, for the first time, a theoretical study of pH-regulated double-barreled nanopores is conducted using three-dimensional Poisson-Nernst-Planck equations and Navier-Stokes equations. Four ionic species and the surface chemistry on the walls of the nanopores are included. The results demonstrate that the properties of the bulk salt solution significantly affect nanopore conductivity and ion transport phenomena in nanopores. There are two ion-enriched zones and two ion-depleted zones in double-barreled nanopores. Due to the symmetry of the double-barreled nanopore structure and surface charge density, there is no ionic rectification effect in double-barreled nanopores. The ion selectivity is similar to that of conventional single pH-regulated nanopores.

Nonlinear electrohydrodynamic ion transport in graphene nanopores

Mechanosensitivity is one of the essential functionalities of biological ion channels. Synthesizing an artificial nanofluidic system to mimic such sensations will not only improve our understanding of these fluidic systems but also inspire applications. In contrast to the electrohydrodynamic ion transport in long nanoslits and nanotubes, coupling hydrodynamical and ion transport at the single-atom thickness remains challenging. Here, we report the pressure-modulated ion conduction in graphene nanopores featuring nonlinear electrohydrodynamic coupling. Increase of ionic conductance, ranging from a few percent to 204.5% induced by the pressure—an effect that was not predicted by the classical linear coupling of molecular streaming to voltage-driven ion transport—was observed experimentally. Computational and theoretical studies reveal that the pressure sensitivity of graphene nanopores arises from the transport of capacitively accumulated ions near the graphene surface. Our findings may help understand the electrohydrodynamic ion transport in nanopores and offer a new ion transport controlling methodology.

Conditions for electroneutrality breakdown in nanopores

It has recently been suggested that a breakdown of electroneutrality occurs in highly confined nanopores that are encompassed by a dielectric material. This work elucidates the conditions for this breakdown. We show that the breakdown within the pore results from the response of the electric field within the dielectric. Namely, we show that this response is highly sensitive to the boundary condition at the dielectric edge. The standard Neumann boundary condition of no-flux predicts that the breakdown does not occur. However, a Dirichlet boundary condition for a zero-potential predicts a breakdown. Within this latter scenario, the breakdown exhibits a dependence on the thickness of the dielectric material. Specifically, infinite thickness dielectrics do not exhibit a breakdown, while dielectrics of finite thickness do exhibit a breakdown. Numerical simulations confirm theoretical predictions. The breakdown outcomes are discussed with regard to single pore systems and multiple pore systems.

Osmotic Transport at the Aqueous Graphene and hBN Interfaces: Scaling Laws from a Unified, First-Principles Description

Osmotic transport in nanoconfined aqueous electrolytes provides alternative venues for water desalination and "blue energy" harvesting. The osmotic response of nanofluidic systems is controlled by the interfacial structure of water and electrolyte solutions in the so-called electrical double layer (EDL), but a molecular-level picture of the EDL is to a large extent still lacking. Particularly, the role of the electronic structure has not been considered in the description of electrolyte/surface interactions. Here, we report enhanced sampling simulations based on ab initio molecular dynamics, aiming at unravelling the free energy of prototypical ions adsorbed at the aqueous graphene and hBN interfaces, and its consequences on nanofluidic osmotic transport. Specifically, we predicted the zeta potential, the diffusio-osmotic mobility, and the diffusio-osmotic conductivity for a wide range of salt concentrations from the ab initio water and ion spatial distributions through an analytical framework based on Stokes equation and a modified Poisson-Boltzmann equation. We observed concentration-dependent scaling laws, together with dramatic differences in osmotic transport between the two interfaces, including diffusio-osmotic flow and current reversal on hBN but not on graphene. We could rationalize the results for the three osmotic responses with a simple model based on characteristic length scales for ion and water adsorption at the surface, which are quite different on graphene and on hBN. Our work provides fundamental insights into the structure and osmotic transport of aqueous electrolytes on 2D materials and explores alternative pathways for efficient water desalination and osmotic energy conversion.

Surface Charge Density Inside a Silicon Nitride Nanopore

Surface charges inside a nanopore determine the zeta potential and ion distributions and play a significant role in affecting ion transport and the sensitivity of detecting biomolecules. It is of great importance to study the fluctuation of surface charges with the salt concentration and pH in various applications of nanopores. Herein, we proposed a theoretical model to predict the surface charge density of a Si₃N₄ nanopore, in which both silanol and amine groups were taken into account. It was demonstrated that the surface charge density in the Si₃N₄ nanopore changes not only with pH but also with the salt concentration. The theoretical model could well predict the experimental results with different salt concentrations, pH values, and pore sizes. The effect of surface functional groups on the isopotential point (pH_iep) of the Si₃N₄ nanopore was also systematically studied. The results indicated that the silanol groups are major determinants of the surface charge, but the influences of the amine groups should not be ignored because the small number of amine groups can change pH_iep dramatically. The pH_iep value of the Si₃N₄ nanopore was measured as 4.1, and the ratio of amine over silanol was ascertained as 0.013.

Electrostatic gating of ion transport in carbon nanotube porins: A modeling study

Carbon nanotube porins (CNTPs) are biomimetic membrane channels that demonstrate excellent biocompatibility and unique water and ion transport properties. Gating transport in CNTPs with external voltage could increase control over ion flow and selectivity. Herein, we used continuum modeling to probe the parameters that enable and further affect CNTP gating efficiency, including the size and composition of the supporting lipid membrane, slip flow in the carbon nanotube, and the intrinsic electronic properties of the nanotube. Our results show that the optimal gated CNTP device consists of a semiconducting CNTP inserted into a small membrane patch containing an internally conductive layer. Moreover, we demonstrate that the ionic transport modulated by gate voltages is controlled by the charge distribution along the CNTP under the external gate electric potential. The theoretical understanding developed in this study offers valuable guidance for the design of gated CNTP devices for nanofluidic studies, novel biomimetic membranes, and cellular interfaces in the future.

Bipolar Nanochannels: A Systematic Approach to Asymmetric Problems

Nanofluidic diodes are capable of rectifying the electrical current by several orders of magnitude. In the current state of affairs, determining the rectification factor is not possible as it depends on many system parameters. In this work, we systematically scan the effects of geometry and excess counterion concentrations (i.e., surface charge effects). We show that the current-voltage response varies between the two extreme behaviors of unipolar and bipolar responses. The exact behavior depends on the geometry and surface charge properties of the system. Here, we have gone beyond the typical setup that only considers the dynamics within the nanochannel itself and we have included the effects of the adjoining microchannels. Systems that include both nanochannels and microchannels exhibit the classical signatures of concentration polarization, such as ionic depletion and enrichment. Here, where we have scanned a wide range of parameters, we show that bipolar and semi-bipolar systems exhibit a wider range of phenomena that are intrinsically more complicated. Our system characterization is for both, the much more investigated case of steady state and the less investigated, but equally interesting, time-transient case. For example, it is common to characterize the system by its steady-state result (current-voltage response, rectification factor, and transport number). Here, we demonstrate that the time-transient behavior of the fluxes can also be used to characterize the system, and that the time-dependent rectification factors and transport numbers are meaningful. The systematic approach taken in this work, and the results presented herein, can be used to further elucidate the complicated behavior of the current-voltage response of nanofluidic diodes and to rationalize experimental results. The insights of this work can be used to enhance and improve the design of all nanofluidic diodes.

From nanotubes to nanoholes: Scaling of selectivity in uniformly charged nanopores through the Dukhin number for 1:1 electrolytes

Scaling of the behavior of a nanodevice means that the device function (selectivity) is a unique smooth and monotonic function of a scaling parameter that is an appropriate combination of the system's parameters. For the uniformly charged cylindrical nanopore studied here, these parameters are the electrolyte concentration, c, voltage, U, the radius and the length of the nanopore, R and H, and the surface charge density on the nanopore's surface, σ. Due to the non-linear dependence of selectivities on these parameters, scaling can only be applied in certain limits. We show that the Dukhin number, Du=|σ|/eRc∼|σ|λ_D ²/eR (λ_D is the Debye length), is an appropriate scaling parameter in the nanotube limit (H → ∞). Decreasing the length of the nanopore, namely, approaching the nanohole limit (H → 0), an alternative scaling parameter has been obtained, which contains the pore length and is called the modified Dukhin number: mDu ∼ Du H/λ_D ∼ |σ|λ_DH/eR. We found that the reason for non-linearity is that the double layers accumulating at the pore wall in the radial dimension correlate with the double layers accumulating at the entrances of the pore near the membrane on the two sides. Our modeling study using the Local Equilibrium Monte Carlo method and the Poisson-Nernst-Planck theory provides concentration, flux, and selectivity profiles that show whether the surface or the volume conduction dominates in a given region of the nanopore for a given combination of the variables. We propose that the inflection point of the scaling curve may be used to characterize the transition point between the surface and volume conductions.

Electric control of ionic transport in sub-nm nanopores

The ion transport behavior through sub-nm nanopores (length (L) ≈ radius (R)) on a film is different from that in nanochannels (L ≫ R), and even more different from the bulk behavior. The many intriguing phenomena in ionic transport are the key to the design and fabrication of solid-state nanofluidic devices. However, ion transport through sub-nm nanopores is not yet clearly understood. We investigate the ionic transport behavior of sub-nm nanopores from the perspective of conductance via molecular dynamics (MD) and experimental methods. Under the action of surface charge, the average ion concentration inside the nanopore is much higher than the bulk value. It is found that 100 mM is the transition point between the surface-charge-governed and the bulk behavior regimes, which is different from the transition point for nanochannels (10 mM). Moreover, by investigating the access, pores, surface charge, electroosmosis and potential leakage conductance, it is found that the conductive properties of the nanopore at low bulk concentration are determined by the surface charge potential leaks into the reservoir. Specifically, there is a huge increase in cation mobility through a cylindrical nanopore, which implies potential applications for the fast charging of supercapacitors and batteries. Sub-nm nanopores also show a strong selectivity toward Na⁺, and a strong repellence toward Cl⁻. These conclusions presented here will be useful not only in understanding the behavior of ion transport, but also in the design of nanofluidic devices.

The ion transport behavior through sub-nm nanopores (length (L) ≈ radius (R)) on a film is different from that in nanochannels (L ≫ R), and even more different from the bulk behavior.