Massively Enhanced Electroosmotic Transport in Nanochannels Grafted with End-Charged Polyelectrolyte Brushes

All-atom molecular dynamics simulations of polymer and polyelectrolyte brushes

Densely grafted polymer and polyelectrolyte (PE) brushes, owing to their significant abilities to functionalize surfaces for a plethora of applications in sensing, diagnostics, current rectification, surface wettability modification, drug delivery, and oil recovery, have attracted significant attention over the past several decades. Unfortunately, most of the attention has primarily focused on understanding the properties of the grafted polymer and the PE chains with little attention devoted to studying the behavior of the brush-supported ions (counterions needed to screen the PE chains) and water molecules. Over the past few years, our group has been at the forefront of addressing this gap: we have employed all-atom molecular dynamics (MD) simulations for studying a wide variety of polymer and PE brush systems with specific attention to unraveling the properties and behavior of the brush-supported water molecules and ions. Our findings have revealed some of the most fascinating properties of such brush-supported ions and water molecules, including the most remarkable control of nanofluidic transport afforded by the specific ion and water responses induced by the PE brushes grafted on the inner walls of the nanochannel. This feature article aims to summarize some of our key contributions associated with such atomistic simulations of polymer and PE brushes and brush-supported water molecules and counterions.

Smart nanochannels: tailoring ion transport properties through variation in nanochannel geometry

This research explores ion transport behavior and functionality in a hybrid nanochannel that consists of two conical and cylindrical parts. The numerical investigation focuses on analyzing the length of each part in the nanochannel. The nanochannels are hybrid cavities embedded in a membrane, where the size of the conical part varies as equal to, larger than, or smaller than the cylindrical part. The nanochannel is coated with a polyelectrolyte layer that exhibits a dense charge density distribution. The charge density of the soft layer is described using the soft step distribution function. We study the electroosmotic flow, ionic current, rectification, and selectivity of the nanochannel versus bulk electrolyte concentration, the charge density of the polyelectrolyte layer, and decay length, while considering the effect of ionic partitioning. The steady-state Poisson-Nernst-Planck and Navier-Stokes equations are solved using the finite element method. The findings reveal that the nanochannel with a more extensive conical section demonstrates increased rectification, with the rectification factor rising from 1.4 to 2 at a bulk concentration of 100 mM. Additionally, the nanochannel with a longer cylindrical part exhibits improved selectivity under negative voltage conditions, while positive voltage introduces a different situation. The nanochannel with equal cylindrical and conical parts significantly affects conductivity by modifying the charge density in the soft layer, resulting in a 3.125-fold increase in conductivity under positive voltage when the charge density in the polyelectrolyte layer is raised from 25 to 100 mol m-3. This research focuses on creating intelligent nanochannels by controlling mass concentration, charge density, and collapse length, improving system performance, and optimizing properties. It also offers valuable insights into ion transport mechanisms in nanochannel systems, advancing our understanding in this field.

Effects of electrostatic coupling and surface polarization on polyelectrolyte brush structure

In this work, we perform molecular dynamics simulations to study a spherical polyelectrolyte brush. We explore the effects of surface polarization and electrostatic coupling on brush size and distribution of counterions. The method of image charges is considered to take into account surface polarization, considering a metallic, an unpolarizable, and a dielectric nano-core. It is observed that, for all cases, a moderate shrinking-swelling effect appears with an increase in the electrostatic coupling parameter. This effect occurs under high Manning ratios. The curves relating the average size of polyelectrolyte brush as a function of coupling show a minimum. The results show that the grafting density of polyelectrolytes on the nano-core surface plays an important role in the polarization effect. We consider a modified Poisson-Boltzmann theory to describe the counterion profiles around the brush in the case of unpolarizable nano-cores and weak electrostatic coupling.

Quantitative model for predicting the electroosmotic flow in dual-pole nanochannels

