Collective dynamics of self-propelled semiflexible filaments

Density and inertia effects on two-dimensional active semiflexible filament suspensions

We examine the influence of density on the transition between chain and spiral structures in planar assemblies of active semiflexible filaments, utilizing detailed numerical simulations. We focus on how increased density, and higher Péclet numbers, affect the activity-induced transition spiral state in a semiflexible, self-avoiding active chain. Our findings show that increasing the density causes the spiral state to break up, reverting to a motile chain-like shape. This results in a density-dependent reentrant phase transition from spirals back to open chains. We attribute this phenomenon to an inertial effect observed at the single polymer level, where increased persistence length due to inertia has been shown in recent three-dimensional studies to cause polymers to open up. Our two-dimensional simulations further reveal that a reduction in the damping coefficient leads to partial unwinding of the spirals, forming longer arms. In suspension, interactions among these extended arms can trigger a complete unwinding of the spirals, driven by the combined effects of density and inertia.

Constrained motion of self-propelling eccentric disks linked by a spring

It has been supposed that the interplay of elasticity and activity plays a key role in triggering the non-equilibrium behaviors in biological systems. However, the experimental model system is missing to investigate the spatiotemporally dynamical phenomena. Here, a model system of an active chain, where active eccentric-disks are linked by a spring, is designed to study the interplay of activity, elasticity, and friction. Individual active chain exhibits longitudinal and transverse motions; however, it starts to self-rotate when pinning one end and self-beat when clamping one end. In addition, our eccentric-disk model can qualitatively reproduce such behaviors and explain the unusual self-rotation of the first disk around its geometric center. Furthermore, the structure and dynamics of long chains were studied via simulations without steric interactions. It was found that a hairpin conformation emerges in free motion, while in the constrained motions, the rotational and beating frequencies scale with the flexure number (the ratio of self-propelling force to bending rigidity), χ, as ∼(χ)4/3. Scaling analysis suggests that it results from the balance between activity and energy dissipation. Our findings show that topological constraints play a vital role in non-equilibrium synergy behaviors.

Inertia and activity: spiral transitions in semi-flexible, self-avoiding polymers

We consider a two-dimensional, tangentially active, semi-flexible, self-avoiding polymer to find a dynamical re-entrant transition between motile open chains and spinning achiral spirals with increasing activity. Utilizing probability distributions of the turning number, we ascertain the comparative stability of the spiral structure and present a detailed phase diagram within the activity inertia plane. The onset of spiral formation at low activity levels is governed by a torque balance and is independent of inertia. At higher activities, however, inertial effects lead to spiral destabilization, an effect absent in the overdamped limit. We further delineate alterations in size and shape by analyzing the end-to-end distance distribution and the radius of gyration tensor. The Kullback-Leibler divergence from equilibrium distributions exhibits a non-monotonic relationship with activity, reaching a peak at the most compact spirals characterized by the most persistent spinning. As inertia increases, this divergence from equilibrium diminishes.

Collective behavior of squirmers in thin films

Bacteria in biofilms form complex structures and can collectively migrate within mobile aggregates, which is referred to as swarming. This behavior is influenced by a combination of various factors, including morphological characteristics and propulsive forces of swimmers, their volume fraction within a confined environment, and hydrodynamic and steric interactions between them. In our study, we employ the squirmer model for microswimmers and the dissipative particle dynamics method for fluid modeling to investigate the collective motion of swimmers in thin films. The film thickness permits a free orientation of non-spherical squirmers, but constraints them to form a two-layered structure at maximum. Structural and dynamic properties of squirmer suspensions confined within the slit are analyzed for different volume fractions of swimmers, motility types (e.g., pusher, neutral squirmer, puller), and the presence of a rotlet dipolar flow field, which mimics the counter-rotating flow generated by flagellated bacteria. Different states are characterized, including a gas-like phase, swarming, and motility-induced phase separation, as a function of increasing volume fraction. Our study highlights the importance of an anisotropic swimmer shape, hydrodynamic interactions between squirmers, and their interaction with the walls for the emergence of different collective behaviors. Interestingly, the formation of collective structures may not be symmetric with respect to the two walls. Furthermore, the presence of a rotlet dipole significantly mitigates differences in the collective behavior between various swimmer types. These results contribute to a better understanding of the formation of bacterial biofilms and the emergence of collective states in confined active matter.

Quantifying gliding forces of filamentous cyanobacteria by self-buckling

Filamentous cyanobacteria are one of the oldest and today still most abundant lifeforms on earth, with manifold implications in ecology and economics. Their flexible filaments, often several hundred cells long, exhibit gliding motility in contact with solid surfaces. The underlying force generating mechanism is not yet understood. Here, we demonstrate that propulsion forces and friction coefficients are strongly coupled in the gliding motility of filamentous cyanobacteria. We directly measure their bending moduli using micropipette force sensors, and quantify propulsion and friction forces by analyzing their self-buckling behavior, complemented with analytical theory and simulations. The results indicate that slime extrusion unlikely generates the gliding forces, but support adhesion-based hypotheses, similar to the better-studied single-celled myxobacteria. The critical self-buckling lengths align well with the peaks of natural length distributions, indicating the importance of self-buckling for the organization of their collective in natural and artificial settings.

Excluded volume effects on tangentially driven active ring polymers

