Defect-mediated phase transitions in active soft matter

Extreme Spontaneous Deformations of Active Crystals

We demonstrate that two-dimensional crystals made of active particles can experience extremely large spontaneous deformations without melting. Using particles mostly interacting via pairwise repulsive forces, we show that such active crystals maintain long-range bond order and algebraically decaying positional order, but with an exponent η not limited by the 1/3 bound given by the (equilibrium) KTHNY theory. We rationalize our findings using linear elastic theory and show the existence of two well-defined effective temperatures quantifying respectively large-scale deformations and bond-order fluctuations. The root of these phenomena lies in the sole time-persistence of the intrinsic axes of particles, and they should thus be observed in many different situations.

Repulsive torques alone trigger crystallization of constant speed active particles

We investigate the possibility for self-propelled particles to crystallize without reducing their intrinsic speed. We illuminate how, in the absence of any force, the competition between self-propulsion and repulsive torques determines the macroscopic phases of constant-speed active particles. This minimal model expands upon existing approaches for an improved understanding of crystallization of active matter.

Unified analysis of topological defects in 2D systems of active and passive disks

We provide a comprehensive quantitative analysis of localized and extended topological defects in the steady state of 2D passive and active repulsive Brownian disk systems. We show that, both in and out-of-equilibrium, the passage from the solid to the hexatic is driven by the unbinding of dislocations, in quantitative agreement with the KTHNY singularity. Instead, extended clusters of defects largely dominate below the solid-hexatic critical line. The latter percolate in the liquid phase very close to the hexatic-liquid transition, both for continuous and discontinuous transitions, in the homogeneous liquid regime. At critical percolation the clusters of defects are fractal with statistical and geometric properties that are independent of the activity and compatible with the universality class of uncorrelated critical percolation. We also characterize the spatial organization of point-like defects and we show that the disclinations are not free, but rather always very near more complex defect structures. At high activity, the bulk of the dense phase generated by Motility-Induced Phase Separation is characterized by a density of point-like defects, and statistics and morphology of defect clusters, set by the amount of activity and not the packing fraction. Hexatic domains within the dense phase are separated by grain-boundaries along which a finite network of topological defects resides, interrupted by gas bubbles in cavitation. This structure is dynamic in the sense that the defect network allows for an unzipping mechanism that leaves free space for gas bubbles to appear, close, and even be released into the dilute phase.

Alignment destabilizes crystal order in active systems

We combine numerical and analytical methods to study two-dimensional active crystals formed by permanently linked swimmers and with two distinct alignment interactions. The system admits a stationary phase with quasi-long-range translational order, as well as a moving phase with quasi-long-range active force director and velocity order. The translational order in the moving phase is significantly influenced by alignment interaction. For Vicsek-like alignment, the translational order is short ranged, whereas the bond-orientational order is quasi-long ranged, implying a moving hexatic phase. For elasticity-based alignment, the translational order is quasi-long ranged parallel to the motion and short ranged in the perpendicular direction, whereas the bond orientational order is long ranged. We also generalize these results to higher dimensions.

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.

Defect Superdiffusion and Unbinding in a 2D XY Model of Self-Driven Rotors

We consider a nonequilibrium extension of the 2D XY model, equivalent to the noisy Kuramoto model of synchronization with short-range coupling, where rotors sitting on a square lattice are self-driven by random intrinsic frequencies. We study the static and dynamic properties of topological defects (vortices) and establish how self-spinning affects the Berezenskii-Kosterlitz-Thouless phase transition scenario. The nonequilibrium drive breaks the quasi-long-range ordered phase of the 2D XY model into a mosaic of ordered domains of controllable size and results in self-propelled vortices that generically unbind at any temperature, featuring superdiffusion ⟨r^{2}(t)⟩∼t^{3/2} with a Gaussian distribution of displacements. Our work provides a simple framework to investigate topological defects in nonequilibrium matter and sheds new light on the problem of synchronization of locally coupled oscillators.

Spatial propagation of electrical signals in circular biofilms: A combined experimental and agent-based fire-diffuse-fire study

