Colloids dragged through a polymer solution: Experiment, theory, and simulation

Nonlinear microrheology with time-dependent forces: Application to recoils in viscoelastic fluids

This work presents a theoretical analysis of the motion of a tracer colloid driven by a time-dependent force through a viscoelastic fluid. The recoil of the colloid after application of a strong force is determined. It provides insights into the elastic forces stored locally in the fluid and their weakening by plastic processes. We generalize the mode-coupling theory of microrheology to include time-dependent forces. After deriving the equations of motion for the tracer correlator and simplifying to a schematic model, we apply the theory to a switch-off force protocol that features the recoiling of the tracer after cessation of the driving. We also include Langevin dynamics simulations to compare to the results of the theory. A nonmonotonic trend of the recoil amplitude is found in the theory and confirmed in the simulations. The linear-response approximation is also verified in the small-force regime. While the overall agreement between simulation and theory is good, simulation shows that the theory predicts a too strong nonmonotonous dependence of the recoil distance on the applied force.

Nonlinear Langevin functionals for a driven probe

When a probe particle immersed in a fluid with nonlinear interactions is subject to strong driving, the cumulants of the stochastic force acting on the probe are nonlinear functionals of the driving protocol. We present a Volterra series for these nonlinear functionals by applying nonlinear response theory in a path integral formalism, where the emerging kernels are shown to be expressed in terms of connected equilibrium correlation functions. The first cumulant is the mean force, the second cumulant characterizes the non-equilibrium force fluctuations (noise), and higher order cumulants quantify non-Gaussian fluctuations. We discuss the interpretation of this formalism in relation to Langevin dynamics. We highlight two example scenarios of this formalism. (i) For a particle driven with the prescribed trajectory, the formalism yields the non-equilibrium statistics of the interaction force with the fluid. (ii) For a particle confined in a moving trapping potential, the formalism yields the non-equilibrium statistics of the trapping force. In simulations of a model of nonlinearly interacting Brownian particles, we find that nonlinear phenomena, such as shear-thinning and oscillating noise covariance, appear in third- or second-order response, respectively.

Lattice Boltzmann approach for acoustic manipulation

We employ a lattice Boltzmann method to compute the acoustic radiation force produced by standing waves on a compressible object for the density matched case. Instead of simulating the fluid mechanics equations directly, the proposed method uses a lattice Boltzmann model that reproduces the wave equation, together with a kernel interpolation scheme, to compute the first-order perturbations of the pressure and velocity fields on the object's surface and, from them, the acoustic radiation force. The procedure reproduces with excellent accuracy the theoretical expressions by Gor'kov and Wei for the sphere as the 3D case and an infinitely long cylinder as the 2D case, respectively, even with a modest number of lattice Boltzmann cells. The proposed method shows to be a promising tool for simulating phenomena where the acoustic radiation force plays a relevant role, like acoustic tweezers or the acoustic manipulation of microswimmers, with applications in medicine and engineering.

Force renormalization for probes immersed in an active bath

Langevin equations or generalized Langevin equations (GLEs) are popular models for describing the motion of a particle in a fluid medium in an effective manner. Here we examine particles immersed in an inherently nonequilibrium fluid, i.e., an active bath, which are subject to an external force. Specifically, we consider two types of forces that are highly relevant for microrheological studies: A harmonic, trapping force and a constant, "drag" force. We study such systems by molecular simulations and use the simulation data to extract an effective GLE description. We find that within this description, in an active bath, the external force in the GLE is not equal to the physical external force, but rather a renormalized external force, which can be significantly smaller. The effect cannot be attributed to the mere temperature renormalization, which is also observed.

Velocity jump process with volume exclusions in a narrow channel

This paper analyzes the impact of collisions in a system of N identical hard-core particles driven according to a velocity jump process. The physical space is essentially a channel in R with a probability of occupants being able to pass each other. The system mimics what nature does, where individuals pass one another in a narrow channel while making incidental contact with those moving in the opposite direction. The passing probability may depend on the particles' size and the channel's width. Starting from the particle level model, we systematically derive a nonlinear transport equation based on an asymptotic expansion. Under low-occupied fractions, numerical solutions of both the kinetic model and the stochastic particle system are compared well during biased and unbiased random velocity changes. Analysis of the subpopulation motility within a large population exhibits the consequences of volume exclusions and channel confinements on the traveling speeds.

