Molecular Friction-Induced Electroosmotic Phenomena in Thin Neutral Nanotubes

From ab initio to continuum: Linking multiple scales using deep-learned forces

We develop a deep learning-based algorithm, called DeepForce, to link ab initio physics with the continuum theory to predict concentration profiles of confined water. We show that the deep-learned forces can be used to predict the structural properties of water confined in a nanochannel with quantum scale accuracy by solving the continuum theory given by Nernst-Planck equation. The DeepForce model has an excellent predictive performance with a relative error less than 7.6% not only for confined water in small channel systems (L < 6 nm) but also for confined water in large channel systems (L = 20 nm) which are computationally inaccessible through the high accuracy ab initio molecular dynamics simulations. Finally, we note that classical Molecular dynamics simulations can be inaccurate in capturing the interfacial physics of water in confinement (L < 4.0 nm) when quantum scale physics are neglected.

Bioinspired Separator with Ion-Selective Nanochannels for Lithium Metal Batteries

The free transport of anions through commercial polyolefin separators used in lithium metal batteries (LMBs) gives rise to concentration polarization and rapid growth of lithium dendrites, leading to poor performance and short circuits. Here, a new poly(ethylene-co-acrylic acid) (EAA) separator with functional active sites (i.e., carboxyl groups) distributing along the pore surface was fabricated, forming bioinspired ion-conducting nanochannels within the separator. As the carboxyl groups effectively desolvated Li⁺ and immobilized anion, the as-prepared EAA separator selectively accelerated the transport of Li⁺ with transference number of Li⁺ (t_Li⁺) up to 0.67, which was further confirmed by molecular dynamics simulations. The battery with the EAA separator can be stably cycled over 500 h at 5 mA cm^-2. The LMBs with the EAA separator have exceptional electrochemical performance of 107 mAh g^-1 at 5 C and a capacity retention of 69% after 200 cycles. This work provides new commercializable separators toward dendrite-free LMBs.

Field-Induced Hydration Shell Reorganization Enables Electro-osmotic Flow in Nanochannels

Electro-osmotic flow is the motion of fluid driven by an applied electric field, for which an electric double layer near a charged surface is deemed essential. Here, we find that electro-osmotic flow can occur in electrically neutral nanochannels in the absence of definable electric double layers through extensive molecular dynamics simulations. An applied electric field is shown to cause an intrinsic channel selectivity between cations and anions, by reorienting the hydration shells of these confined ions. The ion selectivity then results in a net charge density in the channel that induces the unconventional electro-osmotic flow. The flow direction is amenable to manipulation by the field strength and the channel size, which will inform ongoing efforts to develop highly integrated nanofluidic systems capable of complex flow control.

Promoting Electroosmotic Water Flow through a Carbon Nanotube by Weakening the Competition between Cations and Anions in a Lateral Electric Field

Understanding the electroosmotic flow through a nanochannel is essential to the design of novel nanofluidic devices, ranging from desalination to nanometer water pumps. Nonetheless, the competition between cation and anion in electric fields inevitably leads to a limited pumping of water, and thus weakening their competition could be a new avenue for the fundamental control of water transport. In this work, through a series of molecular dynamics simulations, we find a surprising phenomenon in which under the drive of a traditional longitudinal electric field, an additional lateral electric field can significantly weaken the competitive transport of a cation and anion through a carbon nanotube, which spontaneously leads to a massive increase in electroosmotic water flux. Specifically, with the increase in the lateral electric field, the anion flux exhibits an almost linear reduction, and the cation flux is stable and can even be enhanced. As a result, the net water flux along the cation direction increases significantly. The key to this unexpected phenomenon lies in the size and mobility difference between the cation and anion. The anion is larger and has greater mobility and is thus more susceptible to the lateral electric field, which ultimately leads to the reduction of its flux. For different ion types and CNT lengths, we can observe similar electropumping phenomenon, where the friction force induced by the lateral electric field becomes nontrivial for long CNTs. Our results provide a new route to tune the competitive transport of cations and anions and should be useful for the design of novel electroosmotic pumps.

