Tracer particles immersed in suspensions of biological microswimmers such as E. coli or Chlamydomonas display phenomena unseen in conventional equilibrium systems, including strongly enhanced diffusivity relative to the Brownian value and non-Gaussia
n displacement statistics. In dilute, 3-dimensional suspensions, these phenomena have typically been explained by the hydrodynamic advection of point tracers by isolated microswimmers, while, at higher concentrations, correlations between pusher microswimmers such as E. coli can increase the effective diffusivity even further. Anisotropic tracers in active suspensions can be expected to exhibit even more complex behaviour than spherical ones, due to the presence of a nontrivial translation-rotation coupling. Using large-scale lattice Boltzmann simulations of model microswimmers described by extended force dipoles, we study the motion of ellipsoidal point tracers immersed in 3-dimensional microswimmer suspensions. We find that the rotational diffusivity of tracers is much less affected by swimmer-swimmer correlations than the translational diffusivity. We furthermore study the anisotropic translational diffusion in the particle frame and find that, in pusher suspensions, the diffusivity along the ellipsoid major axis is higher than in the direction perpendicular to it, albeit with a smaller ratio than for Brownian diffusion. Thus, we find that far field hydrodynamics cannot account for the anomalous coupling between translation and rotation observed in experiments, as was recently proposed. Finally, we study the probability distributions (PDFs) of translational and rotational displacements. In accordance with experimental observations, for short observation times we observe strongly non-Gaussian PDFs that collapse when rescaled with their variance, which we attribute to the ballistic nature of tracer motion at short times.
The Derjaguin approximation (DA) relates the force between curved surfaces to the interaction free energy between parallel planes. It is typically derived by considering the direct interaction between the bodies involved, thus treating the effect of
an intervening solvent implicitly by a rescaling of the corresponding Hamaker constant. Here, we provide a generalization of DA to the case of a molecular medium between the bodies, as is the case in most applications. The derivation is based on an explicit statistical-mechanical treatment of the contribution to the interaction force from a molecular solvent using a general expression for intermolecular and molecule-surface interactions. Starting from an exact expression for the force, DA is arrived at by a series of well-defined approximations. Our results show that DA remains valid in a molecular solvent as long as (i) the surface-molecule interactions are of much shorter range than the radius R of the sphere and (ii) the density correlation length in the solvent is smaller than R. We then extend our analysis to the case where a phase transition occurs between the surfaces, which cannot easily be covered using a statistical-mechanical formalism due to the discontinuous change in the density of the medium. Instead using a continuum thermodynamic description, we show that this phase transformation induces an attractive force between the bodies, and that the force between curved surfaces can be related to the free energy in the corresponding planar case, in accordance with DA.
Motility-induced phase separation (MIPS) leads to cohesive active matter in the absence of cohesive forces. We present, extend and illustrate a recent generalized thermodynamic formalism which accounts for its binodal curve. Using this formalism, we
identify both a generalized surface tension, that controls finite-size corrections to coexisting densities, and generalized forces, that can be used to construct new thermodynamic ensembles. Our framework is based on a nonequilibrium generalization of the Cahn-Hilliard equation and we discuss its application to active particles interacting either via quorum-sensing interactions or directly through pairwise forces.
In this Letter, we study the collective behaviour of a large number of self-propelled microswimmers immersed in a fluid. Using unprecedently large-scale lattice Boltzmann simulations, we reproduce the transition to bacterial turbulence. We show that,
even well below the transition, swimmers move in a correlated fashion that cannot be described by a mean-field approach. We develop a novel kinetic theory that captures these correlations and is non-perturbative in the swimmer density. To provide an experimentally accessible measure of correlations, we calculate the diffusivity of passive tracers and reveal its non-trivial density dependence. The theory is in quantitative agreement with the lattice Boltzmann simulations and captures the asymmetry between pusher and puller swimmers below the transition to turbulence.
Hydrodynamic interactions in systems comprised of self-propelled particles, such as swimming microorganisms, and passive tracers have a significant impact on the tracer dynamics compared to the equivalent dry sample. However, such interactions are of
ten difficult to take into account in simulations due to their computational cost. Here, we perform a systematic investigation of swimmer-tracer interaction using an efficient force/counter-force based lattice-Boltzmann (LB) algorithm [J. de Graaf~textit{et al.}, J. Chem. Phys.~textbf{144}, 134106 (2016)] in order to validate its ability to capture the relevant low-Reynolds-number physics. We show that the LB algorithm reproduces far-field theoretical results well, both in a system with periodic boundary conditions and in a spherical cavity with no-slip walls, for which we derive expressions here. The force-lattice coupling of the LB algorithm leads to a smearing out of the flow field, which strongly perturbs the tracer trajectories at close swimmer-tracer separations, and we analyze how this effect can be accurately captured using a simple renormalized hydrodynamic theory. Finally, we show that care must be taken when using LB algorithms to simulate systems of self-propelled particles, since its finite momentum transport time can lead to significant deviations from theoretical predictions based on Stokes flow. These insights should prove relevant to the future study of large-scale microswimmer suspensions using these methods.
The interaction between swimming microorganisms or artificial self-propelled colloids and passive (tracer) particles in a fluid leads to enhanced diffusion of the tracers. This enhancement has attracted strong interest, as it could lead to new strate
gies to tackle the difficult problem of mixing on a microfluidic scale. Most of the theoretical work on this topic has focused on hydrodynamic interactions between the tracers and swimmers in a bulk fluid. However, in simulations, periodic boundary conditions (PBCs) are often imposed on the sample and the fluid. Here, we theoretically analyze the effect of PBCs on the hydrodynamic interactions between tracer particles and microswimmers. We formulate an Ewald sum for the leading-order stresslet singularity produced by a swimmer to probe the effect of PBCs on tracer trajectories. We find that introducing periodicity into the system has a surprisingly significant effect, even for relatively small swimmer-tracer separations. We also find that the bulk limit is only reached for very large system sizes, which are challenging to simulate with most hydrodynamic solvers.
We investigate the phase behavior and kinetics of a monodisperse mixture of active (textit{i.e.}, self-propelled) and passive isometric Brownian particles through Brownian dynamics simulations and theory. As in a purely active system, motility of the
active component triggers phase separation into a dense and a dilute phase; in the dense phase we further find active-passive segregation, with rafts of passive particles in a sea of active particles. We find that phase separation from an initially disordered mixture can occur with as little as 15 percent of the particles being active. Finally, we show that a system prepared in a suitable fully segregated initial state reproducibly self-assembles an active corona which triggers crystallization of the passive core by initiating a compression wave. Our findings are relevant to the experimental pursuit of directed self-assembly using active particles.