The conformational and dynamical properties of active ring polymers are studied by numerical simulations. The two-dimensionally confined polymer is modeled as a closed bead-spring chain, driven by tangential forces, put in contact with a heat bath described by the Brownian multiparticle collision dynamics. Both phantom polymers and chains comprising excluded volume interactions are considered for different bending rigidities. The size and shape are found to be dependent on persistence length, driving force, and bead mutual exclusion. The lack of excluded volume interactions is responsible for a shrinkage of active rings when increasing driving force in the flexible limit, while the presence induces a moderate swelling of chains. The internal dynamics of flexible phantom active rings shows activity-enhanced diffusive behavior at large activity values while, in the case of self-avoiding active chains, it is characterized by active ballistic motion not depending on stiffness. The long-time dynamics of active rings is marked by rotational motion whose period scales as the inverse of the applied tangential force, irrespective of persistence length and beads' self-exclusion.

Non-equilibrium structural and dynamic behaviors of active polymers in complex and crowded environments

Active matter systems, which convert internal chemical energy or energy from the environment into directed motion, are ubiquitous in nature and exhibit a range of emerging non-equilibrium behaviors. However, most of the current works on active matter have been devoted to particles, and the study of active polymers has only recently come into the spotlight due to their prevalence within living organisms. The intricate interplay between activity and conformational degrees of freedom gives rise to novel structural and dynamical behaviors of active polymers. Research in active polymers remarkably broadens diverse concepts of polymer physics, such as molecular architecture, dynamics, scaling and so on, which is of significant importance for the development of new polymer materials with unique performance. Furthermore, active polymers are often found in strongly interacting and crowded systems and in complex environments, so that the understanding of this behavior is essential for future developments of novel polymer-based biomaterials. This review thereby focuses on the study of active polymers in complex and crowded environments, and aims to provide insights into the fundamental physics underlying the adaptive and collective behaviors far from equilibrium, as well as the open challenges that the field is currently facing.

Attractive crowding effect on passive and active polymer looping kinetics

Loop formation in complex environments is crucially important to many biological processes in life. In the present work, we adopt three-dimensional Langevin dynamics simulations to investigate passive and active polymer looping kinetics in crowded media featuring polymer-crowder attraction. We find polymers undergo a remarkable coil-globule-coil transition, highlighted by a marked change in the Flory scaling exponent of the gyration radius. Meanwhile, looping time as a function of the crowder's volume fraction demonstrates an apparent non-monotonic alteration. A small number of crowders induce a compact structure, which largely facilitates the looping process. While a large number of crowders heavily impede end-to-end diffusion, looping kinetics is greatly inhibited. For a self-propelled chain, we find that the attractive crowding triggers an unusual activity effect on looping kinetics. Once a globular state is formed, activity takes an effort to open the chain from the compact structure, leading to an unexpected activity-induced inhibition of looping. If the chain maintains a coil state, the dominant role of activity is to enhance diffusivity and, thus, speed up looping kinetics. The novel conformational change and looping kinetics of both passive and active polymers in the presence of attractive crowding highlight a rather distinct scenario that has no analogy in a repulsive crowding counterpart. The underlying mechanism enriches our understanding of the crucial role of attractive interactions in modulating polymer structure and dynamics.

Conformation and dynamics of partially active linear polymers

We perform numerical simulations of isolated, partially active polymers, driven out-of-equilibrium by a fraction of their monomers. We show that, if the active beads are all gathered in a contiguous block, the position of the section along the chain determines the conformational and dynamical properties of the system. Notably, one can modulate the diffusion coefficient of the polymer from active-like to passive-like just by changing the position of the active block. Further, we show that a slight modification of the self-propulsion rule may give rise to an enhancement of diffusion under certain conditions, despite a decrease of the overall polymer activity. Our findings may help in the modelisation of active biophysical systems, such as filamentous bacteria or worms.

Interplay of cell motility and self-secreted extracellular polymeric substance induced depletion effects on spatial patterning in a growing microbial colony

Reproducing bacteria self-organize to develop patterned biofilms in various conditions. Various factors contribute to the shaping of a multicellular bacterial organization. Here we investigate how motility force and self-secreted extracellular polymeric substances (EPS) influence bacterial cell aggregation, leading to phase-separated colonies using a particle-based/individual-based model. Our findings highlight the critical role of the interplay between motility force and depletion effects in regulating phase separation within a growing colony under far-from-equilibrium conditions. We observe that increased motility force hinders depletion-induced cell aggregation and phase segregation, necessitating a higher depletion effect for highly motile bacteria to undergo phase separation within a growing biofilm. We present a phase diagram illustrating the systematic variation of motility force and repulsive mechanical force, shedding light on the combined contributions of these two factors: self-propulsive motion and aggregation due to the depletion effect, resulting in the presence of small to large bacterial aggregates. Furthermore, our study reveals the dynamic nature of clustering, marked by changes in cluster size over time. Additionally, our findings suggest that differential dispersion among the components can lead to the localization of EPS at the periphery of a growing colony. Our study enhances the understanding of the collective dynamics of motile bacterial cells within a growing colony, particularly in the presence of a self-secreted polymer-driven depletion effect.

Active Spaghetti: Collective Organization in Cyanobacteria

Filamentous cyanobacteria can show fascinating examples of nonequilibrium self-organization, which, however, are not well understood from a physical perspective. We investigate the motility and collective organization of colonies of these simple multicellular lifeforms. As their area density increases, linear chains of cells gliding on a substrate show a transition from an isotropic distribution to bundles of filaments arranged in a reticulate pattern. Based on our experimental observations of individual behavior and pairwise interactions, we introduce a nonreciprocal model accounting for the filaments' large aspect ratio, fluctuations in curvature, motility, and nematic interactions. This minimal model of active filaments recapitulates the observations, and rationalizes the appearance of a characteristic length scale in the system, based on the Péclet number of the cyanobacteria filaments.