Bacterial biofilms are a risk to human health, playing critical roles in persistent infections. Recent studies have observed electrical signaling in biofilms and thus biofilms represent a new class of active excitable matter in which cell division is the active process and the spiking of the individual bacterial cells is the excitable process. Electrophysiological models have predominantly been developed to describe eukaryotic systems, but we demonstrate their use in understanding bacterial biofilms. Our agent-based fire-diffuse-fire (ABFDF) model successfully simulates the propagation of both centrifugal (away from the center) and centripetal (toward the center) electrical signals through biofilms of Bacillus subtilis. Furthermore, the ABFDF model allows realistic spatial positioning of the bacteria in two dimensions to be included in the fire-diffuse-fire model and this is the crucial factor that improves agreement with experiments. The speed of propagation is not constant and depends on the radius of the propagating electrical wave front. Centripetal waves are observed to move faster than centrifugal waves, which is a curvature driven effect and is correctly captured by our simulations.

Tuning the self-organization of confined active particles by the steepness of the trap

We consider the collective dynamics of self-propelling particles in two dimensions. They can align themselves according to the direction of propulsion of their neighbours, together with small rotational fluctuations. They also interact with each other via soft, isotropic, repulsive potentials. The particles are confined in a circular trap. The steepness of the trap is tuneable. The average packing fraction of the particles is low. When the trap is steep, particles flock along its boundary. They form a polar cluster that spreads over the boundary. The cluster is not spatially ordered. We show that when the steepness is decreased beyond a threshold value, the cluster becomes round and compact and eventually spatial order (hexagonal) emerges in addition to the pre-established polar order. We investigate the kinetics of such ordering. We find that while rotating around the centre of the trap along its circular boundary, the cluster needs to roll around its centre of mass to be spatially ordered. We have studied the stability of the order when the trap is suddenly switched off. We find that for the particles with velocity alignment interaction, the decay of the spatial order is much slower than the particles without the alignment interaction.

Re-entrant phase separation in nematically aligning active polar particles

We present a numerical study of the phase behavior of repulsively interacting active polar particles that align their active velocities nematically. The amplitude of the active velocity, and the noise in its orientational alignment control the active nature of the system. At high values of orientational noise, the structural fluid undergoes a continuous nematic-isotropic transition in active orientation. This transition is well separated from an active phase separation, characterized by the formation of high density hexatic clusters, observed at lower noise strengths. With increasing activity, the system undergoes a re-entrant fluid - phase separation - fluid transition. The phase coexistence at low activity can be understood in terms of motility induced phase separation. In contrast, the re-melting of hexatic clusters, and the collective motion at low orientational noise are dominated by flocking behavior. At high activity, sliding and jamming of polar sub-clusters, formation of grain boundaries, lane formation, and subsequent fragmentation of the polar patches mediate remelting.

Two types of exclusion interactions for self-propelled objects and collective motion induced by their combination

Exclusive interactions between self-driven objects may play crucial roles in their collective behavior, e.g., in collective migration of living cells. Here, such collective behavior is studied based on a simple but sufficient model taking account the exclusion effects, which incorporate the following two distinct kinds of exclusion interactions in two dimensions: The first is the mechanical exclusion wherein two objects mechanically repel each other when they overlap. The second is the scattering exclusion, wherein the directions along which each object tries to move are modulated to avoid overlapping. We propose a theoretical model based on two principles: (1) Each object maintains its own polarity with a fixed strength and attempts to move into the polarity direction and (2) objects interact with each other through the abovementioned exclusions. Based on this model, we look at the difference of consequences and combinatory effects of these two kinds of exclusions. Furthermore, we calculate the polar order of polarity directions without an external directional bias. Our results suggest that the combination of these two kinds of exclusions leads to effectively inelastic scattering of two objects, which eventually gives rise to global polar ordering. We also find that the traveling band can arise by this mechanism of alignment at the intermediate density, as generally seen in collective motion with polar alignment and investigated in various earlier works. Characteristics of transitions among disordered, traveling band, and homogeneously ordered states of the presented model are investigated, and their similarities and differences with those given by the explicit alignment interaction are discussed.

Thermodynamic phases in two-dimensional active matter

