Generalized Swift-Hohenberg models for dense active suspensions

Self-enhanced mobility enables vortex pattern formation in living matter

Ranging from subcellular organelle biogenesis to embryo development, the formation of self-organized structures is a hallmark of living systems. Whereas the emergence of ordered spatial patterns in biology is often driven by intricate chemical signalling that coordinates cellular behaviour and differentiation1-4, purely physical interactions can drive the formation of regular biological patterns such as crystalline vortex arrays in suspensions of spermatozoa5 and bacteria6. Here we discovered a new route to self-organized pattern formation driven by physical interactions, which creates large-scale regular spatial structures with multiscale ordering. Specifically we found that dense bacterial living matter spontaneously developed a lattice of mesoscale, fast-spinning vortices; these vortices each consisted of around 104-105 motile bacterial cells and were arranged in space at greater than centimetre scale and with apparent hexagonal order, whereas individual cells in the vortices moved in coordinated directions with strong polar and vortical order. Single-cell tracking and numerical simulations suggest that the phenomenon is enabled by self-enhanced mobility in the system-that is, the speed of individual cells increasing with cell-generated collective stresses at a given cell density. Stress-induced mobility enhancement and fluidization is prevalent in dense living matter at various scales of length7-9. Our findings demonstrate that self-enhanced mobility offers a simple physical mechanism for pattern formation in living systems and, more generally, in other active matter systems10 near the boundary of fluid- and solid-like behaviours11-17.

Tracheal Ring Formation

The trachea is a long tube that enables air passage between the larynx and the bronchi. C-shaped cartilage rings on the ventral side stabilise the structure. On its esophagus-facing dorsal side, deformable smooth muscle facilitates the passage of food in the esophagus. While the symmetry break along the dorsal-ventral axis is well understood, the molecular mechanism that results in the periodic Sox9 expression pattern that translates into the cartilage rings has remained elusive. Here, we review the molecular regulatory interactions that have been elucidated, and discuss possible patterning mechanisms. Understanding the principles of self-organisation is important, both to define biomedical interventions and to enable tissue engineering.

Turbulence-induced clustering in compressible active fluids

We study a novel phase of active polar fluids, which is characterized by the continuous creation and destruction of dense clusters due to self-sustained turbulence. This state arises due to the interplay between self-advection of the aligned swimmers and their defect topology. The typical cluster size is determined by the characteristic vortex size. Our results are obtained by investigating a continuum model of compressible polar active fluids, which incorporates typical experimental observations in bacterial suspensions, in particular a non-monotone dependence of speed on density.

Emergence and melting of active vortex crystals

Melting of two-dimensional (2D) equilibrium crystals is a complex phenomenon characterized by the sequential loss of positional and orientational order. In contrast to passive systems, active crystals can self-assemble and melt into an active fluid by virtue of their intrinsic motility and inherent non-equilibrium stresses. Currently, the non-equilibrium physics of active crystallization and melting processes is not well understood. Here, we establish the emergence and investigate the melting of self-organized vortex crystals in 2D active fluids using a generalized Toner-Tu theory. Performing extensive hydrodynamic simulations, we find rich transition scenarios. On small domains, we identify a hysteretic transition as well as a transition featuring temporal coexistence of active vortex lattices and active turbulence, both of which can be controlled by self-propulsion and active stresses. On large domains, an active vortex crystal with solid order forms within the parameter range corresponding to active vortex lattices. The melting of this crystal proceeds through an intermediate hexatic phase. Generally, these results highlight the differences and similarities between crystalline phases in active fluids and their equilibrium counterparts.

Condensate formation and multiscale dynamics in two-dimensional active suspensions

The collective effects of microswimmers in active suspensions result in active turbulence, a spatiotemporally chaotic dynamics at mesoscale, which is characterized by the presence of vortices and jets at scales much larger than the characteristic size of the individual active constituents. To describe this dynamics, Navier-Stokes-based one-fluid models driven by small-scale forces have been proposed. Here, we provide a justification of such models for the case of dense suspensions in two dimensions (2D). We subsequently carry out an in-depth numerical study of the properties of one-fluid models as a function of the active driving in view of possible transition scenarios from active turbulence to large-scale pattern, referred to as condensate, formation induced by the classical inverse energy cascade in Newtonian 2D turbulence. Using a one-fluid model it was recently shown [M. Linkmann et al., Phys. Rev. Lett 122, 214503 (2019)10.1103/PhysRevLett.122.214503] that two-dimensional active suspensions support two nonequilibrium steady states, one with a condensate and one without, which are separated by a subcritical transition. Here, we report further details on this transition such as hysteresis and discuss a low-dimensional model that describes the main features of the transition through nonlocal-in-scale coupling between the small-scale driving and the condensate.

