Search

Magnetotaxis is a well known phenomenon in swimming microorganisms which sense magnetic fields e.g. by incorporating crystalline magnetosomes. In designing artificial active matter with tunable dynamics, external magnetic fields can provide a versatile method for guidance. Here, it is the question what material properties are necessary to elicit a significant response. In this working paper, we document in experiments on self-propelling nematic microdroplets that the weak diamagnetic torques exerted by a sub-Tesla magnetic field are already sufficient to significantly affect the dynamics of these swimmers. We find a rich and nontrivial variety of dynamic modes by varying droplet size, state of confinement and magnetic field strength.

We explore the relaxation dynamics of a tracer in a harmonic trap coupled to a non-equilibrium bath particle in stationary state, finding qualitative differences compared to the well known equilibrium case. These can be attributed to an additional position and time dependent force acting on the tracer, emerging when averaging the bath degree in the non-equilibrium stationary state conditioned to a certain tracer position. Specifically, we provide analytical results for an overdamped tracer coupled linearly to a bath particle in different nonequilibrium scenarios, namely, subjected to a different temperature than the tracer, or to active noise with Gaussian or non-Gaussian fluctuations. For the case of different temperatures, the conditioned tracer-bath force can be as large in magnitude as the force from the trapping potential. For an active bath particle with memory, even the bath noise takes a finite average under conditioning of the tracer. Further, if the noise of the bath particle is non-Gaussian, the relaxation function of the tracer can be non-monotonic as a function of time. We also compute the \emph{pinned relaxation}, proposing that measurement and comparison of conditioned and pinned relaxation allows determination of the non-equilibrium forces in experiments.

We give an overview exploring the role of kinetics in multicomponent mixtures. Compared to the most commonly studied binary (single species plus solvent) case, multicomponent fluids show a rich interplay between kinetics and thermodynamics due to the possibility of fractionation, interdiffusion of mixture components and collective motion. This leads to a competition between multiple timescales that change depending on the underlying kinetics. At high densities, crowding effects are relevant and non-equilibrium structures can become long-lived. We present the main approaches for the study of kinetic effects in multicomponents mixtures, including the role of crowding, and explore their consequences for equilibrium and non-equilibrium scenarios. We conclude by identifying the main challenges in the field.

The precise control of complex reactions is critical for biological processes yet our inability to design for specific outcomes limits the development of synthetic analogues. Here, we leverage differentiable simulators to design nontrivial reaction pathways in colloidal assemblies. By optimizing over external structures, we achieve controlled disassembly and particle release from colloidal shells. Lastly, we characterize the role of configurational entropy in the structure via both forward calculations and optimization, inspiring new parameterizations of designed colloidal reactions.

Latest since the landmark studies of Kovacs and co-workers on the glass transition of polymers, it is clear that thermally induced volume changes are of central importance for the understanding of the nature of the glass transition. Due to the kinetic background of the canonical (thermal) glass transition, it does not seem possible to derive a well-defined glass transition temperature T_g based on susceptibilities such as the thermal volume expansion coefficient, \beta(T), being strongly coupled to the glass transition process. Therefore, in practice, T_g is for example defined via the inflection point of the step-like \beta(T) curve. In this publication, we propose to use a thermo-optical feature, preceding the glass transition in the high-temperature phase, to determine the glass transition temperature T_g of a model polymer in a rather narrow temperature interval., Comment: Preprint

When a colloidal probe is forced through a viscoelastic fluid which is characterized by a long stress-relaxation time, the fluid is excited out of equilibrium. This is leading to a number of interesting effects including a non-trivial recoil of the probe when the driving force is removed. Here, we experimentally and theoretically investigate the transient recoil dynamics of non-spherical particles, i.e., colloidal dumbbells. In addition to a translational recoil of the dumbbells, we also find a pronounced angular reorientation which results from the relaxation of the surrounding fluid. Our findings are in good agreement with a Langevin description based on the symmetries of a director (dumbbell) as well as a microscopic bath-rod model. Remarkably, we find a frustrated state with amplified fluctuations when the dumbbell is oriented perpendicular to the direction of driving. Our results demonstrate the complex behavior of non-spherical objects within a relaxing environment which are of immediate interest for the motion of externally but also self-driven asymmetric objects in viscoelastic fluids., Comment: 10 pages, 6 figures

Topological defects play a central role in the formation and organization of various biological systems. Historically, such nonequilibrium defects have been mainly studied in the context of homogeneous active nematics. Phase-separated systems, in turn, are known to form dense and dynamic nematic bands, but typically lack topological defects. In this paper, we use agent-based simulations of weakly aligning, self-propelled polymers and demonstrate that contrary to the existing paradigm phase-separated active nematics form ${-1/2}$ defects. Moreover, these defects, emerging due to interactions among dense nematic bands, constitute a novel second-order collective state. We investigate the morphology of defects in detail and find that their cores correspond to a strong increase in density, associated with a condensation of nematic fluxes. Unlike their analogs in homogeneous systems, such condensed defects form and decay in a different way and do not involve positively charged partners. We additionally observe and characterize lateral arc-like structures that separate from a band's bulk and move in transverse direction. We show that the key control parameters defining the route from stable bands to the coexistence of dynamic lanes and defects are the total density of particles and their path persistence length. We introduce a hydrodynamic theory that qualitatively recapitulates all the main features of the agent-based model, and use it to show that the emergence of both defects and arcs can be attributed to the same anisotropic active fluxes. Finally, we present a way to artificially engineer and position defects, and speculate about experimental verification of the provided model.

