Dynamically generated patterns in dense suspensions of active filaments

Collective self-caging of active filaments in virtual confinement

Motility coupled to responsive behavior is essential for many microorganisms to seek and establish appropriate habitats. One of the simplest possible responses, reversing the direction of motion, is believed to enable filamentous cyanobacteria to form stable aggregates or accumulate in suitable light conditions. Here, we demonstrate that filamentous morphology in combination with responding to light gradients by reversals has consequences far beyond simple accumulation: Entangled aggregates form at the boundaries of illuminated regions, harnessing the boundary to establish local order. We explore how the light pattern, in particular its boundary curvature, impacts aggregation. A minimal mechanistic model of active flexible filaments resembles the experimental findings, thereby revealing the emergent and generic character of these structures. This phenomenon may enable elongated microorganisms to generate adaptive colony architectures in limited habitats or guide the assembly of biomimetic fibrous materials.

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.

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.

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.

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.

Emergence of diverse patterns driven by molecular motors in the motility assay

Actomyosin contractility originating from interactions between F-actin and myosin motors in the actin cytoskeleton generates mechanical forces and drives a wide range of cellular processes including cell migration and cytokinesis. To probe the interactions between F-actin and myosin motors, the myosin motility assay has been popularly employed, which consists of myosin heads attached to a glass surface and F-actins gliding on the surface via interactions with the heads. Several experiments have shown that F-actins move in a collective fashion due to volume-exclusion effects between neighboring F-actins. Furthermore, Computational models have shown how changes in key parameters lead to diverse pattern formation in motility assay. However, in most of the computational models, myosin motors were implicitly considered by applying a constant propulsion force to filaments to reduce computational cost. This simplification limits the physiological relevance of the insights provided by the models and potentially leads to artifacts. In this study, we employed an agent-based computational model for the motility assay with explicit immobile motors interacting with filaments. We rigorously account for the kinetics of myosin motors including the force-velocity relationship for walking and the binding and unbinding behaviors. We probed the effects of the length, rigidity, and concentration of filaments and repulsive strength on collective movements and pattern formation. It was found that four distinct types of structures-homogeneous networks, flocks, bands, and rings-emerged as a result of collisions between gliding filaments. We further analyzed the frequency and morphology of these structures and the curvature, alignment, and rotational motions of filaments. Our study provides better insights into the origin and properties of patterns formed by gliding filaments beyond what was shown before.

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.

Effects of inertia on conformation and dynamics of tangentially driven active filaments

Active filamentlike systems propelling along their backbone exist across scales ranging from motor-driven biofilaments to worms and robotic chains. In macroscopic active filaments such as a chain of robots, in contrast to their microscopic counterparts, inertial effects on their motion cannot be ignored. Nonetheless, the consequences of the interplay between inertia and flexibility on the shape and dynamics of active filaments remain unexplored. Here we examine inertial effects on a flexible tangentially driven active polymer model pertinent to the above examples and we determine the conditions under which inertia becomes important. Performing Langevin dynamics simulations of active polymers with underdamped and overdamped dynamics for a wide range of contour lengths and activities, we uncover striking inertial effects on conformation and dynamics for high levels of activities. Inertial collisions increase the persistence length of active polymers and remarkably alter their scaling behavior. In stark contrast to passive polymers, inertia leaves its fingerprint at long times by an enhanced diffusion of the center of mass. We rationalize inertia-induced enhanced dynamics by analytical calculations of center-of-mass velocity correlations, applicable to any active polymer model, which reveal significant contributions from active force fluctuations convoluted by inertial relaxation.

A polymer chain with dipolar active forces in connection to spatial organization of chromatin

A living cell is an active environment where the organization and dynamics of chromatin are affected by different forms of activity. Optical experiments report that loci show subdiffusive dynamics and the chromatin fiber is seen to be coherent over micrometer-scale regions. Using a bead-spring polymer chain with dipolar active forces, we study how the subdiffusive motion of the loci generate large-scale coherent motion of the chromatin. We show that in the presence of extensile (contractile) activity, the dynamics of the loci grows faster (slower) and the spatial correlation length increases (decreases) compared to the case with no dipolar forces. Hence, both the dipolar active forces modify the elasticity of the chain. Interestingly in our model, the dynamics and organization of such dipolar active chains largely differ from the passive chain with renormalized elasticity.

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.

Facilitated dynamics of an active polymer in 2D crowded environments with obstacles