Attraction and condensation of driven tracers in a narrow channel

Emergent bath-mediated attraction and condensation arise when multiple particles are simultaneously driven through an equilibrated bath under geometric constraints. While such scenarios are observed in a variety of nonequilibrium phenomena with an abundance of experimental and numerical evidence, little quantitative understanding of how these interactions arise is currently available. Here we approach the problem by studying the behavior of two driven "tracer" particles, propagating through a bath in a 1D lattice with excluded-volume interactions. We apply the mean-field approximation to analytically explore the mechanism responsible for the tracers' emergent interactions and compute the resulting effective attractive potential. This mechanism is then numerically shown to extend to a realistic model of hard driven Brownian disks confined to a narrow 2D channel.

Two step micro-rheological behavior in a viscoelastic fluid

We perform micro-rheological experiments with a colloidal bead driven through a viscoelastic worm-like micellar fluid and observe two distinctive shear thinning regimes, each of them displaying a Newtonian-like plateau. The shear thinning behavior at larger velocities is in qualitative agreement with macroscopic rheological experiments. The second process, observed at Weissenberg numbers as small as a few percent, appears to have no analog in macro-rheological findings. A simple model introduced earlier captured the observed behavior and implied that the two shear thinning processes correspond to two different length scales in the fluid. This model also reproduces oscillations, which have been observed in this system previously. While the system under macro-shear seems to be near equilibrium for shear rates in the regime of the intermediate Newtonian-like plateau, the one under micro-shear is thus still far from it. The analysis suggests the existence of a length scale of a few micrometres, the nature of which remains elusive.

Impact of bacteria motility in the encounter rates with bacteriophage in mucus

Bacteriophages—or phages—are viruses that infect bacteria and are present in large concentrations in the mucosa that cover the internal organs of animals. Immunoglobulin (Ig) domains on the phage surface interact with mucin molecules, and this has been attributed to an increase in the encounter rates of phage with bacteria in mucus. However, the physical mechanism behind this phenomenon remains unclear. A continuous time random walk (CTRW) model simulating the diffusion due to mucin-T4 phage interactions was developed and calibrated to empirical data. A Langevin stochastic method for Escherichia coli (E. coli) run-and-tumble motility was combined with the phage CTRW model to describe phage-bacteria encounter rates in mucus for different mucus concentrations. Contrary to previous theoretical analyses, the emergent subdiffusion of T4 in mucus did not enhance the encounter rate of T4 against bacteria. Instead, for static E. coli, the diffusive T4 mutant lacking Ig domains outperformed the subdiffusive T4 wild type. E. coli’s motility dominated the encounter rates with both phage types in mucus. It is proposed, that the local fluid-flow generated by E. coli’s motility combined with T4 interacting with mucins may be the mechanism for increasing the encounter rates between the T4 phage and E. coli bacteria.

Optical Tweezers Microrheology Maps the Dynamics of Strain-Induced Local Inhomogeneities in Entangled Polymers

Optical tweezers microrheology (OTM) offers a powerful approach to probe the nonlinear response of complex soft matter systems, such as networks of entangled polymers, over wide-ranging spatiotemporal scales. OTM can also uniquely characterize the microstructural dynamics that lead to the intriguing nonlinear rheological properties that these systems exhibit. However, the strain in OTM measurements, applied by optically forcing a microprobe through the material, induces network inhomogeneities in and around the strain path, and the resultant flow field complicates the measured response of the system. Through a robust set of custom-designed OTM protocols, coupled with modeling and analytical calculations, we characterize the time-varying inhomogeneity fields induced by OTM measurements. We show that homogenization following strain does not interfere with the intrinsic stress relaxation dynamics of the system, rather it manifests as an independent component in the stress decay, even in highly nonlinear regimes such as with the microrheological large-amplitude-oscillatory-shear (MLAOS) protocols we introduce. Our specific results show that Rouse-like elastic retraction, rather than disentanglement and disengagement, dominates the nonlinear stress relaxation of entangled polymers at micro- and mesoscales. Thus, our study opens up possibilities of performing precision nonlinear microrheological measurements, such as MLAOS, on a wide range of complex macromolecular systems.

