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

Active transport of cargo-carrying and interconnected chiral particles

Directed motion up a concentration gradient is crucial for the survival and maintenance of numerous biological systems, such as sperms moving towards an egg during fertilization or ciliates moving towards a food source. In these systems, chirality-manifested as a rotational torque-plays a vital role in facilitating directed motion. While systematic studies of active molecules in activity gradients exist, the effect of chirality remains little studied. In this study, we examine the simplest case of a chiral active particle connected to a passive particle in a spatially varying activity field. We demonstrate that this minimal setup can exhibit rich emergent tactic behaviors, with the chiral torque serving as the tuning parameter. Notably, when the chiral torque is sufficiently large, even a small passive particle enables the system to display the desired accumulation behavior. Our results further show that in the dilute limit, this desired accumulation behavior persists despite the presence of excluded volume effects. Additionally, interconnected chiral active particles exhibit emergent chemotaxis beyond a critical chain length, with trimers and longer chains exhibiting strong accumulation at sufficiently high chiral torques. This study provides valuable insights into the design principles of hybrid bio-molecular devices of the future.

Spontaneous generation of angular momentum in chiral active crystals

We study a two-dimensional chiral active crystal composed of underdamped chiral active particles. These particles, characterized by intrinsic handedness and persistence, interact via linear forces derived from harmonic potentials. Chirality plays a pivotal role in shaping the system's behavior: it reduces displacement and velocity fluctuations while inducing cross-spatial correlations among different Cartesian components of velocity. These features distinguish chiral crystals from their non-chiral counterparts, leading to the emergence of net angular momentum, as predicted analytically. This angular momentum, driven by the torque generated by the chiral active force, exhibits a non-monotonic dependence on the degree of chirality. Additionally, it contributes to the entropy production rate, as revealed through a path-integral analysis. We investigate the dynamic properties of the crystal in both Fourier and real space. Chirality induces a non-dispersive peak in the displacement spectrum, which underlies the generation of angular momentum and oscillations in time-dependent autocorrelation functions or mean-square displacement, all of which are analytically predicted.

Fluidization and anomalous density fluctuations in 2D Voronoi cell tissues with pulsating activity

Cells not only can be motile by crawling but are also capable of nonmotility active motions like periodic contraction or pulsation. In this work, based on a Voronoi cell model, we show how this nonmotility activity affects the structure, dynamic, and density fluctuations of cellular monolayers. Our model shows that random cell pulsation fluidizes solid epithelial tissues into a hyperuniform fluid state, while pulsation synchronization inhibits the fluidity and causes a reverse solidification. Our results indicate this solidification is a Berezinskii-Kosterlitz-Thouless-type transition, characterized by strong density/dynamic heterogeneity arising from the annihilation of topological defects in the pulsating phase space. The magnitude and length scale of density heterogeneity diverge with the pulsating period, resulting in an opposite giant density fluctuation or anti-hyperuniformity. We propose a fluctuating hydrodynamic theory that can unify the two opposite anomalous fluctuation phenomena. Our findings can help to understand recent experimental observations in Madin-Darby canine kidney monolayers.

Three-dimensional construction of hyperuniform, nonhyperuniform, and antihyperuniform disordered heterogeneous materials and their transport properties via spectral density functions

Rigorous theories connecting physical properties of a heterogeneous material to its microstructure offer a promising avenue to guide the computational material design and optimization. The spectral density function χ[over ̃]_{_{V}}(k), which can be obtained experimentally from scattering data, enables accurate determination of various transport and wave propagation characteristics, including the time-dependent diffusion spreadability S(t) and effective dynamic dielectric constant ε_{e} for electromagnetic wave propagation. Moreover, χ[over ̃]_{_{V}}(k) determines rigorous upper bounds on the fluid permeability K. Given the importance of χ[over ̃]_{_{V}}(k), we present here an efficient Fourier-space based computational framework to construct three-dimensional (3D) statistically isotropic two-phase heterogeneous materials corresponding to targeted spectral density functions. In particular, we employ a variety of analytical functional forms for χ[over ̃]_{_{V}}(k) that satisfy all known necessary conditions to construct disordered stealthy hyperuniform, standard hyperuniform, nonhyperuniform, and antihyperuniform two-phase heterogeneous material systems at varying phase volume fractions. We show that by tuning the correlations in the system across length scales via the targeted functions, one can generate a rich spectrum of distinct structures within each of the above classes of materials. Importantly, we present the first realization of antihyperuniform two-phase heterogeneous materials in 3D, which are characterized by autocovariance function χ_{_{V}}(r) with a power-law tail, resulting in microstructures that contain clusters of dramatically different sizes and morphologies. We also determine the diffusion spreadability S(t) and estimate the fluid permeability K associated with all of the constructed materials directly from the corresponding spectral densities. Although it is well established that the long-time asymptotic scaling behavior of S(t) only depends on the functional form of χ[over ̃]_{_{V}}(k), with the stealthy hyperuniform and antihyperuniform media, respectively, achieving the most and least efficient transport, we find that varying the length-scale parameter within each class of χ[over ̃]_{_{V}}(k) functions can also lead to orders of magnitude variation of S(t) at intermediate and long time scales. Moreover, we find that increasing the solid volume fraction ϕ_{1} and correlation length a in the constructed media generally leads to a decrease in the dimensionless fluid permeability K/a^{2}, while the antihyperuniform media possess the largest K/a^{2} among the four classes of materials with the same ϕ_{1} and a. These results indicate the feasibility of employing parameterized spectral densities for designing composites with targeted transport properties.

The 2025 motile active matter roadmap

