Active Darcy's Law.

While bacterial swarms can exhibit active turbulence in vacant spaces, they naturally inhabit crowded environments. We numerically show that driving disorderly active fluids through porous media enhances Darcy's law. While purely active flows average to zero flux, hybrid active/driven flows display greater drift than purely pressure-driven flows. This enhancement is nonmonotonic with activity, leading to an optimal activity to maximize flow rate. We incorporate the active contribution into an active Darcy's law, which may serve to help understand anomalous transport of swarming in porous media.

Lock-key microfluidics: simulating nematic colloid advection along wavy-walled channels.

Liquid crystalline media mediate interactions between suspended particles and confining geometries, which not only has potential to guide patterning and bottom-up colloidal assembly, but can also control colloidal migration in microfluidic devices. However, simulating such dynamics is challenging because nemato-elasticity, diffusivity and hydrodynamic interactions must all be accounted for within complex boundaries. We model the advection of colloids dispersed in flowing and fluctuating nematic fluids confined within 2D wavy channels. A lock-key mechanism between homeotropic colloids and troughs is found to be stronger for planar anchoring on the wavy walls compared to homeotropic anchoring on the wavy walls due to the relative location of the colloid-associated defects. Sufficiently large amplitudes result in stick-slip trajectories and even permanent locking of colloids in place. These results demonstrate that wavy walls not only have potential to direct colloids to specific docking sites but also to control site-specific resting duration and intermittent elution.

The type IV pilus chemoreceptor PilJ controls chemotaxis of one bacterial species towards another.

Bacteria live in social communities, where the ability to sense and respond to interspecies and environmental signals is critical for survival. We previously showed the pathogen Pseudomonas aeruginosa detects secreted peptides from bacterial competitors and navigates through interspecies signal gradients using pilus-based motility. Yet, it was unknown whether P. aeruginosa utilizes a designated chemosensory system for this behavior. Here, we performed a systematic genetic analysis of a putative pilus chemosensory system, followed by high-speed live-imaging and single-cell tracking, to reveal behaviors of mutants that retain motility but are blind to interspecies signals. The enzymes predicted to methylate (PilK) and demethylate (ChpB) the putative pilus chemoreceptor, PilJ, are necessary for cells to control the direction of migration. While these findings implicate PilJ as a bona fide chemoreceptor, such function had yet to be experimentally defined, as full-length PilJ is essential for motility. Thus, we constructed systematic genetic modifications of PilJ and found that without the predicted ligand binding domains or predicted methylation sites, cells lose the ability to detect competitor gradients, despite retaining pilus-mediated motility. Chemotaxis trajectory analysis revealed that increased probability and size of P. aeruginosa pilus-mediated steps towards S. aureus peptides, versus steps away, determines motility bias in wild type cells. However, PilJ mutants blind to interspecies signals take less frequent steps towards S. aureus or steps of equal size towards and away. Collectively, this work uncovers the chemosensory nature of PilJ, provides insight into how cell movements are biased during pilus-based chemotaxis, and identifies chemotactic interactions necessary for bacterial survival in polymicrobial communities, revealing putative pathways where therapeutic intervention might disrupt bacterial communication.

Anisotropic run-and-tumble-turn dynamics.

Run-and-tumble processes successfully model several living systems. While studies have typically focused on particles with isotropic tumbles, recent examples exhibit "tumble-turns", in which particles undergo 90° tumbles and so possess explicitly anisotropic dynamics. We study the consequences of such tumble-turn anisotropicity at both short and long-time scales. We model run-and-tumble-turn particles as self-propelled particles subjected to an angular potential that favors directions of movement parallel to Cartesian axes. Using agent-based simulations, we study the effects of the interplay between rotational diffusion and an aligning potential on the particles' trajectories, which leads to the right-angled turns. We demonstrate that the long-time effect is to alter the tumble-turn time, which governs the long-time dynamics. In particular, when normalized by this timescale, trajectories become independent of the underlying details of the potential. As such, we develop a simplified continuum theory, which quantitatively agrees with agent-based simulations. We find that the purely diffusive hydrodynamic limit still exhibits anisotropic features at intermediate times and conclude that the transition to diffusive dynamics precedes the transition to isotropic dynamics. By considering short-range repulsive and alignment particle-particle interactions, we show how the anisotropic features of a single particle are inherited by the global order of the system. We hope this work will shed light on how active systems can extend local anisotropic properties to macroscopic scales, which might be important in biological processes occurring in anisotropic environments.

Complex-tensor theory of simple smectics.

