Erratum: Autonomously Probing Viscoelasticity in Disordered Suspensions [Phys. Rev. Lett. 125, 258002 (2020)].

This corrects the article DOI: 10.1103/PhysRevLett.125.258002.

Magnetic propulsion of colloidal microrollers controlled by electrically modulated friction.

Precise control over the motion of magnetically responsive particles in fluidic chambers is important for probing and manipulating tasks in prospective microrobotic and bio-analytical platforms. We have previously exploited such colloids as shuttles for the microscale manipulation of objects. Here, we study the rolling motion of magnetically driven Janus colloids on solid substrates under the influence of an orthogonal external electric field. Electrically induced attractive interactions were used to tune the load on the Janus colloid and thereby the friction with the underlying substrate, leading to control over the forward velocity of the particle. Our experimental data suggest that the frictional coupling required to achieve translation, transitions from a hydrodynamic regime to one of mixed contact coupling with increasing load force. Based on this insight, we show that our colloidal microrobots can probe the local friction coefficient of various solid surfaces, which makes them potentially useful as tribological microsensors. Lastly, we precisely manipulate porous cargos using our colloidal rollers, a feat that holds promise for bio-analytical applications.

An extensible lattice Boltzmann method for viscoelastic flows: complex and moving boundaries in Oldroyd-B fluids.

 Most biological fluids are viscoelastic, meaning that they have elastic properties in addition to the dissipative properties found in Newtonian fluids. Computational models can help us understand viscoelastic flow, but are often limited in how they deal with complex flow geometries and suspended particles. Here, we present a lattice Boltzmann solver for Oldroyd-B fluids that can handle arbitrarily shaped fixed and moving boundary conditions, which makes it ideally suited for the simulation of confined colloidal suspensions. We validate our method using several standard rheological setups and additionally study a single sedimenting colloid, also finding good agreement with the literature. Our approach can readily be extended to constitutive equations other than Oldroyd-B. This flexibility and the handling of complex boundaries hold promise for the study of microswimmers in viscoelastic fluids.

Impact of surface charge on the motion of light-activated Janus micromotors.

Control over micromotors' motion is of high relevance for lab-on-a-chip and biomedical engineering, wherein such particles encounter complex microenvironments. Here, we introduce an efficient way to influence Janus micromotors' direction of motion and speed by modifying their surface properties and those of their immediate surroundings. We fabricated light-responsive Janus micromotors with positive and negative surface charge, both driven by ionic self-diffusiophoresis. These were used to observe direction-of-motion reversal in proximity to glass substrates for which we varied the surface charge. Quantitative analysis allowed us to extract the dependence of the particle velocity on the surface charge density of the substrate. This constitutes the first quantitative demonstration of the substrate's surface charge on the motility of the light-activated diffusiophoretic motors in water. We provide qualitative understanding of these observations in terms of osmotic flow along the substrate generated through the ions released by the propulsion mechanism. Our results constitute a crucial step in moving toward practical application of self-phoretic artificial micromotors.

Slip Length Dependent Propulsion Speed of Catalytic Colloidal Swimmers near Walls.

Catalytic colloidal swimmers that propel due to self-generated fluid flows exhibit strong affinity for surfaces. Here, we report experimental measurements of a significant dependence of such microswimmers' speed on the nearby substrate material. We find that speeds scale with the solution contact angle θ on the substrate, which relates to the associated hydrodynamic substrate slip length, as V∝(cosθ+1)^{-3/2}. We show that such dependence can be attributed to osmotic coupling between swimmers and substrate. Our work points out that hydrodynamic slip at nearby walls, though often unconsidered, can significantly impact microswimmer self-propulsion.

Autonomously Probing Viscoelasticity in Disordered Suspensions.

Recent experiments show a strong rotational diffusion enhancement for self-propelled microrheological probes in colloidal glasses. Here, we provide microscopic understanding using simulations with a frictional probe-medium coupling that converts active translation into rotation. Diffusive enhancement emerges from the medium's disordered structure and peaks at a second-order transition in the number of contacts. Our results reproduce the salient features of the colloidal glass experiment and support an effective description that is applicable to a broader class of viscoelastic suspensions.

Diffusion-Based Height Analysis Reveals Robust Microswimmer-Wall Separation.