Developing and assessing nanofluidic systems is time-consuming and costly owing to the method's novelty; hence, modeling is essential to determine the optimal areas for implementation and to grasp its workings. In this work, we examined the influence of dual-pole surface and nanopore configuration on ion transfer simultaneously. To achieve this, the two trumpet and cigarette configuration were coated with a dual-pole soft surface so that the negative charge could be positioned in the nanopore's small aperture. Subsequently, the Poisson-Nernst-Planck and Navier-Stokes equations were simultaneously solved under steady-state circumstances using varied values physicochemical properties for the soft surface and electrolyte. The pore's selectivity was S Trumpet > S Cigarette ${S}_{{\rm{Trumpet}}} &gt; {S}_{{\rm{Cigarette}}}$ , and the rectification factor, on the other hand, was R f Cigarette < R f Trumpet ${R}_{{f}_{{\rm{Cigarette}}}} &lt; {R}_{{f}_{{\rm{Trumpet}}}}$ , when the overall concentration was very low. When the ion partitioning effect is taken into account, we clearly show that the rectifying variables for the cigarette configuration and the trumpet configuration can reach values of 45 and 49.2, when the charge density and mass concentration were 100 mol/m³ and 1 mM, respectively. By using dual-pole surfaces, the controllability of nanopores' rectifying behavior may be modified to produce superior separation performance.

Swelling and shrinking of two opposing polyelectrolyte brushes

Salt concentration and confinement effects affect the configuration of polyelectrolyte (PE) brushes due to electrostatic interactions. In this work, we develop a new theoretical model to analyze the electrostatics and swelling-shrinking behavior of two opposing PE brushes. By comparing three length scales, i.e., equilibrium brush height, separation distance, and Debye length, we obtain distinct scaling laws for brush height in different regimes. We provide explanations for the anomalous shrinkage of the PE brush with added salt reported in experiments and simulations, the applicability of the homogeneous brush assumption, and the confinement effect on the brush height. Our model can be used to shed light on the configuration and functionalities of PE-grafted interfaces, which play important roles in ion selective membranes and organism lubrication. We also anticipate that our method will be useful to understand the functionalities of other charged soft matter systems, such as hydrogel swelling and colloidal stability.

Structural and electrostatic properties between pH-responsive polyelectrolyte brushes studied by augmented strong stretching theory

In this paper, we study electrostatic and structural properties between pH-responsive polyelectrolyte brushes by using a strong stretching theory accounting for excluded volume interactions, the density of polyelectrolyte chargeable sites and the Born energy difference between the inside and outside of the brush layer.In a free energy framework, we obtain self-consistent field equations to determine electrostatic properties between two pH-responsive polyelectrolyte brushes. We elucidate that in the region between two pH-responsive polyelectrolyte brushes, electrostatic potential at the centerline and osmotic pressure increase not only with excluded volume interaction, but also with density of chargeable sites on a polyelectrolyte molecule. Importantly, we clarify that when two pH-responsive polyelectrolyte brushes approach each other, the brush thickness becomes short and that a large excluded volume interaction and a large density of chargeable sites yield the enhanced contract of polyelectrolyte brushes. In addition, we also demonstrate how the influence of such quantities as pH, the number of Kuhn monomers, the density of charged sites, the lateral separation between adjacent polyelectrolyte brushes, Kuhn length on the electrostatic and structural properties between the two polyelectrolyte brushes is affected by the exclusion volume interaction. Finally, we investigate the influence of Born energy difference on the thickness of polyelectrolyte brushes and the osmotic pressure between two pH-responsive polyelectrolyte brushes.

Non-monotonic variation of flow strength in nanochannels grafted with end-charged polyelectrolyte layers

The electrokinetic transport of fluids, also called the electroosmotic flow (EOF), in micro/nanoscale devices occurs in promising applications such as electrokinetic energy conversion (EKEC) systems. Recently, EKEC systems grafted with end-charged polyelectrolyte (PE) layers (PELs) have been reported to exhibit higher efficiencies than those of intrinsic systems. Understanding the interplay between the end-charged PELs and electrical double layers (EDLs) on the EOF is crucial for designing highly efficient EKEC systems. The interplay between the end-charged PELs and EDLs on the strength of the EOF (V₀) is studied by explicitly modeling the EOF through nanochannels grafted with end-charged PELs using atomic simulations. The variation of V₀ is examined for nanochannels grafted with PELs at various separations (d = 3.5–0.4 nm) to cover various conformations of PEs, inlcuding mushroom, semi-dilute brushes, and concentrated brushes. We find that V₀ follows a non-monotonic variation as d decreases and this is correlated with the conformation of the PEs. Specifically, as d decreases, V₀ decreases first in the mushroom regime (d = 3.5–2.0 nm), and then V₀ increases in the concentrated brush regime (d = 0.75–0.4 nm). Navigated by the continuum Navier–Stokes–Brinkman model, the above observations are rationalized by the competition between the driving effect from the spatial shift of ions in EDLs and the drag effect from PELs. The insights obtained in this work are important to guide the design of highly efficient EKEC systems by grafting end-charged PELs onto channel surfaces.