Electropumping Phenomenon in Modified Carbon Nanotubes

Controlling the water transport in a given direction is essential to the design of novel nanofluidic devices, which is still a challenge because of thermal fluctuations on the nanoscale. In this work, we find an interesting electropumping phenomenon for charge-modified carbon nanotubes (CNTs) through a series of molecular dynamics simulations. In electric fields, the flowing counterions on the CNT inner surface provide a direct driving force for water conduction. Specifically, the dynamics of cations and anions exhibit distinct behaviors that lead to thoroughly different water dynamics in positively and negatively charged CNTs. Because of the competition between the increased ion number and ion-CNT interaction, the cation flux displays an interesting maximum behavior with the increase in surface charge density; however, the anion flux rises further at higher charge density because it is less attractive to the surface. Thus, the anion flux is always several times larger than cation flux that induces a higher water flux in positive CNTs with nearly 100% pumping efficiency, which highly exceeds the efficiency of pristine CNTs. With the change in charge density, the translocation time, occupancy number, and radial density profiles for water and ions also demonstrate a nontrivial difference for positive and negative CNTs. Furthermore, the ion flux exhibits an excellent linear relationship with the field strength, leading to the same water flux behavior. For the change in salt concentration, the pumping efficiency for positive CNTs is also nearly 100%. Our results provide significant new insight into the ionic transport through modified CNTs and should be helpful for the design of nanometer water pumps.

Effect of temperature on the coupling transport of water and ions through a carbon nanotube in an electric field

Temperature governs the motion of molecules at the nanoscale and thus should play an essential role in determining the transport of water and ions through a nanochannel, which is still poorly understood. This work devotes to revealing the temperature effect on the coupling transport of water and ions through a carbon nanotube by molecular dynamics simulations. A fascinating finding is that the ion flux order changes from cation > anion to anion > cation with the increase in field strength, leading to the same direction change of water flux. The competition between ion hydration strength and mobility should be a partial reason for this ion flux order transition. High temperatures significantly promote the transport of water and ions, stabilize the water flux direction, and enhance the critical field strength. The ion translocation time exhibits an excellent Arrhenius relation with the temperature and a power law relation with the field strength, yielding to the Langevin dynamics. However, because of self-diffusion, the water translocation time displays different behaviors without following the ions. The high temperature also leads to an abnormal maximum behavior of the ion flux, deciphered by the massive increase in water flow that inversely hinders the ion flux, suggesting the coexistence of water-ion coupling transport and competition. Our results shed deep light on the temperature dependence of coupling transport of water and ions, answering a fundamental question on the water flux direction during the ionic transport, and thus should have great implications in the design of high flux nanofluidic devices.

Surfactant-Controlled Mobility of Oil Droplets in Mineral Nanopores

Polymer flooding is one of the widely used enhanced oil recovery (EOR) methods. However, tuning polymer properties to achieve improved performance in porous mineral rocks of diverse oil reservoirs remains one of the challenges of EOR processes. Here, we use molecular dynamics (MD) simulations to examine decane/water mixtures with surfactant additives in calcite and kaolinite mineral nanopores and characterize surfactant properties associated with increased fluid mobility and improved wettability in planar and constricted nanopore geometries. Cetyltrimethylammonium chloride (CTAC) and sodium dodecyl sulfate (SDS) surfactants are found to modulate the contact angles of decane droplets and reduce the decane density on mineral surfaces. CTAC can enhance and unblock the flow of decane droplets through narrowing nanopores with constricted geometries while aiding in decane droplet shape deformation, whereas SDS leads to decane droplets stalling in front of constrictions in nanopores. We hypothesize that the inability of the cationic CTAC headgroup to form hydrogen bonds is one of the key factors leading to enhanced CTAC-coated decane flow through constricted nanopores. The obtained molecular view of equilibrium and dynamic properties of complex fluids typical of oil reservoirs can provide a basis for the future design of new molecules for EOR processes.

