Crystallization of Self-Propelled Hard Discs

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.

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.

Active control of equilibrium, near-equilibrium, and far-from-equilibrium colloidal systems

The development of top-down active control over bottom-up colloidal assembly processes has the potential to produce materials, surfaces, and objects with applications in a wide range of fields spanning from computing to materials science to biomedical engineering. In this review, we summarize recent progress in the field using a taxonomy based on how active control is used to guide assembly. We find there are three distinct scenarios: (1) navigating kinetic pathways to reach a desirable equilibrium state, (2) the creation of a desirable metastable, kinetically trapped, or kinetically arrested state, and (3) the creation of a desirable far-from-equilibrium state through continuous energy input. We review seminal works within this framework, provide a summary of important application areas, and present a brief introduction to the fundamental concepts of control theory that are necessary for the soft materials community to understand this literature. In addition, we outline current and potential future applications of actively-controlled colloidal systems, and we highlight important open questions and future directions.

Polymer-chain configurations in active and passive baths

The configurations taken by polymers embedded in out-of-equilibrium baths may have broad implications in a variety of biological systems. As such, they have attracted considerable interest, particularly in simulation studies. Here we analyze the distribution of configurations taken by a passive flexible chain in a bath of hard, self-propelled, vibrated disks and systematically compare it to that of the same flexible chain in a bath of hard, thermal-like, vibrated disks. We demonstrate experimentally that the mean length and mean radius of gyration of both chains agree with Flory's law. However, the Kuhn length associated with the number of correlated monomers is smaller in the case of the active bath, corresponding to a higher effective temperature. Importantly, the active bath does not just simply map on a hot equilibrium bath. Close examination of the chains' configurations indicates a marked bias, with the chain in the active bath more likely assuming configurations with a single prominent bend.

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.

Disordered Collective Motion in Dense Assemblies of Persistent Particles

We explore the emergence of nonequilibrium collective motion in disordered nonthermal active matter when persistent motion and crowding effects compete, using simulations of a two-dimensional model of size polydisperse self-propelled particles. In stark contrast with monodisperse systems, we find that polydispersity stabilizes a homogeneous active liquid at arbitrary large persistence times, characterized by remarkable velocity correlations and irregular turbulent flows. For all persistence values, the active fluid undergoes a nonequilibrium glass transition at large density. This is accompanied by collective motion, whose nature evolves from near-equilibrium spatially heterogeneous dynamics at small persistence, to a qualitatively different intermittent dynamics when persistence is large. This latter regime involves a complex time evolution of the correlated displacement field.

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.

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.

Particle anisotropy tunes emergent behavior in active colloidal systems

Studies of active particle systems have demonstrated that particle anisotropy can impact the collective behavior of a system, motivating a systematic study. Here, we report a systematic computational investigation of the role of anisotropy in shape and active force director on the collective behavior of a two-dimensional active colloidal system. We find that shape and force anisotropy can combine to produce critical densities both lower and higher than those of disks. We demonstrate that changing particle anisotropy tunes what we define as a "collision efficiency" of inter-particle collisions in leading to motility-induced phase separation (MIPS) of the system. We use this efficiency to determine the relative critical density across systems. Additionally, we observe that local structure in phase-separated clusters is the same as the particle's equilibrium densest packing, suggesting a general connection between equilibrium behavior and non-equilibrium cluster structure of self-propelled anisotropic particles. In engineering applications for active colloidal systems, shape-controlled steric interactions such as those described here may offer a simple route for tailoring emergent behaviors.

Cage Length Controls the Nonmonotonic Dynamics of Active Glassy Matter

Dense active matter is gaining widespread interest due to its remarkable similarity with conventional glass-forming materials. However, active matter is inherently out of equilibrium and even simple models such as active Brownian particles (ABPs) and active Ornstein-Uhlenbeck particles (AOUPs) behave markedly differently from their passive counterparts. Controversially, this difference has been shown to manifest itself via either a speedup, slowdown, or nonmonotonic change of the glassy relaxation dynamics. Here we rationalize these seemingly contrasting views on the departure from equilibrium by identifying the ratio of the short-time length scale to the cage length, i.e., the length scale of local particle caging, as a vital and unifying control parameter for active glassy matter. In particular, we explore the glassy dynamics of both thermal and athermal ABPs and AOUPs upon increasing the persistence time. We find that for all studied systems there is an optimum of the dynamics; this optimum occurs when the cage length coincides with the corresponding short-time length scale of the system, which is either the persistence length for athermal systems or a combination of the persistence length and a diffusive length scale for thermal systems. This new insight, for which we also provide a simple physical argument, allows us to reconcile and explain the manifestly disparate departures from equilibrium reported in many previous studies of dense active materials.