A non-monotonic variation of the strength of electroosmotic flow (V₀) is reported for the electroosmotic flow through nanochannels grafted with end-charged polyelectrolytes at various separation using atomic simulations.

Origin of Ultrahigh Rectification in Polyelectrolyte Bilayers Modified Conical Nanopores

The switching of "ON" and "OFF" states of an ionic diode is investigated by considering a conical nanopore partially functionalized two polyelectrolyte (PE) layers via layer-by-layer deposition. Through observing the inversion of its rectification behavior, we demonstrate the function of the PE bilayers in ionic transport regulation. The ionic diode exhibits an ultrahigh ion rectification at a low level of pH. In an aqueous NaCl solution at pH 2, for example, the ratio of the current at "ON" state and that at "OFF" state can be about 800 and 200 for 1 and 100 mM, respectively. This remarkable gating behavior can be explained by the anion-pump-induced ion accumulation in the neutral region as well as the depletion zone at the interface. Our results further demonstrate the possibility of achieving an ultrahigh rectification in an ionic diode having a unipolar-like configuration.

Friends, Foes, and Favorites: Relative Interactions Determine How Polymer Brushes Absorb Vapors of Binary Solvents

Polymer brushes can absorb vapors from the surrounding atmosphere, which is relevant for many applications such as in sensing and separation technologies. In this article, we report on the absorption of binary mixtures of solvent vapors (A and B) with a thermodynamic mean-field model and with grand-canonical molecular dynamics simulations. Both methods show that the vapor with the strongest vapor-polymer interaction is favored and absorbs preferentially. In addition, the absorption of one vapor (A) influences the absorption of another (B). If the A-B interaction is stronger than the interaction between vapor B and the polymers, the presence of vapor A in the brush can aid the absorption of B: the vapors absorb collaboratively as friends. In contrast, if the A-polymer interaction is stronger than the B-polymer interaction and the brush has reached its maximum sorption capacity, the presence of A can reduce the absorption of B: the vapors absorb competitively as foes.

Efficient electroosmotic mixing in a narrow-fluidic channel: the role of a patterned soft layer

We propose a novel and efficient mixing technique in a soft narrow-fluidic channel under the influence of electrical forcing. We show that a grafted polyelectrolyte layer (PEL) added as a patch to the channel wall modulates the electrical double layer (EDL) so that an applied electric field initiates a local electroosmotic flow (EOF) at the patched section. This EOF develops in the opposite direction to the primary pressure-driven flow. This localized EOF leads to the formation of Lamb vortices at the patched sections through the phenomenon of momentum exchange with the primary stream and promotes the mixing therein. Our study, consistent with the stream-function/vorticity approach, primarily focuses on the numerical analysis of the mixing phenomena. Through a quantitative description, we reveal the effect of different patterns on the underlying mixing phenomena in the convective mixing regime. We also discuss the impact of key parameters on the mixing efficiency, the onset of the recirculation zone, variation in the mixing length, and the shear-driven aggregation kinetics in soft matter systems. Finally, considering the practicability of the present problem, we unveil the values of several design parameters for which the mixing efficiency in the channel reaches the maximum.

Theoretical study on the massively augmented electro-osmotic water transport in polyelectrolyte brush functionalized nanoslits