Microswimmers typically move near walls, which can strongly influence their motion. However, direct experimental measurements of swimmer-wall separation remain elusive to date. Here, we determine this separation for model catalytic microswimmers from the height dependence of the passive component of their mean-squared displacement. We find that swimmers exhibit "ypsotaxis," a tendency to assume a fixed height above the wall for a range of salt concentrations, swimmer surface charges, and swimmer sizes. Our findings indicate that ypsotaxis is activity induced, posing restrictions on future modeling of their still-debated propulsion mechanism.

Regulating the aggregation of colloidal particles in an electro-osmotic micropump.

Unrestricted particle transport through microfluidic channels is of paramount importance to a wide range of applications, including lab-on-a-chip devices. In this article, we study via video microscopy the electro-osmotic aggregation of colloidal particles at the opening of a micrometer-sized silica channel in the presence of a salt gradient. Particle aggregation eventually leads to clogging of the channel, which may be undone by a time-adjusted reversal of the applied electric potential. We numerically model our system via the Stokes-Poisson-Nernst-Planck equations in a geometry that approximates the real sample. This allows us to identify the transport processes induced by the electric field and salt gradient and to provide evidence that a balance thereof leads to aggregation. We further demonstrate experimentally that a net flow of colloids through the channel may be achieved by applying a square-waveform electric potential with an appropriately tuned duty cycle. Our results serve to guide the design of microfluidic and nanofluidic pumps that allow for controlled particle transport and provide new insights for anti-fouling in ultra-filtration.

Self-thermoelectrophoresis at low salinity.

A locally heated Janus colloid can achieve motion in an electrolyte by an effect known as self-thermo(di)electrophoresis. We numerically study the self-propulsion of such a "hot swimmer" in a monovalent electrolyte using the finite-element method and analytic theory. The effect of electrostatic screening for intermediate and large Debye lengths is charted and we report on the fluid flow generated by self-thermoelectrophoresis. We obtain excellent agreement between our analytic theory and numerical calculations in the limit of high salinity, validating our approach. At low salt concentrations, we employ Teubner's integral formalism to arrive at expressions for the speed, which agree semi-quantitatively with our numerical results for conducting swimmers. This lends credibility to the remarkably high swim speed at very low ionic strength, which we numerically obtain for a fully insulating swimmer. We also report on hot swimmers with a mixed electrostatic boundary conditions. Our results should benefit the realization and analysis of further experiments on thermo(di)electrophoretic swimmers.

Hydrodynamic mobility reversal of squirmers near flat and curved surfaces.

Self-propelled particles have been experimentally shown to orbit spherical obstacles and move along surfaces. Here, we theoretically and numerically investigate this behavior for a hydrodynamic squirmer interacting with spherical objects and flat walls using three different methods of approximately solving the Stokes equations: The method of reflections, which is accurate in the far field; lubrication theory, which describes the close-to-contact behavior; and a lattice Boltzmann solver that accurately accounts for near-field flows. The method of reflections predicts three distinct behaviors: orbiting/sliding, scattering, and hovering, with orbiting being favored for lower curvature as in the literature. Surprisingly, it also shows backward orbiting/sliding for sufficiently strong pushers, caused by fluid recirculation in the gap between the squirmer and the obstacle leading to strong forces opposing forward motion. Lubrication theory instead suggests that only hovering is a stable point for the dynamics. We therefore employ lattice Boltzmann to resolve this discrepancy and we qualitatively reproduce the richer far-field predictions. Our results thus provide insight into a possible mechanism of mobility reversal mediated solely through hydrodynamic interactions with a surface.

A lattice Boltzmann model for squirmers.

The squirmer is a simple yet instructive model for microswimmers, which employs an effective slip velocity on the surface of a spherical swimmer to describe its self-propulsion. We solve the hydrodynamic flow problem with the lattice Boltzmann (LB) method, which is well-suited for time-dependent problems involving complex boundary conditions. Incorporating the squirmer into LB is relatively straightforward, but requires an unexpectedly fine grid resolution to capture the physical flow fields and behaviors accurately. We demonstrate this using four basic hydrodynamic tests: two for the far-field flow-accuracy of the hydrodynamic moments and squirmer-squirmer interactions-and two that require the near field to be accurately resolved-a squirmer confined to a tube and one scattering off a spherical obstacle-which LB is capable of doing down to the grid resolution. We find good agreement with (numerical) results obtained using other hydrodynamic solvers in the same geometries and identify a minimum required resolution to achieve this reproduction. We discuss our algorithm in the context of other hydrodynamic solvers and present an outlook on its application to multi-squirmer problems.