Phase Transition to Large Scale Coherent Structures in Two-Dimensional Active Matter Turbulence

The collective motion of microswimmers in suspensions induce patterns of vortices on scales that are much larger than the characteristic size of a microswimmer, attaining a state called bacterial turbulence. Hydrodynamic turbulence acts on even larger scales and is dominated by inertial transport of energy. Using an established modification of the Navier-Stokes equation that accounts for the small-scale forcing of hydrodynamic flow by microswimmers, we study the properties of a dense suspension of microswimmers in two dimensions, where the conservation of enstrophy can drive an inverse cascade through which energy is accumulated on the largest scales. We find that the dynamical and statistical properties of the flow show a sharp transition to the formation of vortices at the largest length scale. The results show that 2D bacterial and hydrodynamic turbulence are separated by a subcritical phase transition.

Insensitivity of active nematic liquid crystal dynamics to topological constraints

Confining a liquid crystal imposes topological constraints on the orientational order, allowing global control of equilibrium systems by manipulation of anchoring boundary conditions. In this article, we investigate whether a similar strategy allows control of active liquid crystals. We study a hydrodynamic model of an extensile active nematic confined in containers, with different anchoring conditions that impose different net topological charges on the nematic director. We show that the dynamics are controlled by a complex interplay between topological defects in the director and their induced vortical flows. We find three distinct states by varying confinement and the strength of the active stress: A topologically minimal state, a circulating defect state, and a turbulent state. In contrast to equilibrium systems, we find that anchoring conditions are screened by the active flow, preserving system behavior across different topological constraints. This observation identifies a fundamental difference between active and equilibrium materials.

Turbulent Dynamics of Epithelial Cell Cultures

We investigate the large length and long time scales collective flows and structural rearrangements within in vitro human bronchial epithelial cell (HBEC) cultures. Activity-driven collective flows result in ensembles of vortices randomly positioned in space. By analyzing a large population of vortices, we show that their area follows an exponential law with a constant mean value and their rotational frequency is size independent, both being characteristic features of the chaotic dynamics of active nematic suspensions. Indeed, we find that HBECs self-organize in nematic domains of several cell lengths. Nematic defects are found at the interface between domains with a total number that remains constant due to the dynamical balance of nucleation and annihilation events. The mean velocity fields in the vicinity of defects are well described by a hydrodynamic theory of extensile active nematics.

Derivation of a hydrodynamic theory for mesoscale dynamics in microswimmer suspensions

In this paper, we systematically derive a fourth-order continuum theory capable of reproducing mesoscale turbulence in a three-dimensional suspension of microswimmers. We start from overdamped Langevin equations for a generic microscopic model (pushers or pullers), which include hydrodynamic interactions on both small length scales (polar alignment of neighboring swimmers) and large length scales, where the solvent flow interacts with the order parameter field. The flow field is determined via the Stokes equation supplemented by an ansatz for the stress tensor. In addition to hydrodynamic interactions, we allow for nematic pair interactions stemming from excluded-volume effects. The results here substantially extend and generalize earlier findings [S. Heidenreich et al., Phys. Rev. E 94, 020601 (2016)2470-004510.1103/PhysRevE.94.020601], in which we derived a two-dimensional hydrodynamic theory. From the corresponding mean-field Fokker-Planck equation combined with a self-consistent closure scheme, we derive nonlinear field equations for the polar and the nematic order parameter, involving gradient terms of up to fourth order. We find that the effective microswimmer dynamics depends on the coupling between solvent flow and orientational order. For very weak coupling corresponding to a high viscosity of the suspension, the dynamics of mesoscale turbulence can be described by a simplified model containing only an effective microswimmer velocity.