We demonstrate that functionalizing nanoslits with pH-responsive polyelectrolyte brushes can lead to extremely fast electro-osmotic (EOS) water transport, where the maximum centreline velocity and the volume flow rate can be an order of magnitude larger than these quantities in identically charged brush-free nanochannels for a wide range of system parameters. Such an enhancement is most remarkable given that the brushes have been known to retard the transport by imparting additional drag on the fluid flow. We argue that this enhancement stems from the localization of the charge density of the brush-induced electric double layer (and, hence, the EOS body force) away from the nanochannel wall (or the location of the wall-induced drag force). This ensures a much larger impact of the EOS body force triggering such fast water transport. Finally, the calculated flux values for the present brush-grafted nanochannels are found to be significantly larger than those for a wide range of nanofluidic membranes and channels, suggesting that the brush functionalization can be considered as a mechanism for enabling such superfast nanofluidic transport.

Strong stretching theory for pH-responsive polyelectrolyte brushes in large salt concentrations

In this paper, we develop a theory for describing the thermodynamics, configuration, and electrostatics of strongly-stretched, pH-responsive polyelectrolyte (PE) brushes in the presence of large salt concentrations.

Counterion-Tunable Thermosensitivity of Strong Polyelectrolyte Brushes

In this work, poly(sodium styrene sulfonate) brushes have been employed as a precursor to prepare thermosensitive strong polyelectrolyte brushes (SPBs) through a counterion exchange strategy. The substitution of hydrophilic Na⁺ counterions by hydrophobic tetraalkylphosphonium counterions leads to a thermoresponsivity of the SPBs. The thermosensitive properties including hydration, stiffness, and surface water wettability of the SPBs can be modulated by the type of the tetraalkylphosphonium counterions. Nevertheless, the wet thickness of the SPBs with tetraalkylphosphonium counterions does not exhibit an obvious temperature dependency due to the high steric barrier in the crowded environment of SPBs generated by the large tetraalkylphosphonium counterions. The mixtures of small Na⁺ counterions and large tetraalkylphosphonium counterions are employed to realize the thermosensitive wet thickness without sacrificing other thermoresponsive properties of the SPBs because the mixed counterions can bring both a certain hydrophobicity and some free space to the brushes. This work opens up the opportunities available for the use of counterions to tune the thermosensitivity of SPBs.

Ionic current in nanochannels grafted with pH-responsive polyelectrolyte brushes modeled using augmented strong stretching theory

In this paper, we provide a theory to quantify the ionic current ( i ion ) in nanochannels grafted with pH-responsive polyelectrolyte (PE) brushes. We consider the PE brushes to be modeled by our recently proposed augmented strong stretching theory (SST) model that improves the existing SST models by incorporating the effects of excluded volume interactions and an extended mass action law. Use of such augmented SST for this problem implies that this is the first study on computing i ion in PE brush-grafted nanochannels accounting for the appropriate coupled configuration-electrostatic description of the PE brushes. i ion is obtained as functions of PE brush grafting density, medium pH and salt concentration ( c ∞ ), and the density of polyelectrolyte chargeable sites (PECS). For large c ∞ , i ion increases linearly with c ∞ (as for such c ∞ , i ion becomes independent of the PE charge and is dominated by the bulk mobility and number density of the electrolyte ions), whereas i ion is independent of c ∞ at small c ∞ (where the electric double layer electrostatics and the total number of ions in the system is dominated by the hydrogen ions). We further witness an enhancement of i ion for smaller pH and larger grafting density at low and moderate c ∞ , while there is little to no effect of the PECS density on the ionic current except for weakly grafted brushes at low c ∞ . We anticipate that this study will serve as a theoretical foundation for a large number of applications that are based on the brush-induced modification of the ionic current in a nanochannel.

Electrokinetic energy conversion in nanochannels grafted with pH-responsive polyelectrolyte brushes modelled using augmented strong stretching theory

