Search

Collective diffusion, characterised by the collective diffusion coefficient $D_\mathrm{coll}$, is a key quantity for describing the macroscopic transport properties of soft matter systems. However, measuring $D_\mathrm{coll}$ is a fundamental experimental and numerical challenge, as it either relies on nonequilibrium techniques that are hard to interpret or on Fourier-based approaches at equilibrium which are fraught with difficulties associated with Fourier transforms. In this work, we present a novel approach to measure collective diffusion properties by analysing the statistics of particle number counts $N(t)$ in virtual observation boxes of an image at equilibrium, a method we term the "Countoscope". By investigating the equilibrium diffusive dynamics of a 2D colloidal suspension experimentally and numerically, we demonstrate this method can accurately measure $D_\mathrm{coll}$. We validate our results against Fourier-based approaches and establish best practices for measuring $D_\mathrm{coll}$ using fluctuating counts. Remarkably, Fourier techniques struggle with long-range collective measurements because of the non-periodic nature of an experimental image, yet counting fully exploits this property by deliberately using finite observation windows. Finally, we discuss the potential of our method to advance our understanding of collective properties in suspensions, particularly the role of hydrodynamic interactions.

Driven by physical questions pertaining to quantifying particle dynamics, microscopy can now resolve complex systems at the single particle level, from cellular organisms to individual ions. Yet, available analysis techniques face challenges reconstructing trajectories in dense and heterogeneous systems where accurately labelling particles is difficult. Furthermore, the inescapable finite field of view of experiments hinders the measurement of collective effects. Inspired by Smoluchowski, we introduce a broadly applicable analysis technique that probes dynamics of interacting particle suspensions based on a remarkably simple principle: counting particles in finite observation boxes. Using colloidal experiments, advanced simulations and theory, we first demonstrate that statistical properties of fluctuating counts can be used to determine self-diffusion coefficients, so alleviating the hurdles associated with trajectory reconstruction. We also provide a recipe for practically extracting the diffusion coefficient from experimental data at variable particle densities, which is sensitive to steric and hydrodynamic interactions. Remarkably, by increasing the observation box size, counting naturally enables the study of collective dynamics in dense suspensions. Using our novel analysis of particle counts, we uncover a surprising enhancement of collective behaviour, as well as a new length scale associated with hyperuniform-like structure. Our counting framework, the Countoscope, thus enables efficient measurements of self and collective dynamics in dense suspensions and opens the way to quantifying dynamics and identifying novel physical mechanisms in diverse complex systems where single particles can be resolved.

Semiflexible slender filaments are ubiquitous in nature and cell biology, including in the cytoskeleton, where reorganization of actin filaments allows the cell to move and divide. Most methods for simulating semiflexible inextensible fibers/polymers are based on discrete (bead-link or blob-link) models, which become prohibitively expensive in the slender limit when hydrodynamics is accounted for. In this paper, we develop a novel coarse-grained approach for simulating fluctuating slender filaments with hydrodynamic interactions. Our approach is tailored to relatively stiff fibers whose persistence length is comparable to or larger than their length, and is based on three major contributions. First, we discretize the filament centerline using a coarse non-uniform Chebyshev grid, on which we formulate a discrete constrained Gibbs-Boltzmann equilibrium distribution and overdamped Langevin equation. Second, we define the hydrodynamic mobility at each point on the filament as an integral of the Rotne-Prager-Yamakawa kernel along the centerline, and apply a spectrally-accurate quadrature to accurately resolve the hydrodynamics. Third, we propose a novel midpoint temporal integrator which can correctly capture the Ito drift terms that arise in the overdamped Langevin equation. We verify that the equilibrium distribution for the Chebyshev grid is a good approximation of the blob-link one, and that our temporal integrator samples the equilibrium distribution for sufficiently small time steps. We also study the dynamics of relaxation of an initially straight filament, and find that as few as 12 Chebyshev nodes provides a good approximation to the dynamics while allowing a time step size two orders of magnitude larger than a resolved blob-link simulation. We conclude by studying how bending fluctuations aid the process of bundling in cross-linked networks of semiflexible fibers.