Coupling effects in electromechanical ion transport in graphene nanochannels

In this work, we use molecular dynamics simulations to study the transport of ions in electromechanical flows in slit-like graphene nanochannels. The variation of ionic currents indicates a nonlinear coupling between pressure-driven and electroosmotic flows, which enhances the ionic currents for electromechanical flows compared with the linear superposition of pressure-driven and electroosmotic flows. The nonlinear coupling is attributed to the reduction of the total potential energy barrier due to the density variations of ions and water molecules in the channel. The numerical results may offer molecular insights into the design of nanofluidic devices for energy conversion.

Phonon-Induced Ratchet Motion of a Water Nanodroplet on a Supported Black Phosphorene

Phonons are not supposed to carry any physical momentum as lattice vibrational modes; thus, it is believed no mass transport could be induced by phonons. In this Letter, we show that a ratchet motion of a water nanodroplet could be induced on a two-dimensional puckered lattice like black phosphorene (BP) by exciting its flexural phonons through a moving substrate. The water nanodroplet exhibits a forward motion along the armchair or a backward motion along the zigzag directions on a BP lattice that is supported on a substrate possessing a relative armchair or zigzag forward motion with BP. Through the analysis of the structure and vibrational density states of BP, it is found that in-plane lattice displacement asymmetry and the in-plane vibration asymmetry are induced by the excited flexural phonons, which determine the water nanodroplet motion as an anisotropic Brownian motor. Simulations of the nanodroplet motion as functions of the substrate relative motion speed and direction and also the substrate coupling strength with BP are performed. Results of the nanodroplet ratchet motion exhibit good agreement with the theoretical predications from calculating the Brownian motor asymmetry. Our findings reveal a promising mass transport strategy and a further understanding of phonon-related interactions in crystalline solids.

Unexpected large impact of small charges on surface frictions with similar wetting properties

Generally, the interface friction on solid surfaces is regarded as consistent with wetting behaviors, characterized by the contact angles. Here using molecular dynamics simulations, we find that even a small charge difference (≤0.36 e) causes a change in the friction coefficient of over an order of magnitude on two-dimensional material and lipid surfaces, despite similar contact angles. This large difference is confirmed by experimentally measuring interfacial friction of graphite and MoS₂ contacting on water, using atomic force microscopy. The large variation in the friction coefficient is attributed to the different fluctuations of localized potential energy under inhomogeneous charge distribution. Our results help to understand the dynamics of two-dimensional materials and biomolecules, generally formed by atoms with small charge, including nanomaterials, such as nitrogen-doped graphene, hydrogen-terminated graphene, or MoS₂, and molecular transport through cell membranes.

Specific ion effects at graphitic interfaces

Improved understanding of aqueous solutions at graphitic interfaces is critical for energy storage and water desalination. However, many mechanistic details remain unclear, including how interfacial structure and response are dictated by intrinsic properties of solvated ions under applied voltage. In this work, we combine hybrid first-principles/continuum simulations with electrochemical measurements to investigate adsorption of several alkali-metal cations at the interface with graphene and within graphene slit-pores. We confirm that adsorption energy increases with ionic radius, while being highly dependent on the pore size. In addition, in contrast with conventional electrochemical models, we find that interfacial charge transfer contributes non-negligibly to this interaction and can be further enhanced by confinement. We conclude that the measured interfacial capacitance trends result from a complex interplay between voltage, confinement, and specific ion effects-including ion hydration and charge transfer.

Understanding aqueous solutions at graphitic interfaces is critical in a wide variety of emerging technologies. Here, the authors unravel specific ion effects at the interface with graphene and within graphene slit-pores by coupling first-principles simulations and electrochemical measurements.

The Role of Interface Ions in the Control of Water Transport through a Carbon Nanotube

