Thermodynamic phases in two-dimensional active matter

Multiple temperatures and melting of a colloidal active crystal

Thermal fluctuations constantly excite all relaxation modes in an equilibrium crystal. As the temperature rises, these fluctuations promote the formation of defects and eventually melting. In active solids, the self-propulsion of "atomic" units provides an additional source of non-equilibrium fluctuations whose effect on the melting scenario is still largely unexplored. Here we show that when a colloidal crystal is activated by a bath of swimming bacteria, solvent temperature and active temperature cooperate to define dynamic and thermodynamic properties. Our system consists of repulsive paramagnetic particles confined in two dimensions and immersed in a bath of light-driven E. coli. The relative balance between fluctuations and interactions can be adjusted in two ways: by changing the strength of the magnetic field and by tuning activity with light. When the persistence time of active fluctuations is short, a single effective temperature controls both the amplitudes of relaxation modes and the melting transition. For more persistent active noise, energy equipartition is broken and multiple temperatures emerge, whereas melting occurs before the Lindemann parameter reaches its equilibrium critical value. We show that this phenomenology is fully confirmed by numerical simulations and framed within a minimal model of a single active particle in a periodic potential.

Extinction and coexistence in a binary mixture of proliferating motile disks

A binary mixture of two-different-size proliferating motile disks is studied. As growth is space limited, we focus on the conditions such that there is a coexistence of both large and small disks, or dominance of the larger disks. The study involves systematically varying some system parameters, such as diffusivities, growth rates, and self-propulsion velocities. In particular, we demonstrate that diffusing faster confers a competitive advantage, so that larger disks can in the long time coexist or even dominate the smaller ones. In the case of self-propelled disks, a coexistence regime is induced by the activity where the two types of disks show the same spatial distribution: both particles are phase separated or both are homogeneously distributed in the whole system.

Motility-induced coexistence of a hot liquid and a cold gas

If two phases exist at the same time, such as a gas and a liquid, they have the same temperature. This fundamental law of equilibrium physics is known to apply even to many non-equilibrium systems. However, recently, there has been much attention in the finding that inertial self-propelled particles like Janus colloids in a plasma or microflyers could self-organize into a hot gas-like phase that coexists with a colder liquid-like phase. Here, we show that a kinetic temperature difference across coexisting phases can occur even in equilibrium systems when adding generic (overdamped) self-propelled particles. In particular, we consider mixtures of overdamped active and inertial passive Brownian particles and show that when they phase separate into a dense and a dilute phase, both phases have different kinetic temperatures. Surprisingly, we find that the dense phase (liquid) cannot only be colder but also hotter than the dilute phase (gas). This effect hinges on correlated motions where active particles collectively push and heat up passive ones primarily within the dense phase. Our results answer the fundamental question if a non-equilibrium gas can be colder than a coexisting liquid and create a route to equip matter with self-organized domains of different kinetic temperatures.

Shaping active matter from crystalline solids to active turbulence

Active matter drives its constituent agents to move autonomously by harnessing free energy, leading to diverse emergent states with relevance to both biological processes and inanimate functionalities. Achieving maximum reconfigurability of active materials with minimal control remains a desirable yet challenging goal. Here, we employ large-scale, agent-resolved simulations to demonstrate that modulating the activity of a wet phoretic medium alone can govern its solid-liquid-gas phase transitions and, subsequently, laminar-turbulent transitions in fluid phases, thereby shaping its emergent pattern. These two progressively emerging transitions, hitherto unreported, bring us closer to perceiving the parallels between active matter and traditional matter. Our work reproduces and reconciles seemingly conflicting experimental observations on chemically active systems, presenting a unified landscape of phoretic collective dynamics. These findings enhance the understanding of long-range, many-body interactions among phoretic agents, offer new insights into their non-equilibrium collective behaviors, and provide potential guidelines for designing reconfigurable materials.

Adsorption of active polymers on attractive nanoparticles

The adsorption of active polymers on an attractive nanoparticle (NP) is studied using Langevin dynamics simulations. The active polymers consist of an active Brownian particle (ABP) at the head and a subsequent passive polymer chain. The ABP experiences an active force of magnitude Fa. The interactions between the active polymer and NP are modeled as Lennard-Jones potential with a strength εpn. We find the critical adsorption point εpn* increases with increasing the active force Fa. The increment of εpn*, denoted as Δεpn*, due to Fa can be expressed approximately as Δεpn* ∝ Fa2.5 for the restricted rotating active polymer (RRAP) where the rotation of the head ABP is restricted and Δεpn* ∝ Fa1.7 for the freely rotating active polymer (FRAP) where the ABP rotates freely. Meanwhile, the conformation of the adsorbed polymer, such as adsorbed trains on NP and the tail near the ABP, are also dependent on Fa. When the tail near the ABP is short, the adsorption is significantly affected by the active force. However, when the tail is long, the whole polymer can be viewed as a long tail stretched by the active force and unperturbed adsorption monomers. Simulation results show that the active force has a direct and significant effect on εpn* and the structure of the adsorbed active polymers.