Active matter has been much studied for its intriguing properties such as collective motion, motility-induced phase separation and giant fluctuations. However, it has remained unclear how the states of active materials connect with the equilibrium phases. For two-dimensional systems, this is also because the understanding of the liquid, hexatic, and solid equilibrium phases and their phase transitions is recent. Here we show that two-dimensional self-propelled point particles with inverse-power-law repulsions moving with a kinetic Monte Carlo algorithm without alignment interactions preserve all equilibrium phases up to very large activities. Furthermore, at high activity within the liquid phase, a critical point opens up a gas–liquid motility-induced phase separation region. In our model, two-step melting and motility-induced phase separation are thus independent phenomena. We discuss the reasons for these findings to be common to a wide class of two-dimensional active systems.

Adapting statistical physics tools to study active systems is challenging due to their non-equilibrium nature. Here the authors use simulations to present a phase diagram of a 2D active system, showing a two-step melting scenario far from equilibrium along with gas-liquid motility-induced phase separation.

Frustration-induced phases in migrating cell clusters

Certain malignant cancer cells form clusters in a chemoattractant gradient, which can spontaneously show three different phases of motion: translational, rotational, and random. Guided by our experiments on the motion of two-dimensional clusters in vitro, we developed an agent-based model in which the cells form a cohesive cluster due to attractive and alignment interactions. We find that when cells at the cluster rim are more motile, all three phases of motion coexist, in agreement with our observations. Using the model, we show that the transitions between different phases are driven by competition between an ordered rim and a disordered core accompanied by the creation and annihilation of topological defects in the velocity field. The model makes specific predictions, which we verify with our experimental data. Our results suggest that heterogeneous behavior of individuals, based on local environment, can lead to novel, experimentally observed phases of collective motion.

Active Atoms and Interstitials in Two-Dimensional Colloidal Crystals

We study experimentally and numerically the motion of a self-phoretic active particle in two-dimensional loosely packed colloidal crystals at fluid interfaces. Two scenarios emerge depending on the interactions between the active particle and the lattice: the active particle either navigates throughout the crystal as an interstitial or is part of the lattice and behaves as an active atom. Active interstitials undergo a run-and-tumble-like motion, with the passive colloids of the crystal acting as tumbling sites. Instead, active atoms exhibit an intermittent motion, stemming from the interplay between the periodic potential landscape of the passive crystal and the particle's self-propulsion. Our results constitute the first step towards the realization of non-close-packed crystalline phases with internal activity.

Spontaneously Flowing Crystal of Self-Propelled Particles

We experimentally and numerically study the structure and dynamics of a monodisperse packing of spontaneously aligning self-propelled hard disks. The packings are such that their equilibrium counterparts form perfectly ordered hexagonal structures. Experimentally, we first form a perfect crystal in a hexagonal arena which respects the same crystalline symmetry. Frustration of the hexagonal order, obtained by removing a few particles, leads to the formation of a rapidly diffusing "droplet." Removing more particles, the whole system spontaneously forms a macroscopic sheared flow, while conserving an overall crystalline structure. This flowing crystalline structure, which we call a "rheocrystal," is made possible by the condensation of shear along localized stacking faults. Numerical simulations very well reproduce the experimental observations and allow us to explore the parameter space. They demonstrate that the rheocrystal is induced neither by frustration nor by noise. They further show that larger systems flow faster while still remaining ordered.

Immotile Active Matter: Activity from Death and Reproduction

Unlike equilibrium atomic solids, biofilms-soft solids composed of bacterial cells-do not experience significant thermal fluctuations at the constituent level. However, living cells stochastically reproduce and die, provoking a mechanical response. We investigate the mechanical consequences of cellular death and reproduction by measuring surface-height fluctuations of biofilms containing two mutually antagonistic strains of Vibrio cholerae that kill one another on contact via the type VI secretion system. While studies of active matter typically focus on activity via constituent mobility, here, activity is mediated by reproduction and death events in otherwise immobilized cells. Biofilm surface topography is measured in the nearly homeostatic limit via white light interferometry. Although biofilms are far from equilibrium systems, measured surface-height fluctuation spectra resemble the spectra of thermal permeable membranes but with an activity-mediated effective temperature, as predicted by Risler, Peilloux, and Prost [Phys. Rev. Lett. 115, 258104 (2015)PRLTAO0031-900710.1103/PhysRevLett.115.258104]. By comparing the activity of killer strains of V. cholerae with that of genetically modified strains that cannot kill each other and validating with individual-based simulations, we demonstrate that extracted effective temperatures increase with the amount of death and reproduction and that death and reproduction can fluidize biofilms. Together, these observations demonstrate the unique physical consequences of activity mediated by death and reproduction events.