Spinning objects which move through air or liquids experience a Magnus force. This effect is commonly exploited in ball sports but also of considerable importance for applications and fundamental science. Opposed to large objects where Magnus forces are strong, they are only weak at small scales and eventually vanish for overdamped micron-sized particles in simple liquids. Here we demonstrate an about one-million-fold enhanced Magnus force of spinning colloids in viscoelastic fluids. Such fluids are characterized by a time-delayed response to external perturbations which causes a deformation of the fluidic network around the moving particle. When the particle additionally spins, the deformation field becomes misaligned relative to the particle's moving direction, leading to a force perpendicular to the direction of travel and the spinning axis. The presence of strongly enhanced memory-induced Magnus forces at microscales opens novel applications for particle sorting and steering, the creation and visualization of anomalous flows and more., Comment: 20 pages, 5 figures

We derive a bound for entropy production in terms of the mean of normalizable path-antisymmetric observables. The optimal observable for this bound is shown to be the signum of entropy production, which is often easier determined or estimated than entropy production itself. It can be preserved under coarse graining by use of a simple path grouping algorithm. We demonstrate this relation and its properties using a driven network on a ring, for which the bound saturates for short times for any driving strength. This work can open a way to systematic coarse graining of entropy production., Comment: 5 pages, 2 figures

The slow viscous flow through a doubly-periodic array of cylinders does not have an analytical solution. However, as a reduced model for the flow within fibrous porous media, this solution is important for many real-world systems. We asymptotically determine the flow around a doubly-periodic array of infinite slender cylinders, by placing doubly-periodic two-dimensional singularity solutions within the cylinder and expanding the no-slip condition on the cylinder's surface in powers of the cylinder radius. The asymptotic solution provides a closed-form estimate for the flow and forces as a function of the radius and the dimensions of the cell. The force is compared to results from lattice-Boltzmann simulations of low-Reynolds-number flows in the same geometry, and the accuracy of the no-slip condition on the surface of the cylinder, predicted by the asymptotic theory, is checked. Finally, the behaviour of the flow, flux, force and effective permeability of the cell is investigated as a function of the geometric parameters. The structure of the asymptotic permeability is consistent with other models for the flow parallel to an array of rods. These models could be used to help understand the flows within porous systems composed of fibres and systems involving periodic arrays such as deterministic lateral displacement.

The motion of a colloidal probe in a viscoelastic fluid is described by friction or mobility, depending on whether the probe is moving with a velocity or feeling a force. While the Einstein relation describes an inverse relationship valid for Newtonian solvents, both concepts are generalized to time-dependent memory kernels in viscoelastic fluids. We theoretically and experimentally investigate their relation by considering two observables: the recoil after releasing a probe that was moved through the fluid and the equilibrium mean squared displacement (MSD). Applying concepts of linear response theory, we generalize Einstein's relation and thereby relate recoil and MSD, which both provide access to the mobility kernel. With increasing concentration, however, MSD and recoil show distinct behaviors, rooted in different behaviors of the two kernels. Using two theoretical models, a linear two-bath particle model and hard spheres treated by mode-coupling theory, we find a Volterra relation between the two kernels, explaining differing timescales in friction and mobility kernels under variation of concentration., Comment: 12 pages, 7 figures

We experimentally investigate the recoil dynamics of a colloidal probe particle after shearing it with constant velocity through a viscoelastic fluid. The recoil displays two distinct timescales which are in excellent agreement with a microscopic model built on a particle being linked to two bath particles by harmonic springs. This model yields analytical expressions which reproduce all experimental protocols, including additional waiting periods before particle release. Notably, two sets of timescales appear, corresponding to reciprocal and nonreciprocal eigenmodes of the model., Comment: 5 pages, 4 figures

We investigate the hopping dynamics of a colloidal particle across a potential barrier and within a viscoelastic, i.e., non-Markovian bath, and report two clearly separated time scales in the corresponding waiting time distributions. While the longer time scale exponentially depends on the barrier height, the shorter one is similar to the relaxation time of the fluid. This short time scale is a signature of the storage and release of elastic energy inside the bath, that strongly increases the hopping rate. Our results are in excellent agreement with numerical simulations of a simple Maxwell model., Comment: 5 pages, 5 figures