Kinetic temperature and pressure of an active Tonks gas

Using computer simulation and analytical theory, we study an active analog of the well-known Tonks gas, where active Brownian particles are confined to a periodic one-dimensional (1D) channel. By introducing the notion of a kinetic temperature, we derive an accurate analytical expression for the pressure and clarify the paradoxical behavior where active Brownian particles confined to 1D exhibit anomalous clustering but no motility-induced phase transition. More generally, this work provides a deeper understanding of pressure in active systems as we uncover a unique link between the kinetic temperature and swim pressure valid for active Brownian particles in higher dimensions.

Entropy production and collective excitations of crystals out of equilibrium: The concept of entropons

We study the collective vibrational excitations of crystals under out-of-equilibrium steady conditions that give rise to entropy production. Their excitation spectrum comprises equilibriumlike phonons of thermal origin and additional collective excitations called entropons because each of them represents a mode of spectral entropy production. Entropons coexist with phonons and dominate them when the system is far from equilibrium while they are negligible in near-equilibrium regimes. The concept of entropons has been recently introduced and verified in a special case of crystals formed by self-propelled particles. Here we show that entropons exist in a broader class of active crystals that are intrinsically out of equilibrium and characterized by the lack of detailed balance. After a general derivation, several explicit examples are discussed, including crystals consisting of particles with alignment interactions and frictional contact forces.

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.

Two-Dimensional Crystals far from Equilibrium

When driven by nonequilibrium fluctuations, particle systems may display phase transitions and physical behavior with no equilibrium counterpart. We study a two-dimensional particle model initially proposed to describe driven non-Brownian suspensions undergoing nonequilibrium absorbing phase transitions. We show that when the transition occurs at large density, the dynamics produces long-range crystalline order. In the ordered phase, long-range translational order is observed because equipartition of energy is lacking, phonons are suppressed, and density fluctuations are hyperuniform. Our study offers an explicit microscopic model where nonequilibrium violations of the Mermin-Wagner theorem stabilize crystalline order in two dimensions.

Entropons as collective excitations in active solids

The vibrational dynamics of solids is described by phonons constituting basic collective excitations in equilibrium crystals. Here, we consider a non-equilibrium active solid, formed by self-propelled particles, which bring the system into a non-equilibrium steady-state. We identify novel vibrational collective excitations of non-equilibrium (active) origin, which coexist with phonons and dominate over them when the system is far from equilibrium. These vibrational excitations are interpreted in the framework of non-equilibrium physics, in particular, stochastic thermodynamics. We call them "entropons" because they are the modes of spectral entropy production (at a given frequency and wave vector). The existence of entropons could be verified in future experiments on dense self-propelled colloidal Janus particles and granular active matter, as well as in living systems, such as dense cell monolayers.

Motility-induced shear thickening in dense colloidal suspensions

Phase transitions and collective dynamics of active colloidal suspensions are fascinating topics in soft matter physics, particularly for out-of-equilibrium systems, which can lead to rich rheological behaviours in the presence of steady shear flow. Here the role of self-propulsion in the rheological response of a dense colloidal suspension is investigated by using particle-resolved Brownian dynamics simulations. First, the combined effect of activity and shear in the solid on the disordering transition of the suspension is analyzed. While both self-propulsion and shear destroy order and melt the system if critical values are exceeded, self-propulsion largely lowers the stress barrier needed to be overcome during the transition. We further explore the rheological response of the active sheared system once a steady state is reached. While passive suspensions show a solid-like behaviour, turning on particle motility fluidises the system. At low self-propulsion, the active suspension behaves in the steady state as a shear-thinning fluid. Increasing the self-propulsion changes the behaviour of the liquid from shear-thinning to shear-thickening. We attribute this to clustering in the sheared suspensions induced by motility. This new phenomenon of motility-induced shear thickening (MIST) can be used to tailor the rheological response of colloidal suspensions.

Detecting Liquid-Liquid Phase Separations Using Molecular Dynamics Simulations and Spectral Clustering

A stringent test of the accuracy of empirical force fields is reproducing the phase diagram of bulk phases and mixtures. Exploring the phase diagram of mixtures requires the detection of phase boundaries and critical points. In contrast to most solid-liquid transitions, in which a global order parameter (average density) can be used to discriminate between two phases, some demixing transitions entail relatively subtle changes in the local environment of each molecule. In such cases, finite sampling errors and finite-size effects make the identification of trends in local order parameters extremely challenging. Here we analyze one such example, namely a methanol/hexane mixture, and compute several local and global structural properties. We simulate the system at various temperatures and study the structural changes associated with demixing. We show that despite a seemingly continuous transformation between mixed and demixed states, the topological properties of the H-bond network change abruptly as the system crosses the demixing line. In particular, by using spectral clustering, we show that the distribution of cluster sizes develops a fat tail (as expected from percolation theory) in the vicinity of the critical point. We illustrate a simple criterion to identify this behavior, which results from the emergence of large system-spanning clusters from a collection of aggregates. We further tested the spectral clustering analysis on a Lennard-Jones system as a standard example of a system with no H-bonds, and also, in this case, we were able to detect the demixing transition.