Activity and autonomous motion are fundamental aspects of many living and engineering systems. Here, the scale of biological agents covers a wide range, from nanomotors, cytoskeleton, and cells, to insects, fish, birds, and people. Inspired by biological active systems, various types of autonomous synthetic nano- and micromachines have been designed, which provide the basis for multifunctional, highly responsive, intelligent active materials. A major challenge for understanding and designing active matter is their inherent non-equilibrium nature due to persistent energy consumption, which invalidates equilibrium concepts such as free energy, detailed balance, and time-reversal symmetry. Furthermore, interactions in ensembles of active agents are often non-additive and non-reciprocal. An important aspect of biological agents is their ability to sense the environment, process this information, and adjust their motion accordingly. It is an important goal for the engineering of micro-robotic systems to achieve similar functionality. Many fundamental properties of motile active matter are by now reasonably well understood and under control. Thus, the ground is now prepared for the study of physical aspects and mechanisms of motion in complex environments, the behavior of systems with new physical features like chirality, the development of novel micromachines and microbots, the emergent collective behavior and swarming of intelligent self-propelled particles, and particular features of microbial systems. The vast complexity of phenomena and mechanisms involved in the self-organization and dynamics of motile active matter poses major challenges, which can only be addressed by a truly interdisciplinary effort involving scientists from biology, chemistry, ecology, engineering, mathematics, and physics. The 2025 motile active matter roadmap of Journal of Physics: Condensed Matter reviews the current state of the art of the field and provides guidance for further progress in this fascinating research area.

Anomalous suppression of large-scale density fluctuations in classical and quantum spin liquids

Classical spin liquids (CSLs) are intriguing states of matter that do not exhibit long-range magnetic order and are characterized by an extensive ground-state degeneracy. Adding quantum fluctuations, which induce dynamics between these different classical ground states, can give rise to quantum spin liquids (QSLs). QSLs are highly entangled quantum phases of matter characterized by fascinating emergent properties, such as fractionalized excitations and topological order. One such exotic quantum liquid is the [Formula: see text] QSL, which can be regarded as a resonating valence bond (RVB) state formed from superpositions of dimer coverings of an underlying lattice. In this work, we unveil a hidden large-scale structural property of archetypal CSLs and QSLs known as hyperuniformity, i.e., normalized infinite-wavelength density fluctuations are completely suppressed in these systems. In particular, we first demonstrate that classical ensembles of close-packed dimers and their corresponding quantum RVB states are perfectly hyperuniform in general. Subsequently, we focus on a ruby-lattice spin liquid that was recently realized in a Rydberg-atom quantum simulator, and show that the QSL remains effectively hyperuniform even in the presence of a finite density of spinon and vison excitations, as long as the dimer constraint is still largely preserved. Moreover, we demonstrate that metrics based on the framework of hyperuniformity can be used to distinguish the QSL from other proximate quantum phases. These metrics can help identify potential QSL candidates, which can then be further analyzed using more advanced, computationally intensive quantum numerics to confirm their status as true QSLs.

Non-equilibrium dynamic hyperuniform states

Disordered hyperuniform structures are an exotic state of matter having suppressed density fluctuations at large length-scale similar to perfect crystals and quasicrystals but without any long range orientational order. In the past decade, an increasing number of non-equilibrium systems were found to have dynamic hyperuniform states, which have emerged as a new research direction coupling both non-equilibrium physics and hyperuniformity. Here we review the recent progress in understanding dynamic hyperuniform states found in various non-equilibrium systems, including the critical hyperuniformity in absorbing phase transitions, non-equilibrium hyperuniform fluids and the hyperuniform structures in phase separating systems via spinodal decomposition.

Emergent mesoscale correlations in active solids with noisy chiral dynamics

We present the linear response theory for an elastic solid composed of active Brownian particles with intrinsic individual chirality, deriving both a normal mode formulation and a continuum elastic formulation. Using this theory, we compute analytically the velocity correlations and energy spectra under different conditions, showing an excellent agreement with simulations. We generate the corresponding phase diagram, identifying chiral and achiral disordered regimes (for high chirality or noise levels), as well as chiral and achiral states with mesoscopic-range order (for low chirality and noise). The chiral ordered states display mesoscopic spatial correlations and oscillating time correlations, but no wave propagation. In the high chirality regime, we find a peak in the elastic energy spectrum that leads to a non-monotonic behavior with increasing noise strength that is consistent with the emergence of the 'hammering state' recently identified in chiral glasses. Finally, we show numerically that our theory, despite its linear response nature, can be applied beyond the idealized homogeneous solid assumed in our derivations. Indeed, by increasing the level of activity, we show that it remains a good approximation of the system dynamics until just below the melting transition. In addition, we show that there is still an excellent agreement between our analytical results and simulations when we extend our results to heterogeneous solids composed of mixtures of active particles with different intrinsic chirality and noise levels. The derived linear response theory is therefore robust and applicable to a broad range of real-world active systems. Our work provides a thorough analytical and numerical description of the emergent states in a densely packed system of chiral self-propelled Brownian disks, thus allowing a detailed understanding of the phases and dynamics identified in a minimal chiral active system.

Unveiling Hidden Hyperuniformity: Radial Turing Pattern Formation of Marangoni-Driven SiO2 Nanoparticles on Liquid Metal Surface

Mastering the self-organization of nanoparticle morphologies is pivotal in soft matter physics and film growth. Silicon dioxide (SiO2) nanoparticles are an archetypical model of nanomotor in soft matter. Here, the emphasis is on the self-organizing behavior of SiO2 nanoparticles under extreme conditions. It is unveiled that manipulating the states of the metal substrate profoundly dictates the motion characteristics of SiO2 nanoparticles. This manipulation triggers the emergence of intricate morphologies and distinctive patterns. Employing a reaction-diffusion model, the fundamental roles played by Brownian motion and Marangoni-driven motion in shaping fractal structures and radial Turing patterns are demonstrated, respectively. Notably, these radial Turing patterns showcase hyperuniform order, challenging conventional notions of film morphology. These discoveries pave the way for crafting non-equilibrium morphological materials, poised with the potential for self-healing, adaptability, and innovative applications.

Enhancing (quasi-)long-range order in a two-dimensional driven crystal