Controlling the water transport toward a given direction is still challenging, particularly due to thermal fluctuations of water motion at the nanoscale. While most of the previous works focus on the symmetric hydrophobic membrane systems, the role of the membrane in affecting the water transport remains largely unexplored. In this work, by using extensive molecular dynamics simulations, we find an interesting electropumping phenomenon, that is, the flowing counterions on an asymmetric hydrophobic-hydrophilic membrane can significantly drive the single-file water transport through a carbon nanotube, suggesting a nanometer water pump in a highly controllable fashion. The ion-water coupling motion in electric fields on the charged surface provides an indirect driving force for this pumping phenomenon. The water dynamics and thermal dynamics demonstrate a unique behavior with the change in electric fields, surface charge density, and even charge species. Particularly, due to the ion flux bifurcation for the positive and negative surfaces, the water dynamics such as the water flow, flux, and translocation time also exhibit similar asymmetry. Surprisingly, the positive surface charge induces an abnormal three-peak dipole distribution for the confined water and subsequent high flipping frequency. This can be attributed to the competition between the surface charge and interface water orientation on it. Our results indicate a new strategy to pump water through a nanochannel, making use of the counterion flowing on an asymmetric charged membrane, which are promising for future studies.

How ions block the single-file water transport through a carbon nanotube

Understanding the blockage of ions for water transport through nanochannels is crucial for the design of desalination nanofluidic devices. In this work, we systematically clarify how ions block the single-file water transport through a (6,6) carbon nanotube (CNT) by using molecular dynamics simulations. We consider various pressure differences and salt concentrations. With the increase of pressure difference, the water flux shows a linear growth that coincides with the Hagen-Poiseuille equation. Interestingly, the dependence of the CNT-ion interaction on the salt concentration results in a distinct ion blockage effect that ultimately leads to water flux bifurcation. The water translocation time shows a power law decay with pressure, depending on the salt concentration. Furthermore, with the increase of salt concentration, the water flux shows a linear decay with a larger slope for higher pressure, while the water translocation time shows an opposite behavior. Therefore, the ions can not only block the water entering but also slow down the water motion inside the CNT. Notably, the probability of cations and anions appearing at the CNT entrance is quite similar, suggesting a similar blockage effect; however, anions show deeper interactions with the CNT because of their larger size. We finally find a unique linear relation between the water flux and occupancy divided by the translocation time. Our results provide insightful information on the ion blockage effect for the single-file water transport, and are thus helpful for the design of novel filtration membranes.

Enantioselective Molecular Transport in Multilayer Graphene Nanopores

Multilayer superstructures based on stacked layered nanomaterials offer the possibility to design three-dimensional (3D) nanopores with highly specific properties analogous to protein channels. In a layer-by-layer design and stacking, analogous to molecular printing, superstructures with lock-and-key molecular nesting and transport characteristics could be prepared. To examine this possibility, we use molecular dynamics simulations to study electric field-driven transport of ions through stacked porous graphene flakes. First, highly selective, tunable, and correlated passage rates of monovalent atomic ions through these superstructures are observed in dependence on the ion type, nanopore type, and relative position and dynamics of neighboring porous flakes. Next, enantioselective molecular transport of ionized l- and d-leucine is observed in graphene stacks with helical nanopores. The outlined approach provides a general scheme for synthesis of functional 3D superstructures.

Asymmetric osmotic water permeation through a vesicle membrane