Understanding the behaviors of a single active chain in complex environments is not only an interesting topic in non-equilibrium physics but also has applicative implications in biological/medical engineering. In this work, by using molecular simulations, we systematically study the dynamical and conformational behaviors of an active polymer in crowded environments, i.e., a single active chain confined in 2D space with randomly arranged obstacles. We found that the competition between the chain's activity and rigidity in the presence of obstacles leads to many interesting dynamical and conformational states, such as the diffusive expanded state, the diffusive collapsed state, and the localized collapsed state. Importantly, we found a counter-intuitive phenomenon, i.e., crowded environments facilitate the diffusion of the active polymer within a large parameter space. As the crowdedness (packing fraction of obstacles) increases, the parameter space in which crowding-enhanced diffusion occurs still remains. This abnormal dynamics is attributed to a structural reason that the obstacles prevent active chains from collapsing. Our findings capture some generic features of active polymers in complex environments and provide insights into the design of novel drug delivery systems.

Mild stratification in drying films of colloidal mixtures

Size stratification of bidisperse colloidal mixtures during vertical drying was investigated using the implicit solvent Langevin dynamics (LD) simulation and the explicit solvent lattice Boltzmann (LB) method. Simulations were performed for the Péclet number (Pe) over a wide range of 1-1000. In the case of a low size ratio of 2, mild stratification was observed in both simulation methods, in contrast to distinct stratification with thick "small-on-top" or "large-on-top" layers. The LD simulations exhibited a "small-on-top" stratification or mixed state. In contrast, the LB simulations exhibited a "large-on-top" or mixed state, according to the variation in Pe. The results demonstrated that the explicit solvent reduced the collective diffusion under moderate Pe conditions. This suppressed the steep concentration gradient of small particles in the packed region of particles near the air-solvent interface. Thus, distinguishable stratification patterns were obtained for the implicit and explicit solvent models.

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.

Structure and dynamics of an active polymer adsorbed on the surface of a cylinder

The structure and dynamics of an active polymer on a smooth cylindrical surface are studied by Brownian dynamics simulations. The effect of an active force on the polymer adsorption behavior and the combined effect of chain mobility, length N, rigidity κ, and cylinder radius, R, on the phase diagrams are systemically investigated. We find that complete adsorption is replaced by the irregular alternative adsorption/desorption process at a large driving force. Three typical (spiral, helix-like, and rod-like) conformations of the active polymer are observed, dependent on N, κ, and R. Dynamically, the polymer shows rotational motion in the spiral state, snake-like motion in the intermediate state, and straight translational motion without turning back in the rod-like state. In the spiral state, we find that the rotation velocity ω and the chain length follow a power-law relation ω ∼ N^-0.42, consistent with the torque-balance theory of general Archimedean spirals. And the polymer shows super-diffusive behavior along the cylinder for a long time in the helix-like and rod-like states. Our results highlight that the mobility, rigidity, and curvature of surface can be used to regulate the polymer behavior.

Absolute negative mobility of active polymer chains in steady laminar flows

We investigate the transport of active polymer chains in steady laminar flows in the presence of thermal noise and an external constant force. In the model, the polymer chain is worm-like and is propelled by active forces along its tangent vectors. Compared with inertial Brownian particles, active polymer chains in steady laminar flows exhibit richer movement patterns due to their specific spatial structures. The simulation results show that the velocity-force relation is strongly dependent on the system parameters such as the chain length, bending rigidity, active force and so on. The polymer chain may move in some preferential movement directions and exhibits absolute negative mobility within appropriate parameter regimes, i.e., the polymer chain can move in a direction opposite to the external constant force. In particular, we can observe giant negative mobility in a broad range of parameter regimes.

Active nematic order and dynamic lane formation of microtubules driven by membrane-bound diffusing motors

Active nematics are ordered liquid crystalline fluids that exhibit spontaneous persistent flows and collective dynamics. The development of motile biopolymer systems inspired by nature has recently attracted considerable attention to out-of-equilibrium soft materials. We report the formation of an active nematic in which microtubules are propelled by kinesin motors coupled to a lipid membrane substrate. The system exhibits apolar order in a globally aligned nematic phase and locally ordered dynamic lanes. Use of a fluid substrate represents a significant advance for active matter as it allows for spatial re-organization of motors, which generate force, in response to the dynamics of the aligning microtubules. This self-organized feedback mechanism may have implications in vivo and for engineering efficient dynamic and reconfigurable materials.