It has been recently shown that 2D systems can exhibit crystalline phases with long-range translational order showcasing a striking violation of the Hohenberg-Mermin-Wagner (HMW) theorem, which is valid at equilibrium. This is made possible by athermal driving mechanisms that inject energy into the system without exciting long wavelength modes of the density field, thereby inducing hyperuniformity. However, as thermal fluctuations are superimposed on the non-equilibrium driving, long-range translational order is inevitably lost. Here, we discuss the possibility of exploiting non-equilibrium effects to suppress arbitrarily large density fluctuations even when a global thermal bath is coupled to the system. We introduce a model of a harmonic crystal driven both by a global thermal bath and by a momentum conserving noise, where the typical observables related to density fluctuations and long-range translational order can be analytically derived and put in relation. This model allows us to rationalize the violation of the HMW theorem observed in previous studies through the prediction of large-wavelength phonons, which thermalize at a vanishing effective temperature when the global bath is switched off. The conceptual framework introduced through this theory is then applied to numerical simulations of a hard-disk solid in contact with a thermal bath and driven out-of-equilibrium by active collisions. Our numerical analysis demonstrates how varying driving and dissipative parameters can lead to an arbitrary enhancement of the quasi-long-range order in the system regardless of the applied global noise amplitude. Finally, we outline a possible experimental procedure to apply our results to a realistic granular system.

Anomalous relaxation and hyperuniform fluctuations in center-of-mass conserving systems with broken time-reversal symmetry

We study the Oslo model, a paradigm for absorbing-phase transition, on a one-dimensional ring of L sites with a fixed global density ρ[over ¯]; we consider the system strictly above critical density ρ_{c}. Notably, microscopic dynamics conserve both mass and center of mass (CoM), but lack time-reversal symmetry. We show that, despite having highly constrained dynamics due to CoM conservation, the system exhibits diffusive relaxation away from criticality and superdiffusive relaxation near criticality. Furthermore, the CoM conservation severely restricts particle movement, causing the mobility-a transport coefficient analogous to the conductivity for charged particles-to vanish exactly. Indeed, the steady-state temporal growth of current fluctuation is qualitatively different from that observed in diffusive systems with a single conservation law. Remarkably, far from criticality where the relative density Δ=ρ[over ¯]-ρ_{c}≫ρ_{c}, the second cumulant, or the variance, 〈Q_{i}^{2}(T,Δ)〉_{c}, of current Q_{i} across the ith bond up to time T in the steady-state saturates as 〈Q_{i}^{2}〉_{c}≃Σ_{Q}^{2}(Δ)-constT^{-1/2}; near criticality, it grows subdiffusively as 〈Q_{i}^{2}〉_{c}∼T^{α}, with 0<α<1/2, and eventually saturates to Σ_{Q}^{2}(Δ). Interestingly, the asymptotic current fluctuation Σ_{Q}^{2}(Δ) is a nonmonotonic function of Δ: It diverges as Σ_{Q}^{2}(Δ)∼Δ^{2} for Δ≫ρ_{c} and Σ_{Q}^{2}(Δ)∼Δ^{-δ}, with δ>0, for Δ→0^{+}. Using a mass-conservation principle, we exactly determine the exponents δ=2(1-1/ν_{⊥})/ν_{⊥} and α=δ/zν_{⊥} via the correlation-length and dynamic exponents, ν_{⊥} and z, respectively. Finally, we show that in the steady state the self-diffusion coefficient D_{s}(ρ[over ¯]) of tagged particles is connected to activity through the relation D_{s}(ρ[over ¯])=a(ρ[over ¯])/ρ[over ¯].

Continuous symmetry breaking of low-dimensional systems driven by inhomogeneous oscillatory driving forces

The driving forces of chiral active particles and deformations of cells are often modeled by spatially inhomogeneous but temporally periodic driving forces. Such inhomogeneous oscillatory driving forces have only recently been proposed in the context of active matter, and their effects on the systems are not yet fully understood. In this work, we theoretically study the impact of spatially inhomogeneous oscillatory driving forces on continuous symmetry breaking. We first analyze the linear model for the soft modes in the ordered phase to derive the lower critical dimension of the model, and then analyze the spherical model to investigate more detailed phase behaviors. Interestingly, our analysis reveals that symmetry breaking occurs even in one and two dimensions, where the Hohenberg-Mermin-Wagner theorem prohibits continuous symmetry breaking in equilibrium. Furthermore, fluctuations of conserved quantities, such as density, are anomalously suppressed in the long-wavelength, i.e., show hyperuniformity.

Spontaneous separation of attractive chiral mixtures

The separation of chiral matter has garnered significant attention due to its wide-ranging applications in biological and chemical processes. In prior researches, particle interactions were predominantly repulsive, but the indiscriminate attraction among particles under attractive interactions makes the separation of mixtures more difficult. The question of whether chiral mixed particles, characterized by attractive effects, can undergo spontaneous separation, remains unresolved. We study a binary mixture of chiral (counterclockwise or clockwise) active particles with attractive interactions. It is demonstrated that attractive chiral particles can undergo spontaneous separation without the aid of any specific strategies. The key factor driving the separation is the attractive interactions, enabling the formation of stable clusters of particles with same chirality. There exist optimal parameters (self-propelled velocity, angular velocity, and packing fraction) at which the separation is optimal. Our results may contribute to a deeper understanding of the mechanisms behind chiral matter separation and potentially catalyze further experimental investigations in this field.

Universal Hyperuniform Organization in Looped Leaf Vein Networks

The leaf vein network is a hierarchical vascular system that transports water and nutrients to the leaf cells. The thick primary veins form a branched network, while the secondary veins can develop closed loops forming a well-defined cellular structure. Through extensive analysis of a variety of distinct leaf species, we discover that the apparently disordered cellular structures of the secondary vein networks exhibit a universal hyperuniform organization and possess a hidden order on large scales. Disorder hyperuniform systems lack conventional long-range order, yet they completely suppress normalized infinite-wavelength density fluctuations like crystals. Specifically, we find that the distributions of the geometric centers associated with the vein network loops possess a vanishing static structure factor in the limit that the wave number k goes to 0, i.e., S(k)∼k^{α}, where α≈0.64±0.021, providing an example of class III hyperuniformity in biology. This hyperuniform organization leads to superior efficiency of diffusive transport, as evidenced by the much faster convergence of the time-dependent spreadability S(t) to its longtime asymptotic limit, compared to that of other uncorrelated or correlated disordered but nonhyperuniform organizations. Our results also have implications for the discovery and design of novel disordered network materials with optimal transport properties.

Collective Hall current in chiral active fluids: Coupling of phase and mass transport through traveling bands