Matter self-assembling into layers generates unique properties, including structures of stacked surfaces, directed transport, and compact area maximization that can be highly functionalized in biology and technology. Smectics represent the paradigm of such lamellar materials - they are a state between fluids and solids, characterized by both orientational and partial positional ordering in one layering direction, making them notoriously difficult to model, particularly in confining geometries. We propose a complex tensor order parameter to describe the local degree of lamellar ordering, layer displacement and orientation of the layers for simple, lamellar smectics. The theory accounts for both dislocations and disclinations, by regularizing singularities within defect cores and so remaining continuous everywhere. The ability to describe disclinations and dislocation allows this theory to simulate arrested configurations and inclusion-induced local ordering. This tensorial theory for simple smectics considerably simplifies numerics, facilitating studies on the mesoscopic structure of topologically complex systems.

Fluctuation-induced dynamics of nematic topological defects.

Topological defects are increasingly being identified in various biological systems, where their characteristic flow fields and stress patterns are associated with continuous active stress generation by biological entities. Here, using numerical simulations of continuum fluctuating nematohydrodynamics, we show that even in the absence of any specific form of active stresses associated with self-propulsion, mesoscopic fluctuations in either orientational alignment or hydrodynamics can independently result in flow patterns around topological defects that resemble the ones observed in active systems. Our simulations further show the possibility of extensile- and contractile-like motion of fluctuation-induced positive half-integer topological defects. Remarkably, isotropic stress fields also reproduce the experimentally measured stress patterns around topological defects in epithelia. Our findings further reveal that extensile- or contractile-like flow and stress patterns around fluctuation-induced defects are governed by passive elastic stresses and flow-aligning behavior of the nematics.

Helical flow states in active nematics.

We show that confining extensile nematics in three-dimensional (3D) channels leads to the emergence of two self-organized flow states with nonzero helicity. The first is a pair of braided antiparallel streams-this double helix occurs when the activity is moderate, anchoring negligible, and reduced temperature high. The second consists of axially aligned counter-rotating vortices-this grinder train arises between spontaneous axial streaming and the vortex lattice. These two unanticipated helical flow states illustrate the potential of active fluids to break symmetries and form complex but organized spatiotemporal structures in 3D fluidic devices.

Mesoscopic simulations of active nematics.

Coarse-grained, mesoscale simulations are invaluable for studying soft condensed matter because of their ability to model systems in which a background solvent plays a substantial role but is not the primary interest. Such methods generally model passive solvents; however, far-from-equilibrium systems may also be composed of complex solutes suspended in an active fluid. Yet, few coarse-grained simulation methods exist to model an active medium. We introduce an algorithm to simulate active nematics, which builds on multiparticle collision dynamics (MPCD) for passive fluctuating nematohydrodynamics by introducing dipolar activity in the local collision operator. Active nematic MPCD (AN-MPCD) simulations not only exhibit the key characteristics of active nematic turbulence but, as a particle-based algorithm, also reproduce crucial attributes of active particle models. Thus, mesoscopic AN-MPCD is an approach that bridges microscopic and continuum descriptions, allowing simulations of composite active-passive systems.

Submersed micropatterned structures control active nematic flow, topology, and concentration.

Coupling between flows and material properties imbues rheological matter with its wide-ranging applicability, hence the excitement for harnessing the rheology of active fluids for which internal structure and continuous energy injection lead to spontaneous flows and complex, out-of-equilibrium dynamics. We propose and demonstrate a convenient, highly tunable method for controlling flow, topology, and composition within active films. Our approach establishes rheological coupling via the indirect presence of fully submersed micropatterned structures within a thin, underlying oil layer. Simulations reveal that micropatterned structures produce effective virtual boundaries within the superjacent active nematic film due to differences in viscous dissipation as a function of depth. This accessible method of applying position-dependent, effective dissipation to the active films presents a nonintrusive pathway for engineering active microfluidic systems.

Driven topological transitions in active nematic films.

The topological properties of many materials are central to their behavior. In intrinsically out-of-equilibrium active materials, the dynamics of topological defects can be particularly important. In this paper, local manipulation of the order, dynamics, and topological properties of microtubule-based active nematic films is demonstrated in a joint experimental and simulation study. Hydrodynamic stresses created by magnetically actuated rotation of disk-shaped colloids in proximity to the films compete with internal stresses in the active nematic, influencing the local motion of +1/2 charge topological defects that are intrinsic to the nematic order in the spontaneously turbulent active films. Sufficiently large applied stresses drive the formation of +1 charge topological vortices through the merger of two +1/2 defects. The directed motion of the defects is accompanied by ordering of the vorticity and velocity of the active flows within the film that is qualitatively unlike the response of passive viscous films. Many features of the film's response to the stress are captured by lattice Boltzmann simulations, providing insight into the anomalous viscoelastic nature of the active nematic. The topological vortex formation is accompanied by a rheological instability in the film that leads to significant increase in the flow velocities. Comparison of the velocity profile in vicinity of the vortex with fluid-dynamics calculations provides an estimate of the film viscosity.