Inertia changes evolution of motility-induced phase separation in active matter across particle activity

The effects of inertia in active matter and motility-induced phase separation (MIPS) have attracted growing interest but still remain poorly studied. We studied MIPS behavior in the Langevin dynamics across a broad range of particle activity and damping rate values with molecular dynamic simulations. Here we show that the MIPS stability region across particle activity values consists of several domains separated by discontinuous or sharp changes in susceptibility of mean kinetic energy. These domain boundaries have fingerprints in the system's kinetic energy fluctuations and characteristics of gas, liquid, and solid subphases, such as the number of particles, densities, or the power of energy release due to activity. The observed domain cascade is most stable at intermediate damping rates but loses its distinctness in the Brownian limit or vanishes along with phase separation at lower damping values.

Alignment rule and geometric confinement lead to stability of a vortex in active flow

Vortices are hallmarks of a wide range of nonequilibrium phenomena in fluids at multiple length scales. In this work, we numerically study the whirling motion of self-propelled soft point particles confined in circular domain, and aim at addressing the stability issue of the coherent vortex structure. By the combination of dynamical and statistical analysis at the individual particle level, we reveal the persistence of the whirling motion resulting from the subtle competition of activity and geometric confinement. In the stable whirling motion, the scenario of the coexistence of the irregular microscopic motions of individual particles and the regular global whirling motion is fundamentally different from the motion of a vortex in passive fluid. Possible orientational order coexisting with the whirling are further explored. This work shows the stability mechanism of vortical dynamics in active media under the alignment rule in confined space and may have implications in creating and harnessing macroscale coherent dynamical states by tuning the confining geometry.

Collective behavior of soft self-propelled disks with rotational inertia

We investigate collective properties of a large system of soft self-propelled inertial disks with active Langevin dynamics simulation in two dimensions. Rotational inertia of the disks is found to favor motility induced phase separation (MIPS), due to increased effective persistence of the disks. The MIPS phase diagram in the parameter space of rotational inertia and disk softness is reported over a range of values of translation inertia and self-propulsion strength of the disks. Our analytical prediction of the phase boundary between the homogeneous (no-MIPS) and MIPS state in the limit of small and large rotational inertia is found to agree with the numerical data over a large range of translational inertia. Shape of the high density MIPS phase is found to change from circular to rectangular one as the system moves away from the phase boundary. Structural and dynamical properties of the system, measured by several physical quantities, are found to be invariant in the central region of the high density MIPS phase, whereas they are found to vary gradually near the peripheral region of the high density phase. Importantly, the width of the peripheral region near the phase boundary is much larger compared to the narrow peripheral region far away from the phase boundary. Rich dynamics of the disks inside the high density MIPS phase is addressed. Spatial correlation of velocity of the disks is found to increase with rotational inertia and disk hardness. However, temporal correlation of the disks' velocity is found to be a function of rotational inertia, while it is independent of disk softness.

Biomolecular condensate phase diagrams with a combinatorial microdroplet platform

The assembly of biomolecules into condensates is a fundamental process underlying the organisation of the intracellular space and the regulation of many cellular functions. Mapping and characterising phase behaviour of biomolecules is essential to understand the mechanisms of condensate assembly, and to develop therapeutic strategies targeting biomolecular condensate systems. A central concept for characterising phase-separating systems is the phase diagram. Phase diagrams are typically built from numerous individual measurements sampling different parts of the parameter space. However, even when performed in microwell plate format, this process is slow, low throughput and requires significant sample consumption. To address this challenge, we present here a combinatorial droplet microfluidic platform, termed PhaseScan, for rapid and high-resolution acquisition of multidimensional biomolecular phase diagrams. Using this platform, we characterise the phase behaviour of a wide range of systems under a variety of conditions and demonstrate that this approach allows the quantitative characterisation of the effect of small molecules on biomolecular phase transitions.

Liquid-hexatic-solid phases in active and passive Brownian particles determined by stochastic birth and death events

We study the effects of stochastic birth and death processes on the structural phases of systems of active and passive Brownian particles subject to volume exclusion. The total number of particles in the system is a fluctuating quantity, determined by the birth and death parameters and on the activity of the particles. As the birth and death parameters are varied, we find liquid, hexatic, and solid phases. For passive particles, these phases are found to be spatially homogeneous. For active particles, motility-induced phase separation (coexisting hexatic and liquid phases) occurs for large activity and sufficiently small birth rates. We also observe a reentrant transition to the hexatic phase when the birth rate is increased. This results from a balance of an increasing number of particles filling the system, and a larger number of defects resulting from the birth and death dynamics.

Active Refrigerators Powered by Inertia

We present the operational principle for a refrigerator that uses inertial effects in active Brownian particles to locally reduce their (kinetic) temperature by 2 orders of magnitude below the environmental temperature. This principle exploits the peculiar but so-far unknown shape of the phase diagram of inertial active Brownian particles to initiate motility-induced phase separation in the targeted cooling regime only. Remarkably, active refrigerators operate without requiring isolating walls opening the route toward using them to systematically absorb and trap, e.g., toxic substances from the environment.