Active fluids composed of constituents that are constantly driven away from thermal equilibrium can support spontaneous currents and can be engineered to have unconventional transport properties. Here, we report the emergence of (meta)stable traveling bands in computer simulations of aligning circle swimmers. These bands are different from polar flocks and, through coupling phase with mass transport, induce a bulk particle current with a component perpendicular to the propagation direction, thus giving rise to a collective Hall (or Magnus) effect. Traveling bands require sufficiently small orbits and undergo a discontinuous transition into a synchronized state with transient polar clusters for large orbital radii. Within a minimal hydrodynamic theory, we show that the bands can be understood as nondispersive soliton solutions fully accounting for the numerically observed properties.

Nonequilibrium hyperuniform states in active turbulence

We demonstrate that the complex spatiotemporal structure in active fluids can feature characteristics of hyperuniformity. Using a hydrodynamic model, we show that the transition from hyperuniformity to nonhyperuniformity and antihyperuniformity depends on the strength of active forcing and can be related to features of active turbulence without and with scaling characteristics of inertial turbulence. Combined with identified signatures of Levy walks and nonuniversal diffusion in these systems, this allows for a biological interpretation and the speculation of nonequilibrium hyperuniform states in active fluids as optimal states with respect to robustness and strategies of evasion and foraging.

Rotational inertia-induced glassy transition in chiral particle systems

The dense active matter exhibits characteristics reminiscent of traditional glassy phenomena, yet the role of rotational inertia in glass dynamics remains elusive. In this study, we investigate the glass dynamics of chiral active particles influenced by rotational inertia. Rotational inertia endows exponential memory to particle orientation, restricting its alteration and amplifying the effective persistence time. At lower spinning frequencies, the diffusion coefficient exhibits a peak function relative to rotational inertia for shorter persistence times, while it steadily increases with rotational inertia for longer persistence times. In the realm of high-frequency spinning, the impact of rotational inertia on diffusion behavior becomes more pronounced, resulting in a nonmonotonic and intricate relationship between the diffusion coefficient and rotational inertia. Consequently, the introduction of rotational inertia significantly alters the glassy dynamics of chiral active particles, allowing for the control over transitions between fluid and glassy states by modulating rotational inertia. Moreover, our findings indicate that at a specific spinning temperature, there exists an optimal spinning frequency at which the diffusion coefficient attains its maximum value.

Phase separation kinetics and cluster dynamics in two-dimensional active dumbbell systems

Molecular dynamics simulations were employed to investigate the phase separation process of a two-dimensional active Brownian dumbbell model. We evaluated the time dependence of the typical size of the dense component using the scaling properties of the structure factor, along with the averaged number of clusters and their radii of gyration. The growth observed is faster than in active disk models, and this effect is further enhanced under stronger activity. Next, we focused on studying the hexatic order of the clusters. The length associated with the orientational order increases algebraically with time and faster than for spherical active Brownian particles. Under weak active forces, most clusters exhibit a uniform internal orientational order. However, under strong forces, large clusters consist of domains with different orientational orders. We demonstrated that the latter configurations are not stable, and given sufficient time to evolve, they eventually achieve homogeneous configurations as well. No gas bubbles are formed within the clusters, even when there are patches of different hexatic order. Finally, attention was directed towards the geometry and motion of the clusters themselves. By employing a tracking algorithm, we showed that clusters smaller than the typical size at the observation time exhibit regular shapes, while larger ones display fractal characteristics. In between collisions or break-ups, the clusters behave as solid bodies. Their centers of mass undergo circular motion, with radii increasing with the cluster size. The angular velocity of the center of mass equals that of the constituents with respect to their center of mass. These observations were rationalised with a simple mechanical model.

Vibrational properties of disordered stealthy hyperuniform 1D atomic chains

Disorder hyperuniformity is a recently discovered exotic state of many-body systems that possess a hidden order in between that of a perfect crystal and a completely disordered system. Recently, this novel disordered state has been observed in a number of quantum materials including amorphous 2D graphene and silica, which are endowed with unexpected electronic transport properties. Here, we numerically investigate 1D atomic chain models, including perfect crystalline, disordered stealthy hyperuniform (SHU) as well as randomly perturbed atom packing configurations to obtain a quantitative understanding of how the unique SHU disorder affects the vibrational properties of these low-dimensional materials. We find that the disordered SHU chains possess lower cohesive energies compared to the randomly perturbed chains, implying their potential reliability in experiments. Our inverse partition ratio (IPR) calculations indicate that the SHU chains can support fully delocalized states just like perfect crystalline chains over a wide range of frequencies, i.e.ω∈(0,100)cm-1, suggesting superior phonon transport behaviors within these frequencies, which was traditionally considered impossible in disordered systems. Interestingly, we observe the emergence of a group of highly localized states associated withω∼200cm-1, which is characterized by a significant peak in the IPR and a peak in phonon density of states at the corresponding frequency, and is potentially useful for decoupling electron and phonon degrees of freedom. These unique properties of disordered SHU chains have implications in the design and engineering of novel quantum materials for thermal and phononic applications.

Emergent Chirality and Hyperuniformity in an Active Mixture with Nonreciprocal Interactions

We investigate collective dynamics in a binary mixture of programmable robots in experiments and simulations. While robots of the same species align their motion direction, interaction between species is distinctly nonreciprocal: species A aligns with B and species B antialigns with A. This nonreciprocal interaction gives rise to the emergence of collective chiral motion that can be stabilized by limiting the robot angular speed to be below a threshold. Within the chiral phase, increasing the robot density or extending the range of local repulsive interactions can drive the system through an absorbing-active transition. At the transition point, the robots exhibit a remarkable capacity for self-organization, forming disordered hyperuniform states.

Chiral active particles are sensitive reporters to environmental geometry

Chiral active particles (CAPs) are self-propelling particles that break time-reversal symmetry by orbiting or spinning, leading to intriguing behaviors. Here, we examined the dynamics of CAPs moving in 2D lattices of disk obstacles through active Brownian dynamics simulations and granular experiments with grass seeds. We find that the effective diffusivity of the CAPs is sensitive to the structure of the obstacle lattice, a feature absent in achiral active particles. We further studied the transport of CAPs in obstacle arrays under an external field and found a reentrant directional locking effect, which can be used to sort CAPs with different activities. Finally, we demonstrated that parallelogram lattices of obstacles without mirror symmetry can separate clockwise and counter-clockwise CAPs. The mechanisms of the above three novel phenomena are qualitatively explained. As such, our work provides a basis for designing chirality-based tools for single-cell diagnosis and separation, and active particle-based environmental sensors.