Hydrodynamics strongly affect the dynamics of colloidal gelation but not gel structure.

Colloidal particles with strong, short-ranged attractions can form a gel. We simulate this process without and with hydrodynamic interactions (HI), using the lattice-Boltzmann method to account for presence of a thermalized solvent. We show that HI speed up and slow down gelation at low and high volume fractions, respectively. The transition between these two regimes is linked to the existence of a percolating cluster shortly after quenching the system. However, when we compare gels at matched 'structural age', we find nearly indistinguishable structures with and without HI. Our result explains longstanding, unresolved conflicts in the literature.

Ionic screening and dissociation are crucial for understanding chemical self-propulsion in polar solvents.

Polar solvents like water support the bulk dissociation of themselves and their solutes into ions, and the re-association of these ions into neutral molecules in a dynamic equilibrium, e.g., H2O2 ⇌ H+ + HO2-. Using continuum theory, we study the influence of these association-dissociation reactions on the self-propulsion of colloids driven by surface chemical reactions (chemical swimmers). We find that association-dissociation reactions should have a strong influence on swimmers' behaviour, and therefore should be included in future modelling. In particular, such bulk reactions should permit charged swimmers to propel electrophoretically even if all species involved in the surface reactions are neutral. The bulk reactions also significantly modify the predicted speed of chemical swimmers propelled by ionic currents, by up to an order of magnitude. For swimmers whose surface reactions produce both anions and cations (ionic self-diffusiophoresis), the bulk reactions produce an additional reactive screening length, analogous to the Debye length in electrostatics. This in turn leads to an inverse relationship between swimmer radius and swimming speed, which could provide an alternative explanation for recent experimental observations on Pt-polystyrene Janus swimmers [S. Ebbens et al., Phys. Rev. E: Stat., Nonlinear, Soft Matter Phys., 2012, 85, 020401]. We also use our continuum theory to investigate the effect of the Debye screening length itself, going beyond the infinitely-thin-screening-length approximation used by previous analytical theories. We identify significant departures from this limiting behaviour for micron-sized swimmers under typical experimental conditions and find that the approximation fails entirely for nanoscale swimmers.

Microfluidic pumping by micromolar salt concentrations.

An ion-exchange-resin-based microfluidic pump is introduced that utilizes trace amounts of ions to generate fluid flows. We show experimentally that our pump operates in almost deionized water for periods exceeding 24 h and induces fluid flows of µm s-1 over hundreds of µm. This flow displays a far-field, power-law decay which is characteristic of two-dimensional (2D) flow when the system is strongly confined and of three-dimensional (3D) flow when it is not. Using theory and numerical calculations we demonstrate that our observations are consistent with electroosmotic pumping driven by µmol L-1 ion concentrations in the sample cell that serve as 'fuel' to the pump. Our study thus reveals that trace amounts of charge carriers can produce surprisingly strong fluid flows; an insight that should benefit the design of a new class of microfluidic pumps that operate at very low fuel concentrations.

Selective Trapping of DNA Using Glass Microcapillaries.

We show experimentally that an inexpensive glass microcapillary can accumulate λ-phage DNA at its tip and deliver the DNA into the capillary using a combination of electro-osmotic flow, pressure-driven flow, and electrophoresis. We develop an efficient simulation model based on the electrokinetic equations and the finite-element method to explain this phenomenon. As a proof of concept for the generality of this trapping mechanism we use our numerical model to explore the effect of the salt concentration, the capillary surface charge, the applied voltage, the pressure difference, and the mobility of the analyte molecules. Our results indicate that the simple microcapillary system has the potential to capture a wide range of analyte molecules based on their electrophoretic mobility that extends well beyond our experimental example of λ-phage DNA. Our method for separation and preconcentration of analytes therefore has implications for the development of low-cost lab-on-a-chip devices.

Understanding the onset of oscillatory swimming in microchannels.

Self-propelled colloids (swimmers) in confining geometries follow trajectories determined by hydrodynamic interactions with the bounding surfaces. However, typically these interactions are ignored or truncated to the lowest order. We demonstrate that higher-order hydrodynamic moments cause rod-like swimmers to follow oscillatory trajectories in quiescent fluid between two parallel plates, using a combination of lattice-Boltzmann simulations and far-field calculations. This behavior occurs even far from the confining walls and does not require lubrication results. We show that a swimmer's hydrodynamic quadrupole moment is crucial to the onset of the oscillatory trajectories. This insight allows us to develop a simple model for the dynamics near the channel center based on these higher hydrodynamic moments, and suggests opportunities for trajectory-based experimental characterization of swimmers' hydrodynamic properties.