Worm blobs as entangled living polymers: from topological active matter to flexible soft robot collectives

Recently, the study of long, slender living worms has gained attention due to their unique ability to form highly entangled physical structures, exhibiting emergent behaviors. These organisms can assemble into an active three-dimensional soft entity referred to as the "blob", which exhibits both solid-like and liquid-like properties. This blob can respond to external stimuli such as light, to move or change shape. In this perspective article, we acknowledge the extensive and rich history of polymer physics, while illustrating how these living worms provide a fascinating experimental platform for investigating the physics of active, polymer-like entities. The combination of activity, long aspect ratio, and entanglement in these worms gives rise to a diverse range of emergent behaviors. By understanding the intricate dynamics of the worm blob, we could potentially stimulate further research into the behavior of entangled active polymers, and guide the advancement of synthetic topological active matter and bioinspired tangling soft robot collectives.

Conformation and dynamics of an active filament in crowded media

The structural and dynamical properties of active filamentous objects under macromolecular crowding have a great relevance in biology. By means of Brownian dynamics simulations, we perform a comparative study for the conformational change and diffusion dynamics of an active chain in pure solvents and in crowded media. Our result shows a robust compaction-to-swelling conformational change with the augment of the Péclet number. The presence of crowding facilitates self-trapping of monomers and, thus, reinforces the activity mediated compaction. In addition, the efficient collisions between the self-propelled monomers and crowders induce a coil-to-globulelike transition, indicated by a marked change of the Flory scaling exponent of the gyration radius. Moreover, the diffusion dynamics of the active chain in crowded solutions demonstrates activity-enhanced subdiffusion. The center of mass diffusion manifests rather new scaling relations with respect to both the chain length and Péclet number. The interplay of chain activity and medium crowding provides a new mechanism to understand the non-trivial properties of active filaments in complex environments.

A mechanistic understanding of microcolony morphogenesis: coexistence of mobile and sessile aggregates

Most bacteria in the natural environment self-organize into collective phases such as cell clusters, swarms, patterned colonies, or biofilms. Several intrinsic and extrinsic factors, such as growth, motion, and physicochemical interactions, govern the occurrence of different phases and their coexistence. Hence, predicting the conditions under which a collective phase emerges due to individual-level interactions is crucial. Here we develop a particle-based biophysical model of bacterial cells and self-secreted extracellular polymeric substances (EPS) to decipher the interplay of growth, motility-mediated dispersal, and mechanical interactions during microcolony morphogenesis. We show that the microcolony dynamics and architecture significantly vary depending upon the heterogeneous EPS production. In particular, microcolony shows the coexistence of both motile and sessile aggregates rendering a transition towards biofilm formation. We identified that the interplay of differential dispersion and the mechanical interactions among the components of the colony determines the fate of the colony morphology. Our results provide a significant understanding of the mechano-self-regulation during biofilm morphogenesis and open up possibilities of designing experiments to test the predictions.

Collective vortical motion and vorticity reversals of self-propelled particles on circularly patterned substrates

The collective behavior of self-propelled particles (SPPs) under the combined effects of a circularly patterned substrate and circular confinement is investigated through coarse-grained molecular dynamics simulations of polarized and disjoint ring polymers. The study is performed over a wide range of values of the SPPs packing fraction ϕ[over ¯], motility force F_{D}, and area fraction of the patterned region. At low packing fractions, the SPPs are excluded from the system's center and exhibit a vortical motion that is dominated by the substrate at intermediate values of F_{D}. This exclusion zone is due to the coupling between the driving force and torque induced by the substrate, which induces an outward spiral motion of the SPPs. For high values of F_{D}, the SPPs exclusion from the center is dominated by the confining boundary. At high values of ϕ[over ¯], the substrate pattern leads to reversals in the vorticity, which become quasiperiodic with increasing ϕ[over ¯]. We also found that the substrate pattern is able to separate SPPs based on their motilities.

Durability of Aligned Microtubules Dependent on Persistence Length Determines Phase Transition and Pattern Formation in Collective Motion

Collective motion is a ubiquitous phenomenon in nature. The collective motion of cytoskeleton filaments results mainly from dynamic collisions and alignments; however, the detailed mechanism of pattern formation still needs to be clarified. In particular, the influence of persistence length, which is a measure of filament flexibility, on collective motion is still unclear and lacks experimental verifications although it is likely to directly affect the orientational flexibility of filaments. Here, we investigated the collective motion of microtubules with different persistence lengths using a microtubule-kinesin motility system. We showed that local interactions between microtubules significantly vary depending on their persistence length. We demonstrated that the bundling of microtubules is enhanced by more durable alignment, rather than by greater likelihood of alignment. An agent-based computational model confirmed that the rigidity-dependent durability of microtubule alignment dominates their collective behavior.

Obstacle-induced giant jammed aggregation of active semiflexible filaments

