No Arabic abstract
We develop efficient numerical methods for performing many-body Brownian dynamics simulations of a recently-observed fingering instability in an active suspension of colloidal rollers sedimented above a wall [M. Driscoll, B. Delmotte, M. Youssef, S. Sacanna, A. Donev and P. Chaikin, Nature Physics, 2016, doi:10.1038/nphys3970]. We present a stochastic Adams-Bashforth integrator for the equations of Brownian dynamics, which has the same cost as but is more accurate than the widely-used Euler-Maruyama scheme, and uses a random finite difference to capture the stochastic drift proportional to the divergence of the configuration-dependent mobility matrix. We generate the Brownian increments using a Krylov method, and show that for particles confined to remain in the vicinity of a no-slip wall by gravity or active flows the number of iterations is independent of the number of particles. Our numerical experiments with active rollers show that the thermal fluctuations set the characteristic height of the colloids above the wall, both in the initial condition and the subsequent evolution dominated by active flows. The characteristic height in turn controls the timescale and wavelength for the development of the fingering instability.
We perform detailed computational and experimental measurements of the driven dynamics of a dense, uniform suspension of sedimented microrollers driven by a magnetic field rotating around an axis parallel to the floor. We develop a lubrication-corrected Brownian Dynamics method for dense suspensions of driven colloids sedimented above a bottom wall. The numerical method adds lubrication friction between nearby pairs of particles, as well as particles and the bottom wall, to a minimally-resolved model of the far-field hydrodynamic interactions. Our experiments combine fluorescent labeling with particle tracking to trace the trajectories of individual particles in a dense suspension, and to measure their propulsion velocities. Previous computational studies [B. Sprinkle et al., J. Chem. Phys., 147, 244103, 2017] predicted that at sufficiently high densities a uniform suspension of microrollers separates into two layers, a slow monolayer right above the wall, and a fast layer on top of the bottom layer. Here we verify this prediction, showing good quantitative agreement between the bimodal distribution of particle velocities predicted by the lubrication-corrected Brownian Dynamics and those measured in the experiments. The computational method accurately predicts the rate at which particles are observed to switch between the slow and fast layers in the experiments. We also use our numerical method to demonstrate the important role that pairwise lubrication plays in motility-induced phase separation in dense monolayers of colloidal microrollers, as recently suggested for suspensions of Quincke rollers [D. Geyer et al., Physical Review X, 9(3), 031043, 2019].
We introduce methods for large scale Brownian Dynamics (BD) simulation of many rigid particles of arbitrary shape suspended in a fluctuating fluid. Our method adds Brownian motion to the rigid multiblob method at a cost comparable to the cost of deterministic simulations. We demonstrate that we can efficiently generate deterministic and random displacements for many particles using preconditioned Krylov iterative methods, if kernel methods to efficiently compute the action of the Rotne-Prager-Yamakawa (RPY) mobility matrix and it square root are available for the given boundary conditions. We address a major challenge in large-scale BD simulations, capturing the stochastic drift term that arises because of the configuration-dependent mobility. Unlike the widely-used Fixman midpoint scheme, our methods utilize random finite differences and do not require the solution of resistance problems or the computation of the action of the inverse square root of the RPY mobility matrix. We construct two temporal schemes which are viable for large scale simulations, an Euler-Maruyama traction scheme and a Trapezoidal Slip scheme, which minimize the number of mobility solves per time step while capturing the required stochastic drift terms. We validate and compare these schemes numerically by modeling suspensions of boomerang shaped particles sedimented near a bottom wall. Using the trapezoidal scheme, we investigate the steady-state active motion in a dense suspensions of confined microrollers, whose height above the wall is set by a combination of thermal noise and active flows. We find the existence of two populations of active particles, slower ones closer to the bottom and faster ones above them, and demonstrate that our method provides quantitative accuracy even with relatively coarse resolutions of the particle geometry.
We derive an analytic expression for the mechanical pressure of a generic one-dimensional model of confined active Brownian particles (ABPs) that is valid for all values of Peclet number Pe and all confining scenarios. Our model reproduces the known scaling of bulk pressure with Pe^2 while in strong confinement pressure scales with Pe. Our analytic results are very well reproduced by simulations of ABPs in 2D. We use the pressure formula to calculate both the work performed by an active engine and its efficiency. In particular, efficiency is maximized for work cycles with finite period and not in the limit of infinitely slow cycles as in thermodynamic engines.
Self-propelled colloids constitute an important class of intrinsically non-equilibrium matter. Typically, such a particle moves ballistically at short times, but eventually changes its orientation, and displays random-walk behavior in the long-time limit. Theory predicts that if the velocity of non-interacting swimmers varies spatially in 1D, $v(x)$, then their density $rho(x)$ satisfies $rho(x) = rho(0)v(0)/v(x)$, where $x = 0$ is an arbitrary reference point. Such a dependence of steady-state $rho(x)$ on the particle dynamics, which was the qualitative basis of recent work demonstrating how to `paint with bacteria, is forbidden in thermal equilibrium. We verify this prediction quantitatively by constructing bacteria that swim with an intensity-dependent speed when illuminated. A spatial light pattern therefore creates a speed profile, along which we find that, indeed, $rho(x)v(x) = mathrm{constant}$, provided that steady state is reached.
We derive equations of motion for the mean-squared displacement (MSD) of an active Brownian particle (ABP) in a crowded environment modeled by a dense system of passive Brownian particles, and of a passive tracer particle in a dense active-Brownian particle system, using a projection-operator scheme. The interaction of the tracer particle with the dense host environment gives rise to strong memory effects. Evaluating these approximately in the framework of a recently developed mode-coupling theory for the glass transition in active Brownian particles (ABP-MCT), we discuss the various regimes of activity-induced super-diffusive motion and density-induced sub-diffusive motion. The predictions of the theory are shown to be in good agreement with results from an event-driven Brownian dynamics simulation scheme for the dynamics of two-dimensional active Brownian hard disks.