Collective Dynamics of Model Pili-Based Twitcher-Mode Bacilliforms.

Pseudomonas aeruginosa, like many bacilliforms, are not limited only to swimming motility but rather possess many motility strategies. In particular, twitching-mode motility employs hair-like pili to transverse moist surfaces with a jittery irregular crawl. Twitching motility plays a critical role in redistributing cells on surfaces prior to and during colony formation. We combine molecular dynamics and rule-based simulations to study twitching-mode motility of model bacilliforms and show that there is a critical surface coverage fraction at which collective effects arise. Our simulations demonstrate dynamic clustering of twitcher-type bacteria with polydomains of local alignment that exhibit spontaneous correlated motions, similar to rafts in many bacterial communities.

Morphology of depletant-induced erythrocyte aggregates.

Red blood cells suspended in quiescent plasma tend to aggregate into multicellular assemblages, including linearly stacked columnar rouleaux, which can reversibly form more complex clusters or branching networks. While these aggregates play an essential role in establishing hemorheological and pathological properties, the biophysics behind their self-assembly into dynamic mesoscopic structures remains under-explored. We employ coarse-grained molecular simulations to model low-hematocrit erythrocytes subject to short-range implicit depletion forces, and demonstrate not only that depletion interactions are sufficient to account for a sudden dispersion-aggregate transition, but also that the volume fraction of depletant macromolecules controls small aggregate morphology. We observe a sudden transition from a dispersion to a linear column rouleau, followed by a slow emergence of disorderly amorphous clusters of many short rouleaux at larger volume fractions. This work demonstrates how discocyte topology and short-range, non-specific, physical interactions are sufficient to self-assemble erythrocytes into various aggregate structures, with markedly different morphologies and biomedical consequences.

Twist-induced crossover from two-dimensional to three-dimensional turbulence in active nematics.

While studies of active nematics in two dimensions have shed light on various aspects of the flow regimes and topology of active matter, three-dimensional properties of topological defects and chaotic flows remain unexplored. By confining a film of active nematics between two parallel plates, we use continuum simulations and analytical arguments to demonstrate that the crossover from quasi-two-dimensional (quasi-2D) to three-dimensional (3D) chaotic flows is controlled by the morphology of the disclination lines. For small plate separations, the active nematic behaves as a quasi-2D material, with straight topological disclination lines spanning the height of the channel and exhibiting effectively 2D active turbulence. Upon increasing channel height, we find a crossover to 3D chaotic flows due to the contortion of disclinations above a critical activity. Above this critical activity highly contorted disclination lines and disclination loops are formed. We further show that these contortions are engendered by twist perturbations producing a sharp change in the curvature of disclinations.

Enhanced gel formation in binary mixtures of nanocolloids with short-range attraction.

Colloidal suspensions transform between fluid and disordered solid states as parameters such as the colloid volume fraction and the strength and nature of the colloidal interactions are varied. Seemingly subtle changes in the characteristics of the colloids can markedly alter the mechanical rigidity and flow behavior of these soft composite materials. This sensitivity creates both a scientific challenge and an opportunity for designing suspensions for specific applications. In this paper, we report a novel mechanism of gel formation in mixtures of weakly attractive nanocolloids with modest size ratio. Employing a combination of x-ray photon correlation spectroscopy, rheometry, and molecular dynamics simulations, we find that gels are stable at remarkably weaker attraction in mixtures with size ratio near two than in the corresponding monodisperse suspensions. In contrast with depletion-driven gelation at larger size ratio, gel formation in the mixtures is triggered by microphase demixing of the species into dense regions of immobile smaller colloids surrounded by clusters of mobile larger colloids that is not predicted by mean-field thermodynamic considerations. These results point to a new route for tailoring nanostructured colloidal solids through judicious combination of interparticle interaction and size distribution.

Onset of meso-scale turbulence in active nematics.

Meso-scale turbulence is an innate phenomenon, distinct from inertial turbulence, that spontaneously occurs at low Reynolds number in fluidized biological systems. This spatiotemporal disordered flow radically changes nutrient and molecular transport in living fluids and can strongly affect the collective behaviour in prominent biological processes, including biofilm formation, morphogenesis and cancer invasion. Despite its crucial role in such physiological processes, understanding meso-scale turbulence and any relation to classical inertial turbulence remains obscure. Here we show how the motion of active matter along a micro-channel transitions to meso-scale turbulence through the evolution of locally disordered patches (active puffs) from an ordered vortex-lattice flow state. We demonstrate that the stationary critical exponents of this transition to meso-scale turbulence in a channel coincide with the directed percolation universality class. This finding bridges our understanding of the onset of low-Reynolds-number meso-scale turbulence and traditional scale-invariant turbulence in confinement.