Dynamic lane formation and long-range active nematic alignment are reported using a geometry in which kinesin motors are directly coupled to a lipid bilayer, allowing for in-plane motor diffusion during microtubule gliding. We use fluorescence microscopy to image protein distributions in and below the dense two-dimensional microtubule layer, revealing evidence of diffusion-enabled kinesin restructuring within the fluid membrane substrate as microtubules collectively glide above. We find that the lipid membrane acts to promote filament–filament alignment within the gliding layer, enhancing the formation of a globally aligned active nematic state. We also report the emergence of an intermediate, locally ordered state in which apolar dynamic lanes of nematically aligned microtubules migrate across the substrate. To understand this emergent behavior, we implement a continuum model obtained from coarse graining a collection of self-propelled rods, with propulsion set by the local motor kinetics. Tuning the microtubule and kinesin concentrations as well as active propulsion in these simulations reveals that increasing motor activity promotes dynamic nematic lane formation. Simulations and experiments show that, following fluid bilayer substrate mediated spatial motor restructuring, the total motor concentration becomes enriched below the microtubule lanes that they drive, with the feedback leading to more dynamic lanes. Our results have implications for membrane-coupled active nematics in vivo as well as for engineering dynamic and reconfigurable materials where the structural elements and power sources can dynamically colocalize, enabling efficient mechanical work.

Cluster and conquer: the morphodynamics of invasion of a compliant substrate by active rods

The colonisation of a soft passive material by motile cells such as bacteria is common in biology. The resulting colonies of the invading cells are often observed to exhibit intricate patterns whose morphology and dynamics can depend on a number of factors, particularly the mechanical properties of the substrate and the motility of the individual cells. We use simulations of a minimal 2D model of self-propelled rods moving through a passive compliant medium consisting of particles that offer elastic resistance before being plastically displaced from their equilibrium positions. It is observed that the clustering of active (self-propelled) particles is crucial for understanding the morphodynamics of colonisation. Clustering enables motile colonies to spread faster than they would have as isolated particles. The colonisation rate depends non-monotonically on substrate stiffness with a distinct maximum at a non-zero value of substrate stiffness. This is observed to be due to a change in the morphology of clusters. Furrow networks created by the active particles have a fractal-like structure whose dimension varies systematically with substrate stiffness but is less sensitive to particle activity. The power-law growth exponent of the furrowed area is smaller than unity, suggesting that, to sustain such extensive furrow networks, colonies must regulate their overall growth rate.

A Review on the Some Issues of Multiphase Flow with Self-Driven Particles

Multiphase flow with self-driven particles is ubiquitous and complex. Exploring the flow properties has both important academic meaning and engineering value. This review emphasizes some recent studies on multiphase flow with self-driven particles: the hydrodynamic interactions between self-propelled/self-rotary particles and passive particles; the aggregation, phase separation and sedimentation of squirmers; the influence of rheological properties on its motion; and the kinematic characteristics of axisymmetric squirmers. Finally, some open problems, challenges, and future directions are highlighted.

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.

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.

Periodic oscillations in a string of camphor infused disks

The rhythmic beating motion of autonomously motile filaments has many practical applications. Here, we present an experimental study on a filament made of camphor infused paper disks, stitched together adjacent to each other using nylon thread. The filament displays spontaneous translatory motion when it is placed on the surface of water due to the surface tension gradients created by camphor molecules on the water surface. When this filament is clamped on one end, we obtain regular oscillatory motion instead of translation. The filament shows qualitatively different dynamics at different activity levels, which is controlled by the amount of camphor infused into the paper disks. For a better physical understanding of the filament dynamics, we develop a minimal numerical model involving a semi-flexible filament made of active polar disks, where the polarity is coupled to the instantaneous velocity of the particle. This model qualitatively reproduces different oscillatory modes of the filament. Moreover, our model reveals a rich dynamical state diagram of the system, as a function of filament activity and the coupling strength.

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.

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).

Matrix feedback enables diverse higher-order patterning of the extracellular matrix