The motion of Brownian particles in nonlinear baths, such as, e.g., viscoelastic fluids, is of great interest. We theoretically study a simple model for such bath, where two particles are coupled via a sinusoidal potential. This model, which is an extension of the famous Prandtl Tomlinson model, has been found to reproduce some aspects of recent experiments, such as shear-thinning and position oscillations [J. Chem. Phys. {\bf 154}, 184904 (2021)]. Analyzing this model in detail, we show that the predicted behavior of position oscillations agrees qualitatively with experimentally observed trends; (i) oscillations appear only in a certain regime of velocity and trap stiffness of the confining potential, and (ii), the amplitude and frequency of oscillations increase with driving velocity, the latter in a linear fashion. Increasing the potential barrier height of the model yields a rupture transition as a function of driving velocity, where the system abruptly changes from a mildly driven state to a strongly driven state. The frequency of oscillations scales as $(v_0-v_0^*)^{1/2}$ near the rupture velocity $v_0^*$, found for infinite trap stiffness. Investigating the (micro-)viscosity for different parameter ranges, we note that position oscillations leave their signature by an additional (mild) plateau in the flow curves, suggesting that oscillations influence the micro-viscosity. For a time-modulated driving, the mean friction force of the driven particle shows a pronounced resonance behavior, i.e, it changes strongly as a function of driving frequency. The model has two known limits: For infinite trap stiffness, it can be mapped to diffusion in a tilted periodic potential. For infinite bath friction, the original Prandtl Tomlinson model is recovered. We find that the flow curve of the model (roughly) crosses over between these two limiting cases., Comment: 15 pages, 12 figures

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 captures the observed behavior, and implies 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., Comment: 7 pages, 5 figures

The onset of polar flocking in active matter is discontinuous, akin to gas-liquid phase transitions, except that the steady state exhibits microphase separation into polar clusters. While these features have been observed in theoretical models and experiments, little is known about the underlying mesoscopic processes at the cluster level. Here we show that emergence and maintenance of polar order are governed by the interplay between the assembly and disassembly dynamics of clusters with varying size and degree of polar order. Using agent-based simulations of propelled filaments in a parameter regime relevant for actomyosin motility assays, we monitor the temporal evolution of cluster statistics and the transport processes of filaments between clusters. We find that, over a broad parameter range, the emergence of order is determined by nucleation and growth of polar clusters, where the nucleation threshold depends not only on the cluster size but also on its polar moment. Growth involves cluster self-replication, and polar order is established by cluster growth and fragmentation. Maintenance of the microphase-separated, polar-ordered state results from a cyclic dynamics in cluster size and order, driven by an interplay between cluster nucleation, coagulation, fragmentation and evaporation of single filaments. These findings are corroborated by a kinetic model for the cluster dynamics that includes these elementary cluster-level processes. It consistently reproduces the cluster statistics as well as the cyclic turnover from disordered to ordered clusters and back. Such cyclic kinetic processes could represent a general mechanism for the maintenance of order in active matter systems.

Perturbations of fluid media can give rise to non-equilibrium dynamics, which may in turn cause motion of immersed inclusions. We consider perturbations ("activations") that are local in space and time, of a fluid density which is conserved, and study the resulting diffusiophoretic phenomena that emerge at a large distance. Specifically, we consider cases where the perturbations propagate diffusively, providing examples from passive and active matter for which this is expected to be the case. Activations can, for instance, be realized by sudden and local changes in interaction potentials of the medium, or by local changes of its activity. Various analytical results are provided for the case of confinement by two parallel walls. We investigate the possibility of extracting work from inclusions which are moving through the activated fluid. Further, we show that a time-dependent density profile, created via suitable activation protocols, allows for conveyance of inclusions along controlled and stable trajectories. In contrast, in states with a steady density, inclusions cannot be held at stable positions, reminiscent of Earnshaw's theorem of electrostatics. We expect these findings to be applicable in a range of experimental systems., Comment: 14 pages (12+2), 14 Figures

A colloidal particle is a prominent example of a stochastic system, and, if suspended in a simple viscous liquid, very closely resembles the case of an ideal random walker. A variety of new phenomena have been observed when such colloid is suspended in a viscoelastic fluid instead, for example pronounced nonlinear responses when the viscoelastic bath is driven out of equilibrium. Here, using a micron-sized particle in a micellar solution, we investigate in detail, how these nonlinear bath properties leave their fingerprints already in equilibrium measurements, for the cases where the particle is unconfined or trapped in a harmonic potential. We find that the coefficients in an effective linear (generalized) Langevin equation show intriguing inter-dependencies, which can be shown to arise only in nonlinear baths: For example, the friction memory can depend on the external potential that acts only on the colloidal particle (as recently noted in simulations of molecular tracers in water in Phys. Rev. X 7, 041065 (2017)), it can depend on the mass of the colloid, or, in an overdamped setting, on its bare diffusivity. These inter-dependencies, caused by so-called fluctuation renormalizations, are seen in an exact small time expansion of the friction memory based on microscopic starting points. Using linear response theory, they can be interpreted in terms of microrheological modes of force-controlled or velocity-controlled driving. The mentioned nonlinear markers are observed in our experiments, which are astonishingly well reproduced by a stochastic Prandtl-Tomlinson model mimicking the nonlinear viscoelastic bath. The pronounced nonlinearities seen in our experiments together with the good understanding in a simple theoretical model make this system a promising candidate for exploration of colloidal motion in nonlinear stochastic environments., Comment: 14 pages, 12 figures changes v2: minor changes regarding abstract, and Fig. 11