We develop a linearly-scaling variant of the Force Coupling Method [K. Yeo and M. R. Maxey, J. Fluid Mech. 649, 205-231 (2010)] for computing hydrodynamic interactions among particles confined to a doubly-periodic geometry with either a single bottom wall or two walls (slit channel) in the aperiodic direction. Our spectrally-accurate Stokes solver uses the Fast Fourier Transform (FFT) in the periodic $xy$ plane and Chebyshev polynomials in the aperiodic $z$ direction normal to the wall(s). We decompose the problem into two problems. The first is a doubly-periodic subproblem in the presence of particles (source terms) with free-space boundary conditions in the $z$ direction, which we solve by borrowing ideas from a recent method for rapid evaluation of electrostatic interactions in doubly-periodic geometries [O. Maxian, R. P. Pel\'aez, L. Greengard and A. Donev, J. Chem. Phys. 154, 204107 (2021)]. The second is a correction subproblem to impose the boundary conditions on the wall(s). Instead of the traditional Gaussian kernel, we use the exponential of a semicircle kernel to model the source terms (body force) due to the presence of particles, and provide optimum values for the kernel parameters that ensure a given hydrodynamic radius with at least two digits of accuracy and rotational and translational invariance. The computation time of our solver, which is implemented in graphical processing units, scales linearly with the number of particles, and allows computations with about a million particles in less than a second for a sedimented layer of colloidal microrollers. We find that in a slit channel, a driven dense suspension of microrollers maintains the same two-layer structure as above a single wall, but moves at a substantially lower collective speed due to increased confinement.

In swimming microorganisms and the cell cytoskeleton, inextensible fibers resist bending and twisting, and interact with the surrounding fluid to cause or resist large-scale fluid motion. In this paper, we develop a novel numerical method for the simulation of cylindrical fibers by extending our previous work on inextensible bending fibers [Maxian et al., Phys. Rev. Fluids 6 (1), 014102] to fibers with twist elasticity. In our "Euler" model, twist is a scalar function that measures the deviation of the fiber cross section relative to a twist-free frame, the fiber exerts only torque parallel to the centerline on the fluid, and the perpendicular components of the rotational fluid velocity are discarded in favor of the translational velocity. In the first part of this paper, we justify this model by comparing it to another commonly-used "Kirchhoff" formulation where the fiber exerts both perpendicular and parallel torque on the fluid, and the perpendicular angular fluid velocity is required to be consistent with the translational fluid velocity. We then develop a spectral numerical method for the hydrodynamics of the Euler model. We define hydrodynamic mobility operators using integrals of the Rotne-Prager-Yamakawa tensor, and evaluate these integrals through a novel slender-body quadrature, which requires on the order of 10 points along the fiber to obtain several digits of accuracy. We demonstrate that this choice of mobility removes the unphysical negative eigenvalues in the translation-translation mobility associated with asymptotic slender body theories, and ensures strong convergence of the fiber velocity and weak convergence of the fiber constraint forces. We pair the spatial discretization with a semi-implicit temporal integrator to confirm the negligible contribution of twist elasticity to the relaxation dynamics of a bent fiber and study the instability of a twirling fiber., Comment: 76 pages, 13 figures. See https://github.com/stochasticHydroTools/SlenderBody for MATLAB codes that accompany this paper and https://cims.nyu.edu/~om759/TwistVideos/index.html for supplementary animations

We study the driven collective dynamics of a colloidal monolayer sedimentating down an inclined plane. The action of the gravity force parallel to the bottom wall creates a flow around each colloid, and the hydrodynamic interactions among the colloids accelerate the sedimentation as the local density increases. This leads to the creation of a universal "triangular" inhomogeneous density profile, with a traveling density shock at the leading front moving in the downhill direction. Unlike density shocks in a colloidal monolayer driven by applied torques rather than forces [Phys. Rev. Fluids, 2(9):092301, 2017], the density front during sedimentation remains stable over long periods of time even though it develops a roughness on the order of tens of particle diameters. Through experimental measurements and particle-based computer simulations, we find that the Burgers equation can model the density profile along the sedimentation direction as a function of time remarkably well, with a modest improvement if the nonlinear conservation law accounts for the sub-linear dependence of the collective sedimentation velocity on density.

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]., Comment: Submitted to Soft Matter. See https://cims.nyu.edu/~donev/Publications.html for supplementary material animations

We introduce a Rigid-Body Fluctuating Immersed Boundary (RB-FIB) method to perform large-scale Brownian dynamics simulations of suspensions of rigid particles in fully confined domains, without any need to explicitly construct Green's functions or mobility operators. In the RB-FIB approach, discretized fluctuating Stokes equations are solved with prescribed boundary conditions in conjunction with a rigid-body immersed boundary method to discretize arbitrarily-shaped colloidal particles with no-slip or active-slip prescribed on their surface. We design a specialized Split--Euler--Maruyama temporal integrator that uses a combination of random finite differences to capture the stochastic drift appearing in the overdamped Langevin equation. The RB-FIB method presented in this work only solves mobility problems in each time step using a preconditioned iterative solver, and has a computational complexity that scales linearly in the number of particles and fluid grid cells. We demonstrate that the RB-FIB method correctly reproduces the Gibbs-Boltzmann equilibrium distribution, and use the method to examine the time correlation functions for two spheres tightly confined in a cuboid. We model a quasi--two-dimensional colloidal crystal confined in a narrow microchannel and hydrodynamically driven across a commensurate periodic substrate potential mimicking the effect of a corrugated wall. We observe partial and full depinning of the colloidal monolayer from the substrate potential above a certain wall speed, consistent with a transition from static to kinetic friction through propagating kink solitons. Unexpectedly, we find that particles nearest the boundaries of the domain are the first to be displaced, followed by particles in the middle of the domain.