In this paper, we develop a theory to quantify the electrokinetic energy conversion in electrolyte-filled nanochannels grafted with pH-responsive polyelectrolyte (PE) brushes. A pressure-driven flow drives the mobile electrolyte ions of the electric double layer (EDL) supported by the charged PE brushes leading to the generation of a streaming current, a streaming electric field and eventually an electrical energy. The salient feature of this study is that the brushes are described using our recently developed augmented Strong Stretching Theory (SST) model. In all the previous theoretical studies on liquid transport in PE-brush-grafted nanochannels, the brushes have either been assumed to be of constant height (independent of salt concentration or pH) or modelled using the Alexander-de-Gennes model that considers uniform monomer distribution along the brush height. Such simplifications have meant that the salt and the pH dependence of the brush height, the monomer distribution, and the resulting electrostatics have not been appropriately accounted for in the transport calculations. This paper addresses these limitations and provides a much more detailed description of the brushes while capturing the corresponding electrokinetic energy conversion. The results establish that the presence of the PE brushes ensures a localization of the average EDL charge density away from the grafting surface, thereby enabling the migration of the EDL ions with a larger background flow velocity; as a consequence, there is an enhancement of the streaming current, streaming electric field, and the resulting electrical energy generation under certain grafting densities of the PE brushes.

Revisiting the strong stretching theory for pH-responsive polyelectrolyte brushes: effects of consideration of excluded volume interactions and an expanded form of the mass action law

In this paper, we develop a theory to account for the effect of excluded volume (EV) interactions in the strong stretching theory (SST) based description of pH-responsive polyelectrolyte (PE) brushes. The existing studies have considered the PE brushes to be present in a θ-solvent and hence have neglected the EV interactions; however, such a consideration cannot describe the situations where the pH-responsive brushes are in a "good" solvent. Secondly, we consider a more expanded form of the mass action law, governing the pH-dependent ionization of the PE molecules, in the SST description of the PE brushes. This expanded form of the mass action law considers different values of γa³ (γ is the density of chargeable sites on the PE molecule and a is the PE Kuhn length) and therefore is an improvement over the existing SST models of PE brushes as well as other theories involving pH-responsive PE molecules that always consider γa³ = 1. Our results demonstrate that the EV effects enhance the brush height by inducing additional PE inter-segmental repulsion. Similarly, the consideration of the expanded form of the mass action law would lead to a reduced (enhanced) brush height for γa³ < 1 (γa³ > 1). We also quantify variables such as the monomer density distribution, the distribution of the ends of the PE brush, and the EDL electrostatic potential and explain their differences with respect to those obtained with no EV interactions or γa³ = 1.

Charge polarization, local electroneutrality breakdown and eddy formation due to electroosmosis in varying-section channels

We characterize the dynamics of an electrolyte embedded in a varying-section channel under the action of a constant external electrostatic field. By means of molecular dynamics simulations we determine the stationary density, charge and velocity profiles of the electrolyte. Our results show that when the Debye length is comparable to the width of the channel bottlenecks a concentration polarization along with two eddies sets inside the channel. Interestingly, upon increasing the external field, local electroneutrality breaks down and charge polarization sets leading to the onset of net dipolar field. This novel scenario, that cannot be captured by the standard approaches based on local electroneutrality, opens the route for the realization of novel micro and nano-fluidic devices.

Overlimiting current due to electro-diffusive amplification of the second Wien effect at a cation-anion bipolar membrane junction

Numerical simulations are presented for the transient and steady-state response of a model electrodiffusive cell with a bipolar ion-selective membrane under electric current. The model uses a continuum Poisson-Nernst-Planck theory including source terms to account for the catalytic second Wien effect between ionogenic groups in the membranes and resolves the Debye layers at interfaces. The resulting electric field at the membrane junction is increased by as much as four orders of magnitude in comparison to the field external to the membrane. This leads to a significant amplification of the second Wien effect, creating an increased ionic flux due to the catalytic decomposition of water. The effect also induces an exaltation effect wherein the salt ion flux undergoes a concomitant increase as well. The interplay of effects results in a unique over-limiting current mechanism due to concentration polarization internal, rather than external, to the membranes. In addition to the case of two equal but oppositely charged membranes under the standard simplifying assumption of equal ionic diffusivities, two variations on this model are studied. Asymmetric diffusivities, representative of the actual mobility difference in dissociated water ions, and the effect of the membrane charge density ratio were also considered. The latter elucidates an overlimiting current shift mechanism for DNA adsorption on anion-selective membranes proposed by Slouka et al. [Langmuir 29, 8275 (2013)]. The former provides more realistic picture of multi-ion transport and demonstrates a surprising steady-state effect due to the asymmetry in the diffusivity of hydroxide and hydronium.