Driven spheres, ellipsoids and rods in explicitly modeled polymer solutions

Understanding the transport of driven nano- and micro-particles in complex fluids is of relevance for many biological and technological applications. Here we perform hydrodynamic multiparticle collision dynamics simulations of spherical and elongated particles driven through polymeric fluids containing different concentrations of polymers. We determine the mean particle velocities which are larger than expected from Stokes law for all particle shapes and polymer densities. Furthermore we measure the fluid flow fields and local polymer density and polymer conformation around the particles. We find that polymer-depleted regions close to the particles are responsible for an apparent tangential slip velocity which accounts for the measured flow fields and transport velocities. A simple two-layer fluid model gives a good match to the simulation results.

A mesoscale study of creep in a microgel using the acoustic radiation force

We study the motion of a sphere of diameter 330 μm embedded in a Carbopol microgel under the effect of the acoustic radiation pressure exerted by a focused ultrasonic field. The sphere motion within the microgel is tracked using videomicroscopy and compared to conventional creep and recovery measurements performed with a rheometer. We find that under moderate ultrasonic intensities, the sphere creeps as a power law of time with an exponent α ≃ 0.2 that is significantly smaller than the one inferred from global creep measurements below the yield stress of the microgel (α ≃ 0.4). Moreover, the sphere relaxation motion after creep and the global recovery are respectively consistent with these two different exponents. By allowing a rheological characterization at the scale of the sphere with forces of the order of micronewtons, the present experiments pave the way for acoustic "mesorheology" which probes volumes and forces an intermediate between standard macroscopic rheology and classical microrheology. They also open new questions about the effects of the geometry of the deformation field and of the sphere size and surface properties on the creep behaviour of soft materials.

Two-point active microrheology in a viscous medium exploiting a motional resonance excited in dual-trap optical tweezers

Two-point microrheology measurements from widely separated colloidal particles approach the bulk viscosity of the host medium more reliably than corresponding single-point measurements. In addition, active microrheology offers the advantage of enhanced signal to noise over passive techniques. Recently, we reported the observation of a motional resonance induced in a probe particle in dual-trap optical tweezers when the control particle was driven externally [Paul et al., Phys. Rev. E 96, 050102(R) (2017)2470-004510.1103/PhysRevE.96.050102]. We now demonstrate that the amplitude and phase characteristics of the motional resonance can be used as a sensitive tool for active two-point microrheology to measure the viscosity of a viscous fluid. Thus, we measure the viscosity of viscous liquids from both the amplitude and phase response of the resonance, and demonstrate that the zero crossing of the phase response of the probe particle with respect to the external drive is superior compared to the amplitude response in measuring viscosity at large particle separations. We compare our viscosity measurements with those using a commercial rheometer and obtain an agreement ∼1%. The method can be extended to viscoelastic material where the frequency dependence of the resonance may provide further accuracy for active microrheological measurements.

Oscillating modes of driven colloids in overdamped systems

Microscopic colloidal particles suspended in liquids are a prominent example of an overdamped system where viscous forces dominate over inertial effects. Frequently, colloids are used as sensitive probes, e.g., in biophysical applications from which molecular forces are inferred. The interpretation of such experiments rests on the assumption that, even when the particles are driven, the liquid remains in equilibrium. Here we experimentally demonstrate that this is not valid for particles in viscoelastic fluids. Even at small driving forces, we observe particle oscillations with several tens of seconds. They are attributed to non-equilibrium fluctuations of the fluid being excited by the particle’s motion. The oscillatory dynamics is in quantitative agreement with an overdamped Langevin equation with negative friction-memory term being equivalent to a stochastically driven underdamped oscillator. Such oscillatory modes are expected to widen the use of colloids as model systems but must also be considered in colloidal probe experiments.

The motion of microparticles suspended in liquids is assumed to be dominated by viscous forces. Here, Berner et al. challenge this consensus by observing underdamped particle oscillations in a viscoelastic fluid and attributing it to the non-equilibrium fluctuations of liquid excited by particles.

Nonequilibrium depletion interactions in active microrheology