Motility-induced phase separation of self-propelled soft inertial disks

The phase diagram of the phenomenon of motility-induced phase separation (MIPS) for a collection of self-propelled interacting disks over a large inertial range is explored using active Langevin dynamics simulation with particular emphasis on disk softness and effective size. It is shown that the parabola-like phase boundary between the homogeneous and MIPS states in the semi-log space of disk softness and effective size moves towards the hard disk limit with increase in inertia, before complete disappearance in the limit of large inertia. With increase in effective size of the disks, re-entrant phase separation, that is the system phase-separating from a homogeneous phase and eventually re-entering the homogeneous phase, is reported. The structural and the dynamical properties of the different phases are investigated in the considered inertial range. The particular shape of the phase boundary and the re-entrant behavior is explained based on several qualitative and quantitative results. Unlike most of the earlier studies on MIPS, which consider hard particle limits, our findings may be directly applicable to soft active matter for a range of physical and biological systems.

Reentrant phase separation of a sparse collection of nonreciprocally aligning self-propelled disks

We study a model of aligning self-propelled disks that nonreciprocally reorient the self-propulsion directions along the interparticle separation and towards the other disks. In the limit of small inertia and large softness, where conventional motility-induced phase separation is absent, we demonstrate that the homogeneous system at a small area fraction phase-separates into clusters and a low-density phase that, eventually, reenters the homogeneous phase with a monotonic increase in alignment strength. The disks inside the clusters move with a finite space-dependent speed, constantly shuttling between clusters through the surrounding low-density homogeneous phase while maintaining the hexatic structure properties within the clusters. The area fraction gradually increases from the periphery towards the center of the clusters with a negligible correlation of the velocity and propulsion direction inside the clusters. The novel collective behavior of reentrant phase separation is found to follow from both the limits of hard disks and extremely small inertia, tending towards the overdamped limit. However, important differences in the structural and dynamical properties are shown in the limit of hard disks and extremely small inertia, as compared to that for soft disks at finite inertia. We show that the cluster phase is associated with an effective temperature for a wide range of values of alignment strength, whereas an effective temperature is associated with the specific range of alignment in the low-density phase. We believe that the reentrant phase behavior in the limit of small area fraction and the remarkable properties of the clusters should be useful in understanding a wide range of physics issues, ranging from clogging and unclogging to information exchange and transport, in biological and synthetic self-propelled systems.

Annealing and melting of active two-dimensional soliton lattices in chiral nematic films

In this work, thousands of electrically driven dissipative solitons, called directrons, are generated in a chiral nematic liquid crystal. The directrons start with random motions but soon synchronize their motions and self-organize into a two-dimensional hexagonal lattice. The directron lattice moves collectively and forms a hexatic phase. By increasing the applied voltage, the lattice exhibits a first-order hexatic-to-liquid phase transition.

Nonequilibrium thermodynamics of Janus particle self-assembly

We compute the energetic cost of formation of Janus particle structures. Using an approach that couples particle dynamics to the evolution of fuel concentration in the medium, which we consider to be initially inhomogeneous, we show the different types of emerging structures. The energy dissipated in the formation of such structures is obtained from the entropy production rate, which is a non-monotonic function of the fraction of assembled particles and, thus, different in each self-assembly regime. An analysis of the free energy of these particles allows us to establish a thermodynamic criterion of structure formation based on the behavior of chemical potential as a function of the fraction of assembled particles.

Hydrodynamic effects on the liquid-hexatic transition of active colloids

We study numerically the role of hydrodynamics in the liquid-hexatic transition of active colloids at intermediate activity, where motility induced phase separation (MIPS) does not occur. We show that in the case of active Brownian particles (ABP), the critical density of the transition decreases upon increasing the particle's mass, enhancing ordering, while self-propulsion has the opposite effect in the activity regime considered. Active hydrodynamic particles (AHP), instead, undergo the liquid-hexatic transition at higher values of packing fraction [Formula: see text] than the corresponding ABP, suggesting that hydrodynamics have the net effect of disordering the system. At increasing densities, close to the hexatic-liquid transition, we found in the case of AHP the appearance of self-sustained organized motion with clusters of particles moving coherently.

Geometry-controlled phase transition in vibrated granular media

We report experiments on the dynamics of vibrated particles constrained in a two-dimensional vertical container, motivated by the following question: how to get the most out of a given external vibration to maximize internal disorder (e.g. to blend particles) and agitation (e.g. to absorb vibrations)? Granular media are analogs to classical thermodynamic systems, where the injection of energy can be achieved by shaking them: fluidization arises by tuning either the amplitude or the frequency of the oscillations. Alternatively, we explore what happens when another feature, the container geometry, is modified while keeping constant the energy injection. Our method consists in modifying the container base into a V-shape to break the symmetries of the inner particulate arrangement. The lattice contains a compact hexagonal solid-like crystalline phase coexisting with a loose amorphous fluid-like phase, at any thermal agitation. We show that both the solid-to-fluid volume fraction and the granular temperature depend not only on the external vibration but also on the number of topological defects triggered by the asymmetry of the container. The former relies on the statistics of the energy fluctuations and the latter is consistent with a two-dimensional melting transition described by the KTHNY theory.