Understanding the water permeation through a cell membrane is of primary importance for biological activities and a key step to capture its shape transformation in salt solution. In this work, we reveal the dynamical behaviors of osmotically driven transport of water molecules across a vesicle membrane by molecular dynamics simulations. Of particular interest is that the water transport in and out of vesicles is highly distinguishable given the osmotic force are the same, suggesting an asymmetric osmotic transportation. This asymmetric phenomenon exists in a broad range of parameter space such as the salt concentration, temperature, and vesicle size and can be ascribed to the similar asymmetric potential energy of lipid-ion, lipid-water, lipid-solution, lipid-lipid, and the lipid-lipid energy fluctuation. Specifically, the water flux has a linear increase with the salt concentration, similar to the prediction by Nernst-Planck equation or Fick's first law. Furthermore, due to the Arrhenius relation between the membrane permeability and temperature, the water flux also exhibits excellent Arrhenius dependence on the temperature. Meanwhile, the water flux shows a linear increase with the vesicle surface area since the flux amount across a unit membrane area should be a constant. Finally, we also present the anonymous diffusion behaviors for the vesicle itself, where transitions from normal diffusion at short times to subdiffusion at long times are identified. Our results provide significant new physical insights for the osmotic water permeation through a vesicle membrane and are helpful for future experimental studies.

Correlated Rectification Transport in Ultranarrow Charged Nanocones

Using molecular dynamics simulations, we reveal ion rectification in charged nanocones with exit diameters of 1-2 nm. The simulations exhibit an opposite rectification current direction than experiments performed in conical channels with exit diameters larger than 5 nm. This can be understood by the fact that in ultranarrow charged cones screening ions are trapped close to the cone tip at both field directions, which necessitates them to be released from the cone in a correlated multi-ion fashion. Electroosmosis induced by a unidirectional ion flow is also observed.

Role of Electroosmosis in the Permeation of Neutral Molecules: CymA and Cyclodextrin as an Example

To quantify the flow of small uncharged molecules into and across nanopores, one often uses ion currents. The respective ion-current fluctuations caused by the presence of the analyte make it possible to draw some conclusions about the direction and magnitude of the analyte flow. However, often this flow appears to be asymmetric with respect to the applied voltage. As a possible reason for this asymmetry, we identified the electroosmotic flow (EOF), which is the water transport associated with ions driven by the external transmembrane voltage. As an example, we quantify the contribution of the EOF through a nanopore by investigating the permeation of α-cyclodextrin through CymA, a cyclodextrin-specific channel from Klebsiella oxytoca. To understand the results from electrophysiology on a molecular level, all-atom molecular dynamics simulations are used to detail the effect of the EOF on substrate entry to and exit from a CymA channel in which the N-terminus has been deleted. The combined experimental and computational results strongly suggest that one needs to account for the significant contribution of the EOF when analyzing the penetration of cyclodextrins through the CymA pore. This example study at the same time points to the more general finding that the EOF needs to be considered in translocation studies of neutral molecules and, at least in many cases, should be able to help in discriminating between translocation and binding events.

Transport behavior of a single Ca(2+), K(+), and Na(+) in a water-filled transmembrane cyclic peptide nanotube

Molecular dynamics simulations have been performed to investigate the transport properties of a single Ca(2+), K(+), and Na(+) in a water-filled transmembrane cyclic peptide nanotube (CPNT). Two transmembrane CPNTs, i.e., 8×(WL)n=4,5/POPE (with uniform lengths but various radii), were applied to clarify the dependence of ionic transport properties on the channel radius. A huge energy barrier keeps Ca(2+) out of the octa-CPNT, while Na(+) and K(+) can be trapped in two CPNTs. The dominant electrostatic interaction of a cation with water molecules leads to a high distribution of channel water around the cation and D-defects in the first and last gaps, and significantly reduces the axial diffusion of channel water. Water-bridged interactions were mostly found between the artificially introduced Ca(2+) and the framework of the octa-CPNT, and direct coordinations with the tube wall mostly occur for K(+) in the octa-CPNT. A cation may drift rapidly or behave lazily in a CPNT. K(+) behaves most actively and can visit the whole deca-CPNT quickly. The first solvation shells of Ca(2+) and Na(+) are basically saturated in two CPNTs, while the hydration of K(+) is incomplete in the octa-CPNT. The solvation structure of Ca(2+) in the octa-CPNT is most stable, while that of K(+) in the deca-CPNT is most labile. Increasing the channel radius induces numerous interchange attempts between the first-shell water molecules of a cation and the ones in the outer region, especially for the K(+) system.