Topographic Control of Order in Quasi-2D Granular Phase Transitions

We experimentally investigate the nature of 2D phase transitions in a quasi-2D granular fluid. Using a surface decorated with periodically spaced dimples we observe interfacial tension between coexisting granular liquid and crystal phases. Measurements of the orientational and translational order parameters and associated susceptibilities indicate that the surface topography alters the order of the phase transition from a two-step continuous one to a first-order liquid-solid one. The interplay of boundary inelasticity and geometry, either order promoting or inhibiting, controls whether it is the granular crystal or the granular fluid which makes contact with the edge. This order induced wetting has important consequences, determining how coexisting phases separate spatially.

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.

Emergent Field-Driven Robot Swarm States

We present an ecology-inspired form of active matter consisting of a robot swarm. Each robot moves over a planar dynamic resource environment represented by a large light-emitting diode array in search of maximum light intensity; the robots deplete (dim) locally by their presence the local light intensity and seek maximum light intensity. Their movement is directed along the steepest local light intensity gradient; we call this emergent symmetry breaking motion "field drive." We show there emerge dynamic and spatial transitions similar to gas, crystalline, liquid, glass, and jammed states as a function of robot density, resource consumption rates, and resource recovery rates. Paradoxically the nongas states emerge from smooth, flat resource landscapes, not rough ones, and each state can directly move to a glassy state if the resource recovery rate is slow enough, at any robot density.

FastTrack: An open-source software for tracking varying numbers of deformable objects

Analyzing the dynamical properties of mobile objects requires to extract trajectories from recordings, which is often done by tracking movies. We compiled a database of two-dimensional movies for very different biological and physical systems spanning a wide range of length scales and developed a general-purpose, optimized, open-source, cross-platform, easy to install and use, self-updating software called FastTrack. It can handle a changing number of deformable objects in a region of interest, and is particularly suitable for animal and cell tracking in two-dimensions. Furthermore, we introduce the probability of incursions as a new measure of a movie’s trackability that doesn’t require the knowledge of ground truth trajectories, since it is resilient to small amounts of errors and can be computed on the basis of an ad hoc tracking. We also leveraged the versatility and speed of FastTrack to implement an iterative algorithm determining a set of nearly-optimized tracking parameters—yet further reducing the amount of human intervention—and demonstrate that FastTrack can be used to explore the space of tracking parameters to optimize the number of swaps for a batch of similar movies. A benchmark shows that FastTrack is orders of magnitude faster than state-of-the-art tracking algorithms, with a comparable tracking accuracy. The source code is available under the GNU GPLv3 at https://github.com/FastTrackOrg/FastTrack and pre-compiled binaries for Windows, Mac and Linux are available at http://www.fasttrack.sh.

Many researchers and engineers face the challenge of tracking objects from very different systems across several fields of research. We observed that despite this diversity the core of the tracking task is very general and can be formalized. We thus introduce the notion of incursions—i.e. to what extent an object can enter a neighbor’s space—which can be defined on a statistical basis and captures the interplay between the acquisition rate, the objects’ dynamics and the geometrical characteristics of the scene, including density. To validate this approach, we compiled a dataset from various fields of Physics, Biology and human activities to serve as a benchmark for general-purpose tracking softwares. This dataset is open and accepts new submissions. We also developped a software called FastTrack that is able to track most of the movies in the dataset by proposing standard image processing tools and state-of-the-art implementation of the matching algorithm, which is at the core of the tracking task. Besides, it is open-source, simple to install and use and has an ergonomic interface to obtain fast and reliable results. FastTrack is particularly convenient for small-scale research projects, typically when the development of a dedicated software is overkill.

Pattern detection in colloidal assembly: A mosaic of analysis techniques