How does a hyperuniform fluid freeze?

All phase transitions can be categorized into two different types: continuous and discontinuous phase transitions. Discontinuous phase transitions are normally accompanied with significant structural changes, and nearly all of them have the kinetic pathway of nucleation and growth, if the system does not suffer from glassy dynamics. Here, in a system of barrier-controlled reactive particles, we find that the discontinuous freezing transition of a nonequilibrium hyperuniform fluid into an absorbing state does not have the kinetic pathway of nucleation and growth, and the transition is triggered by long-wavelength fluctuations. The transition rate decreases with increasing the system size, which suggests that the metastable hyperuniform fluid could be kinetically stable in an infinitely large system. This challenges the common understanding of metastability in discontinuous phase transitions. Moreover, we find that the "metastable yet kinetically stable" hyperuniform fluid features a scaling in the structure factor [Formula: see text] in 2D, which is the third dynamic hyperuniform state in addition to the critical hyperuniform state with [Formula: see text] and the nonequilibrium hyperuniform fluid with [Formula: see text].

Computational design of anisotropic stealthy hyperuniform composites with engineered directional scattering properties

Disordered hyperuniform materials are an emerging class of exotic amorphous states of matter that endow them with singular physical properties, including large isotropic photonic band gaps, superior resistance to fracture, and nearly optimal electrical and thermal transport properties, to name but a few. Here we generalize the Fourier-space-based numerical construction procedure for designing and generating digital realizations of isotropic disordered hyperuniform two-phase heterogeneous materials (i.e., composites) developed by Chen and Torquato [Acta Mater. 142, 152 (2018)1359-645410.1016/j.actamat.2017.09.053] to anisotropic microstructures with targeted spectral densities. Our generalized construction procedure explicitly incorporates the vector-dependent spectral density function χ[over ̃]_{_{V}}(k) of arbitrary form that is realizable. We demonstrate the utility of the procedure by generating a wide spectrum of anisotropic stealthy hyperuniform microstructures with χ[over ̃]_{_{V}}(k)=0 for k∈Ω, i.e., complete suppression of scattering in an "exclusion" region Ω around the origin in Fourier space. We show how different exclusion-region shapes with various discrete symmetries, including circular-disk, elliptical-disk, square, rectangular, butterfly-shaped, and lemniscate-shaped regions of varying size, affect the resulting statistically anisotropic microstructures as a function of the phase volume fraction. The latter two cases of Ω lead to directionally hyperuniform composites, which are stealthy hyperuniform only along certain directions and are nonhyperuniform along others. We find that while the circular-disk exclusion regions give rise to isotropic hyperuniform composite microstructures, the directional hyperuniform behaviors imposed by the shape asymmetry (or anisotropy) of certain exclusion regions give rise to distinct anisotropic structures and degree of uniformity in the distribution of the phases on intermediate and large length scales along different directions. Moreover, while the anisotropic exclusion regions impose strong constraints on the global symmetry of the resulting media, they can still possess structures at a local level that are nearly isotropic. Both the isotropic and anisotropic hyperuniform microstructures associated with the elliptical-disk, square, and rectangular Ω possess phase-inversion symmetry over certain range of volume fractions and a percolation threshold ϕ_{c}≈0.5. On the other hand, the directionally hyperuniform microstructures associated with the butterfly-shaped and lemniscate-shaped Ω do not possess phase-inversion symmetry and percolate along certain directions at much lower volume fractions. We also apply our general procedure to construct stealthy nonhyperuniform systems. Our construction algorithm enables one to control the statistical anisotropy of composite microstructures via the shape, size, and symmetries of Ω, which is crucial to engineering directional optical, transport, and mechanical properties of two-phase composite media.

Chiral active matter in external potentials

We investigate the interplay between chirality and confinement induced by the presence of an external potential. For potentials having radial symmetry, the circular character of the trajectories induced by the chiral motion reduces the spatial fluctuations of the particle, thus providing an extra effective confining mechanism, that can be interpreted as a lowering of the effective temperature. In the case of non-radial potentials, for instance, with an elliptic shape, chirality displays a richer scenario. Indeed, the chirality can break the parity symmetry of the potential that is always fulfilled in the non-chiral system. The probability distribution displays a strong non-Maxwell-Boltzmann shape that emerges in cross-correlations between the two Cartesian components of the position, that vanishes in the absence of chirality or when radial symmetry of the potential is restored. These results are obtained by considering two popular models in active matter, i.e. chiral Active Brownian particles and chiral active Ornstein-Uhlenbeck particles.

Exploiting compositional disorder in collectives of light-driven circle walkers

Emergent behavior in collectives of "robotic" units with limited capabilities that is robust and programmable is a promising route to perform tasks on the micro and nanoscale that are otherwise difficult to realize. However, a comprehensive theoretical understanding of the physical principles, in particular steric interactions in crowded environments, is still largely missing. Here, we study simple light-driven walkers propelled through internal vibrations. We demonstrate that their dynamics is well captured by the model of active Brownian particles, albeit with an angular speed that differs between individual units. Transferring to a numerical model, we show that this polydispersity of angular speeds gives rise to specific collective behavior: self-sorting under confinement and enhancement of translational diffusion. Our results show that, while naively perceived as imperfection, disorder of individual properties can provide another route to realize programmable active matter.

Collective behavior of chiral active particles with anisotropic interactions in a confined space

Extensive studies so far have indicated that chirality, anisotropic interactions and spatial confinement play important roles in collective dynamics in active matter systems. However, how the overall interplay of these crucial factors affects the novel phases and macroscopic properties remains less explored. Here, using Langevin dynamics simulations, we investigate the self-organization of a chiral active system composed of amphiphilic Janus particles, where the embedded anisotropic interaction orientation is assumed to be either the same or just opposite to the direction of active force. A wealth of dynamic phases are observed including formation of phase separation, clustering state, homogeneous state, spiral vortex flow, swarm and spatiotemporal oscillation. By tuning self-propelled angular speed and anisotropic interaction strength, we identify the non-equilibrium phase diagrams, and reveal the very non-trivial modulation of both vortex and swarm patterns. Intriguingly, we find that strong chirality-alignment-confinement coupling yields a self-driven spatial and temporal organization periodically oscillating between a counterclockwise vortex and a clockwise one. Our work provides a new understanding of the novel self-assembly arising in such a confined system and enables new strategies for achieving ordered dynamic structures.