Entropic depletion forces arise between mesoscopic bodies that are immersed in a suspension of macromolecules, such as colloid-polymer mixtures. Here we consider the case of a driven colloidal probe in the presence of another, passive colloidal particle, both solvated in an ideal bath of small spherical particles. We calculate the nonequilibrium forces mediated by the depletants on the two colloidal particles within a dynamical superposition approximation (DSA) scheme. In order to assess the quality of this approximation, and to obtain the colloidal microstructure around the driven probe, we corroborate our theoretical results with Brownian dynamics simulations.

Driven tracers in narrow channels

Steady-state properties of a driven tracer moving in a narrow two-dimensional (2D) channel of quiescent medium are studied. The tracer drives the system out of equilibrium, perturbs the density and pressure fields, and gives the bath particles a nonzero average velocity, creating a current in the channel. Three models in which the confining effect of the channel is probed are analyzed and compared in this study: the first is the simple symmetric exclusion process (SSEP), for which the stationary density profile and the pressure on the walls in the frame of the tracer are computed. We show that the tracer acts like a dipolar source in an average velocity field. The spatial structure of this 2D strip is then simplified to a one-dimensional (1D) SSEP, in which exchanges of position between the tracer and the bath particles are allowed. Using a combination of mean-field theory and exact solution in the limit where no exchange is allowed gives good predictions of the velocity of the tracer and the density field. Finally, we show that results obtained for the 1D SSEP with exchanges also apply to a gas of overdamped hard disks in a narrow channel. The correspondence between the parameters of the SSEP and of the gas of hard disks is systematic and follows from simple intuitive arguments. Our analytical results are checked numerically.

Nonlinear response and emerging nonequilibrium microstructures for biased diffusion in confined crowded environments

We study analytically the dynamics and the microstructural changes of a host medium caused by a driven tracer particle moving in a confined, quiescent molecular crowding environment. Imitating typical settings of active microrheology experiments, we consider here a minimal model comprising a geometrically confined lattice system (a two-dimensional striplike or a three-dimensional capillary-like system) populated by two types of hard-core particles with stochastic dynamics (a tracer particle driven by a constant external force and bath particles moving completely at random). Resorting to a decoupling scheme, which permits us to go beyond the linear-response approximation (Stokes regime) for arbitrary densities of the lattice gas particles, we determine the force-velocity relation for the tracer particle and the stationary density profiles of the host medium particles around it. These results are validated a posteriori by extensive numerical simulations for a wide range of parameters. Our theoretical analysis reveals two striking features: (a) We show that, under certain conditions, the terminal velocity of the driven tracer particle is a nonmonotonic function of the force, so in some parameter range the differential mobility becomes negative, and (b) the biased particle drives the whole system into a nonequilibrium steady state with a stationary particle density profile past the tracer, which decays exponentially, in sharp contrast with the behavior observed for unbounded lattices, where an algebraic decay is known to take place.

Nonlinear microrheology reveals entanglement-driven molecular-level viscoelasticity of concentrated DNA

We optically drive a trapped microscale probe through entangled DNA at rates up to 100× the disentanglement rate (Wi≈100), then remove the trap and track subsequent probe recoil motion. We identify a unique crossover to the nonlinear regime at Wi≈20. Recoil dynamics display rate-dependent dilation and complex power-law healing of the reptation tube. The force response during strain exhibits key nonlinear features such as shear thinning and yielding with power-law rate dependence. Our results, distinctly nonclassical and in accord with recent theoretical predictions, reveal molecular dynamics governed by individual stress-dependent entanglements rather than chain stretching.

Microrheology of colloidal systems