Reducing spurious flow in simulations of electrokinetic phenomena.

Electrokinetic transport phenomena can strongly influence the behaviour of macromolecules and colloidal particles in solution, with applications in, e.g., DNA translocation through nanopores, electro-osmotic flow in nanocapillaries, and electrophoresis of charged macromolecules. Numerical simulations are an important tool to investigate these electrokinetic phenomena, but are often plagued by spurious fluxes and spurious flows that can easily exceed physical fluxes and flows. Here, we present a method that reduces one of these spurious currents, spurious flow, by several orders of magnitude. We demonstrate the effectiveness and generality of our method for both the electrokinetic lattice-Boltzmann and finite-element-method based algorithms by simulating a charged sphere in an electrolyte solution and flow through a nanopore. We also show that previous attempts to suppress these spurious currents introduce other sources of error.

The efficiency of self-phoretic propulsion mechanisms with surface reaction heterogeneity.

We consider the efficiency of self-phoretic colloidal particles (swimmers) as a function of the heterogeneity in the surface reaction rate. The set of fluid, species, and electrostatic continuity equations is solved analytically using a linearization and numerically using a finite-element method. To compare spherical swimmers of different size and with heterogeneous catalytic conversion rates, a "swimmer efficiency" functional η is introduced. It is proven that in order to obtain maximum swimmer efficiency, the reactivity has to be localized at the pole(s). Our results also shed light on the sensitivity of the propulsion speed to details of the surface reactivity, a property that is notoriously hard to measure. This insight can be utilized in the design of new self-phoretic swimmers.

Surface roughness stabilizes the clustering of self-propelled triangles.

Self-propelled particles can spontaneously form dense phases from a dilute suspension in a process referred to as motility-induced phase separation. The properties of the out-of-equilibrium structures that are formed are governed by the specifics of the particle interactions and the strength of the activity. Thus far, most studies into the formation of these structures have focused on spherical colloids, dumbbells, and rod-like particles endowed with various interaction potentials. Only a few studies have examined the collective behavior of more complex particle shapes. Here, we increase the geometric complexity and use molecular dynamics simulations to consider the structures formed by triangular self-propelled particles with surface roughness. These triangles either move towards their apex or towards their base, i.e., they possess a polarity. We find that apex-directed triangles cluster more readily, more stably, and have a smoother cluster interface than their base-directed counterparts. A difference between the two polarities is in line with the results of Wensink et al. [Phys. Rev. E 89, 010302 (2014)]; however, we obtain the reversed result when it comes to clustering, namely, that apex-directed triangles cluster more successfully. We further show that reducing the surface roughness negatively impacts the stability of the base-directed structures, suggesting that their formation is in large part due to surface roughness. Our results lay a solid foundation for future experimental and computational studies into the effect of roughness on the collective dynamics of swimmers.

Moving charged particles in lattice Boltzmann-based electrokinetics.

The motion of ionic solutes and charged particles under the influence of an electric field and the ensuing hydrodynamic flow of the underlying solvent is ubiquitous in aqueous colloidal suspensions. The physics of such systems is described by a coupled set of differential equations, along with boundary conditions, collectively referred to as the electrokinetic equations. Capuani et al. [J. Chem. Phys. 121, 973 (2004)] introduced a lattice-based method for solving this system of equations, which builds upon the lattice Boltzmann algorithm for the simulation of hydrodynamic flow and exploits computational locality. However, thus far, a description of how to incorporate moving boundary conditions into the Capuani scheme has been lacking. Moving boundary conditions are needed to simulate multiple arbitrarily moving colloids. In this paper, we detail how to introduce such a particle coupling scheme, based on an analogue to the moving boundary method for the pure lattice Boltzmann solver. The key ingredients in our method are mass and charge conservation for the solute species and a partial-volume smoothing of the solute fluxes to minimize discretization artifacts. We demonstrate our algorithm's effectiveness by simulating the electrophoresis of charged spheres in an external field; for a single sphere we compare to the equivalent electro-osmotic (co-moving) problem. Our method's efficiency and ease of implementation should prove beneficial to future simulations of the dynamics in a wide range of complex nanoscopic and colloidal systems that were previously inaccessible to lattice-based continuum algorithms.