Characterization of the morphology, identification of patterns and quantification of order encountered in colloidal assemblies is essential for several reasons. First of all, it is useful to compare different self-assembly methods and assess the influence of different process parameters on the final colloidal pattern. In addition, casting light on the structures formed by colloidal particles can help to get better insight into colloidal interactions and understand phase transitions. Finally, the growing interest in colloidal assemblies in materials science for practical applications going from optoelectronics to biosensing imposes a thorough characterization of the morphology of colloidal assemblies because of the intimate relationship between morphology and physical properties (e.g. optical and mechanical) of a material. Several image analysis techniques developed to investigate images (acquired via scanning electron microscopy, digital video microscopy and other imaging methods) provide variegated and complementary information on the colloidal structures under scrutiny. However, understanding how to use such image analysis tools to get information on the characteristics of the colloidal assemblies may represent a non-trivial task, because it requires the combination of approaches drawn from diverse disciplines such as image processing, computational geometry and computational topology and their application to a primarily physico-chemical process. Moreover, the lack of a systematic description of such analysis tools makes it difficult to select the ones more suitable for the features of the colloidal assembly under examination. In this review we provide a methodical and extensive description of real-space image analysis tools by explaining their principles and their application to the investigation of two-dimensional colloidal assemblies with different morphological characteristics.

Competing active and passive interactions drive amoebalike crystallites and ordered bands in active colloids

Swimmers and self-propelled particles are physical models for the collective behavior and motility of a wide variety of living systems, such as bacteria colonies, bird flocks, and fish schools. Such artificial active materials are amenable to physical models which reveal the microscopic mechanisms underlying the collective behavior. Here we study colloids in a dc electric field. Our quasi-two-dimensional system of electrically driven particles exhibits a rich and exotic phase behavior exhibiting passive crystallites, motile crystallites, an active gas, and banding. Amongst these are two mesophases, reminiscent of systems with competing interactions. At low field strengths activity suppresses demixing, leading to motile crystallites. Meanwhile, at high field strengths, activity drives partial demixing to traveling bands. We parametrize a particulate simulation model which reproduces the experimentally observed phases.

Phase Separation by Entanglement of Active Polymerlike Worms

We investigate the aggregation and phase separation of thin, living T. tubifex worms that behave as active polymers. Randomly dispersed active worms spontaneously aggregate to form compact, highly entangled blobs, a process similar to polymer phase separation, and for which we observe power-law growth kinetics. We find that the phase separation of active polymerlike worms does not occur through Ostwald ripening, but through active motion and coalescence of the phase domains. Interestingly, the growth mechanism differs from conventional growth by droplet coalescence: the diffusion constant characterizing the random motion of a worm blob is independent of its size, a phenomenon that can be explained from the fact that the active random motion arises from the worms at the surface of the blob. This leads to a fundamentally different phase-separation mechanism that may be unique to active polymers.

Dense active matter model of motion patterns in confluent cell monolayers

Epithelial cell monolayers show remarkable displacement and velocity correlations over distances of ten or more cell sizes that are reminiscent of supercooled liquids and active nematics. We show that many observed features can be described within the framework of dense active matter, and argue that persistent uncoordinated cell motility coupled to the collective elastic modes of the cell sheet is sufficient to produce swirl-like correlations. We obtain this result using both continuum active linear elasticity and a normal modes formalism, and validate analytical predictions with numerical simulations of two agent-based cell models, soft elastic particles and the self-propelled Voronoi model together with in-vitro experiments of confluent corneal epithelial cell sheets. Simulations and normal mode analysis perfectly match when tissue-level reorganisation occurs on times longer than the persistence time of cell motility. Our analytical model quantitatively matches measured velocity correlation functions over more than a decade with a single fitting parameter.

Activity-modulated phase transition in a two-dimensional mixture of active and passive colloids

We study a two-dimensional binary mixture of active and passive colloids as an idealized model of an hybrid aggregate of living cells and inert particles. We perform molecular dynamics simulations of this system using two different thermostats, and we systematically investigate the effect of varying these two effective temperatures on the system behavior, as characterized by its density, structure and thermoelastic properties. Our results indicate that the presence of active colloids shifts the mixture towards the liquid state and renders it more deformable. Such system softening and melting effects due to the addition of active particles are larger than expected from a linear combination of temperatures of the active and passive components. This heightened effect becomes more pronounced as the effective temperature difference between the two components becomes larger. The binary mixture remains homogeneous for moderate colloidal activity, but segregation arises for large effective temperature difference. Our results provide insights to guide future experimental hybrid aggregate studies with promising biomedical applications.

Nonmonotonic behavior in dense assemblies of active colloids