Efficient electrochemomechanical energy conversion in nanochannels grafted with end-charged polyelectrolyte brushes at medium and high salt concentration

We develop a theory to study the generation of the streaming potential and the resulting electrochemomechanical energy conversion (ECMEC) in the presence of pressure-driven transport in nanochannels grafted with end-charged polyelectrolyte (PE) brushes. Our theory gives a thermodynamically self-consistent coupled description of the PE-brush and the electrostatics of the electric double layer (EDL) induced by the PE charges. The end-charged brushes localize the maximum EDL charge density away from the wall, thereby enabling a larger magnitude of pressure-driven transport to stream the ions downstream. This effect is retarded by the drag force imparted by the brushes as well as by the enhanced electroosmotic transport in a direction opposite to the pressure-driven transport. An interplay of these three issues leads to highly non-trivial electrohydrodynamic transport that eventually allows us to converge on appropriate properties of the brushes (e.g., grafting density and the number of monomers) that lead to the generation of a significantly larger streaming potential and a much improved efficiency of the ECMEC as compared to the brush-free nanochannels particularly at medium and high salt concentrations.

Softness Induced Enhancement in Net Throughput of Non-Linear Bio-Fluids in Nanofluidic Channel under EDL Phenomenon

In this article, we describe the electro-hydrodynamics of non-Newtonian fluid in narrow fluidic channel with solvent permeable and ion-penetrable polyelectrolyte layer (PEL) grafted on channel surface with an interaction of non-overlapping electric double layer (EDL) phenomenon. In this analysis, we integrate power-law model in the momentum equation for describing the non-Newtonian rheology. The complex interplay between the non-Newtonian rheology and interfacial electrochemistry in presence of PEL on the walls leads to non-intuitive variations in the underlying flow dynamics in the channels. As such, we bring out the variations in flow dynamics and their implications on the net throughput in the channel in terms of different parameters like power-law index (n), drag parameter (α), PEL thickness (d) and Debye length ratio (κ/κ_PEL) are discussed. We show, in this analysis, a relative enhancement in the net throughput through a soft nanofluidic channel for both the shear-thinning and shear-thickening fluids, attributed to the stronger electrical body forces stemming from ionic interactions between polyelectrolyte layer and electrolyte layer. Also, we illustrate that higher apparent viscosity inherent with the class of shear-thickening fluid weakens the softness induced enhancement in the volumetric flow rate for the shear-thickening fluids, since the viscous drag offered to the f low f ield becomes higher for the transport of shear-thickening fluid.

Polyelectrolyte brush bilayers in weak interpenetration regime: Scaling theory and molecular dynamics simulations

We employ molecular dynamics (MD) simulations and develop scaling theories to quantify the equilibrium behavior of polyelectrolyte (PE) brush bilayers (BBLs) in the weakly interpenetrated regime, which is characterized by d_{0}<d_{g}<2d_{0}, where d_{g} is the gap between the opposing plates where the PE brushes are grafted and d_{0} is the unperturbed height of a PE brush grafted at a single plate. Scaling predictions establish that, for the weakly interpenetrated osmotic PE BBLs δ∼N^{1/2}(2-d_{g}/d_{0})^{1/2} (where δ is the interpenetration length and N is the number of Kuhn segments in PE brush). MD simulations excellently recover this dependence of δ on N and the extent of interpenetration (quantified by d_{g}/d_{0}). These predictions, unlike the existing studies, establish a finite interpenetration for all values of d_{g}/d_{0} as long as d_{g}<2d_{0}. Finally, we quantify the monomer and counterion concentration distributions and point out that these two distributions may quantitatively deviate from each other at locations very close to the channel centerline, where the interpenetration-induced monomer concentration can be significantly low.

Highly enhanced liquid flows via thermoosmotic effects in soft and charged nanochannels

This paper proposes a massively augmented thermoosmotic transport in nanochannels grafted with end-charged polyelectrolyte brushes.