Filaments driven by bound motor proteins and chains of self-propelled colloidal particles are a typical example of active polymers (APs). Due to deformability, APs exhibit very rich dynamic behaviors and collective assembling structures. Here, we are concerned with a basic question: how APs behave near a single obstacle? We find that, in the presence of a big single obstacle, the assembly of APs becomes a two-state system, i.e. APs either gather nearly completely together into a giant jammed aggregate (GJA) on the surface of the obstacle or distribute freely in space. No partial aggregation is observed. Such a complete aggregation/collection is unexpected since it happens on a smooth convex surface instead of, e.g., a concave wedge. We find that the formation of a GJA experiences a process of nucleation and the curves of the transition between the GJA and the non-aggregate state form hysteresis-like loops. Statistical analysis of massive data on the growing time, chirality and angular velocity of both the GJAs and the corresponding nuclei shows the strong random nature of the phenomenon. Our results provide new insights into the behavior of APs in contact with porous media and also a reference for the design and application of polymeric active materials.

Emergent collective dynamics of pusher and puller squirmer rods: swarming, clustering, and turbulence

We study the interplay of steric and hydrodynamic interactions in suspensions of elongated microswimmers by simulating the full hydrodynamics of squirmer rods in the quasi two-dimensional geometry of a Hele-Shaw cell. To create pusher or puller-type squirmer rods, we concentrate the surface slip-velocity field more to the back or to the front of the rod and thereby are able to tune the rod's force-dipole strength. We study a wide range of aspect ratios and area fractions and provide corresponding state diagrams. The flow field of pusher-type squirmer rods destabilizes ordered structures and favors the disordered state at small area fractions and aspect ratios. Only when steric interactions become relevant, we observe a turbulent and dynamic cluster state, while for large aspect ratios a single swarm and jammed cluster occurs. The power spectrum of the turbulent state shows two distinct energy cascades at small and large wave numbers with power-law scaling and non-universal exponents. Pullers show a strong tendency to form swarms instead of the disordered state found for neutral and pusher rods. At large area fractions a dynamic cluster is observed and at larger aspect ratio a single swarm or jammed cluster occurs.

Emergent conformational properties of end-tailored transversely propelling polymers

We study the dynamics and conformations of a single active semiflexible polymer whose monomers experience a propulsion force perpendicular to the local tangent, with the end beads being different from the inner beads ("end-tailored"). Using Langevin simulations, we demonstrate that, apart from sideways motion, the relative propulsion strength between the end beads and the polymer backbone significantly changes the conformational properties of the polymers as a function of bending stiffness, end-tailoring and propulsion force. Expectedly, for slower ends the polymer curves away from the moving direction, while faster ends lead to opposite curving, in both cases slightly reducing the center of mass velocity compared to a straight fiber. Interestingly, for faster end beads there is a rich and dynamic morphology diagram: the polymer ends may get folded together to 2D loops or hairpin-like conformations that rotate due to their asymmetry in shape and periodic flapping motion around a rather straight state during full propulsion is also possible. We rationalize the simulations using scaling and kinematic arguments and present the state diagram of the conformations. Sideways propelled fibers comprise a rather unexplored and versatile class of self-propellers, and their study will open novel ways for designing, e.g. motile actuators or mixers in soft robotics.

Effects of alignment activity on the collapse kinetics of a flexible polymer

The dynamics of various biological filaments can be understood within the framework of active polymer models. Here we consider a bead-spring model for a flexible polymer chain in which the active interaction among the beads is introduced via an alignment rule adapted from the Vicsek model. Following quenching from the high-temperature coil phase to a low-temperature state point, we study the coarsening kinetics via molecular dynamics (MD) simulations using the Langevin thermostat. For the passive polymer case the low-temperature equilibrium state is a compact globule. The results from our MD simulations reveal that though the globular state is also the typical final state in the active case, the nonequilibrium pathways to arrive at such a state differ from the picture for the passive case due to the alignment interaction among the beads. We notice that deviations from the intermediate "pearl-necklace"-like arrangement, which is observed in the passive case, and the formation of more elongated dumbbell-like structures increase with increasing activity. Furthermore, it appears that while a small active force on the beads certainly makes the coarsening process much faster, there exists a nonmonotonic dependence of the collapse time on the strength of active interaction. We quantify these observations by comparing the scaling laws for the collapse time and growth of pearls with the passive case.

Collective motion of cells modeled as ring polymers

In this article, we use a coarse-grained model of disjoint semi-flexible ring polymers to investigate computationally the spatiotemporal collective behavior of cell colonies. A ring polymer in this model is self-propelled by a motility force along the cell's polarity, which depends on its historical kinetics. Despite the repulsive interaction between the cells, a collective behavior sets in as a result of cells pushing against each other. This cooperative motion emerges as the amplitude of the motility force is increased and/or their areal density is increased. The degree of collectivity, characterized by the average cluster size, the velocity field order parameter, and the polarity field nematic order parameter, is found to increase with increasing the amplitude of the motility force and area coverage of the cells. Furthermore, the degree of alignment exhibited by the cell velocity field within a cluster is found to be stronger than that exhibited by the cell polarity. Comparison between the collective behavior of elongated cells and that of circular cells, at the same area coverage and motility force, shows that elongated cells exhibit a stronger collective behavior than circular cells, in agreement with earlier studies of self-propelled anisotropic particles. An investigation of two-cell collisions shows that while two clustered cells move in tandem, their polarities are misaligned. As such the cells push against each other while moving coherently.

Vesicle shape transformations driven by confined active filaments