Percolation transition in phase-separating active fluid

The motility-induced phase separation exhibited by active particles with repulsive interactions is well known. We show that the interaction softness of active particles destabilizes the highly ordered dense phase, leading to the formation of a porous cluster which spans the system. This soft limit can also be achieved if the particle motility is increased beyond a critical value, at which the system clearly exhibits all the characteristics of a standard percolation transition. We also show that in the athermal limit, active particles exhibit similar transitions even at low motility. With these additional new phases, the phase diagram of repulsive active particles is revealed to be richer than what was previously conceived.

Critical behavior of active Brownian particles: Connection to field theories

We explore the relation between active Brownian particles, a minimal particle-based model for active matter, and scalar field theories. Both show a liquid-gas-like phase transition toward stable coexistence of a dense liquid with a dilute active gas that terminates in a critical point. However, a comprehensive mapping between the particle-based model parameters and the effective coefficients governing the field theories has not been established yet. We discuss conflicting recent numerical results for the critical exponents of active Brownian particles in two dimensions. Starting from the intermediate effective hydrodynamic equations, we then present a construction for a scalar order parameter for active Brownian particles that yields the active model B+. We argue that a crucial ingredient is the coupling between density and polarization in the particle current. The renormalization flow close to two dimensions exhibits a pair of perturbative fixed points that limit the attractive basin of the Wilson-Fisher fixed point, with the perspective that the critical behavior of active Brownian particles in two dimensions is governed by a strong-coupling fixed point different from Wilson-Fisher and not necessarily corresponding to Ising universality.

Biased-angle effect on diffusion dynamics and phase separation in anisotropic active particle system

A deep understanding for collective behavior in an active matter system with complex interactions has far-reaching impact in biology. In the present work, we adopt Langevin dynamics simulations to investigate diffusion dynamics and phase separation in an anisotropic active particle system with a tunable biased angle α defined as the deviation between the active force direction and anisotropic orientation. Our results demonstrate that the biased angle can induce super-rotational diffusion dynamics characterized by a power-law relationship between the mean square angle displacement (MSAD) and the time interval Δ t in the form of MSAD ∼ Δ t β with β > 1 and also result in non-trivial phase separation kinetics. As activity is dominant, nucleation time shows a non-monotonic dependence on the biased angle. Moreover, there arises a distinct transition of phase separation, from spinodal decomposition without apparent nucleation time to binodal decomposition with prominent nucleation delay. A significant inhibition effect occurs at right and obtuse angles, where the remarkable super-rotational diffusion prevents particle aggregation, leading to a slow nucleation process. As active force is competitive to anisotropic interactions, the system is almost homogeneous, while, intriguingly, we observe a re-entrant phase separation as a small acute angle is introduced. The prominent super-rotational diffusion under small angles provides an optimum condition for particle adsorption and cluster growth and, thus, accounts for the re-entrance of phase separation. A consistent scenario for the physical mechanism of our observations is achieved by properly considering the modulation of the biased angle on the interplay between activity and anisotropic interactions.

Active Brownian motion of strongly coupled charged grains driven by laser radiation in plasma

The systems of active Brownian grains can be considered as open systems, in which there is an exchange of energy and matter with the environment. The collective phenomena of active Brownian grains can demonstrate analogies with ordinary phase transitions. We study the active Brownian motion of light-absorbing and strongly interacting grains far from equilibrium suspended in gas discharge under laser irradiation when the nature and intensity of the active motion depend on the effect of radiation. Active Brownian motion is caused by photophoresis, i.e., absorption of laser radiation at the metal-coated surface of the grain creates radiometric force, which in turn drives the grains. We experimentally observed the active Brownian motion of charged grains in the transition of the grain monolayer from the solid to liquid state. An analysis of the character of motion, including the mean-square and linear displacement and persistence length at various values of the randomization (coupling parameter) of the grain structure, was presented.

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.

Electrically tunable collective motion of dissipative solitons in chiral nematic films

From the motion of fish and birds, to migrating herds of ungulates, collective motion has attracted people for centuries. Active soft matter exhibits a plethora of emergent dynamic behaviors that mimic those of biological systems. Here we introduce an active system composed of dynamic dissipative solitons, i.e. directrons, which mimics the collective motion of living systems. Although the directrons are inanimate, artificial particle-like solitonic field configurations, they locally align their motions like their biological counterparts. Driven by external electric fields, hundreds of directrons are generated in a chiral nematic film. They start with random motions but self-organize into flocks and synchronize their motions. The directron flocks exhibit rich dynamic behaviors and induce population density fluctuations far larger than those in thermal equilibrium systems. They exhibit “turbulent” swimming patterns manifested by transient vortices and jets. They even distinguish topological defects, heading towards defects of positive topological strength and avoiding negative ones.