Phase Coexistence in Two-Dimensional Passive and Active Dumbbell Systems

We demonstrate that there is a macroscopic coexistence between regions with hexatic order and regions in the liquid or gas phase over a finite interval of packing fractions in active dumbbell systems with repulsive power-law interactions in two dimensions. In the passive limit, this interval remains finite, similar to what has been found in two-dimensional systems of hard and soft disks. We did not find discontinuous behavior upon increasing activity from the passive limit.

Control of active turbulence through addressable soft interfaces

We present an experimental study of a kinesin/tubulin active nematic formed at different oil interfaces. By tuning the interfacial rheology of the contacting oil, we have been able to condition and control the seemingly chaotic motion that characterizes the self-sustained active flows in our preparations. The active nematic is inherently unstable and spontaneously develops defects from an initial homogeneous state. We show that the steady state and, in particular, the density and dynamics of the defects strongly depends on the rheology of the contacting oil. Using a smectic-A thermotropic liquid crystal as the oil phase, we pattern the interface thanks to the anisotropy of the shear viscosity in this material. The geometry of the active nematic adapts to the boundary conditions at the interface by changing from the so-called active turbulent regime to laminar flows along the easy flow directions. The latter can be either a lattice of self-assembled circular paths or reconfigurable homogeneous orientations that can be addressed by means of an external magnetic field. We show that, under all confinement conditions, the spatiotemporal modes exhibited by the active liquid are consistent with a single intrinsic length scale, which can be tuned by the material parameters, and obey basic topological requirements imposed on the defects that drive the active flows. Future control strategies, including a tunable depleting agent, are discussed.

Taming active turbulence with patterned soft interfaces

Active matter embraces systems that self-organize at different length and time scales, often exhibiting turbulent flows apparently deprived of spatiotemporal coherence. Here, we use a layer of a tubulin-based active gel to demonstrate that the geometry of active flows is determined by a single length scale, which we reveal in the exponential distribution of vortex sizes of active turbulence. Our experiments demonstrate that the same length scale reemerges as a cutoff for a scale-free power law distribution of swirling laminar flows when the material evolves in contact with a lattice of circular domains. The observed prevalence of this active length scale can be understood by considering the role of the topological defects that form during the spontaneous folding of microtubule bundles. These results demonstrate an unexpected strategy for active systems to adapt to external stimuli, and provide with a handle to probe the existence of intrinsic length and time scales.

Active nematics consist of self-driven components that develop orientational order and turbulent flow. Here Guillamat et al. investigate an active nematic constrained in a quasi-2D geometrical setup and show that there exists an intrinsic length scale that determines the geometry in all forcing regimes.

Critical mingling and universal correlations in model binary active liquids

Ensembles of driven or motile bodies moving along opposite directions are generically reported to self-organize into strongly anisotropic lanes. Here, building on a minimal model of self-propelled bodies targeting opposite directions, we first evidence a critical phase transition between a mingled state and a phase-separated lane state specific to active particles. We then demonstrate that the mingled state displays algebraic structural correlations also found in driven binary mixtures. Finally, constructing a hydrodynamic theory, we single out the physical mechanisms responsible for these universal long-range correlations typical of ensembles of oppositely moving bodies.

Self-driven particles can drive the system out of equilibrium without the help of external forces. Here Bain et al. develop a model for two groups of motile bodies targeting opposite directions and show that the system exhibits critical behaviour as a function of particle density.

Gaussian theory for spatially distributed self-propelled particles