Glassy Dynamics in Chiral Fluids

Chiral active matter is enjoying a rapid increase of interest, spurred by the rich variety of asymmetries that can be attained in, e.g., the shape or self-propulsion mechanism of active particles. Though this has already led to the observance of so-called chiral crystals, active chiral glasses remain largely unexplored. A possible reason for this could be the naive expectation that interactions dominate the glassy dynamics and the details of the active motion become increasingly less relevant. Here, we show that quite the opposite is true by studying the glassy dynamics of interacting chiral active Brownian particles. We demonstrate that when our chiral fluid is pushed to glassy conditions, it exhibits highly nontrivial dynamics, especially compared to a standard linear active fluid such as common active Brownian particles. Despite the added complexity, we are still able to present a full rationalization for all identified dynamical regimes. Most notably, we introduce a new "hammering" mechanism, unique to rapidly spinning particles in high-density conditions, that can fluidize a chiral active solid.

Tunable collective dynamics of ellipsoidal Quincke particles

Collective behaviors in active systems become dramatically complicated in the presence of chirality. In this study, we show that ellipsoidal Quincke particles driven by an electric field exhibit flexible and tunable chirality because of the tilting of the spinning axis. As the tilting torque decreases with the increase of angular speed, the motion of individual particles transforms from localized circle motion to global rolling. However, because of the anisotropic shape and the resulting anisotropic polar interactions, it is dynamically easier for ellipsoids to bind and form rotating structures rather than to align their velocities. In dense systems, the suppression of velocity aligning produces transient dense clusters which produce dynamic heterogeneity. The formation and dissociation of dense clusters prohibit the emergence of large-scale collective motions and limit the amplitude of density fluctuations. These findings demonstrate that collective dynamics and thus the scale of density fluctuations in active systems with tunable chirality can be well controlled. This has potential applications in exploring disordered hyperuniform states.

Equilibrium states corresponding to targeted hyperuniform nonequilibrium pair statistics

The Zhang-Torquato conjecture [G. Zhang and S. Torquato, Phys. Rev. E, 2020, 101, 032124.] states that any realizable pair correlation function g₂(r) or structure factor S(k) of a translationally invariant nonequilibrium system can be attained by an equilibrium ensemble involving only (up to) effective two-body interactions. To further test and study this conjecture, we consider two singular nonequilibrium models of recent interest that also have the exotic hyperuniformity property: a 2D "perfect glass" and a 3D critical absorbing-state model. We find that each nonequilibrium target can be achieved accurately by equilibrium states with effective one- and two-body potentials, lending further support to the conjecture. To characterize the structural degeneracy of such a nonequilibrium-equilibrium correspondence, we compute higher-order statistics for both models, as well as those for a hyperuniform 3D uniformly randomized lattice (URL), whose higher-order statistics can be very precisely ascertained. Interestingly, we find that the differences in the higher-order statistics between nonequilibrium and equilibrium systems with matching pair statistics, as measured by the "hole" probability distribution, provide measures of the degree to which a system is out of equilibrium. We show that all three systems studied possess the bounded-hole property and that holes near the maximum hole size in the URL are much rarer than those in the underlying simple cubic lattice. Remarkably, upon quenching, the effective potentials for all three systems possess local energy minima (i.e., inherent structures) with stronger forms of hyperuniformity compared to their target counterparts. Our methods are expected to facilitate the self-assembly of tunable hyperuniform soft-matter systems.

Hydrodynamic Synchronization of Chiral Microswimmers

We study synchronization in bulk suspensions of spherical microswimmers with chiral trajectories using large scale numerics. The model is generic. It corresponds to the lowest order solution of a general model for self-propulsion at low Reynolds numbers, consisting of a nonaxisymmetric rotating source dipole. We show that both purely circular and helical swimmers can spontaneously synchronize their rotation. The synchronized state corresponds to velocity alignment with high orientational order in both the polar and azimuthal directions. Finally, we consider a racemic mixture of helical swimmers where intraspecies synchronization is observed while the system remains as a spatially uniform fluid. Our results demonstrate hydrodynamic synchronization as a natural collective phenomenon for microswimmers with chiral trajectories.

Collective behavior of passive and active circle swimming particle mixtures

We present a numerical study on a binary mixture of passive and circle swimming, self-propelling particles which interact via the Lennard-Jones (LJ) potential in two dimensions. Using Brownian Dynamics (BD) simulations, we present state diagrams using the control parameters such as attraction strength, angular velocity, self-propulsion velocity and composition. In a symmetric mixture, the system undergoes a transition from a mixed gel to a rotating passive cluster state and finally to a homogeneous fluid state as translational activity increases. The formation of the rotating cluster of passive particles surrounded by active and passive monomers is attributed to the combined effect of composition, activity and strength of attraction of the active particles. Different phases are characterized using radial distribution functions, bond order parameters, cluster fraction and probability distribution of local volume fractions. The present study addresses comprehensively the intricate role of activity, angular velocity, inter-particle interaction and compositional variation on the phase behavior. The predictions presented in the study can be experimentally realized in synthetic colloidal swimmers and motile bacterial suspensions.

Directed Autonomous Motion and Chiral Separation of Self-Propelled Janus Particles in Convection Roll Arrays

Self-propelled Janus particles exhibit autonomous motion thanks to engines of their own. However, due to the randomly changing direction of such motion they are of little use for emerging nanotechnological and biomedical applications. Here, we numerically show that the motion of chiral active Janus particles can be directed, subjecting them to a linear array of convection rolls. The rectification power of self-propulsion motion here can be made to be more than 60%, which is much larger than earlier reports. We show that rectification of a chiral Janus particle's motion leads to conspicuous segregation of dextrogyre and levogyre active particles from a racemic binary mixture. Further, we demonstrate how efficiently the rectification effect can be exploited to separate dextrogyre and levogyre particles when their intrinsic torques are distributed with Gaussian statistics.