Understanding of the collective motion seen in biological systems is crucial for development of autonomous robots and swarm computing. The authors introduce an experimental platform with liquid crystal driven by external electric field, that mimics the collective motion of living systems.

Intrinsic structure perspective for MIPS interfaces in two-dimensional systems of active Brownian particles

Suspensions of active Brownian particles (ABPs) undergo motility-induced phase separation (MIPS) over a wide range of mean density and activity strength, which implies the spontaneous aggregation of particles due to the persistence of their direction of motion, even in the absence of an explicit attraction. Both similarities and qualitative differences have been obtained when the MIPS is analysed in the same terms as a liquid-gas phase coexistence in an equilibrium attractive system. Negative values of the mechanical surface tension have been reported, from the total forces across the interface, while stable fluctuations of the interfacial line could be interpreted as a positive capillary surface tension; in equilibrium liquid surfaces, these two magnitudes are equal. We present here the analysis of 2D-ABP interfaces in terms of the intrinsic density and force profiles, calculated with the particle distance to the instantaneous interfacial line. Our results provide new insight into the origin of MIPS from the local rectification of the random active force on the particles near the interface. As has been reported, this effect acts as an external potential that produces a pressure gradient across the interface, such that the mechanical surface tension of the MIPS cannot be described as that of equilibrium coexisting phases; however, our analysis shows that most of that effect comes from the tightly caged particles at the dense (inner) side of the MIPS interface, rather than from the free moving particles at the outer side that collide with the dense cluster. Moreover, a clear correlation appears between the decay of the hexatic order parameter at the dense slab and the end of the MIPS as the strength of the active force is lowered. We show that, using the strong active forces required for MIPS, the interfacial structure and properties are very similar for ABPs with purely repulsive interactions (the Weeks-Chandler-Andersen-Lennard-Jones (WCA-LJ) model truncated at its minimum) and when the interaction includes a range of the LJ attractive force.

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.

Rapid Determination of Phase Diagrams for Biomolecular Liquid-Liquid Phase Separation with Microfluidics

Biomolecular phase separation is currently emerging in both the medical and life science fields. Meanwhile, the application of liquid-liquid phase separation has been extended to many fields including drug discovery, fibrous material fabrication, 3D printing, and polymer design. Although more than 8600 proteins and other synthetic macromolecules are capable of phase separation as recently reported, there is still a lack of a high-throughput approach to quantitatively characterize its phase behaviors. To meet this requirement, here, we proposed fast and high-resolution acquisition of biomolecular phase diagrams using microfluidic chips. Using this platform, we demonstrated the phase behavior of polyU/RRASLRRASLRRASL in a quantitative manner. Up to 1750 concentration conditions can be generated in 140 min. The detection limitation of our device to capture the saturation concentration for phase separation is about 5 times lower than that of the traditional turbidity method. Thus, our results provide a basis for the rapid acquisition of phase diagrams with high-throughput and pave the way for its wide application.

Dynamic overlap concentration scale of active colloids

By introducing the notion of a dynamic overlap concentration scale, we identify additional universal features of the mechanical properties of active colloids. We codify these features by recognizing that the characteristic length scale of an active particle's trajectory, the run length, introduces a concentration scale ϕ^{*}. Large-scale simulations of repulsive active Brownian particles (ABPs) confirm that this run-length dependent concentration, the trajectory-space analog of the overlap concentration in polymer solutions, delineates distinct concentration regimes in which interparticle collisions alter particle trajectories. Using ϕ^{*} and concentration scales associated with colloidal jamming, the mechanical equation of state for ABPs collapses onto a set of principal curves that contain several overlooked features. The inclusion of these features qualitatively alters previous predictions of the behavior for active colloids, as we demonstrate by computing the spinodal for a suspension of purely repulsive ABPs. Our findings suggest that dynamic overlap concentration scales should help unravel the behavior of active and driven systems.

Data-driven criterion for the solid-liquid transition of two-dimensional self-propelled colloidal particles far from equilibrium

We establish an explicit data-driven criterion for identifying the solid-liquid transition of two-dimensional self-propelled colloidal particles in the far from equilibrium parameter regime, where the transition points predicted by different conventional empirical criteria for melting and freezing diverge. This is achieved by applying a hybrid machine learning approach that combines unsupervised learning with supervised learning to analyze a huge amount of the system's configurations in the nonequilibrium parameter regime on an equal footing. Furthermore, we establish a generic data-driven evaluation function, according to which the performance of different empirical criteria can be systematically evaluated and improved. In particular, by applying this evaluation function, we identify a new nonequilibrium threshold value for the long-time diffusion coefficient, based on which the predictions of the corresponding empirical criterion are greatly improved in the far from equilibrium parameter regime. These data-driven approaches provide a generic tool for investigating phase transitions in complex systems where conventional empirical ones face difficulties.

Critical yielding rheology: from externally deformed glasses to active systems