We study experimentally a sediment of self-propelled Brownian particles with densities ranging from dilute to ergodic supercooled to nonergodic glass to nonergodic polycrystal. In a companion paper, we observe a nonmonotonic response to activity of relaxation of the nonergodic glass state: a dramatic slowdown when particles become weakly self-propelled, followed by a speedup at higher activities. Here we map ergodic supercooled states to standard passive glassy physics, provided a monotonic shift of the glass packing fraction and the replacement of the ambient temperature by the effective temperature. However, we show that this mapping fails beyond glass transition. This failure is responsible for the nonmonotonic response. Furthermore, we generalize our finding by examining the dynamical response of another class of nonergodic systems: polycrystals. We observe the same nonmonotonic response to activity. To explain this phenomenon, we measure the size of domains where particles move in the same direction. This size also shows a nonmonotonic response, with small lengths corresponding to slow relaxation. This suggests that the failure of the mapping of nonergodic active states to a passive situation is general and is linked to anisotropic relaxation mechanisms specific to active matter.

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

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

Inertial delay of self-propelled particles

The motion of self-propelled massive particles through a gaseous medium is dominated by inertial effects. Examples include vibrated granulates, activated complex plasmas and flying insects. However, inertia is usually neglected in standard models. Here, we experimentally demonstrate the significance of inertia on macroscopic self-propelled particles. We observe a distinct inertial delay between orientation and velocity of particles, originating from the finite relaxation times in the system. This effect is fully explained by an underdamped generalisation of the Langevin model of active Brownian motion. In stark contrast to passive systems, the inertial delay profoundly influences the long-time dynamics and enables new fundamental strategies for controlling self-propulsion in active matter.

Full Phase Diagram of Active Brownian Disks: From Melting to Motility-Induced Phase Separation

We establish the complete phase diagram of self-propelled hard disks in two spatial dimensions from the analysis of the equation of state and the statistics of local order parameters. The equilibrium melting scenario is maintained at small activities, with coexistence between active liquid and hexatic order, followed by a proper hexatic phase, and a further transition to an active solid. As activity increases, the emergence of hexatic and solid order is shifted towards higher densities. Above a critical activity and for a certain range of packing fractions, the system undergoes motility-induced phase separation and demixes into low and high density phases; the latter can be either disordered (liquid) or ordered (hexatic or solid) depending on the activity.

Active Atoms and Interstitials in Two-Dimensional Colloidal Crystals

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

Spontaneously Flowing Crystal of Self-Propelled Particles

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

Active crystals on a sphere

Two-dimensional crystals on curved manifolds exhibit nontrivial defect structures. Here we consider "active crystals" on a sphere, which are composed of self-propelled colloidal particles. Our work is based on a phase-field-crystal-type model that involves a density and a polarization field on the sphere. Depending on the strength of the self-propulsion, three different types of crystals are found: a static crystal, a self-spinning "vortex-vortex" crystal containing two vortical poles of the local velocity, and a self-translating "source-sink" crystal with a source pole where crystallization occurs and a sink pole where the active crystal melts. These different crystalline states as well as their defects are studied theoretically here and can in principle be confirmed in experiments.

Spontaneous membrane formation and self-encapsulation of active rods in an inhomogeneous motility field

We study the collective dynamics of self-propelled rods in an inhomogeneous motility field. At the interface between two regions of constant but different motility, a smectic rod layer is spontaneously created through aligning interactions between the active rods, reminiscent of an artificial, semipermeable membrane. This "active membrane" engulfes rods which are locally trapped in low-motility regions and thereby further enhances the trapping efficiency by self-organization, an effect which we call "self-encapsulation." Our results are gained by computer simulations of self-propelled rod models confined on a two-dimensional planar or spherical surface with a stepwise constant motility field, but the phenomenon should be observable in any geometry with sufficiently large spatial inhomogeneity. We also discuss possibilities to verify our predictions of active-membrane formation in experiments of self-propelled colloidal rods and vibrated granular matter.

Criticality of the zero-temperature jamming transition probed by self-propelled particles