Microfluidic technologies are commonly used for the manipulation of red blood cell (RBC) suspensions and analyses of flow-mediated biomechanics. To enhance the performance of microfluidic devices, understanding the dynamics of the suspensions processed within is crucial. We report novel aspects of the spatio-temporal dynamics of RBC suspensions flowing through a typical microchannel at low Reynolds number. Through experiments with dilute RBC suspensions, we find an off-centre two-peak (OCTP) profile of cells contrary to the centralised distribution commonly reported for low-inertia flows. This is reminiscent of the well-known "tubular pinch effect" which arises from inertial effects. However, given the conditions of negligible inertia in our experiments, an alternative explanation is needed for this OCTP profile. Our massively-parallel simulations of RBC flow in real-size microfluidic dimensions using the immersed-boundary-lattice-Boltzmann method (IB-LBM) confirm the experimental findings and elucidate the underlying mechanism for the counterintuitive RBC pattern. By analysing the RBC migration and cell-free layer (CFL) development within a high-aspect-ratio channel, we show that such a distribution is co-determined by the spatial decay of hydrodynamic lift and the global deficiency of cell dispersion in dilute suspensions. We find a CFL development length greater than 46 and 28 hydraulic diameters in the experiment and simulation, respectively, exceeding typical lengths of microfluidic designs. Our work highlights the key role of transient cell distribution in dilute suspensions, which may negatively affect the reliability of experimental results if not taken into account., Comment: 19+15 pages, 7+18 figures, including supplementary material

We present a thermodynamically consistent model of a ternary fluid interacting with elastic membranes. Following a free-energy modelling approach and taking into account the thermodynamics laws, we derive the equations governing the ternary fluid flow and dynamics of the membranes. We also provide the numerical framework for simulating such fluid-structure interaction problems. It is based on the lattice Boltzmann method, employed for resolving the evolution equations of the ternary fluid in an Eulerian description, coupled to the immersed boundary method, allowing for the membrane equations of motion to be solved in a Lagrangian system. The configuration of an elastic capsule placed at a fluid-fluid interface is considered for validation purposes. Systematic simulations are performed for a detailed comparison with reference numerical results obtained by Surface Evolver, and the Galilean invariance of the proposed model is also proven. The proposed approach is versatile, and a wide range of geometries can be simulated. To demonstrate this, the problem of a capillary bridge formed between two deformable capsules is investigated here.

Soft particles at fluid interfaces play an important role in many aspects of our daily life, such as the food industry, paints and coatings, and medical applications. Analytical methods are not capable of describing the emergent effects of the complex dynamics of suspensions of many soft particles, whereas experiments typically either only capture bulk properties or require invasive methods. Computational methods are therefore a great tool to complement experimental work. However, an efficient and versatile numerical method is needed to model dense suspensions of many soft particles. In this article we propose a method to simulate soft particles in a multi-component fluid, both at and near fluid-fluid interfaces, based on the lattice Boltzmann method, and characterize the error stemming from the fluid-structure coupling for the particle equilibrium shape when adsorbed onto a fluid-fluid interface. Furthermore, we characterize the influence of the preferential contact angle of the particle surface and the particle softness on the vertical displacement of the center of mass relative to the fluid interface. Finally, we demonstrate the capability of our model by simulating a soft capsule adsorbing onto a fluid-fluid interface with a shear flow parallel to the interface, and the covering of a droplet suspended in another fluid by soft particles with different wettability., Comment: 16 pages, 14 figures

Liquid shells (e.g. double emulsions, vesicles etc.) are susceptible to interfacial instability and rupturing when driven out of mechanical equilibrium. This poses a significant challenge for the design of liquid shell based micro-machines, where the goal is to maintain stability and dynamical control in combination with motility. Here we present our solution to this problem with controllable self-propelling liquid shells, which we have stabilized using the soft topological constraints imposed by a nematogen oil. We demonstrate, through experiments and simulations, that anisotropic elasticity can counterbalance the destabilizing effect of viscous drag induced by shell motility, and inhibit rupturing. We analyze their propulsion dynamics, and identify a peculiar meandering behavior driven by a combination of topological and chemical spontaneously broken symmetries. Based on our understanding of these symmetry breaking mechanisms, we provide routes to control shell motion via topology, chemical signaling and hydrodynamic interactions.

Analytic relations are derived for finite volume integrals over the radial distribution function of a fluid, so-called Kirkwood-Buff integrals. Closed form expressions are obtained for cubes and cuboids, the system shapes commonly employed in molecular simulations. When finite volume Kirkwood-Buff integrals are expanded over inverse system size, the leading term depends on shape only through the surface area to volume ratio. This conjecture is proved for arbitrary shapes and a general expression for the leading term is derived. From this, a new extrapolation to the infinite volume limit is proposed, which converges much faster with system size than previous approximations and thus significantly simplifies the numerical computations.