We use extensive computer simulations to study the yielding transition under two different loading schemes: standard simple shear dynamics and self-propelled dense active systems. In the active systems, a yielding transition toward an out-of-equilibrium flowing state known as the liquid phase is observed when self-propulsion is increased. The range of self-propulsions in which this pure liquid regime exists appears to vanish upon approaching the so-called 'jamming point' at which the solidity of soft-sphere packings is lost. Such an 'active yielding' transition shares similarities with the generic yielding transition for shear flows. A Herschel-Bulkley law is observed along the liquid regime in both loading scenarios, with a clear difference in the critical scaling exponents between the two, suggesting the existence of different universality classes for the yielding transition under different driving conditions. In addition, we present the direct measurements of growing length and time scales for both driving scenarios. A comparison with theoretical predictions from the recent literature reveals poor agreement with our numerical results.

Kinetic Monte Carlo Algorithms for Active Matter Systems

We study kinetic Monte Carlo (KMC) descriptions of active particles. We show that, when they rely on purely persistent, active steps, their continuous-time limit is ill-defined, leading to the vanishing of trademark behaviors of active matter such as the motility-induced phase separation, ratchet effects, as well as to a diverging mechanical pressure. We then show how, under an appropriate scaling, mixing passive steps with active ones leads to a well-defined continuous-time limit that however differs from standard active dynamics. Finally, we propose new KMC algorithms whose continuous-time limits lead to the dynamics of active Ornstein-Uhlenbeck, active Brownian, and run-and-tumble particles.

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.

Machine learning for phase behavior in active matter systems

We demonstrate that deep learning techniques can be used to predict motility-induced phase separation (MIPS) in suspensions of active Brownian particles (ABPs) by creating a notion of phase at the particle level. Using a fully connected network in conjunction with a graph neural network we use individual particle features to predict to which phase a particle belongs. From this, we are able to compute the fraction of dilute particles to determine if the system is in the homogeneous dilute, dense, or coexistence region. Our predictions are compared against the MIPS binodal computed from simulation. The strong agreement between the two suggests that machine learning provides an effective way to determine the phase behavior of ABPs and could prove useful for determining more complex phase diagrams.

Influence of interaction softness on phase separation of active particles

Using a minimal model of active Brownian particles, we study the effect of a crucial parameter, namely the softness of the interparticle repulsion, on motility-induced phase separation. We show that an increase in particle softness reduces the ability of the system to phase separate and the system exhibits a delayed transition. After phase separation, the system state properties can be explained by a single relevant length scale, the effective interparticle distance. We estimate this length scale analytically and use it to rescale the state properties at dense phase for systems with different interaction softness. Using this length scale, we provide a scaling relation for the time taken to phase separate which shows a high sensitivity to the interaction softness.

Phase Diagram of Active Brownian Spheres: Crystallization and the Metastability of Motility-Induced Phase Separation

Motility-induced phase separation (MIPS), the phenomenon in which purely repulsive active particles undergo a liquid-gas phase separation, is among the simplest and most widely studied examples of a nonequilibrium phase transition. Here, we show that states of MIPS coexistence are in fact only metastable for three-dimensional active Brownian particles over a very broad range of conditions, decaying at long times through an ordering transition we call active crystallization. At an activity just above the MIPS critical point, the liquid-gas binodal is superseded by the crystal-fluid coexistence curve, with solid, liquid, and gas all coexisting at the triple point where the two curves intersect. Nucleating an active crystal from a disordered fluid, however, requires a rare fluctuation that exhibits the nearly close-packed density of the solid phase. The corresponding barrier to crystallization is surmountable on a feasible timescale only at high activity, and only at fluid densities near maximal packing. The glassiness expected for such dense liquids at equilibrium is strongly mitigated by active forces, so that the lifetime of liquid-gas coexistence declines steadily with increasing activity, manifesting in simulations as a facile spontaneous crystallization at extremely high activity.

Active particles with polar alignment in ring-shaped confinement

We study steady-state properties of active, nonchiral and chiral Brownian particles with polar alignment and steric interactions confined within a ring-shaped confinement (annulus) in two dimensions. Exploring possible interplays between polar interparticle alignment, geometric confinement and the surface curvature, being incorporated here on minimal levels, we report a surface-population reversal effect, whereby active particles migrate from the outer concave boundary of the annulus to accumulate on its inner convex boundary. This contrasts the conventional picture, implying stronger accumulation of active particles on concave boundaries relative to the convex ones. The population reversal is caused by both particle alignment and surface curvature, disappearing when either of these factors is absent. We explore the ensuing consequences for the chirality-induced current and swim pressure of active particles and analyze possible roles of system parameters, such as the mean number density of particles and particle self-propulsion, chirality, and alignment strengths.

Wetting of a solid surface by active matter