Resonant Diffusion of a Gravitactic Circle Swimmer

We investigate the dynamics of a single chiral active particle subject to an external torque due to the presence of a gravitational field. Our computer simulations reveal an arbitrarily strong increase of the long-time diffusivity of the gravitactic agent when the external torque approaches the intrinsic angular drift. We provide analytic expressions for the mean-square displacement in terms of eigenfunctions and eigenvalues of the noisy-driven-pendulum problem. The pronounced maximum in the diffusivity is then rationalized by the vanishing of the lowest eigenvalues of the Fokker-Planck equation for the angular motion as the rotational diffusion decreases and the underlying classical bifurcation is approached. A simple harmonic-oscillator picture for the barrier-dominated motion provides a quantitative description for the onset of the resonance while its range of validity is determined by the crossover to a critical-fluctuation-dominated regime.

Emergent synchronization and flocking in purely repulsive self-navigating particles

Inspired by groups of animals and robots, we study the collective dynamics of large numbers of active particles, each one trying to get to its own randomly placed target, while avoiding collisions with each other. The particles we study are repulsive homing active Brownian particles, self-propelled particles whose orientation relaxes at a finite rate towards an absorbing target in two-dimensional continuous space. For a wide range of parameters, these particles form synchronized system-wide chiral flocks, in spite of the absence of explicit alignment interactions. We show that this dramatic behavior obtains for different system sizes and density, that it is robust against the addition of noise, polydispersity, and bounding walls, and that it can exhibit dynamical topological defects. We develop an analogy to an off-lattice, ferromagnetic XY model, which allows us to interpret the different phases, as well as the topological defects.

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.

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.

Hyperuniform Active Chiral Fluids with Tunable Internal Structure

Large density fluctuations observed in active systems and hyperuniformity are two seemingly incompatible phenomena. However, the formation of hyperuniform states has been recently predicted in nonequilibrium fluids formed by chiral particles performing circular motion with the same handedness. Here we report evidence of hyperuniformity realized in a chiral active fluid comprised of pear-shaped Quincke rollers of arbitrary handedness. We show that hyperuniformity and large density fluctuations, triggered by dynamic clustering, coexist in this system at different length scales. The system loses its hyperuniformity as the curvature of particles' motion increases, transforming them into localized spinners. Our results experimentally demonstrate a novel hyperuniform active fluid and provide new insights into an interplay between chirality, activity, and hyperuniformity.

Analytical approach to chiral active systems: suppressed phase separation of interacting Brownian circle swimmers

We consider chirality in active systems by exemplarily studying the phase behavior of planar systems of interacting Brownian circle swimmers with a spherical shape. For this purpose, we derive a predictive field theory that is able to describe the collective dynamics of circle swimmers. The theory yields a mapping between circle swimmers and noncircling active Brownian particles and predicts that the angular propulsion of the particles leads to a suppression of their motility-induced phase separation, being in line with recent simulation results. In addition, the theory provides analytical expressions for the spinodal corresponding to the onset of motility-induced phase separation and the associated critical point as well as for their dependence on the angular propulsion of the circle swimmers. We confirm our findings by Brownian dynamics simulations. The agreement between results from theory and simulations is found to be good.

Migration of an active colloidal cell in inhomogeneous environments

Living cells on a substrate with mechanical inhomogeneities often migrate along or against the mechanical gradient, i.e., mechanotaxis, which inspires us to ask how biomimetic cells without biochemical signaling processes respond to environmental inhomogeneity. Here, we perform computer simulations to study the migration of a 2D active colloidal cell (ACC), which consists of active particles enclosed by a passive vesicle, in a heterogeneous environment composed of two adjoining uniform regions with different attributes (influencing the persistent length of the active particle). We find that the ACC can migrate unidirectionally across the interface separating the heterogeneous region and behave tactically. Interestingly, the tactic motion of the ACC is qualitatively different from that of the constituent active particles themselves. In addition, the ACC may also experience a directed drift along the interface of the heterogeneous environment. The tactic behavior of the ACC can be explained by analyzing the pressure distribution on the cell membrane exerted by the enclosed active particles. The findings provide insights into understanding the taxis of biological cells and designing biomimetic cells with environment-sensitive capabilities.

Ultraefficient reconstruction of effectively hyperuniform disordered biphase materials via non-Gaussian random fields

Disordered hyperuniform systems are statistically isotropic and possess no Bragg peaks like liquids and glasses, yet they suppress large-scale density fluctuations in a similar manner as in perfect crystals. The unique hyperuniform long-range order in these systems endow them with nearly optimal transport, electronic, and mechanical properties. The concept of hyperuniformity was originally introduced for many-particle systems and has subsequently been generalized to biphase heterogeneous materials such as porous media, composites, polymers, and biological tissues for unconventional property discovery. Existing methods for rendering realizations of disordered hyperuniform biphase materials reconstruction typically employ stochastic optimization such as the simulated annealing approach, which requires many iterations. Here, we propose an explicit ultraefficient method for reconstructing effectively hyperuniform biphase materials, based on the second-order non-Gaussian random fields where no additional tuning step or iteration is needed. Both the effectively hyperuniform microstructure and the latent material property field can be simultaneously generated in a single reconstruction. Moreover, our method can also incorporate hierarchical uncertainties in the heterogeneous materials, including both uncertainties in the disordered material microstructure and material property variation within each phase. The efficiency and feasibility of the proposed reconstruction method are demonstrated via a wide spectrum of examples spanning from isotropic to anisotropic, effectively hyperuniform to nonhyperuniform, and antihyperuniform systems. Our ultraefficient reconstruction method can be readily incorporated into material design, probabilistic analysis, optimization, and discovery of novel disordered hyperuniform heterogeneous materials.

Dynamical clustering interrupts motility-induced phase separation in chiral active Brownian particles

