Nonequilibrium strongly hyperuniform fluids of circle active particles with large local density fluctuations

Circular swimming motility and disordered hyperuniform state in an algae system

Active matter comprises individually driven units that convert locally stored energy into mechanical motion. Interactions between driven units lead to a variety of nonequilibrium collective phenomena in active matter. One of such phenomena is anomalously large density fluctuations, which have been observed in both experiments and theories. Here we show that, on the contrary, density fluctuations in active matter can also be greatly suppressed. Our experiments are carried out with marine algae ([Formula: see text]), which swim in circles at the air-liquid interfaces with two different eukaryotic flagella. Cell swimming generates fluid flow that leads to effective repulsions between cells in the far field. The long-range nature of such repulsive interactions suppresses density fluctuations and generates disordered hyperuniform states under a wide range of density conditions. Emergence of hyperuniformity and associated scaling exponent are quantitatively reproduced in a numerical model whose main ingredients are effective hydrodynamic interactions and uncorrelated random cell motion. Our results demonstrate the existence of disordered hyperuniform states in active matter and suggest the possibility of using hydrodynamic flow for self-assembly in active matter.

Emergent stereoselective interactions and self-recognition in polar chiral active ellipsoids

In many active matter systems, particle trajectories have a well-defined handedness or chirality. Whether such chiral activity can introduce stereoselective interactions between particles is not known. Here, we developed a strategy to tune the nature of chiral activity of three-dimensionally printed granular ellipsoids without altering their shape or size. In vertically agitated monolayers of these particles, we observed two types of dimers form depending on the chirality of the pairing monomers. Heterochiral dimers moved collectively as a single achiral active unit, while homochiral ones formed a translationally immobile spinner. In active racemic mixtures, the former was more abundant than the latter, indicating that interactions were stereoselective. Through dimer lifetime measurements, we further provide evidence for chiral self-recognition in mixtures of particles with different chiral activities. We lastly show that, at fixed particle number density, changing the net chirality of a dense active liquid fundamentally alters the nature of collective relaxation.

Connecting relaxation time to a dynamical length scale in athermal active glass formers

Supercooled liquids display dynamics that are inherently heterogeneous in space. This essentially means that at temperatures below the melting point, particle dynamics in certain regions of the liquid can be orders of magnitude faster than other regions. Often dubbed dynamical heterogeneity, this behavior has fascinated researchers involved in the study of glass transition for over two decades. A fundamentally important question in all glass transition studies is whether one can connect the growing relaxation time to a concomitantly growing length scale. In this paper, we go beyond the realm of ordinary glass forming liquids and study the origin of a growing dynamical length scale ξ in a self-propelled "active" glass former. This length scale, which is constructed using structural correlations, agrees well with the average size of the clusters of slow-moving particles that are formed as the liquid becomes spatially heterogeneous. We further report that the concomitantly growing α-relaxation time exhibits a simple scaling law, τ_{α}∼exp(μξ/T_{eff}), with μ as an effective chemical potential, T_{eff} as the effective temperature, and μξ as the growing free energy barrier for cluster rearrangements. The findings of our study are valid over four decades of persistence times, and hence they could be very useful in understanding the slow dynamics of a generic active liquid such as an active colloidal suspension, or a self-propelled granular medium.

Stone-Wales defects preserve hyperuniformity in amorphous two-dimensional networks

Disordered hyperuniformity (DHU) is a recently discovered novel state of many-body systems that possesses vanishing normalized infinite-wavelength density fluctuations similar to a perfect crystal and an amorphous structure like a liquid or glass. Here, we discover a hyperuniformity-preserving topological transformation in two-dimensional (2D) network structures that involves continuous introduction of Stone-Wales (SW) defects. Specifically, the static structure factor [Formula: see text] of the resulting defected networks possesses the scaling [Formula: see text] for small wave number k, where [Formula: see text] monotonically decreases as the SW defect concentration p increases, reaches [Formula: see text] at [Formula: see text], and remains almost flat beyond this p. Our findings have important implications for amorphous 2D materials since the SW defects are well known to capture the salient feature of disorder in these materials. Verified by recently synthesized single-layer amorphous graphene, our network models reveal unique electronic transport mechanisms and mechanical behaviors associated with distinct classes of disorder in 2D materials.

Dynamical Crystallites of Active Chiral Particles

One of the intrinsic characteristics of far-from-equilibrium systems is the nonrelaxational nature of the system dynamics, which leads to novel properties that cannot be understood and described by conventional pathways based on thermodynamic potentials. Of particular interest are the formation and evolution of ordered patterns composed of active particles that exhibit collective behavior. Here we examine such a type of nonpotential active system, focusing on effects of coupling and competition between chiral particle self-propulsion and self-spinning. It leads to the transition between three bulk dynamical regimes dominated by collective translative motion, spinning-induced structural arrest, and dynamical frustration. In addition, a persistently dynamical state of self-rotating crystallites is identified as a result of a localized-delocalized transition induced by the crystal-melt interface. The mechanism for the breaking of localized bulk states can also be utilized to achieve self-shearing or self-flow of active crystalline layers.