We perform simulations of athermal systems of self-propelled particles (SPPs) interacting via harmonic repulsion in the vicinity of the zero-temperature jamming transition at point J. Every particle is propelled by a constant force f pointing to a randomly assigned and fixed direction. When f is smaller than the yield force f_y, the system is statically jammed. We find that f_y increases with packing fraction and exhibits finite size scaling, implying the criticality of point J. When f > f_y, SPPs flow forever and their velocities satisfy the k-Gamma distribution. Velocity distributions at various packing fractions and f collapse when the particle velocity is scaled by the average velocity v[combining macron], suggesting that v[combining macron] is a reasonable quantity to characterize the response to f. We thus define a response function R(ϕ,f) = v[combining macron](ϕ,f)/f. The function exhibits critical scaling nicely, implying again the criticality of point J. Our analysis and results indicate that systems of SPPs behave analogically to sheared systems, although their driving mechanisms are apparently distinct.

Cooperative standing-horizontal-standing reentrant transition for numerous solid particles under external vibration

We report the collective behavior of numerous plastic bolt-like particles exhibiting one of two distinct states, either standing stationary or horizontal accompanied by tumbling motion, when placed on a horizontal plate undergoing sinusoidal vertical vibration. Experimentally, we prepared an initial state in which all of the particles were standing except for a single particle that was placed at the center of the plate. Under continuous vertical vibration, the initially horizontal particle triggers neighboring particles to fall over into a horizontal state through tumbling-induced collision, and this effect gradually spreads to all of the particles, i.e., the number of horizontal particles is increased. Interestingly, within a certain range of vibration intensity, almost all of the horizontal particles revert back to standing in association with the formation of apparent 2D hexagonal dense-packing. Thus, phase segregation between high and low densities, or crystalline and disperse domains, of standing particles is generated as a result of the reentrant transition. The essential features of such cooperative dynamics through the reentrant transition are elucidated with a simple kinetic model. We also demonstrate that an excitable wave with the reentrant transition is observed when particles are situated in a quasi-one-dimensional confinement on a vibrating plate.

Driving dynamic colloidal assembly using eccentric self-propelled colloids

Designing protocols to dynamically direct the self-assembly of colloidal particles has become an important direction in soft matter physics because of promising applications in the fabrication of dynamic responsive functional materials. Here, using computer simulations, we found that in the mixture of passive colloids and eccentric self-propelled active particles, when the eccentricity and self-propulsion of active particles are high enough, the eccentric active particles can push passive colloids to form a large dense dynamic cluster, and the system undergoes a novel dynamic demixing transition. Our simulations show that the dynamic demixing occurs when the eccentric active particles move much faster than the passive particles such that the dynamic trajectories of different active particles can overlap each other while passive particles are depleted from the dynamic trajectories of active particles. Our results suggest that this is in analogy to the entropy-driven demixing in colloid-polymer mixtures, in which polymer random coils can overlap with each other while depleting the colloids. More interestingly, we find that by fixing the passive colloid composition at a certain value with increasing density, the system undergoes an intriguing re-entrant mixing, and the demixing only occurs within a certain intermediate density range. This suggests a new way of designing active matter to drive the self-assembly of passive colloids and fabricate dynamic responsive materials.

Noise and diffusion of a vibrated self-propelled granular particle

Granular materials are an important physical realization of active matter. In vibration-fluidized granular matter, both diffusion and self-propulsion derive from the same collisional forcing, unlike many other active systems where there is a clean separation between the origin of single-particle mobility and the coupling to noise. Here we present experimental studies of single-particle motion in a vibrated granular monolayer, along with theoretical analysis that compares grain motion at short and long time scales to the assumptions and predictions, respectively, of the active Brownian particle (ABP) model. Our results show that despite the unique relation between noise and propulsion, a variety of granular particles are correctly described by the ABP model. Additionally, our scheme of analysis for validating the inputs and outputs of the model can be applied to other granular and non-granular active systems.

Brownian motion of a circle swimmer in a harmonic trap

We study the dynamics of a Brownian circle swimmer with a time-dependent self-propulsion velocity in an external temporally varying harmonic potential. For several situations, the noise-free swimming paths, the noise-averaged mean trajectories, and the mean-square displacements are calculated analytically or by computer simulation. Based on our results, we discuss optimal swimming strategies in order to explore a maximum spatial range around the trap center. In particular, we find a resonance situation for the maximum escape distance as a function of the various frequencies in the system. Moreover, the influence of the Brownian noise is analyzed by comparing noise-free trajectories at zero temperature with the corresponding noise-averaged trajectories at finite temperature. The latter reveal various complex self-similar spiral or rosette-like patterns. Our predictions can be tested in experiments on artificial and biological microswimmers under dynamical external confinement.