Tuning chemotactic and diffusiophoretic spreading via hydrodynamic flows

Diffusioosmotic flow reversals due to ion-ion electrostatic correlations

Existing theories of diffusioosmosis have neglected ion-ion electrostatic correlations, which are important in concentrated electrolytes. Here, we develop a mathematical model to numerically compute the diffusioosmotic mobilities of binary symmetric electrolytes across low to high concentrations in a charged parallel-plate channel. We use the modified Poisson equation to model the ion-ion electrostatic correlations and the Bikerman model to account for the finite size of ions. We report two key findings. First, ion-ion electrostatic correlations can cause a unique reversal in the direction of diffusioosmosis. Such a reversal is not captured by existing theories, occurs at ≈ 0.4 M for a monovalent electrolyte, and at a much lower concentration of ≈ 0.003 M for a divalent electrolyte in a channel with the same surface charge. This highlights that diffusioosmosis of a concentrated electrolyte can be qualitatively different from that of a dilute electrolyte, not just in its magnitude but also its direction. Second, we predict a separate diffusioosmotic flow reversal, which is not due to electrostatic correlations but the competition between the underlying chemiosmosis and electroosmosis. This reversal can be achieved by varying the magnitude of the channel surface charge without changing its sign. However, electrostatic correlations can radically change how this flow reversal depends on the channel surface charge and ion diffusivity between a concentrated and a dilute electrolyte. The mathematical model developed here can be used to design diffusioosmosis of dilute and concentrated electrolytes, which is central to applications such as species mixing and separation, enhanced oil recovery, and reverse electrodialysis.

Asymmetric rectified electric and concentration fields in multicomponent electrolytes with surface reactions

Recent experimental studies have utilized AC electric fields and electrochemical reactions in multicomponent electrolyte solutions to control colloidal assembly. However, theoretical investigations have thus far been limited to binary electrolytes and have overlooked the impact of electrochemical reactions. In this study, we address these limitations by analyzing a system with multicomponent electrolytes, while also relaxing the assumption of ideally blocking electrodes to capture the effect of surface electrochemical reactions. Through a regular perturbation analysis in the low-applied-potential regime, we solve the Poisson-Nernst-Planck equations and obtain effective equations for electrical potential and ion concentrations. By employing a combination of numerical and analytical calculations, our analysis reveals a significant finding: electrochemical reactions alone can generate asymmetric rectified electric fields (AREFs), i.e., time-averaged, long-range electric fields, even when the diffusivities of the ionic species are equal. This finding expands our understanding beyond the conventional notion that AREFs arise solely from diffusivity contrast. Furthermore, we demonstrate that AREFs induced by electrochemical reactions can be stronger than those resulting from asymmetric diffusivities. Additionally, we report the emergence of asymmetric rectified concentration fields (ARCFs), i.e., time-averaged, long-range concentration fields, which supports the electrodiffusiophoresis mechanism of colloidal assembly observed in experiments. We also derive analytical expressions for AREFs and ARCFs, emphasizing the role of imbalances in ionic strength and charge density, respectively, as the driving forces behind their formation. The results presented in this article advance the field of colloidal assembly and also have implications for improved understanding of electrolyte transport in electrochemical devices.

Unidirectional drying of a suspension of diffusiophoretic colloids under gravity