Optimized Large Hyperuniform Binary Colloidal Suspensions in Two Dimensions

The creation of disordered hyperuniform materials with extraordinary optical properties (e.g., large complete photonic band gaps) requires a capacity to synthesize large samples that are effectively hyperuniform down to the nanoscale. Motivated by this challenge, we propose a feasible equilibrium fabrication protocol using binary paramagnetic colloidal particles confined in a 2D plane. The strong and long-ranged dipolar interaction induced by a tunable magnetic field is free from screening effects that attenuate long-ranged electrostatic interactions in charged colloidal systems. Specifically, we numerically find a family of optimal size ratios that makes the two-phase system effectively hyperuniform. We show that hyperuniformity is a general consequence of low isothermal compressibilities, which makes our protocol suitable to treat more general systems with other long-ranged interactions, dimensionalities, and/or polydispersity. Our methodology paves the way to synthesize large photonic hyperuniform materials that function in the visible to infrared range and hence may accelerate the discovery of novel photonic materials.

Traveling bands, clouds, and vortices of chiral active matter

We consider stochastic dynamics of self-propelled particles with nonlocal normalized alignment interactions subject to phase lag. The role of the lag is to indirectly generate chirality into particle motion. To understand large-scale behavior, we derive a continuum description of an active Brownian particle flow with macroscopic scaling in the form of a partial differential equation for a one-particle probability density function. Due to indirect chirality, we find a spatially homogeneous nonstationary analytic solution for this class of equations. Our development of kinetic and hydrodynamic theories towards such a solution reveals the existence of a wide variety of spatially nonhomogeneous patterns reminiscent of traveling bands, clouds, and vortical structures of linear active matter. Our model may thereby serve as the basis for understanding the nature of chiral active media and designing multiagent swarms with designated behavior.

Transport of self-propelled particles across a porous medium: trapping, clogging, and the Matthew effect

We study the transport of self-propelled particles from one free chamber to another across two stripe-like areas of dense porous medium. The medium is mimicked by arrays of obstacles. We find that active motion could greatly speed up the transport of particles. However, more and more particles become trapped in the obstacle arrays with the enhancement of activity. At high persistence (low rotational diffusion rate) and moderate particle concentration, we observe the Matthew effect in the aggregation of particles in the two obstacle arrays. This effect is weakened by introduction of randomness or deformability into the obstacle arrays. Moreover, the dependence on deformability shows the characteristics of first-order phase transition. In rare situations, the system could be stuck in a dynamic unstable state, e.g. the particles alternatively gather more in one of the two obstacle arrays, exhibiting oscillation of particle number between the arrays. Our results reveal new features in the transport of active objects in a complex medium and have implications for manipulating their collective assembly.

Realizable hyperuniform and nonhyperuniform particle configurations with targeted spectral functions via effective pair interactions

The capacity to identify realizable many-body configurations associated with targeted functional forms for the pair correlation function g_{2}(r) or its corresponding structure factor S(k) is of great fundamental and practical importance. While there are obvious necessary conditions that a prescribed structure factor at number density ρ must satisfy to be configurationally realizable, sufficient conditions are generally not known due to the infinite degeneracy of configurations with different higher-order correlation functions. A major aim of this paper is to expand our theoretical knowledge of the class of pair correlation functions or structure factors that are realizable by classical disordered ensembles of particle configurations, including exotic "hyperuniform" varieties. We first introduce a theoretical formalism that provides a means to draw classical particle configurations from canonical ensembles with certain pairwise-additive potentials that could correspond to targeted analytical functional forms for the structure factor. This formulation enables us to devise an improved algorithm to construct systematically canonical-ensemble particle configurations with such targeted pair statistics, whenever realizable. As a proof of concept, we test the algorithm by targeting several different structure factors across dimensions that are known to be realizable and one hyperuniform target that is known to be nontrivially unrealizable. Our algorithm succeeds for all realizable targets and appropriately fails for the unrealizable target, demonstrating the accuracy and power of the method to numerically investigate the realizability problem. Subsequently, we also target several families of structure-factor functions that meet the known necessary realizability conditions but are not known to be realizable by disordered hyperuniform point configurations, including d-dimensional Gaussian structure factors, d-dimensional generalizations of the two-dimensional one-component plasma (OCP), and the d-dimensional Fourier duals of the previous OCP cases. Moreover, we also explore unusual nonhyperuniform targets, including "hyposurficial" and "antihyperuniform" examples. In all of these instances, the targeted structure factors are achieved with high accuracy, suggesting that they are indeed realizable by equilibrium configurations with pairwise interactions at positive temperatures. Remarkably, we also show that the structure factor of nonequilibrium perfect glass, specified by two-, three-, and four-body interactions, can also be realized by equilibrium pair interactions at positive temperatures. Our findings lead us to the conjecture that any realizable structure factor corresponding to either a translationally invariant equilibrium or nonequilibrium system can be attained by an equilibrium ensemble involving only effective pair interactions. Our investigation not only broadens our knowledge of analytical functional forms for g_{2}(r) and S(k) associated with disordered point configurations across dimensions but also deepens our understanding of many-body physics. Moreover, our work can be applied to the design of materials with desirable physical properties that can be tuned by their pair statistics.