In active matter systems, deformable boundaries provide a mechanism to organize internal active stresses. To study a minimal model of such a system, we perform particle-based simulations of an elastic vesicle containing a collection of polar active filaments. The interplay between the active stress organization due to interparticle interactions and that due to the deformability of the confinement leads to a variety of filament spatiotemporal organizations that have not been observed in bulk systems or under rigid confinement, including highly-aligned rings and caps. In turn, these filament assemblies drive dramatic and tunable transformations of the vesicle shape and its dynamics. We present simple scaling models that reveal the mechanisms underlying these emergent behaviors and yield design principles for engineering active materials with targeted shape dynamics.

Dependency of active pressure and equation of state on stiffness of wall

Autonomous motion and motility are hallmarks of active matter. Active agents, such as biological cells and synthetic colloidal particles, consume internal energy or extract energy from the environment to generate self-propulsion and locomotion. These systems are persistently out of equilibrium due to continuous energy consumption. It is known that pressure is not always a state function for generic active matter. Torque interaction between active constituents and confinement renders the pressure of the system a boundary-dependent property. The mechanical pressure of anisotropic active particles depends on their microscopic interactions with a solid wall. Using self-propelled dumbbells confined by solid walls as a model system, we perform numerical simulations to explore how variations in the wall stiffness influence the mechanical pressure of dry active matter. In contrast to previous findings, we find that mechanical pressure can be independent of the interaction of anisotropic active particles with walls, even in the presence of intrinsic torque interaction. Particularly, the dependency of pressure on the wall stiffness vanishes when the stiffness is above a critical level. In such a limit, the dynamics of dumbbells near the walls are randomized due to the large torque experienced by the dumbbells, leading to the recovery of pressure as a state variable of density.

Multi-particle collision dynamics with a non-ideal equation of state. II. Collective dynamics of elongated squirmer rods

Simulations of flow fields around microscopic objects typically require methods that both solve the Navier-Stokes equations and also include thermal fluctuations. One such method popular in the field of soft-matter physics is the particle-based simulation method of multi-particle collision dynamics (MPCD). However, in contrast to the typically incompressible real fluid, the fluid of the traditional MPCD methods obeys the ideal-gas equation of state. This can be problematic because most fluid properties strongly depend on the fluid density. In a recent article, we proposed an extended MPCD algorithm and derived its non-ideal equation of state and an expression for the viscosity. In the present work, we demonstrate its accuracy and efficiency for the simulations of the flow fields of single squirmers and of the collective dynamics of squirmer rods. We use two exemplary squirmer-rod systems for which we compare the outcome of the extended MPCD method to the well-established MPCD version with an Andersen thermostat. First, we explicitly demonstrate the reduced compressibility of the MPCD fluid in a cluster of squirmer rods. Second, for shorter rods, we show the interesting result that in simulations with the extended MPCD method, dynamic swarms are more pronounced and have a higher polar order. Finally, we present a thorough study of the state diagram of squirmer rods moving in the center plane of a Hele-Shaw geometry. From a small to large aspect ratio and density, we observe a disordered state, dynamic swarms, a single swarm, and a jammed cluster, which we characterize accordingly.

A single active ring model with velocity self-alignment

Cellular tissue behavior is a multiscale problem. At the cell level, out of equilibrium, biochemical reactions drive physical cell-cell interactions in a typical active matter process. Cell modeling computer simulations are a robust tool to explore countless possibilities and test hypotheses. Here, we introduce a two-dimensional, extended active matter model for biological cells. A ring of interconnected self-propelled particles represents the cell. Neighboring particles are subject to harmonic and bending potentials. Within a characteristic time, each particle's self-velocity tends to align with its scattering velocity after an interaction. Translational modes, rotational modes, and mixtures of these appear as collective states. Using analytical results derived from active Brownian particles, we identify effective characteristic time scales for ballistic and diffusive movements. Finite-size scale investigation shows that the ring diffusion increases linearly with its size when in collective movement. A study on the ring shape reveals that all collective states are present even when bending forces are weak. In that case, when in a translational mode, the collective velocity aligns with the largest ring's direction in a spontaneous polarization emergence.

Pattern formation and polarity sorting of driven actin filaments on lipid membranes

Pattern formation processes in active systems give rise to a plethora of collective structures. Predicting how the emergent structures depend on the microscopic interactions between the moving agents remains a challenge. By introducing a high-density actin gliding assay on a fluid membrane, we demonstrate the emergence of polar structures in a regime of nematic binary interactions dominated by steric repulsion. The transition from a microscopic nematic symmetry to a macroscopic polar structure is linked to microscopic polarity sorting mechanisms, including accumulation in wedge-like topological defects. Our results should be instrumental for a better understanding of pattern formation and polarity sorting processes in active matter.

Collective motion of active matter is ubiquitously observed, ranging from propelled colloids to flocks of bird, and often features the formation of complex structures composed of agents moving coherently. However, it remains extremely challenging to predict emergent patterns from the binary interaction between agents, especially as only a limited number of interaction regimes have been experimentally observed so far. Here, we introduce an actin gliding assay coupled to a supported lipid bilayer, whose fluidity forces the interaction between self-propelled filaments to be dominated by steric repulsion. This results in filaments stopping upon binary collisions and eventually aligning nematically. Such a binary interaction rule results at high densities in the emergence of dynamic collectively moving structures including clusters, vortices, and streams of filaments. Despite the microscopic interaction having a nematic symmetry, the emergent structures are found to be polar, with filaments collectively moving in the same direction. This is due to polar biases introduced by the stopping upon collision, both on the individual filaments scale as well as on the scale of collective structures. In this context, positive half-charged topological defects turn out to be a most efficient trapping and polarity sorting conformation.