Microscopic particles suspended in liquids are the prime example of an overdamped system because viscous forces dominate over inertial effects. Apart from their use as model systems, they receive considerable attention as sensitive probes from which forces on molecular scales can be inferred. The interpretation of such experiments rests on the assumption, that, even if the particles are driven, the liquid remains in equilibrium, and all modes are overdamped. Here, we experimentally demonstrate that this is no longer valid when a particle is forced through a viscoelastic fluid. Even at small driving velocities where Stokes law remains valid, we observe particle oscillations with periods up to several tens of seconds. We attribute these to non-equilibrium fluctuations of the fluid, which are excited by the particle's motion. The observed oscillatory dynamics is in quantitative agreement with an overdamped Langevin equation with negative friction-memory term and which is equivalent to the motion of a stochastically driven underdamped oscillator. This fundamentally new oscillatory mode will largely expand the variety of model systems but has also considerable implications on how molecular forces are determined by colloidal probe particles under natural viscoelastic conditions., Comment: Accepted with Nat. Comm. (originally submitted version, complying with Nature policies). 10 pages, 8 figures

Starting from the stochastic equation for the density operator, we formulate the exact (instantaneous) stress tensor for interacting Brownian particles, whose average value agrees with expressions derived previously. We analyze the relation between the stress tensor and forces on external potentials, and observe that, out of equilibrium, particle currents give rise to extra forces. Next, we derive the stress tensor for a Landau-Ginzburg theory in non-equilibrium situations, finding an expression analogous to that of the exact microscopic stress tensor, and discuss the computation of out-of-equilibrium (classical) Casimir forces. We use these relations to study the spatio-temporal correlations of the stress tensor in a Brownian fluid, which we derive exactly to leading order in the interaction potential strength. We observe that, after integration over time, the spatial correlations generally decay as power laws in space. These are expected to be of importance for driven confined systems. We also show that divergence-free parts of the stress tensor do not contribute in the Green-Kubo relation for the viscosity., Comment: 12 pages

The deformability of soft condensed matter often requires modelling of hydrodynamical aspects to gain quantitative understanding. This, however, requires specialised methods that can resolve the multiscale nature of soft matter systems. We review a number of the most popular simulation methods that have emerged, such as Langevin dynamics, dissipative particle dynamics, multi-particle collision dynamics, sometimes also referred to as stochastic rotation dynamics, and the lattice-Boltzmann method. We conclude this review with a short glance at current compute architectures for high-performance computing and community codes for soft matter simulation., Comment: 16 pages, 6 figures

Here we investigate how the positions of a condensed phase can be controlled by using concentration gradients of a regulator that influences phase separation. We consider a mean field model of a ternary mixture where a concentration gradient of a regulator is imposed by an external potential. We show that novel first order phase transition exists at which the position of the condensed phase switches in a discontinuous manner. This mechanism could have implications for the spatial organization of biological cells and provides a control mechanism for droplets in microfluidic systems.

Assuming an effective quadratic Hamiltonian, we derive an approximate, linear stochastic equation of motion for the density-fluctuations in liquids, composed of overdamped Brownian particles. From this approach, time dependent two point correlation functions (such as the intermediate scattering function) are derived. We show that this correlation function is exact at short times, for any interaction and, in particular, for arbitrary external potentials so that it applies to confined systems. Furthermore, we discuss the relation of this approach to previous ones, such as dynamical density functional theory as well as the formally exact treatment. This approach, inspired by the well known Landau-Ginzburg Hamiltonians, and the corresponding "Model B" equation of motion, may be seen as its microscopic version, containing information about the details on the particle level., Comment: 8 pages

In this work we consider the following question: given a mechanical microswimming mechanism, does increased deformability of the swimmer body hinder or promote the motility of the swimmer? To answer this we study a microswimmer model composed of deformable beads connected with springs. We determine the velocity of the swimmer analytically, starting from the forces driving the motion and assuming that the oscillations in the effective radii of the beads are known and are much smaller than the radii themselves. We find that to the lowest order, only the driving frequency mode of the surface oscillations contributes to the swimming velocity, and that this velocity may both rise and fall with the deformability of the beads depending on the spring constant. To test these results, we run immersed boundary lattice Boltzmann simulations of the swimmer, and show that they reproduce both the velocity-promoting and velocity-hindering effects of bead deformability correctly in the predicted parameter ranges. Our results mean that for a general swimmer, its elasticity determines whether passive deformations of the swimmer body, induced by the fluid flow, aid or oppose the motion.