Motility-Induced Temperature Difference in Coexisting Phases

Unlike in thermodynamic equilibrium where coexisting phases always have the same temperature, here we show that systems comprising "active" self-propelled particles can self-organize into two coexisting phases at different kinetic temperatures, which are separated from each other by a sharp and persistent temperature gradient. Contrasting previous studies that have focused on overdamped descriptions of active particles, we show that a "hot-cold coexistence" occurs if and only if accounting for inertia, which is significant, e.g., in activated dusty plasmas, microflyers, whirling fruits, or beetles at interfaces. Our results exemplify a route to use active particles to create a self-sustained temperature gradient across coexisting phases. This phenomenon is fundamentally beyond equilibrium physics and is accompanied by a slow coarsening law with an exponent significantly smaller than the universal 1/3 exponent seen in both equilibrium systems and overdamped active Brownian particles.

Hydrodynamics of random-organizing hyperuniform fluids

Disordered hyperuniform structures are locally random while uniform like crystals at large length scales. Recently, an exotic hyperuniform fluid state was found in several nonequilibrium systems, while the underlying physics remains unknown. In this work, we propose a nonequilibrium (driven-dissipative) hard-sphere model and formulate a hydrodynamic theory based on Navier-Stokes equations to uncover the general mechanism of the fluidic hyperuniformity (HU). At a fixed density, this model system undergoes a smooth transition from an absorbing state to an active hyperuniform fluid and then, to the equilibrium fluid by changing the dissipation strength. We study the criticality of the absorbing-phase transition. We find that the origin of fluidic HU can be understood as the damping of a stochastic harmonic oscillator in q space, which indicates that the suppressed long-wavelength density fluctuation in the hyperuniform fluid can exhibit as either acoustic (resonance) mode or diffusive (overdamped) mode. Importantly, our theory reveals that the damping dissipation and active reciprocal interaction (driving) are the two ingredients for fluidic HU. Based on this principle, we further demonstrate how to realize the fluidic HU in an experimentally accessible active spinner system and discuss the possible realization in other systems.

Active Patchy Colloids with Shape-Tunable Dynamics

Controlling the complex dynamics of active colloids-the autonomous locomotion of colloidal particles and their spontaneous assembly-is challenging yet crucial for creating functional, out-of-equilibrium colloidal systems potentially useful for nano- and micromachines. Herein, by introducing the synthesis of active "patchy" colloids of various low-symmetry shapes, we demonstrate that the dynamics of such systems can be precisely tuned. The low-symmetry patchy colloids are made in bulk via a cluster-encapsulation-dewetting method. They carry essential information encoded in their shapes (particle geometry, number, size, and configurations of surface patches, etc.) that programs their locomotive and assembling behaviors. Under AC electric field, we show that the velocity of particle propulsion and the ability to brake and steer can be modulated by having two asymmetrical patches with various bending angles. The assembly of monopatch particles leads to the formation of dynamic and reconfigurable structures such as spinners and "cooperative swimmers" depending on the particle's aspect ratios. A particle with two patches of different sizes allows for "directional bonding", a concept popular in static assemblies but rare in dynamic ones. With the capability to make tunable and complex shapes, we anticipate the discovery of a diverse range of new dynamics and structures when other external stimuli (e.g., magnetic, optical, chemical, etc.) are employed and spark synergy with shapes.

Hyperuniform states generated by a critical friction field

Hyperuniform states are an efficient way to fill up space for disordered systems. In these states the particle distribution is disordered at the short scale but becomes increasingly uniform when looked at large scales. Hyperuniformity appears in several systems, in static or quasistatic regimes, as well as close to transitions to absorbing states. Here, we show that a vibrated granular layer, at the critical point of the liquid-to-solid transition, displays dynamic hyperuniformity. Prior to the transition, patches of the solid phase form, with length scales and mean lifetimes that diverge critically at the transition point. When reducing the wave number, density fluctuations encounter increasingly more patches that block their propagation, resulting in a static structure factor that tends to zero for small wave numbers at the critical point, which is a signature of hyperuniformity. A simple model demonstrates that this coupling of a density field to a highly fluctuating scalar friction field gives rise to dynamic hyperuniform states. Finally, we show that the structure factor detects better the emergence of hyperuniformity, compared to the particle number variance.