Chiral self-sorting of active semiflexible filaments with intrinsic curvature

Many-body interactions in systems of active matter can cause particles to move collectively and self-organize into dynamic structures with long-range order. In cells, the self-assembly of cytoskeletal filaments is critical for cellular motility, structure, intracellular transport, and division. Semiflexible cytoskeletal filaments driven by polymerization or motor-protein interactions on a two-dimensional substrate, such as the cell cortex, can induce filament bending and curvature leading to interesting collective behavior. For example, the bacterial cell-division filament FtsZ is known to have intrinsic curvature that causes it to self-organize into rings and vortices, and recent experiments reconstituting the collective motion of microtubules driven by motor proteins on a surface have observed chiral symmetry breaking of the collective behavior due to motor-induced curvature of the filaments. Previous work on the self-organization of driven filament systems have not studied the effects of curvature and filament structure on collective behavior. In this work, we present Brownian dynamics simulation results of driven semiflexible filaments with intrinsic curvature and investigate how the interplay between filament rigidity and radius of curvature can tune the self-organization behavior in homochiral systems and heterochiral mixtures. We find a curvature-induced reorganization from polar flocks to self-sorted chiral clusters, which is modified by filament flexibility. This transition changes filament transport from ballistic to diffusive at long timescales.

A semiflexible polymer in a gliding assay: reentrant transition, role of turnover and activity

We consider a model of an extensible semiflexible filament moving in two dimensions on a motility assay of motor proteins represented explicitly as active harmonic linkers. Their heads bind stochastically to polymer segments within a capture radius, and extend along the filament in a directed fashion before detaching. Both the extension and detachment rates are load-dependent and generate an active drive on the filament. The filament undergoes a first order phase transition from the open chain to spiral conformation and shows a reentrant behavior in both the active extension and the turnover, defined as the ratio of attachment-detachment rates. Associated with the phase transition, the size and shape of the polymer change non-monotonically, and the relevant autocorrelation functions display a double-exponential decay. The corresponding correlation times show a maximum signifying the dominance of spirals. The orientational dynamics captures the rotation of spirals, and its correlation time decays with activity as a power law.

Heterogeneous bacterial swarms with mixed lengths

Heterogeneous systems of active matter exhibit a range of complex emergent dynamical patterns. In particular, it is difficult to predict the properties of the mixed system based on its constituents. These considerations are particularly significant for understanding realistic bacterial swarms, which typically develop heterogeneities even when grown from a single cell. Here, mixed swarms of cells with different aspect ratios are studied both experimentally and in simulations. In contrast with previous theory, there is no macroscopic phase segregation. However, locally, long cells act as nucleation cites, around which aggregates of short, rapidly moving cells can form, resulting in enhanced swarming speeds. On the other hand, high fractions of long cells form a bottleneck for efficient swarming. Our results suggest a physical advantage for the spontaneous heterogeneity of bacterial swarm populations.

Synchronized oscillations, traveling waves, and jammed clusters induced by steric interactions in active filament arrays

Autonomous active, elastic filaments that interact with each other to achieve cooperation and synchrony underlie many critical functions in biology. The mechanisms underlying this collective response and the essential ingredients for stable synchronization remain a mystery. Inspired by how these biological entities integrate elasticity with molecular motor activity to generate sustained oscillations, a number of synthetic active filament systems have been developed that mimic oscillations of these biological active filaments. Here, we describe the collective dynamics and stable spatiotemporal patterns that emerge in such biomimetic multi-filament arrays, under conditions where steric interactions may impact or dominate the collective dynamics. To focus on the role of steric interactions, we study the system using Brownian dynamics, without considering long-ranged hydrodynamic interactions. The simulations treat each filament as a connected chain of self-propelling colloids. We demonstrate that short-range steric inter-filament interactions and filament roughness are sufficient - even in the absence of inter-filament hydrodynamic interactions - to generate a rich variety of collective spatiotemporal oscillatory, traveling and static patterns. We first analyze the collective dynamics of two- and three-filament clusters and identify parameter ranges in which steric interactions lead to synchronized oscillations and strongly occluded states. Generalizing these results to large one-dimensional arrays, we find rich emergent behaviors, including traveling metachronal waves, and modulated wavetrains that are controlled by the interplay between the array geometry, filament activity, and filament elasticity. Interestingly, the existence of metachronal waves is non-monotonic with respect to the inter-filament spacing. We also find that the degree of filament roughness significantly affects the dynamics - specifically, filament roughness generates a locking-mechanism that transforms traveling wave patterns into statically stuck and jammed configurations. Taken together, simulations suggest that short-ranged steric inter-filament interactions could combine with complementary hydrodynamic interactions to control the development and regulation of oscillatory collective patterns. Furthermore, roughness and steric interactions may be critical to the development of jammed spatially periodic states; a spatiotemporal feature not observed in purely hydrodynamically interacting systems.

Collective motion of driven semiflexible filaments tuned by soft repulsion and stiffness