A lattice model is used to study repulsive active particles at a planar surface. A rejection-free Kinetic Monte Carlo method is employed to characterize the wetting behaviour. The model predicts a motility-induced phase separation of active particles, and the bulk coexistence of dense liquid-like and dilute vapour-like steady states is determined. An "ensemble", with a varying number of particles, analogous to a grand canonical ensemble in equilibrium, is introduced. The formation and growth of the liquid film between the solid surface and the vapour phase is investigated. At constant activity, as the system is brought towards coexistence from the vapour side, the thickness of the adsorbed film exhibits a divergent behaviour regardless of the activity. This suggests a complete wetting scenario along the full coexistence curve.

Statistical mechanics of active Ornstein-Uhlenbeck particles

We study the statistical properties of active Ornstein-Uhlenbeck particles (AOUPs). In this simplest of models, the Gaussian white noise of overdamped Brownian colloids is replaced by a Gaussian colored noise. This suffices to grant this system the hallmark properties of active matter, while still allowing for analytical progress. We study in detail the steady-state distribution of AOUPs in the small persistence time limit and for spatially varying activity. At the collective level, we show AOUPs to experience motility-induced phase separation both in the presence of pairwise forces or due to quorum-sensing interactions. We characterize both the instability mechanism leading to phase separation and the resulting phase coexistence. We probe how, in the stationary state, AOUPs depart from their thermal equilibrium limit by investigating the emergence of ratchet currents and entropy production. In the small persistence time limit, we show how fluctuation-dissipation relations are recovered. Finally, we discuss how the emerging properties of AOUPs can be characterized from the dynamics of their collective modes.

Pair correlation of dilute active Brownian particles: From low-activity dipolar correction to high-activity algebraic depletion wings

We study the pair correlation of active Brownian particles at low density using numerical simulations and analytical calculations. We observe a winged pair correlation: While particles accumulate in front of an active particle as expected, the depletion wake consists of two depletion wings. In the limit of soft particles, we obtain a closed equation for the pair correlation, allowing us to characterize the depletion wings. In particular, we unveil two regimes at high activity, where the wings adopt a self-similar profile and decay algebraically. We also perform experiments of self-propelled Janus particles and indeed observe the depletion wings.

Phase Separation and Multibody Effects in Three-Dimensional Active Brownian Particles

Simulation studies of the phase diagram of repulsive active Brownian particles in three dimensions reveal that the region of motility-induced phase separation between a high and low density phase is enclosed by a region of gas-crystal phase separation. Near-critical loci and structural crossovers can additionally be identified in analogy with simple fluids. Motivated by the striking similarity to the behavior of equilibrium fluids with short-ranged pairwise attractions, we show that a direct mapping to pair potentials in the dilute limit implies interactions that are insufficiently attractive to engender phase separation. Instead, this is driven by the emergence of multibody effects associated with particle caging that occurs at sufficiently high number density. We quantify these effects via information-theoretical measures of n-body effective interactions extracted from the configurational structure.

The "isothermal" compressibility of active matter

We demonstrate that the mechanically defined "isothermal" compressibility behaves as a thermodynamic-like response function for suspensions of active Brownian particles. The compressibility computed from the active pressure-a combination of the collision and unique swim pressures-is capable of predicting the critical point for motility induced phase separation, as expected from the mechanical stability criterion. We relate this mechanical definition to the static structure factor via an active form of the thermodynamic compressibility equation and find the two to be equivalent, as would be the case for equilibrium systems. This equivalence indicates that compressibility behaves like a thermodynamic response function, even when activity is large. Finally, we discuss the importance of the phase interface when defining an active chemical potential. Previous definitions of the active chemical potential are shown to be accurate above the critical point but breakdown in the coexistence region. Inclusion of the swim pressure in the mechanical compressibility definition suggests that the interface is essential for determining phase behavior.

Motility-Induced Microphase and Macrophase Separation in a Two-Dimensional Active Brownian Particle System

As a result of nonequilibrium forces, purely repulsive self-propelled particles undergo macrophase separation between a dense and a dilute phase. We present a thorough study of the ordering kinetics of such motility-induced phase separation (MIPS) in active Brownian particles in two dimensions, and we show that it is generically accompanied by microphase separation. The growth of the dense phase follows a law akin to the one of liquid-gas phase separation. However, it is made of a mosaic of hexatic microdomains whose size does not coarsen indefinitely, leaving behind a network of extended topological defects from which microscopic dilute bubbles arise. The characteristic length of these finite-size structures increases with activity, independently of the choice of initial conditions.

Self-Organized Critical Coexistence Phase in Repulsive Active Particles

We revisit motility-induced phase separation in two models of active particles interacting by pairwise repulsion and uncover new qualitative features: the resulting dense phase contains gas bubbles distributed algebraically up to a typically extremely large cutoff scale. At large enough system size and/or global density, all the gas may be contained inside the bubbles, at which point the system is microphase separated with a finite cutoff bubble scale. We further observe that the ordering is clearly anomalous, with different dynamics for the coarsening of the dense phase and of the gas bubbles. This self-organized critical phenomenology is reproduced by a "reduced bubble model" that implements the basic idea of reverse Ostwald ripening put forward in Tjhung et al. [Phys. Rev. X 8, 031080 (2018)PRXHAE2160-330810.1103/PhysRevX.8.031080].