Microrheology was proposed almost twenty years ago as a technique to obtain rheological properties in soft matter from the microscopic motion of colloidal tracers used as probes, either freely diffusing in the host medium, or subjected to external forces. The former case is known as passive microrheology, and is based on generalizations of the Stokes-Einstein relation between the friction experienced by the probe and the host-fluid viscosity. The latter is termed active microrheology, and extends the measurement of the friction coefficient to the nonlinear-response regime of strongly driven probes. In this review article, we discuss theoretical models available in the literature for both passive and active microrheology, focusing on the case of single-probe motion in model colloidal host media. A brief overview of the theory of passive microrheology is given, starting from the work of Mason and Weitz. Further developments include refined models of the host suspension beyond that of a Newtonian-fluid continuum, and the investigation of probe-size effects. Active microrheology is described starting from microscopic equations of motion for the whole system including both the host-fluid particles and the tracer; the many-body Smoluchowski equation for the case of colloidal suspensions. At low fluid densities, this can be simplified to a two-particle equation that allows the calculation of the friction coefficient with the input of the density distribution around the tracer, as shown by Brady and coworkers. The results need to be upscaled to agree with simulations at moderate density, in both the case of pulling the tracer with a constant force or dragging it at a constant velocity. The full many-particle equation has been tackled by Fuchs and coworkers, using a mode-coupling approximation and the scheme of integration through transients, valid at high densities. A localization transition is predicted for a probe embedded in a glass-forming host suspension. The nonlinear probe-friction coefficient is calculated from the tracer's position correlation function. Computer simulations show qualitative agreement with the theory, but also some unexpected features, such as superdiffusive motion of the probe related to the breaking of nearest-neighbor cages. We conclude with some perspectives and future directions of theoretical models of microrheology.

Anomalous field-induced growth of fluctuations in dynamics of a biased intruder moving in a quiescent medium

We present exact results on the dynamics of a biased, by an external force F, intruder (BI) in a two-dimensional lattice gas of unbiased, randomly moving hard-core particles. Going beyond the usual analysis of the force-velocity relation, we study the probability distribution P(R(n)) of the BI displacement R(n) at time n. We show that despite the fact that the BI drives the gas to a nonequilibrium steady state, P(R(n)) converges to a Gaussian distribution as n→∞. We find that the variance σ(x)(2) of P(R(n)) along F exhibits a weakly superdiffusive growth σ(x)(2)~ν(1)nln(n), and a usual diffusive growth, σ(y)(2)~ν(2)n, in the perpendicular direction. We determine ν(1) and ν(2) exactly for arbitrary bias, in the lowest order in the density of vacancies, and show that ν(1)~|F|(2) for small bias, which signifies that superdiffusive behavior emerges beyond the linear-response approximation. We also present analytical arguments predicting a striking field-induced superdiffusive behavior σ(x)(2)~n(3/2) for two-dimensional stripes and three-dimensional capillaries, which is confirmed by Monte Carlo simulations.

Micro-macro-discrepancies in nonlinear microrheology: I. Quantifying mechanisms in a suspension of Brownian ellipsoids

Active and nonlinear microrheology experiments involve a colloidal probe that is forced to move within a material, with the goal of recovering the nonlinear rheological response properties of the material. Various mechanisms cause discrepancies between the nonlinear rheology measured microrheologically and macroscopically, including direct probe-bath collisions, the Lagrangian unsteadiness experienced by the material elements, and the spatially inhomogeneous and rheologically mixed strain field set up around the probe. Here, we perform computational nonlinear microrheology experiments, in which a colloidal probe translates through a dilute suspension of Brownian ellipsoids, whose results we compare against analogous computational experiments on the macroscopic shear rheology of the same model material. The quantitative impact of each of the mechanisms for micro-macro-discrepancy can thus be computed directly, with additional computational experiments performed where the processes in question are 'turned off'. We show that all three discrepancy mechanisms impact the microrheological measurement quantitatively, and that none can be neglected. This motivates a search for microrheological probes whose geometry or forcing is optimized to minimize these impacts, which we present in a companion article.

Hydrodynamic interactions in active colloidal crystal microrheology

In dense colloids it is commonly assumed that hydrodynamic interactions do not play a role. However, a found theoretical quantification is often missing. We present computer simulations that are motivated by experiments where a large colloidal particle is dragged through a colloidal crystal. To qualify the influence of long-ranged hydrodynamics, we model the setup by conventional Langevin dynamics simulations and by an improved scheme with limited hydrodynamic interactions. This scheme significantly improves our results and allows to show that hydrodynamics strongly impacts the development of defects, the crystal regeneration, as well as the jamming behavior.

Micro-rheology on (polymer-grafted) colloids using optical tweezers