In active matter systems, self-propelled particles can self-organize to undergo collective motion, leading to persistent dynamical behavior out of equilibrium. In cells, cytoskeletal filaments and motor proteins form complex structures important for cell mechanics, motility, and division. Collective dynamics of cytoskeletal systems can be reconstituted using filament gliding experiments, in which cytoskeletal filaments are propelled by surface-bound motor proteins. These experiments have observed diverse dynamical states, including flocks, polar streams, swirling vortices, and single-filament spirals. Recent experiments with microtubules and kinesin motor proteins found that the collective behavior of gliding filaments can be tuned by altering the concentration of the crowding macromolecule methylcellulose in solution. Increasing the methylcellulose concentration reduced filament crossing, promoted alignment, and led to a transition from active, isotropically oriented filaments to locally aligned polar streams. This emergence of collective motion is typically explained as an increase in alignment interactions by Vicsek-type models of active polar particles. However, it is not yet understood how steric interactions and bending stiffness modify the collective behavior of active semiflexible filaments. Here we use simulations of driven filaments with tunable soft repulsion and rigidity in order to better understand how the interplay between filament flexibility and steric effects can lead to different active dynamic states. We find that increasing filament stiffness decreases the probability of filament alignment, yet increases collective motion and long-range order, in contrast to the assumptions of a Vicsek-type model. We identify swirling flocks, polar streams, buckling bands, and spirals, and describe the physics that govern transitions between these states. In addition to repulsion and driving, tuning filament stiffness can promote collective behavior, and controls the transition between active isotropic filaments, locally aligned flocks, and polar streams.

Squirmer rods as elongated microswimmers: flow fields and confinement

Microswimmers or active elements, such as bacteria and active filaments, have an elongated shape, which determines their individual and collective dynamics. There is still a need to identify what role long-range hydrodynamic interactions play in their fascinating dynamic structure formation. We construct rods of different aspect ratios using several spherical squirmer model swimmers. With the help of the mesoscale simulation method of multi-particle collision dynamics we analyze the flow fields of these squirmer rods both in a bulk fluid and in Hele-Shaw geometries of different slab widths. Based on the hydrodynamic multipole expansion either for bulk or confinement between two parallel plates, we categorize the different multipole contributions of neutral as well as pusher-type squirmer rods. We demonstrate how confinement alters the radial decay of the flow fields for a given force or source multipole moment compared to the bulk fluid.

Rheology of Entangled Active Polymer-Like T. Tubifex Worms

We experimentally study the rheology of long, slender, and entangled living worms (Tubifex Tubifex). Their level of activity can be controlled by changing the temperature or by adding small amounts of alcohol to make the worms temporarily inactive. Performing classical rheology experiments on this entangled polymer-like system, we find that the rheology is qualitatively similar to that of usual polymers, but, quantitatively, (i) shear thinning is reduced by activity, (ii) the characteristic shear rate for the onset of shear-thinning is given by the time scale of the activity, and (iii) the low shear viscosity as a function of concentration shows a very different scaling from that of regular polymers. Our study paves the way towards a new experimental research field of active "polymer-like worms."

How a local active force modifies the structural properties of polymers

We study the dynamics of a polymer, described as a variant of a Rouse chain, driven by an active terminal monomer (head). The local active force induces a transition from a globule-like to an elongated state, as revealed by the study of the end-to-end distance, the variance of which is analytically predicted under suitable approximations. The change in the relaxation times of the Rouse-modes produced by the local self-propulsion is consistent with the transition from globule to elongated conformations. Moreover, also the bond-bond spatial correlation for the chain head are affected by the self-propulsion and a gradient of over-stretched bonds along the chain is observed. We compare our numerical results both with the phenomenological stiff-polymer theory and several analytical predictions in the Rouse-chain approximation.

Conformation and dynamics of a self-avoiding active flexible polymer

We investigate conformations and dynamics of a polymer considering its monomers to be active Brownian particles. This active polymer shows very intriguing physical behavior which is absent in an active Rouse chain. The chain initially shrinks with active force, which starts swelling on further increase in force. The shrinkage followed by swelling is attributed purely to excluded-volume interactions among the monomers. In the swelling regime, the chain shows a crossover from the self-avoiding behavior to the Rouse behavior with scaling exponent ν_{a}≈1/2 for end-to-end distance. The nonmonotonicity in the structure is analyzed through various physical quantities; specifically, radial distribution function of monomers, scattering time, as well as various energy calculations. The chain relaxes faster than the Rouse chain in the intermediate force regime, with a crossover in variation of relaxation time at large active force as given by a power law τ_{r}∼Pe^{-4/3} (Pe is Péclet number).

Dynamics of Active Filaments in Porous Media

The motion of active polymers in a two-dimensional porous medium is shown to depend critically on flexibility, activity, and degree of polymerization. For a given Péclet number, we observe a transition from localization to diffusion as the stiffness of the chains is increased. Whereas stiff chains move almost unhindered through the porous medium, flexible ones spiral and get stuck. Their motion can be accounted for by the model of a continuous time random walk with a renewal process corresponding to unspiraling. The waiting time distribution is shown to develop heavy tails for decreasing stiffness, resulting in subdiffusive and ultimately caged behavior.

Morphological transitions of axially-driven microfilaments

The interactions of microtubules with motor proteins are ubiquitous in cellular and sub-cellular processes that involve motility and cargo transport. In vitro motility assays have demonstrated that motor-driven microtubules exhibit rich dynamical behaviors from straight to curved configurations. Here, we theoretically investigate the dynamic instabilities of elastic filaments, with free-ends, driven by single follower forces that emulate the action of molecular motors. Using the resistive force theory at low Reynolds number, and a combination of numerical techniques with linear stability analysis, we show the existence of four distinct regimes of filament behavior, including a novel buckled state with locked curvature. These successive instabilities recapitulate the full range of experimentally-observed microtubule behavior, implying that neither structural nor actuation asymmetry are needed to elicit this rich repertoire of motion.