Chemotaxis and auto-chemotaxis play an important role in many essential biological processes. We present a self-propelling artificial swimmer system which exhibits chemotaxis as well as negative auto-chemotaxis. Oil droplets in an aqueous surfactant solution are driven by interfacial Marangoni flows induced by micellar solubilization of the oil phase. We demonstrate that chemotaxis along micellar surfactant gradients can guide these swimmers through a microfluidic maze. Similarly, a depletion of empty micelles in the wake of a droplet swimmer causes negative autochemotaxis and thereby trail avoidance. We have studied autochemotaxis quantitatively in a microfluidic device of bifurcating channels: Branch choices of consecutive swimmers are anticorrelated, an effect decaying over time due to trail dispersion. We have modeled this process by a simple one-dimensional diffusion process and stochastic Langevin dynamics. Our results are consistent with a linear surfactant gradient force and diffusion constants appropriate for micellar diffusion, and provide a measure of autochemotactic feedback strength versus stochastic forces.

We study the second order response functions of a colloidal particle being subjected to an anharmonic potential. Contrary to typical response measurements which require an external perturbation, here we experimentally confirm a recently developed approach where the system's susceptibilities up to second order are obtained from the particle's equilibrium trajectory [PCCP $\mathrm{{\bf 17}}$, 6653 (2015)]. The measured susceptibilities are in quantitative agreement with those obtained from the response to an external perturbation., Comment: 4 figures

We compute fluctuation-induced (Casimir) forces for classical systems after a temperature quench. Using a generic coarse-grained model for fluctuations of a conserved density, we find that transient forces arise even if the initial and final states are force-free. In setups reminiscent of Casimir (planar walls) and van der Waals (small inclusions) interactions, we find comparable exact universal expressions for the force. Dynamical details only scale the time axis of transient force curves. We propose that such quenches can be achieved, for instance, in experiments on active matter, employing tunable activity or interaction protocols., Comment: 5 pages, 2 figures, added equation references/comments

When the motion of a probe strongly disturbs the thermal equilibrium of the solvent or bath, the nonlinear response of the latter must enter the probe's effective evolution equation. We derive that induced stochastic dynamics using second order response around the bath thermal equilibrium. We discuss the nature of the new term in the evolution equation which is not longer purely dissipative, and the appearance of a novel time-scale for the probe related to changes in the dynamical activity of the bath. A major application for the obtained nonlinear generalized Langevin equation is in the study of colloid motion in a visco-elastic medium., Comment: 17 pages in preprint format, no figures

We report curling self-propulsion in aqueous emulsions of common mesogenic compounds. Nematic liquid crystal droplets self-propel in a surfactant solution with concentrations above the critical micelle concentration while undergoing micellar solubilization. We analyzed trajectories both in a Hele-Shaw geometry and in a 3D setup at variable buoyancy. The coupling between the nematic director field and the convective flow inside the droplet leads to a second symmetry breaking which gives rise to curling motion in 2D. This is demonstrated through a reversible transition to non-helical persistent swimming by heating to the isotropic phase. Furthermore, auto-chemotaxis can spontaneously break the inversion symmetry, leading to helical trajectories.

We present a new ternary free energy lattice Boltzmann model. The distinguishing feature of our model is that we are able to analytically derive and independently vary all fluid-fluid surface tensions and the solid surface contact angles. We carry out a number of benchmark tests: (i) double emulsions and liquid lenses to validate the surface tensions, (ii) ternary fluids in contact with a square well to compare the contact angles against analytical predictions, and (iii) ternary phase separation to verify that the multicomponent fluid dynamics is accurately captured. Additionally we also describe how the model here presented here can be extended to include an arbitrary number of fluid components.

The effect of tube diameter $D$ and capillary number $Ca$ on platelet margination in blood flow at $\approx 37\%$ tube haematocrit is investigated. The system is modelled as three-dimensional suspension of deformable red blood cells and nearly rigid platelets using a combination of the lattice-Boltzmann, immersed boundary and finite element methods. Results show that margination is facilitated by a non-diffusive radial platelet transport. This effect is important near the edge of the cell-free layer, but it is only observed for $Ca > 0.2$, when red blood cells are tank-treading rather than tumbling. It is also shown that platelet trapping in the cell-free layer is reversible for $Ca \leq 0.2$. Only for the smallest investigated tube ($D = 10 \mu\text{m}$) margination is essentially independent of $Ca$. Once platelets have reached the cell-free layer, they tend to slide rather than tumble. The tumbling rate is essentially independent of $Ca$ but increases with $D$. Tumbling is suppressed by the strong confinement due to the relatively small cell-free layer thickness at $\approx 37\%$ tube haematocrit., Comment: 16 pages, 10 figures

We investigate the effect of particle volume fraction on the efficiency of deterministic lateral displacement (DLD) devices. DLD is a popular passive sorting technique for microfluidic applications. Yet, it has been designed for treating dilute suspensions, and its efficiency for denser samples is not well known. We perform 3D simulations based on the immersed-boundary, lattice-Boltzmann and finite-element methods to model the flow of red blood cells (RBCs) in different DLD devices. We quantify the DLD efficiency in terms of appropriate "failure" probabilities and RBC counts in designated device outlets. Our main result is that the displacement mode breaks down upon an increase of RBC volume fraction, while the zigzag mode remains relatively robust. This suggests that the separation of larger particles (such as white blood cells) from a dense RBC background is simpler than separating smaller particles (such as platelets) from the same background. The observed breakdown stems from non-deterministic particle collisions interfering with the designed deterministic nature of DLD devices. Therefore, we postulate that dense suspension effects generally hamper efficient particle separation in devices based on deterministic principles., Comment: 21 pages, 9 figures