Can a Large Packing be Assembled from Smaller Ones?

We consider zero temperature packings of soft spheres that undergo a jamming to unjamming transition as a function of packing fraction. We compare differences in the structure, as measured from the contact statistics, of a finite subsystem of a large packing to a whole packing with periodic boundaries of an equivalent size and pressure. We find that the fluctuations of the ensemble of whole packings are smaller than those of the ensemble of subsystems. Convergence of these two quantities appears to occur at very large systems, which are usually not attainable in numerical simulations. Finding differences between packings in two dimensions and three dimensions, we also consider four dimensions and mean-field models, and find that they show similar system size dependence. Mean-field critical exponents appear to be consistent with the 3D and 4D packings, suggesting they are above the upper critical dimension. We also find that the convergence as a function of system size to the thermodynamic limit is characterized by two different length scales. We argue that this is the result of the system being above the upper critical dimension.

Simultaneous phase separation and pattern formation in chiral active mixtures

Chiral active particles, or self-propelled circle swimmers, from sperm cells to asymmetric Janus colloids, form a rich set of patterns, which are different from those seen in linear swimmers. Such patterns have mainly been explored for identical circle swimmers, while real-world circle swimmers typically possess a frequency distribution. Here we show that even the simplest mixture of (velocity-aligning) circle swimmers with two different frequencies hosts a complex world of superstructures: The most remarkable example comprises a microflock pattern, formed in one species, while the other species phase separates and forms a macrocluster, coexisting with a gas phase. Here one species microphase separates and selects a characteristic length scale, whereas the other one macrophase separates and selects a density. A second notable example, here occurring in an isotropic system, are patterns comprising two different characteristic length scales, which are controllable via frequency and swimming speed of the individual particles.

Assembly structures and dynamics of active colloidal cells

Many types of active matter are deformable, such as epithelial cells and bacteria. To mimic the feature of deformability, we built a model called an active colloidal cell (ACC), i.e. a vesicle enclosed with self-propelled particles (SPPs), which as a whole can move actively. Based on the model, we then study the role of deformability in the assembly structures and dynamics of ACCs by Langevin dynamics simulation. We find that deformability weakens the self-trapping effect and hence suppresses the clustering and phase separation of the deformable soft ACCs (sACCs). Instead of forming a large compact cluster like ordinary SPPs, sACCs pack into a loose network or porous structure in the phase-separation region. The condensed phase is liquid-like, in which sACCs are strongly compressed and deformed but still keep high motility. The interface between the gas and the condensed phases is blurry and unstable, and the effective interfacial energy is very low. Our work gives new insights into the role of deformability in the assembly of active matter and also provides a reference for further studies on different types of deformable active matter.

Active particles in noninertial frames: How to self-propel on a carousel

Typically the motion of self-propelled active particles is described in a quiescent environment establishing an inertial frame of reference. Here we assume that friction, self-propulsion, and fluctuations occur relative to a noninertial frame and thereby the active Brownian motion model is generalized to noninertial frames. First, analytical solutions are presented for the overdamped case, both for linear swimmers and for circle swimmers. For a frame rotating with constant angular velocity ("carousel"), the resulting noise-free trajectories in the static laboratory frame are trochoids if these are circles in the rotating frame. For systems governed by inertia, such as vibrated granulates or active complex plasmas, centrifugal and Coriolis forces become relevant. For both linear and circling self-propulsion, these forces lead to out-spiraling trajectories which for long times approach a spira mirabilis. This implies that a self-propelled particle will typically leave a rotating carousel. A navigation strategy is proposed to avoid this expulsion, by adjusting the self-propulsion direction at will. For a particle, initially quiescent in the rotating frame, it is shown that this strategy only works if the initial distance to the rotation center is smaller than a critical radius R_{c} which scales with the self-propulsion velocity. Possible experiments to verify the theoretical predictions are discussed.

Rectification of chiral active particles driven by transversal temperature difference

Rectification of chiral active particles driven by transversal temperature difference is investigated in a two-dimensional periodic channel. Chiral active particles can be rectified by transversal temperature difference. Transport behaviors are qualitatively different for different wall boundary conditions. For the sliding boundary condition, the direction of transport completely depends on the chirality of particles. The average velocity is a peaked function of angular velocity or temperature difference. The average velocity increases linearly with the self-propulsion speed, while it decreases monotonically with the increase in the packing fraction. For randomized boundary condition, the transport behaviors become complex. When self-propulsion speed is small, in contrast with the sliding boundary condition, particles move in the opposite direction. However, for large self-propulsion speed, current reversals can occur by continuously changing the system parameters (angular velocity, temperature difference, packing fraction, and width of the channel).