Dancing disclinations in confined active nematics.

The spontaneous emergence of collective flows is a generic property of active fluids and often leads to chaotic flow patterns characterised by swirls, jets, and topological disclinations in their orientation field. However, the ability to achieve structured flows and ordered disclinations is of particular importance in the design and control of active systems. By confining an active nematic fluid within a channel, we find a regular motion of disclinations, in conjunction with a well defined and dynamic vortex lattice. As pairs of moving disclinations travel through the channel, they continually exchange partners producing a dynamic ordered state, reminiscent of Ceilidh dancing. We anticipate that this biomimetic ability to self-assemble organised topological disclinations and dynamically structured flow fields in engineered geometries will pave the road towards establishing new active topological microfluidic devices.

Rapid dynamics of cell-shape recovery in response to local deformations.

It is vital that cells respond rapidly to mechanical cues within their microenvironment through changes in cell shape and volume, which rely upon the mechanical properties of cells' highly interconnected cytoskeletal networks and intracellular fluid redistributions. While previous research has largely investigated deformation mechanics, we now focus on the immediate cell-shape recovery response following mechanical perturbation by inducing large, local, and reproducible cellular deformations using AFM. By continuous imaging within the plane of deformation, we characterize the membrane and cortical response of HeLa cells to unloading, and model the recovery via overdamped viscoelastic dynamics. Importantly, the majority (90%) of HeLa cells recover their cell shape in <1 s. Despite actin remodelling on this time scale, we show that cell-shape recovery time is not affected by load duration, nor magnitude for untreated cells. To further explore this rapid recovery response, we expose cells to cytoskeletal destabilizers and osmotic shock conditions, which uncovers the interplay between actin and osmotic pressure. We show that the rapid dynamics of recovery depend crucially on intracellular pressure, and provide strong evidence that cortical actin is the key regulator in the cell-shape recovery processes, in both cancerous and non-cancerous epithelial cells.

Active micromachines: Microfluidics powered by mesoscale turbulence.

Dense active matter, from bacterial suspensions and microtubule bundles driven by motor proteins to cellular monolayers and synthetic Janus particles, is characterized by mesoscale turbulence, which is the emergence of chaotic flow structures. By immersing an ordered array of symmetric rotors in an active fluid, we introduce a microfluidic system that exploits spontaneous symmetry breaking in mesoscale turbulence to generate work. The lattice of rotors self-organizes into a spin state where neighboring discs continuously rotate in permanent alternating directions due to combined hydrodynamic and elastic effects. Our virtual prototype demonstrates a new research direction for the design of micromachines powered by the nematohydrodynamic properties of active turbulence.

Encapsulation-free controlled release: Electrostatic adsorption eliminates the need for protein encapsulation in PLGA nanoparticles.

Encapsulation of therapeutic molecules within polymer particles is a well-established method for achieving controlled release, yet challenges such as low loading, poor encapsulation efficiency, and loss of protein activity limit clinical translation. Despite this, the paradigm for the use of polymer particles in drug delivery has remained essentially unchanged for several decades. By taking advantage of the adsorption of protein therapeutics to poly(lactic-co-glycolic acid) (PLGA) nanoparticles, we demonstrate controlled release without encapsulation. In fact, we obtain identical, burst-free, extended-release profiles for three different protein therapeutics with and without encapsulation in PLGA nanoparticles embedded within a hydrogel. Using both positively and negatively charged proteins, we show that short-range electrostatic interactions between the proteins and the PLGA nanoparticles are the underlying mechanism for controlled release. Moreover, we demonstrate tunable release by modifying nanoparticle concentration, nanoparticle size, or environmental pH. These new insights obviate the need for encapsulation and offer promising, translatable strategies for a more effective delivery of therapeutic biomolecules.

Understanding the onset of oscillatory swimming in microchannels.

Self-propelled colloids (swimmers) in confining geometries follow trajectories determined by hydrodynamic interactions with the bounding surfaces. However, typically these interactions are ignored or truncated to the lowest order. We demonstrate that higher-order hydrodynamic moments cause rod-like swimmers to follow oscillatory trajectories in quiescent fluid between two parallel plates, using a combination of lattice-Boltzmann simulations and far-field calculations. This behavior occurs even far from the confining walls and does not require lubrication results. We show that a swimmer's hydrodynamic quadrupole moment is crucial to the onset of the oscillatory trajectories. This insight allows us to develop a simple model for the dynamics near the channel center based on these higher hydrodynamic moments, and suggests opportunities for trajectory-based experimental characterization of swimmers' hydrodynamic properties.