Active Brownian filaments with hydrodynamic interactions: conformations and dynamics

The conformational and dynamical properties of active self-propelled filaments/polymers are investigated in the presence of hydrodynamic interactions by both, Brownian dynamics simulations and analytical theory. Numerically, a discrete linear chain composed of active Brownian particles is considered, analytically, a continuous linear semiflexible polymer with active velocities changing diffusively. The force-free nature of active monomers is accounted for-no Stokeslet fluid flow induced by active forces-and higher order hydrodynamic multipole moments are neglected. Hence, fluid-mediated interactions are assumed to arise solely due to intramolecular forces. The hydrodynamic interactions (HI) are taken into account analytically by the preaveraged Oseen tensor, and numerically by the Rotne-Prager-Yamakawa tensor. The nonequilibrium character of the active process implies a dependence of the stationary-state properties on HI via the polymer relaxation times. In particular, at moderate activities, HI lead to a substantial shrinkage of flexible and semiflexible polymers to an extent far beyond shrinkage of comparable free-draining polymers; even flexible HI-polymers shrink, while active free-draining polymers swell monotonically. Large activities imply a reswelling, however, to a less extent than for non-HI polymers, caused by the shorter polymer relaxation times due to hydrodynamic interactions. The polymer mean square displacement is enhanced, and an activity-determined ballistic regime appears. Over a wide range of time scales, flexible active polymers exhibit a hydrodynamically governed subdiffusive regime, with an exponent significantly smaller than that of the Rouse and Zimm models of passive polymers. Compared to simulations, the analytical approach predicts a weaker hydrodynamic effect. Overall, hydrodynamic interactions modify the conformational and dynamical properties of active polymers substantially.

Morphological and dynamical properties of semiflexible filaments driven by molecular motors

We consider an explicit model of a semiflexible filament moving in two dimensions on a gliding assay of motor proteins, which attach to and detach from filament segments stochastically, with a detachment rate that depends on the local load experienced. Attached motor proteins move along the filament to one of its ends with a velocity that varies nonlinearly with the motor protein extension. The resultant force on the filament drives it out of equilibrium. The distance from equilibrium is reflected in the end-to-end distribution, modified bending stiffness, and a transition to spiral morphology of the polymer. The local stress dependence of activity results in correlated fluctuations in the speed and direction of the center of mass leading to a series of ballistic-diffusive crossovers in its dynamics.

Emergent topological phenomena in active polymeric fluids

Polymeric fluids show a wealth of topological phenomena, from entanglement and reptation at microscales to orientational ordering and defect production at macroscales, which can be explained by statistical-mechanical theories. In the presence of activity, the latter must be augmented by forces that cause spontaneous chain motion and fluid flow. Here, using such augmented Langevin equations, we study active polymeric solutions and melts composed of chains of hydrodynamically interacting stresslets. In a spherical volume, contractile chains are unstable and self-knot into entangled melts at both low and high densities. Extensile chains in the same geometry form an unentangled reptating state at low densities and an entangled, coherently moving, non-reptating state at high densities. On a spherical surface, contractile chains show transitions, with increasing areal density, between isotropic, orientationally ordered and micro-phase separated states. Extensile chains in the same geometry show a transition between isotropic and nematic states. In both cases, defects in orientationally ordered states are produced athermally and without conserving topological charge. Our work reproduces the phenomenology of several recent experiments, highlights the importance of hydrodynamic interactions in active polymer fluids, and suggests non-equilibrium kinetic routes to topological structures that are otherwise difficult to obtain in equilibrium.

Chronology of motor-mediated microtubule streaming

We introduce a filament-based simulation model for coarse-grained, effective motor-mediated interaction between microtubule pairs to study the time-scales that compose cytoplasmic streaming. We characterise microtubule dynamics in two-dimensional systems by chronologically arranging five distinct processes of varying duration that make up streaming, from microtubule pairs to collective dynamics. The structures found were polarity sorted due to the propulsion of antialigned microtubules. This also gave rise to the formation of large polar-aligned domains, and streaming at the domain boundaries. Correlation functions, mean squared displacements, and velocity distributions reveal a cascade of processes ultimately leading to microtubule streaming and advection, spanning multiple microtubule lengths. The characteristic times for the processes extend over three orders of magnitude from fast single-microtubule processes to slow collective processes. Our approach can be used to directly test the importance of molecular components, such as motors and crosslinking proteins between microtubules, on the collective dynamics at cellular scale.

Active Brownian Filamentous Polymers under Shear Flow

The conformational and rheological properties of active filaments/polymers exposed to shear flow are studied analytically. Using the continuous Gaussian semiflexible polymer model extended by the activity, we derive analytical expressions for the dependence of the deformation, orientation, relaxation times, and viscosity on the persistence length, shear rate, and activity. The model yields a Weissenberg-number dependent shear-induced deformation, alignment, and shear thinning behavior, similarly to the passive counterpart. Thereby, the model shows an intimate coupling between activity and shear flow. As a consequence, activity enhances the shear-induced polymer deformation for flexible polymers. For semiflexible polymers/filaments, a nonmonotonic deformation is obtained because of the activity-induced shrinkage at moderate and swelling at large activities. Independent of stiffness, activity-induced swelling facilitates and enhances alignment and shear thinning compared to a passive polymer. In the asymptotic limit of large activities, a polymer length- and stiffness-independent behavior is obtained, with universal shear-rate dependencies for the conformations, dynamics, and rheology.