An official website of the United States government
The role of particle shape in evaporation-induced auto-stratification in polydisperse colloidal suspensions is explored with molecular dynamics simulations of mixtures of spheres and aspherical particles. A unified framework based on the competition between diffusion and diffusiophoresis is proposed to understand the effects of shape and size dispersity. In general, particles diffusing more slowly (e.g., larger particles) tend to accumulate more strongly at the evaporation front. However, larger particles have larger surface areas and therefore greater diffusiophoretic mobility. Hence, they are more likely to be driven away from the evaporation front via diffusiophoresis. For a rapidly dried bidisperse suspension containing small and large spheres, the competition leads to “small-on-top” stratification. Here, we employ a computational model in which the diffusion coefficient is inversely proportional to particle mass. For a mixture of spheres and aspherical particles with similar mass, the diffusion contrast is reduced, and the spheres are always enriched at the evaporation front as they have the smallest surface area for a given mass and, therefore, the lowest diffusiophoretic mobility. For a mixture of solid and hollow spheres that have the same outer radius and thus the same surface area, the diffusiophoretic contrast is suppressed, and the system is dominated by diffusion. Consequently, the solid spheres, which have a larger mass and diffuse more slowly, accumulate on top of the hollow spheres. Finally, for a mixture of thin disks and long rods that differ significantly in shape but have similar mass and surface area, both diffusion and diffusiophoresis contrasts are suppressed, and the mixture does not stratify.
more »
« less
Diffusiophoresis in one-dimensional solute gradients
The diffusiophoretic motion of suspended colloidal particles under one-dimensional solute gradients is solved using numerical and analytical techniques. Similarity solutions are developed for the injection and withdrawal dynamics of particles into semi-infinite pores. Furthermore, a method of characteristics formulation of the diffusion-free particle transport model is presented and integrated to realize particle trajectories. Analytical solutions are presented for the limit of small particle diffusiophoretic mobility Γ p relative to the solute diffusivity D s for particle motions in both semi-infinite and finite domains. Results confirm the build up of local maxima and minima in the propagating particle front dynamics. The method of characteristics is shown to successfully predict particle motions and the position of the particle front, although it fails to accurately predict suspended particle concentrations in the vicinity of sharp gradients, such as at the particle front peak seen in some injection cases, where particle diffusion inevitably plays an important role. Results inform the design of applications in which the use of applied solute gradients can greatly enhance particle injection into and withdrawal from pores.
more »
« less
- Award ID(s):
- 1702693
- PAR ID:
- 10056144
- Date Published:
- Journal Name:
- Soft Matter
- Volume:
- 13
- Issue:
- 47
- ISSN:
- 1744-683X
- Page Range / eLocation ID:
- 9015 to 9023
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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.more » « less
-
Lipid vesicles immersed in solute gradients are predicted to migrate from regions of high to low solute concentration due to osmotic flows induced across their semipermeable membranes. This process─known as osmophoresis─is potentially relevant to biological processes such as vesicle trafficking and cell migration; however, there exist significant discrepancies (several orders of magnitude) between experimental observations and theoretical predictions for the vesicle speed. Here, we seek to reconcile predictions of osmophoresis with observations of vesicle motion in osmotic gradients. We prepare quasi-steady solute gradients in a microfluidic chamber using density-matched solutions of sucrose and glucose to eliminate buoyancy-driven flows. We quantify the motions of giant DLPC vesicles and Brownian tracer particles in such gradients using Bayesian analysis of particle tracking data. Despite efforts to mitigate convective flows, we observe directed motion of both lipid vesicles and tracer particles in a common direction at comparable speeds of order 10 nm/s. These observations are not inconsistent with models of osmophoresis, which predict slower motion at ca. 1 nm/s; however, experimental uncertainty and the confounding effects of fluid convection prohibit a quantitative comparison. In contrast to previous reports, we find no evidence for anomalously fast osmophoresis of lipid vesicles when fluid convection is mitigated and quantified. We discuss strategies for enhancing the speed of osmophoresis using high permeability membranes and geometric confinement.more » « less
-
For isotropic swimming particles driven by self-diffusiophoresis at zero Reynolds number (where particle velocity responds instantaneously to applied force), the diffusive timescale of emitted solute can produce an emergent quasi-inertial behavior. These particles can orbit in a central potential and reorient under second-order dynamics, not the first-order dynamics of classical zero-Reynolds motion. They are described by a simple effective model that embeds their history-dependent behavior as an effective inertia, this being the most primitive expression of memory. The system can be parameterized with dynamic quantities such as particle size and swimming speed, without detailed knowledge of the diffusiophoretic mechanism. Published by the American Physical Society2024more » « less
-
We present a new method for introducing stable non-equilibrium concentration gradients in molecular dynamics simulations of mixtures. This method extends earlier Reverse Non-Equilibrium Molecular Dynamics (RNEMD) methods which use kinetic energy scaling moves to create temperature or velocity gradients. In the new scaled particle flux (SPF-RNEMD) algorithm, energies and forces are computed simultaneously for a molecule existing in two non-adjacent regions of a simulation box, and the system evolves under a linear combination of these interactions. A continuously increasing particle scaling variable is responsible for migration of the molecule between the regions as the simulation progresses, allowing for simulations under an applied particle flux. To test the method, we investigate diffusivity in mixtures of identical, but distinguishable particles, and in a simple mixture of multiple Lennard-Jones particles. The resulting concentration gradients provide Fick diffusion constants for mixtures. We also discuss using the new method to obtain coupled transport properties using simultaneous particle \textit{and} thermal fluxes to compute the temperature dependence of the diffusion coefficient and activation energies for diffusion from a single simulation. Lastly, we demonstrate the use of this new method in interfacial systems by computing the diffusive permeability for a molecular fluid moving through a nanoporous graphene membrane.more » « less