The higher-order patterning of extra-cellular matrix in normal and pathological tissues has profound consequences on tissue function. Whilst studies have documented both how fibroblasts create and maintain individual matrix fibers and how cell migration is altered by the fibers they interact with, a model unifying these two aspects of tissue organization is lacking. Here we use computational modelling to understand the effect of this interconnectivity between fibroblasts and matrix at the mesoscale level. We created a unique adaptation to the Vicsek flocking model to include feedback from a second layer representing the matrix, and use experimentation to parameterize our model and validate model-driven hypotheses. Our two-layer model demonstrates that feedback between fibroblasts and matrix increases matrix diversity creating higher-order patterns. The model can quantitatively recapitulate matrix patterns of tissues in vivo. Cells follow matrix fibers irrespective of when the matrix fibers were deposited, resulting in feedback with the matrix acting as temporal 'memory' to collective behaviour, which creates diversity in topology. We also establish conditions under which matrix can be remodelled from one pattern to another. Our model elucidates how simple rules defining fibroblast-matrix interactions are sufficient to generate complex tissue patterns.

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.

The interplay between activity and filament flexibility determines the emergent properties of active nematics

Active nematics are microscopically driven liquid crystals that exhibit dynamical steady states characterized by the creation and annihilation of topological defects. Motivated by differences between previous simulations of active nematics based on rigid rods and experimental realizations based on semiflexible biopolymer filaments, we describe a large-scale simulation study of a particle-based computational model that explicitly incorporates filament semiflexibility. We find that energy injected into the system at the particle scale preferentially excites bend deformations, reducing the apparent filament bend modulus. The emergent characteristics of the active nematic depend on activity and flexibility only through this activity-renormalized bend 'modulus', demonstrating that apparent values of material parameters, such as the Frank 'constants', depend on activity. Thus, phenomenological parameters within continuum hydrodynamic descriptions of active nematics must account for this dependence. Further, we present a systematic way to estimate these parameters from observations of deformation fields and defect shapes in experimental or simulation data.

Inertial delay of self-propelled particles

The motion of self-propelled massive particles through a gaseous medium is dominated by inertial effects. Examples include vibrated granulates, activated complex plasmas and flying insects. However, inertia is usually neglected in standard models. Here, we experimentally demonstrate the significance of inertia on macroscopic self-propelled particles. We observe a distinct inertial delay between orientation and velocity of particles, originating from the finite relaxation times in the system. This effect is fully explained by an underdamped generalisation of the Langevin model of active Brownian motion. In stark contrast to passive systems, the inertial delay profoundly influences the long-time dynamics and enables new fundamental strategies for controlling self-propulsion in active matter.

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.

Collective dynamics of self-propelled semiflexible filaments

The collective behavior of active semiflexible filaments is studied with a model of tangentially driven self-propelled worm-like chains. The combination of excluded-volume interactions and self-propulsion leads to several distinct dynamic phases as a function of bending rigidity, activity, and aspect ratio of individual filaments. We consider first the case of intermediate filament density. For high-aspect-ratio filaments, we identify a transition with increasing propulsion from a state of free-swimming filaments to a state of spiraled filaments with nearly frozen translational motion. For lower aspect ratios, this gas-of-spirals phase is suppressed with growing density due to filament collisions; instead, filaments form clusters similar to self-propelled rods. As activity increases, finite bending rigidity strongly effects the dynamics and phase behavior. Flexible filaments form small and transient clusters, while stiffer filaments organize into giant clusters, similarly to self-propelled rods, but with a reentrant phase behavior from giant to smaller clusters as activity becomes large enough to bend the filaments. For high filament densities, we identify a nearly frozen jamming state at low activities, a nematic laning state at intermediate activities, and an active-turbulence state at high activities. The latter state is characterized by a power-law decay of the energy spectrum as a function of wave number. The resulting phase diagrams encapsulate tunable non-equilibrium steady states that can be used in the organization of living matter.

Dynamical patterns in nematic active matter on a sphere

Using simulations of self-propelled agents with short-range repulsion and nematic alignment, we explore the dynamical phases of a dense active nematic confined to the surface of a sphere. We map the nonequilibrium phase diagram as a function of curvature, alignment strength, and activity. Our model reproduces several phases seen in recent experiments on active microtubule bundles confined the surfaces of vesicles. At low driving, we recover the equilibrium nematic ground state with four +1/2 defects. As the driving is increased, geodesic forces drive the transition to a polar band wrapping around an equator, with large empty spherical caps corresponding to two +1 defects at the poles. Upon further increasing activity, the bands fold onto themselves, and the system eventually transitions to a turbulent state marked by the proliferation of pairs of topological defects. We highlight the key role of the nematic persistence length in controlling pattern formation in these confined systems with positive Gaussian curvature.