The viscosity of fluids is generally understood in terms of kinetic mechanisms, i.e., particle collisions, or thermodynamic ones as imposed through structural distortions upon e.g. applying shear. Often the former is less relevant, and (damped) Brownian particles are considered good fluid model systems. We formulate a general theoretical approach for rheology in confinement, based on the many particle diffusion equation, evaluated via classical density functional theory. We discuss the viscosity for the situation of two parallel walls in relative motion as a function of wall-to-wall distance., Comment: 5 pages, 3 figures

It has long been known that some microswimmers seem to swim counter-intuitively faster when the viscosity of the surrounding fluid is increased, whereas others slow down. This conflicting dependence of the swimming velocity on the viscosity is poorly understood theoretically. Here we explain that any mechanical microswimmer with an elastic degree of freedom in a simple Newtonian fluid can exhibit both kinds of response to an increase in the fluid viscosity for different viscosity ranges, if the driving is weak. The velocity response is controlled by a single parameter $\varGamma$, the ratio of the relaxation time of the elastic component of the swimmer in the viscous fluid and the swimming stroke period. This defines two velocity-viscosity regimes, which we characterize using the bead-spring microswimmer model and analyzing the different forces acting on the parts of this swimmer. The analytical calculations are supported by lattice-Boltzmann simulations, which accurately reproduce the two velocity regimes for the predicted values of $\varGamma$.

The fabrication of novel soft materials is an important scientific and technological challenge. We investigate the response of magnetic ellipsoidal particles adsorbed at fluid-fluid interfaces to external magnetic fields. By exploiting previously discovered first-order orientation phase transitions we show how to switch on and off dipolar capillary interactions between particles, leading to the formation of distinctive self-assembled structures and allowing dynamic control of the bottom-up fabrication of reconfigurable novel-structured materials, Comment: 5 pages, 5 Figures

The energy required to detach a single particle from a fluid-fluid interface is an important parameter for designing certain soft materials, for example, emulsions stabilised by colloidal particles, colloidosomes designed for targeted drug delivery, and bio-sensors composed of magnetic particles adsorbed at interfaces. For a fixed particle volume, prolate and oblate spheroids attach more strongly to interfaces because they have larger particle-interface areas. Calculating the detachment energy of spheroids necessitates the difficult measurement of particle-liquid surface tensions, in contrast with spheres, where the contact angle suffices. We develop a simplified detachment energy model for spheroids which depends only on the particle aspect ratio and the height of the particle centre of mass above the fluid-fluid interface. We use lattice Boltzmann simulations to validate the model and provide quantitative evidence that the approach can be applied to simulate particle-stabilized emulsions, and highlight the experimental implications of this validation., Comment: 9 pages, 8 figures

We show, via three-dimensional immersed-boundary-finite-element-lattice-Boltzmann simulations, that deformability-based red blood cell (RBC) separation in deterministic lateral displacement (DLD) devices is possible. This is due to the deformability-dependent lateral extension of RBCs and enables us to predict a priori which RBCs will be displaced in a given DLD geometry. Several diseases affect the deformability of human cells. Malaria-infected RBCs or sickle cells, for example, tend to become stiffer than their healthy counterparts. It is therefore desirable to design microfluidic devices which can detect those diseases based on the cells' deformability fingerprint, rather than preparing samples using expensive and time-consuming biochemical preparation steps. Our findings should be helpful in the development of new methods for sorting cells and particles by deformability., Comment: 10 pages, 10 figures

Manufacturing new soft materials with specific optical, mechanical and magnetic properties is a significant challenge. Assembling and manipulating colloidal particles at fluid interfaces is a promising way to make such materials. We use lattice-Boltzmann simulations to investigate the response of magnetic ellipsoidal particles adsorbed at liquid-liquid interfaces to external magnetic fields. We provide further evidence for the first-order orientation phase transition predicted by Bresme and Faraudo [Journal of Physics: Condensed Matter 19 (2007), 375110]. We show that capillary interface deformations around the ellipsoidal particle significantly affect the tilt-angle of the particle for a given dipole-field strength, altering the properties of the orientation transition. We propose scaling laws governing this transition, and suggest how to use these deformations to facilitate particle assembly at fluid-fluid interfaces., Comment: 7 pages, 8 figures