Optical tweezers are experimental tools with extraordinary resolution in positioning (± 1 nm) a micron-sized colloid and in the measurement of forces (± 50 fN) acting on it-without any mechanical contact. This enables one to carry out a multitude of novel experiments in nano- and microfluidics, of which the following will be presented in this review: (i) forces within single pairs of colloids in media of varying concentration and valency of the surrounding ionic solution, (ii) measurements of the electrophoretic mobility of single colloids in different solvents (concentration, valency of the ionic solution and pH), (iii) similar experiments as in (i) with DNA-grafted colloids, (iv) the nonlinear response of single DNA-grafted colloids in shear flow and (v) the drag force on single colloids pulled through a polymer solution. The experiments will be described in detail and their analysis discussed.

Dynamics of colloids in confined geometries

We discuss the Brownian dynamics of colloids in confinement with special emphasis on the influence of the solvent dynamics. We review the derivation of a dynamic density functional theory (DDFT) including some aspects of hydrodynamic interactions and its application to the micro-rheology of suspensions. In particular we discuss the failure of Stokes' law in suspensions and non-equilibrium solvent structure mediated interactions. With regard to hydrodynamic chromatography we also discuss the stationary transport of particles in narrow channels, and the reasons for the failure of DDFT in this situation.

Small-world rheology: an introduction to probe-based active microrheology

We introduce active, probe-based microrheological techniques for measuring the flow and deformation of complex fluids. These techniques are ideal for mechanical characterization either when little sample is available, or when samples show significant spatial heterogeneity. We review recent results, paying particular attention to comparing and contrasting rheological parameters obtained from micro- and macro-rheological techniques.

Active microrheology: a proposed technique to measure normal stress coefficients of complex fluids

We propose a microrheological technique to measure normal stress coefficients (NSCs) of complex fluids, which would represent the first quantitatively accurate measurement of a nonlinear rheological property by microrheology. Specifically, the mechanical response of almost all complex fluids to "weakly nonlinear" deformations is described by the second-order fluid model. Two microrheological probes pulled with equal velocities through a second-order fluid experience a relative force that is linear in the first and second NSCs of the complex fluid. We compute the coupling matrix between NSCs and relative forces for probes translating parallel and perpendicular to their line of centers, which can be inverted to yield NSCs from measured relative forces. There exists an optimum probe separation for inversion of the coupling matrix and, hence, experimental recovery of NSCs.

DDFT for Brownian particles and hydrodynamics

Classical dynamic density functional theory (DDFT) has developed into a versatile tool for describing the dynamics of overdamped Brownian particles. The underlying equilibrium density functional is a quantitative model for colloidal suspensions. Non-equilibrium properties of these systems are strongly influenced by the solvent hydrodynamics. In this paper we discuss ways to introduce certain aspects of solvent hydrodynamics into the DDFT, in particular advection by a flowing solvent and hydrodynamic interactions among suspended particles as well as with container walls.

Diffusion of a sphere in a dilute solution of polymer coils

We calculate the short time and the long time diffusion coefficients of a spherical tracer particle in a polymer solution in the low density limit by solving the Smoluchowski equation for a two-particle system and applying a generalized Einstein relation (fluctuation dissipation theorem). The tracer particle as well as the polymer coils are idealized as hard spheres with a no-slip boundary condition for the solvent but the hydrodynamic radius of the polymer coils is allowed to be smaller than the direct-interaction radius. We take hydrodynamic interactions up to 11th order in the particle distance into account. For the limit of small polymers, the expected generalized Stokes-Einstein relation is found. The long time diffusion coefficient also roughly obeys the generalized Stokes-Einstein relation for larger polymers whereas the short time coefficient does not. We find good qualitative and quantitative agreement to experiments.

Selectivity in binary fluid mixtures: static and dynamical properties

Selectivity of particles in a region of space can be achieved by applying external potentials to influence the particles in that region. We investigate static and dynamical properties of size selectivity in binary fluid mixtures of two particles sizes. We find that by applying an external potential that is attractive to both kinds of particles, due to crowding effects, this can lead to one species of particles being expelled from that region, while the other species is attracted into the region where the potential is applied. This selectivity of one species of particle over the other in a localized region of space depends on the density and composition of the fluid mixture. Applying an external potential that repels both kinds of particles leads to selectivity of the opposite species of particles to the selectivity with attractive potentials. We use equilibrium and dynamical density-functional theory to describe and understand the static and dynamical properties of this striking phenomenon. Selectivity by some ion channels is believed to be due to this effect.