Recent experiments (K. Inoue and S. Inasawa, RSC Adv., 2020, 10, 15763-15768) and simulations (J.-B. Salmon and F. Doumenc, Phys. Rev. Fluids, 2020, 5, 024201) demonstrated the significant impact of gravity on unidirectional drying of a colloidal suspension. However, under gravity, the role of colloid transport induced by an electrolyte concentration gradient, a mechanism known as diffusiophoresis, is unexplored to date. In this work, we employ direct numerical simulations and develop a macrotransport theory to analyze the advective-diffusive transport of an electrolyte-colloid suspension in a unidirectional drying cell under the influence of gravity and diffusiophoresis. We report three key findings. First, drying a suspension of solute-attracted diffusiophoretic colloids causes the strongest phase separation and generates the thinnest colloidal layer compared to non-diffusiophoretic or solute-repelled colloids. Second, when colloids are strongly solute-repelled, diffusiophoresis prevents the formation of colloid concentration gradient and hence gravity has a negligible effect on colloidal layer formation. Third, our macrotransport theory predicts new scalings for the growth of the colloidal layer. The scalings match with direct numerical simulations and indicate that the colloidal layer produced by solute-repelled diffusiophoretic colloids could be an order of magnitude thicker compared to non-diffusiophoretic or solute-attracted colloids. Our results enable tailoring the separation of colloid-electrolyte suspensions by tuning the interactions between the solvent, electrolyte, and colloids under Earth's or microgravity, which is central to ground-based and in-space applications.

Diffusiophoresis of a spherical particle in porous media

Recent experiments by Doan et al. (Nano Lett., 2021, 21, 7625-7630) demonstrated and measured colloid diffusiophoresis in porous media but existing theories cannot predict the observed colloid motion. Here, using regular perturbation method, we develop a mathematical model that can predict the diffusiophoretic motion of a charged colloidal particle driven by a binary monovalent electrolyte concentration gradient in a porous medium. The porous medium is modeled as a Brinkman medium with a constant Darcy permeability. The linearized Poisson-Boltzmann equation is employed to model the equilibrium electric potential distribution that is driven out-of-equilibrium under diffusiophoresis. We report three key findings. First, we demonstrate that colloid diffusiophoresis could be drastically hindered in a porous medium due to the additional hydrodynamic drag compared to diffusiophoresis in a free electrolyte solution. Second, we show that the variation of the diffusiophoretic motion with respect to a change in the electrolyte concentration in a porous medium could be qualitatively different from that in a free electrolyte solution. Third, our results match quantitatively with experimental measurements, highlighting the predictive power of the present model. The mathematical model developed here could be employed to design diffusiophoretic colloid transport in porous media, which are central to applications such as nanoparticle drug delivery and enhanced oil recovery.

Two-dimensional diffusiophoretic colloidal banding: optimizing the spatial and temporal design of solute sinks and sources

Diffusiophoresis refers to the phenomenon where colloidal particles move in response to solute concentration gradients. Existing studies on diffusiophoresis, both experimental and theoretical, primarily focus on the movement of colloidal particles in response to one-dimensional solute gradients. In this work, we numerically investigate the impact of two-dimensional solute gradients on the distribution of colloidal particles, i.e., colloidal banding, induced via diffusiophoresis. The solute gradients are generated by spatially arranged sources and sinks that emit/absorb a time-dependent solute molar rate. First we study a dipole system, i.e., one source and one sink, and discover that interdipole diffusion and molar rate decay timescales dictate colloidal banding. At timescales shorter than the interdipole diffusion timescale, we observe a rapid enhancement in particle enrichment around the source due to repulsion from the sink. However, at timescales longer than the interdipole diffusion timescale, the source and sink screen each other, leading to a slower enhancement. If the solute molar rate decays at the timescale of interdipole diffusion, an optimal separation distance is obtained such that particle enrichment is maximized. We find that the partition coefficient of solute at the interface between the source and bulk strongly impacts the optimal separation distance. Surprisingly, the diffusivity ratio of solute in the source and bulk has a much weaker impact on the optimal dipole separation distance. We also examine an octupole configuration, i.e., four sinks and four sources, arranged in a circle, and demonstrate that the geometric arrangement that maximizes enrichment depends on the radius of the circle. If the radius of the circle is small, it is preferred to have sources and sinks arranged in an alternating fashion. However, if the radius of the circle is large, a consecutive arrangement of sources and sinks is optimal. Our numerical framework introduces a novel method for spatially and temporally designing the banded structure of colloidal particles in two dimensions using diffusiophoresis and opens up new avenues in a field that has primarily focused on one-dimensional solute gradients.