We discuss two topics that we have encountered in our lattice-Boltzmann simulations of complex fluids: the sizes of droplets in particle-stabilised emulsions and deformable particles in fluid flow. The common factor in these seemingly disparate subjects is that both represent an opportunity for a novel application of the Hoshen-Kopelman algorithm. This algorithm is based on detecting connected clusters on a lattice and labelling the involved lattice sites such that all sites that are connected share the same label. The assumption of the presence of a lattice makes it a convenient algorithm to use in combination with lattice-Boltzmann simulations. In order to apply this algorithm on the fly during massively parallel simulations, it needs to be parallelised as well. We present our parallel implementation, which is tailored to a common parallelisation scheme for the lattice-Boltzmann method, and compare it to previous work. We then briefly discuss some examples of results obtained using this procedure., Comment: 11 pages, 7 figures

We analyze fluctuations of particle displacements and stresses in a sheared athermal suspension of elastic capsules (red blood cells). Upon variation of the volume fraction from the dilute up to the highly concentrated regime, our numerical simulations reveal different characteristic power-law regimes of the fluctuation variances and relaxation times. In the jammed phase and at high shear rates, anomalous scaling exponents are found that deviate from pure dimensional predictions. The observed behavior is rationalized via kinetic arguments and a dissipation balance model that takes into account the local fluid flows between the particles. Our findings support the view that the rheology of dense suspensions is essentially governed by the non-affine displacements., Comment: 6 pages, 6 figures

We theoretically investigate general properties of driven (sheared) colloidal suspensions in confinement, based on methods of classical density functional theory. In the absence of an exact closed (Smoluchowski-) equation for the one-particle density under shear, we formulate a set of general conditions for approximations, and show that a simple closure fulfills them. The exact microscopic stress tensor is identified. Exemplifying the situation near a wall (oriented parallel to the direction of shear), we note that the microscopic shear stress is not necessarily homogeneous. Formulating a second equation additional to the Smoluchowski equation, we achieve a homogeneous shear stress, and thereby compute the local flow velocity near the wall. This finally leads to a slip length of the complex fluid at the wall., Comment: 11 pages, 8 figures

We study the shearing rheology of dense suspensions of elastic capsules, taking aggregation-free red blood cells as a physiologically relevant example. Particles are non-Brownian and interact only via hydrodynamics and short-range repulsive forces. An analysis of the different stress mechanisms in the suspension shows that the viscosity is governed by the shear elasticity of the capsules, whereas the repulsive forces are subdominant. Evidence for a dynamic yield stress above a critical volume fraction is provided and related to the elastic properties of the capsules. The shear stress is found to follow a critical jamming scenario and is rather insensitive to the tumbling-to-tank-treading transition. The particle pressure and normal stress differences display some sensitivity to the dynamical state of the cells and exhibit a characteristic scaling, following the behavior of a single particle, in the tank-treading regime. The behavior of the viscosity in the fluid phase is rationalized in terms of effective medium models. Furthermore, the role of confinement effects, which increase the overall magnitude and enhance the shear-thinning of the viscosity, is discussed., Comment: 14 pages, 11 figures; as published. Soft Matter, 2014

The interplay of inertia and deformability has a substantial impact on the transport of soft particles suspended in a fluid. However, to date a thorough understanding of these systems is still missing and only a limited number of experimental and theoretical studies is available. We combine the finite-element, immersed-boundary and lattice-Boltzmann methods to simulate three-dimensional suspensions of soft particles subjected to planar Poiseuille flow at finite Reynolds numbers. Our findings confirm that the particle deformation and inclination increase when inertia is present. We observe that the Segr\'e-Silberberg effect is suppressed with respect to the particle deformability. Depending on the deformability and strength of inertial effects, inward or outward lateral migration of the particles takes place. In particular, for increasing Reynolds numbers and strongly deformable particles, a hitherto unreported distinct flow focusing effect emerges which is accompanied by a non-monotonic behaviour of the apparent suspension viscosity and thickness of the particle-free layer close to the channel walls. This effect can be explained by the behaviour of a single particle and the change of the particle collision mechanism when both deformability and inertia effects are relevant., Comment: 20 pages, 9 figures

Via computer simulations, we provide evidence that the shear rate induced red blood cell tumbling-to-tank-treading transition also occurs at quite high volume fractions, where collective effects are important. The transition takes place as the ratio of suspension stress to the characteristic cell membrane stress exceeds a certain value, independent of volume fraction and cell deformability. This value coincides with that for a transition from an orientationally less ordered to a highly ordered phase. The average cell deformation does not show any signature of the transition, but rather follows a simple scaling law independent of volume fraction., Comment: 8 pages, 9 figures, as published

The history dependence of the glasses formed from flow-melted steady states by a sudden cessation of the shear rate $\dot\gamma$ is studied in colloidal suspensions, by molecular dynamics simulations, and mode-coupling theory. In an ideal glass, stresses relax only partially, leaving behind a finite persistent residual stress. For intermediate times, relaxation curves scale as a function of $\dot\gamma t$, even though no flow is present. The macroscopic stress evolution is connected to a length scale of residual liquefaction displayed by microscopic mean-squared displacements. The theory describes this history dependence of glasses sharing the same thermodynamic state variables, but differing static properties., Comment: submitted to Physical Review