Self-generated electrokinetic flows from active-charged boundary patterns

We develop a hydrodynamic description of self-generated electrolyte flow in capillaries whose bounding walls feature both non-uniform distributions of charge and non-uniform active ionic fluxes. The hydrodynamic velocity arising in such a system has components that are forbidden by symmetry in the absence of charge and fluxes. However, when these two boundary mechanisms are simultaneously present, they can lead to a symmetry broken state where steady flows with both unidirectional and circulatory components emerge. We show that these flow states arise when modulated boundary patterns of charge and fluxes are offset by a flux-charge phase difference, which is associated with the separation between sites of their peak densities on the wall. Mismatch in diffusivity of cationic and anionic species can modify the flow states and becomes an enhancing factor when fluxes of both ion species are being produced together at the same site. We demonstrate that this mechanism can be realized with a microfluidic generator which is powered by enzyme-coated patches that catalyzes reactants in the solution to produce fluxes of ions. The local ionic elevation or depletion that disrupts a non-uniform double layer, promotes self-induced gradients yielding persistent body forces to generate bulk fluid motion. Our work quantifies a boundary-driven mechanism behind self-sustained electrolyte flow in confined environments that exists without any external bulk-imposed fields or gradients. It provides a theoretical framework for understanding the combined effect of active and charged boundaries that are relevant in biological or soft matter systems, and can be utilized in electrofluidic and iontronic applications.

💡 Research Summary

In this paper the authors develop a comprehensive theoretical framework for a novel class of electrokinetic flows that arise spontaneously in confined electrolyte systems when the bounding walls possess simultaneously non‑uniform surface charge and spatially varying active ionic fluxes. Unlike conventional electro‑osmotic or diffusio‑osmotic flows, which require externally imposed electric fields or bulk concentration gradients, the flows described here are generated entirely by the interplay of the two patterned boundary conditions, without any external bulk driving forces.

The study begins with a symmetry analysis of a long cylindrical capillary of radius R containing a symmetric 1:1 electrolyte. The surface charge is prescribed as σ(z)=σ₀ cos(kz+β) and the normal ionic fluxes of cations and anions as j±(z)=j₀± cos(kz+α±), where k=2π/l is the modulation wavenumber and β, α± are phase offsets. By examining how the axial component of the velocity field v_z transforms under spatial inversion (r,z)→(−r,−z) and k→−k, the authors show that only two types of terms are allowed in v_z when both charge and flux modulations are present: (i) a term proportional to the product of two flux amplitudes, γ_{m,m′} j₀m j₀m′, which is odd in k and gives rise to a circulatory flow consisting of alternating vortices along the channel, and (ii) a mixed term γ_{m,σ} j₀m σ₀, also odd in k, which produces a non‑zero spatial average (a “zero‑mode”) and drives a net unidirectional flow along the capillary axis. The coefficients γ are derived analytically later and depend on the Debye length λ, the channel geometry, and the wavenumber k.

To obtain quantitative predictions, the authors solve the coupled Poisson‑Nernst‑Planck‑Navier‑Stokes (PNP‑NS) equations under steady‑state conditions, assuming small amplitude perturbations of the ion concentrations and electric potential about the bulk values. By nondimensionalising with the Debye length and linearising the equations, they obtain analytic expressions for the leading‑order charge density perturbation δρ(r,z) and potential perturbation δϕ(r,z) in terms of modified Bessel functions I₀ and I₁. The boundary conditions at r=R enforce the prescribed surface charge (−ε ∂ϕ/∂r=σ(z)) and the prescribed normal ionic fluxes (J±·n=j±(z)). The resulting fields contain cosine and sine components with the same wavenumber k, and the relative phase between the charge and flux modulations, Δ=β−α, determines whether the electric body force e(c⁺−c⁻)E is in phase (producing a single‑sign axial force) or out of phase (producing alternating sign forces that generate vortex pairs).

The analytical results are corroborated by finite‑element simulations of the full nonlinear PNP‑NS system. The simulations illustrate three representative flow regimes: (a) pure circulatory flow when Δ≈π, where adjacent vortices of opposite rotation fill the channel; (b) dominant unidirectional flow when Δ≈π/2, where the axial velocity profile is essentially uniform and the net flux is maximal; and (c) mixed states for intermediate Δ (e.g., 3π/4), where both a net drift and a vortex lattice coexist. The magnitude of the net drift scales linearly with the product j₀σ₀ and inversely with the Debye length, reaching velocities of order 10–50 µm s⁻¹ for realistic parameters (R≈10–100 λ, l≈200 λ, C∞≈10⁻²–10² mM).

A further key insight concerns the role of unequal diffusivities of the two ionic species (D⁺≠D⁻). When the diffusivity mismatch is present, the symmetry‑allowed mixed term γ_{m,σ} acquires an additional contribution that can either amplify or reverse the net flow, depending on the sign of (D⁺−D⁻) and on whether both cations and anions are generated at the same location. The authors illustrate this by comparing two enzyme‑catalyzed reactions that produce different ion pairs with distinct diffusivities; the case with a larger diffusivity contrast yields a substantially higher flow speed, confirming that diffusivity mismatch is an enhancing factor.

To demonstrate practical feasibility, the paper proposes a microfluidic generator based on enzyme‑coated patches. A urease‑coated region on the wall catalyzes the hydrolysis of urea, releasing NH₄⁺ and OH⁻ locally, thereby creating a prescribed inward ionic flux. A separate patch with a fixed surface charge (e.g., silanized silica) is placed at a controlled distance, establishing a well‑defined phase shift Δ between charge and flux patterns. Numerical simulations show that varying the separation distance tunes both the magnitude and direction of the net flow, offering a simple design knob for micro‑pumps without external electrodes.

Overall, the work provides a rigorous theoretical foundation for a boundary‑driven, self‑sustained electrokinetic mechanism that can operate in biological contexts (e.g., coupled ion‑pump and channel activity on cell membranes) and in engineered iontronic devices (e.g., neuromorphic fluidic circuits, autonomous mixers, and pumps). By highlighting the essential role of the flux‑charge phase difference and the diffusivity mismatch, the authors open new avenues for designing soft‑matter systems that harness internal chemical activity to generate controlled fluid motion.