Obtaining a reduced description with particle and momentum flux densities outgoing from the microscopic equations of motion of the particles requires approximations. The usual method, we refer to as truncation method, is to zero Fourier modes of the orientation distribution starting from a given number. Here we propose another method to derive continuum equations for interacting self-propelled particles. The derivation is based on a Gaussian approximation (GA) of the distribution of the direction of particles. First, by means of simulation of the microscopic model, we justify that the distribution of individual directions fits well to a wrapped Gaussian distribution. Second, we numerically integrate the continuum equations derived in the GA in order to compare with results of simulations. We obtain that the global polarization in the GA exhibits a hysteresis in dependence on the noise intensity. It shows qualitatively the same behavior as we find in particles simulations. Moreover, both global polarizations agree perfectly for low noise intensities. The spatiotemporal structures of the GA are also in agreement with simulations. We conclude that the GA shows qualitative agreement for a wide range of noise intensities. In particular, for low noise intensities the agreement with simulations is better as other approximations, making the GA to an acceptable candidates of describing spatially distributed self-propelled particles.

Probing the shear viscosity of an active nematic film

In vitro reconstituted active systems, such as the adenosine triphosphate (ATP)-driven microtubule bundle suspension developed by the Dogic group [T. Sanchez, D. T. Chen, S. J. DeCamp, M. Heymann, and Z. Dogic, Nature (London) 491, 431 (2012)10.1038/nature11591], provide a fertile testing ground for elucidating the phenomenology of active liquid crystalline states. Controlling such novel phases of matter crucially depends on our knowledge of their material and physical properties. In this Rapid Communication, we show that the shear viscosity of an active nematic film can be probed by varying its hydrodynamic coupling to a bounding oil layer. Using the motion of disclinations as intrinsic tracers of the flow field and a hydrodynamic model, we obtain an estimate for the shear viscosity of the nematic film. Knowing this now provides us with an additional handle for robust and precision tunable control of the emergent dynamics of active fluids.

Control of active liquid crystals with a magnetic field

Living cells sense the mechanical features of their environment and adapt to it by actively remodeling their peripheral network of filamentary proteins, known as cortical cytoskeleton. By mimicking this principle, we demonstrate an effective control strategy for a microtubule-based active nematic in contact with a hydrophobic thermotropic liquid crystal. By using well-established protocols for the orientation of liquid crystals with a uniform magnetic field, and through the mediation of anisotropic shear stresses, the active nematic reversibly self-assembles with aligned flows and textures that feature orientational order at the millimeter scale. The turbulent flow, characteristic of active nematics, is in this way regularized into a laminar flow with periodic velocity oscillations. Once patterned, the microtubule assembly reveals its intrinsic length and time scales, which we correlate with the activity of motor proteins, as predicted by existing theories of active nematics. The demonstrated commanding strategy should be compatible with other viable active biomaterials at interfaces, and we envision its use to probe the mechanics of the intracellular matrix.

Topological defects in liquid crystals as templates for molecular self-assembly

Topological defects in liquid crystals (LCs) have been widely used to organize colloidal dispersions and template polymerization, leading to a range of assemblies, elastomers and gels. However, little is understood about molecular-level assembly processes within defects. Here, we report that nanoscopic environments defined by LC topological defects can selectively trigger processes of molecular self-assembly. By using fluorescence microscopy, cryogenic transmission electron microscopy and super-resolution optical microscopy, we observed signatures of molecular self-assembly of amphiphilic molecules in topological defects, including cooperativity, reversibility and controlled growth. We also show that nanoscopic o-rings synthesized from Saturn-ring disclinations and other molecular assemblies templated by defects can be preserved by using photocrosslinkable amphiphiles. Our results reveal that, in analogy to other classes of macromolecular templates such as polymer-surfactant complexes, topological defects in LCs are a versatile class of three-dimensional, dynamic and reconfigurable templates that can direct processes of molecular self-assembly.

Orientational order of motile defects in active nematics

The study of liquid crystals at equilibrium has led to fundamental insights into the nature of ordered materials, as well as to practical applications such as display technologies. Active nematics are a fundamentally different class of liquid crystals, driven away from equilibrium by the autonomous motion of their constituent rod-like particles. This internally generated activity powers the continuous creation and annihilation of topological defects, which leads to complex streaming flows whose chaotic dynamics seem to destroy long-range order. Here, we study these dynamics in experimental and computational realizations of active nematics. By tracking thousands of defects over centimetre-scale distances in microtubule-based active nematics, we identify a non-equilibrium phase characterized by a system-spanning orientational order of defects. This emergent order persists over hours despite defect lifetimes of only seconds. Similar dynamical structures are observed in coarse-grained simulations, suggesting that defect-ordered phases are a generic feature of active nematics.