One of the most intriguing phenomena in active matter has been the gas-liquid-like motility-induced phase separation (MIPS) observed in repulsive active particles. However, experimentally, no particle can be a perfect sphere, and the asymmetric shape, mass distribution, or catalysis coating can induce an active torque on the particle, which makes it a chiral active particle. Here, using computer simulations and dynamic mean-field theory, we demonstrate that the large enough torque of circle active Brownian particles in two dimensions generates a dynamical clustering state interrupting the conventional MIPS. Multiple clusters arise from the combination of the conventional MIPS cohesion, and the circulating current caused disintegration. The nonvanishing current in non-equilibrium steady states microscopically originates from the motility "relieved" by automatic rotation, which breaks the detailed balance at the continuum level. This suggests that no equilibrium-like phase separation theory can be constructed for chiral active colloids even with tiny active torque, in which no visible collective motion exists. This mechanism also sheds light on the understanding of dynamic clusters observed in a variety of active matter systems.

Learning self-driven collective dynamics with graph networks

Despite decades of theoretical research, the nature of the self-driven collective motion remains indigestible and controversial, while the phase transition process of its dynamic is a major research issue. Recent methods propose to infer the phase transition process from various artificially extracted features using machine learning. In this thesis, we propose a new order parameter by using machine learning to quantify the synchronization degree of the self-driven collective system from the perspective of the number of clusters. Furthermore, we construct a powerful model based on the graph network to determine the long-term evolution of the self-driven collective system from the initial position of the particles, without any manual features. Results show that this method has strong predictive power, and is suitable for various noises. Our method can provide reference for the research of other physical systems with local interactions.

Polyhedral Micromotors of Metal-Organic Frameworks: Symmetry Breaking and Propulsion

Colloidal micromotors can autonomously propel due to their broken symmetry that leads to unbalanced local mechanical forces. Most commonly, micromotors are synthesized to possess a Janus structure or its variants, having two components distinct in shape, composition, or surface joined together on opposite sides. Here, we report on an alternative approach for creating micromotors, where microcrystals of metal-organic frameworks (MOFs) with various polyhedral shapes are propelled under an AC electric field. In these cases, symmetry breaking is realized by orienting the polyhedral particles in a unique direction to generate uneven electrohydrodynamic flow. The particle orientations are controlled by a delicate competition between the electric and gravitational forces exerted on the particle, which we rationalize using experiments and a theoretical model. Furthermore, by leveraging the MOF types and shapes, or surface properties, we show that the propulsion of MOF motors can be tuned or reversed. Because of the flexibility in designing MOFs and their one-step scalable synthesis, our strategy is simple yet versatile for making well-defined functional micromotors.

Geometrical properties of mechanically annealed systems near the jamming transition

Geometrical properties of two-dimensional mixtures near the jamming transition point are numerically investigated using harmonic particles under mechanical training. The configurations generated by the quasi-static compression and oscillatory shear deformations exhibit anomalous suppression of the density fluctuations, known as hyperuniformity, below and above the jamming transition. For the jammed system trained by compression above the transition point, the hyperuniformity exponent increases. For the system below the transition point under oscillatory shear, the hyperuniformity exponent also increases until the shear amplitude reaches the threshold value. The threshold value matches with the transition point from the point-reversible phase where the particles experience no collision to the loop-reversible phase where the particles' displacements are non-affine during a shear cycle before coming back to an original position. The results demonstrated in this paper are explained in terms of neither of universal criticality of the jamming transition nor the nonequilibrium phase transitions.

Self-Organization of Active Droplets into Vortex-like Structures

Natural or artificial active objects can demonstrate mirror asymmetry of collective motion when they are moving coherently in a vortex. The majority of known cases related to the emergence of collective dynamical chirality are referred to as active objects with individual structure chirality and/or dynamical chirality. Here, we demonstrate that dynamically and structurally achiral active droplets can self-organize into vortex-like structures. Octane droplets dispersed in the aqueous solution of an anionic surfactant are activated with ammonia addition. The motion of droplets is due to the Marangoni flow emerging at the interfaces of the droplets. We found out that different modes of vortex motion of droplets in the emulsion can arise depending on the size of the region that confines the motion of the droplets and their number density and velocity.

Emergent vortices and phase separation in systems of chiral active particles with dipolar interactions

Using Brownian dynamics (BD) simulations we investigate the self-organization of a monolayer of chiral active particles with dipolar interactions. Each particle is driven by both, translational and rotational self-propulsion, and carries a permanent point dipole moment at its center. The direction of the translational propulsion for each particle is chosen to be parallel to its dipole moment. Simulations are performed at high dipolar coupling strength and a density below that related to motility-induced phase separation in simple active Brownian particles. Despite this restriction, we observe a wealth of phenomena including formation of two types of vortices, phase separation, and flocking transitions. To understand the appearance and disappearance of vortices in the many-particle system, we further investigate the dynamics of simple ring structures under the impact of self-propulsion.

Dynamic evolution of hyperuniformity in a driven dissipative colloidal system

Hyperuniformity is evolving to become a unifying concept that can help classify and characterize equilibrium and nonequilibrium states of matter. Therefore, understanding the extent of hyperuniformity in dissipative systems is critical. Here, we study the dynamic evolution of hyperuniformity in a driven dissipative colloidal system. We experimentally show and numerically verify that the hyperuniformity of a colloidal crystal is robust against various lattice imperfections and environmental perturbations. This robustness even manifests during crystal disassembly as the system switches between strong (class I), logarithmic (class II), weak (class III), and non-hyperuniform states. To aid analyses, we developed a comprehensive computational toolbox, enabling real-time characterization of hyperuniformity in real- and reciprocal-spaces together with the evolution of several order metric features, and measurements showing the effect of external perturbations on the spatiotemporal distribution of the particles. Our findings provide a new framework to understand the basic principles that drive a dissipative system to a hyperuniform state.

Barrier-controlled nonequilibrium criticality in reactive particle systems

Nonequilibrium critical phenomena generally exist in many dynamic systems, like chemical reactions and some driven-dissipative reactive particle systems. Here, by using computer simulation and theoretical analysis, we demonstrate the crucial role of the activation barrier on the criticality of dynamic phase transitions in a minimal reactive hard-sphere model. We find that at zero thermal noise, with increasing the activation barrier, the type of transition changes from a continuous conserved directed percolation into a discontinuous dynamic transition by crossing a tricritical point. A mean-field theory combined with field simulation is proposed to explain this phenomenon. The possibility of Ising-type criticality in the nonequilibrium system at finite thermal noise is also discussed.