Inertial effects of self-propelled particles: From active Brownian to active Langevin motion

Direct measurement of self-diffusiophoretic force generated by active colloids of different patch coverage using optical tweezers

HYPOTHESIS

Synthetic micro/nanomotors are gaining extensive attention for various biomedical applications (especially in drug delivery) due to their ability to mimic the motion of biological micro/nanoscale swimmers. The feasibility of these applications relies on tight control of propulsion speed, direction, and type of motion (translation, circular, etc.) along with the exerted self-propulsive force. We propose to exploit the variation of both self-propulsion speed and force of active colloids with different patch coverages (with and without supporting layer) for engineering diffusiophoretic micro/nanomotors.

EXPERIMENTS

The microswimmers were designed at various patch coverages (10°, 30°, and 90°) with (Ti/Pt) and without (Pt) an adhesion layer for the catalytic patch through glancing angle metal deposition (GLAD) technique. Mean-square displacement (MSD) analysis was performed to obtain the self-propulsion parameters like speed and angular speed. Using optical tweezers (OT), the self-propulsive force was measured from the force power spectral density.

FINDINGS

The findings of our experiments suggest the non-requirement of any adhesion layer preceding the catalyst deposition since the Pt 10° colloidal batch had the maximal self-propulsion speed (4.61±0.3μm/s) and force (345±57fN) for 5% w/v H2O2 fuel concentration. Moreover, the self-propulsion speed and force decreased with increasing patch size, contrary to theoretical estimates. Also, the self-propulsive force obtained from MSD is 2 to 4 times lower in magnitude than the OT based force values. We believe that the self-propelling motion of the micromotors is possibly hindered due to interactions with the surface of the quartz cuvette during the optical microscopic analysis. Further, the MSD is limited to the self-propulsive motion in two dimensions. On the other hand, OT based force measurement involve trapping the particles in the bulk of the solution entirely avoiding the particle-substrate interactions. Hence, OT based force measurements are better than the propulsion velocity based stokes drag force estimates. We believe that this study can lay the foundation in designing efficient micro/nanomotors for translational biomedical applications.

Particle-wall alignment interaction and active Brownian diffusion through narrow channels

We numerically examine the impacts of particle-wall alignment interactions on active species diffusion through a structureless narrow two-dimensional channel. We consider particle-wall interaction to depend on the self-propulsion velocity direction whereby some specific particle's alignments with respect to the boundary walls are stabilized more. Further, the alignment interaction is meaningful as long as particles are close to the confining boundaries. Unbiased diffusion of active particles for various possible stable velocity alignments against the walls has been examined. We show that for the most stable configuration leading to the self-propulsion velocity direction perpendicular to the wall, diffusivity becomes inversely proportional to the square of the alignment interaction torque. On the other hand, when the self-propulsion velocity direction making an acute angle to the channel walls is the most stable configuration, diffusion exponentially grows with strengthening alignment interaction. Hence, particle-wall interaction plays a pivotal role in the transport control of active particles through narrow channels. Moreover, the impacts of the alignment interactions on diffusion largely depend on the particle's self-propulsion properties and its chirality. Our simulation results can potentially be used to understand unbiased diffusion of artificial or living micro/nano-objects (such as virus, bacteria, Janus particles, etc.) though narrow confined structures.

Active Brownian information engine: Self-propulsion induced colossal performance

The information engine is a feedback mechanism that extorts work from a single heat bath using the mutual information earned during the measurement. We consider an overdamped active Ornstein-Uhlenbeck particle trapped in a 1D harmonic oscillator. The particle experiences fluctuations from an inherent thermal bath with a diffusion coefficient (D) and an active reservoir, with characteristic correlation time (τa) and strength (Da). We design a feedback-driven active Brownian information engine (ABIE) and analyze its best performance criteria. The optimal functioning criteria, the information gained during measurement, and the excess output work are reliant on the dispersion of the steady-state distribution of the particle's position. The extent of enhanced performance of such ABIE depends on the relative values of two underlying time scales of the process, namely, thermal relaxation time (τr) and the characteristic correlation time (τa). In the limit of τa/τr → 0, one can achieve the upper bound on colossal work extraction as ∼0.202γ(D+Da) (γ is the friction coefficient). The excess amount of extracted work reduces and converges to its passive counterpart (∼0.202γD) in the limit of τa/τr → high. Interestingly, when τa/τr = 1, half the upper bound of excess work is achieved irrespective of the strength of either reservoirs, thermal or active. Finally, we look into the average displacement of active Brownian particles in each feedback cycle, which surpasses its thermal analog due to the broader marginal probability distribution.

Optimal Control of Underdamped Systems: An Analytic Approach

Optimal control theory deals with finding protocols to steer a system between assigned initial and final states, such that a trajectory-dependent cost function is minimized. The application of optimal control to stochastic systems is an open and challenging research frontier, with a spectrum of applications ranging from stochastic thermodynamics to biophysics and data science. Among these, the design of nanoscale electronic components motivates the study of underdamped dynamics, leading to practical and conceptual difficulties. In this work, we develop analytic techniques to determine protocols steering finite time transitions at a minimum thermodynamic cost for stochastic underdamped dynamics. As cost functions, we consider two paradigmatic thermodynamic indicators. The first is the Kullback-Leibler divergence between the probability measure of the controlled process and that of a reference process. The corresponding optimization problem is the underdamped version of the Schrödinger diffusion problem that has been widely studied in the overdamped regime. The second is the mean entropy production during the transition, corresponding to the second law of modern stochastic thermodynamics. For transitions between Gaussian states, we show that optimal protocols satisfy a Lyapunov equation, a central tool in stability analysis of dynamical systems. For transitions between states described by general Maxwell-Boltzmann distributions, we introduce an infinite-dimensional version of the Poincaré-Lindstedt multiscale perturbation theory around the overdamped limit. This technique fundamentally improves the standard multiscale expansion. Indeed, it enables the explicit computation of momentum cumulants, whose variation in time is a distinctive trait of underdamped dynamics and is directly accessible to experimental observation. Our results allow us to numerically study cost asymmetries in expansion and compression processes and make predictions for inertial corrections to optimal protocols in the Landauer erasure problem at the nanoscale.

Density and inertia effects on two-dimensional active semiflexible filament suspensions

We examine the influence of density on the transition between chain and spiral structures in planar assemblies of active semiflexible filaments, utilizing detailed numerical simulations. We focus on how increased density, and higher Péclet numbers, affect the activity-induced transition spiral state in a semiflexible, self-avoiding active chain. Our findings show that increasing the density causes the spiral state to break up, reverting to a motile chain-like shape. This results in a density-dependent reentrant phase transition from spirals back to open chains. We attribute this phenomenon to an inertial effect observed at the single polymer level, where increased persistence length due to inertia has been shown in recent three-dimensional studies to cause polymers to open up. Our two-dimensional simulations further reveal that a reduction in the damping coefficient leads to partial unwinding of the spirals, forming longer arms. In suspension, interactions among these extended arms can trigger a complete unwinding of the spirals, driven by the combined effects of density and inertia.

Tunable colloidal spinners: Active chirality and hydrodynamic interactions governed by rotating external electric fields

The rotational dynamics of microparticles in liquids have a wide range of applications, including chemical microreactors, biotechnologies, microfluidic devices, tunable heat and mass transfer, and fundamental understanding of chiral active soft matter which refers to systems composed of particles that exhibit a handedness in their rotation, breaking mirror symmetry at the microscopic level. Here, we report on the study of two effects in colloids in rotating electric fields: (i) the rotation of individual colloidal particles in rotating electric field and related to that (ii) precession of pairs of particles. We show that the mechanism responsible for the rotation of individual particles is related to the time lag between the external field applied to the particle and the particle polarization. Using numerical simulations and experiments with silica particles in a water-based solvent, we prove that the observed rotation of particle pairs and triplets is governed by the tunable rotation of individual particles and can be explained and described by the action of hydrodynamic forces. Our findings demonstrate that colloidal suspensions in rotating electric fields, under some conditions, represent a novel class of chiral soft active matter-tunable colloidal spinners. The experiments and the corresponding theoretical framework we developed open novel prospects for future studies of these systems and for their potential applications.

Langevin dynamics simulations for the critical adsorption of end-grafted active polymers

The critical adsorption of end-grafted active polymer chains on an attractive surface is studied using Langevin dynamics simulations. The active polymers are composed of an active Langevin particle located at the head and a sequential passive chain. Results show that the active force exerted by the active head pulls the active polymer away from the surface. Consequently, the adsorption of the active polymer is hindered, and the critical surface attraction strength, , increases proportionally to the square of the active force, Fa2. The increase in depends on the rotation behavior of the active head. Specifically, for the restricted rotating active polymer (RRAP) chain with a longer rotational persistence time as the rotation of the active head is restricted, increases significantly with Fa. On the other hand, for the freely rotating active polymer (FRAP) chain with a shorter rotational persistence time as the rotation of the active head is free, shows a weak dependence on Fa. The results show that the active force has a significantly stronger pulling effect on the RRAP chain than on the FRAP chain. Furthermore, knotted conformations are observed for the adsorbed RRAP chain at large Fa. These knots reduce the adsorption of monomers near the grafted end. In contrast, no knotted conformations are observed for the FRAP chains due to the comparatively weaker pulling effect of the active force.

Self-aligning active agents with inertia and active torque

We extend the study of the inertial effects on the two-dimensional dynamics of active agents to the case where self-alignment is present. In contrast with the most common models of active particles, we find that self-alignment, which couples the rotational dynamics to the translational one, produces unexpected and nontrivial dynamics, already at the deterministic level. Examining first the motion of a free particle, we contrast the role of inertia depending on the sign of the self-aligning torque. When positive, inertia does not alter the steady-state linear motion of an a-chiral self-propelled particle. On the contrary, for a negative self-aligning torque, inertia leads to the destabilization of the linear motion into a spontaneously broken chiral symmetry orbiting dynamics. Adding an active torque, or bias, to the angular dynamics, the bifurcation becomes imperfect in favor of the chiral orientation selected by the bias. In the case of a positive self-alignment, the interplay of the active torque and inertia leads to the emergence, out of a saddle-node bifurcation, of solutions which coexist with the simply biased linear motion. In the context of a free particle, the rotational inertia leaves unchanged the families of steady-state solutions but sets their stability properties. The situation is radically different when considering the case of a collision with a wall, where a very singular oscillating dynamics takes place which can only be captured if both translational and rotational inertia are present.

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.

Discovering dynamic laws from observations: The case of self-propelled, interacting colloids

Active matter spans a wide range of time and length scales, from groups of cells and synthetic self-propelled colloids to schools of fish and flocks of birds. The theoretical framework describing these systems has shown tremendous success in finding universal phenomenology. However, further progress is often burdened by the difficulty of determining forces controlling the dynamics of individual elements within each system. Accessing this local information is pivotal for the understanding of the physics governing an ensemble of active particles and for the creation of numerical models capable of explaining the observed collective phenomena. In this work, we present ActiveNet, a machine-learning tool consisting of a graph neural network that uses the collective motion of particles to learn active and two-body forces controlling their individual dynamics. We verify our approach using numerical simulations of active Brownian particles, active particles undergoing underdamped Langevin dynamics, and chiral active Brownian particles considering different interaction potentials and values of activity. Interestingly, ActiveNet can equally learn conservative or nonconservative forces as well as torques. Moreover, ActiveNet has proven to be a useful tool to learn the stochastic contribution to the forces, enabling the estimation of the diffusion coefficients. Therefore, all coefficients of the equation of motion of Active Brownian Particles are captured. Finally, we apply ActiveNet to experiments of electrophoretic Janus particles, extracting the active and two-body forces controlling colloids' dynamics. On the one side, we have learned that the active force depends on the electric field and area fraction. On the other side, we have also discovered a dependence of the two-body interaction with the electric field that leads us to propose that the dominant force between active colloids is a screened electrostatic interaction with a constant length scale. We believe that the proposed methodological tool, ActiveNet, might open a new avenue for the study and modeling of experimental suspensions of active particles.

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.

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

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

Inertia suppresses signatures of activity of active Brownian particles in a harmonic potential

A harmonically trapped active Brownian particle exhibits two types of positional distributions-one has a single peak and the other has a single well-that signify steady-state dynamics with low and high activity, respectively. Adding inertia to the translational motion preserves this strict classification of either single-peak or single-well densities but shifts the dividing boundary between the states in the parameter space. We characterize this shift for the dynamics in one spatial dimension using the static Fokker-Planck equation for the full joint distribution of the state space. We derive local results analytically with a perturbation method for a small rotational velocity and then extend them globally with a numerical approach.

Force renormalization for probes immersed in an active bath

Langevin equations or generalized Langevin equations (GLEs) are popular models for describing the motion of a particle in a fluid medium in an effective manner. Here we examine particles immersed in an inherently nonequilibrium fluid, i.e., an active bath, which are subject to an external force. Specifically, we consider two types of forces that are highly relevant for microrheological studies: A harmonic, trapping force and a constant, "drag" force. We study such systems by molecular simulations and use the simulation data to extract an effective GLE description. We find that within this description, in an active bath, the external force in the GLE is not equal to the physical external force, but rather a renormalized external force, which can be significantly smaller. The effect cannot be attributed to the mere temperature renormalization, which is also observed.

When an active bath behaves as an equilibrium one

Active scalar baths consisting of active Brownian particles are characterized by a non-Gaussian velocity distribution, a kinetic temperature, and a diffusion coefficient that scale with the square of the active velocity v_{0}. While these results hold in overdamped active systems, inertial effects lead to normal velocity distributions, with kinetic temperature and diffusion coefficient increasing as ∼v_{0}^{α} with 1<α<2. Remarkably, the late-time diffusivity and mobility decrease with mass. Moreover, we show that the equilibrium Einstein relation is asymptotically recovered with inertia. In summary, the inertial mass restores an equilibriumlike behavior.

Long-range ordering of velocity-aligned active polymers

In this work, we study the effect of covalent bonding on the behavior of non-equilibrium systems with the active force acting on particles along their velocity. Self-ordering of single particles does not occur in this model. However, starting from some critical polymerization degree, the ordered state is observed. It is homogeneous and exhibits no phase separation. In the ordered state, the chains prefer a near-two-dimensional configuration and all move in one direction. Importantly, the self-ordering is obtained only at intermediate active force magnitudes. At high magnitudes, the transition from the disordered to ordered state is suppressed by the swelling of the chains during the transition, as we show by the transition kinetics analysis. We demonstrate the bistable behavior of the system in a particular range of polymerization degrees, amplitudes of active force, densities, and thermostat temperatures. Overall, we show that covalent bonding greatly aids the self-ordering in this active particle model, in contrast to active Brownian particles.

Adsorption of active polymers on attractive nanoparticles

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

Structure of the active Fokker-Planck equation

This paper solves in one and two dimensions the steady noninteractive active Fokker-Planck (FP) equation and finds that its velocity distribution admits, under limiting cases, a dual behavior. Briefly, when the inertial relaxation time is smaller than the orientation time, the active FP equation admits a bimodal shape, whereas the inverse condition is seen to admit a Gaussian one. Once the velocity distribution functions are available, they are used to find their effect on the system's transport properties, such as its mean-square speed. In the process, a useful mathematical identity for the first kind Bessel function as a sum of bimodal exponential functions is spotted.

Pair-distribution function of active Brownian spheres in three spatial dimensions: simulation results and analytical representation

The pair-distribution function, which provides information about correlations in a system of interacting particles, is one of the key objects of theoretical soft matter physics. In particular, it allows for microscopic insights into the phase behavior of active particles. While this function is by now well studied for two-dimensional active matter systems, the more complex and more realistic case of three-dimensional systems is not well understood by now. In this work, we analyze the full pair-distribution function of spherical active Brownian particles interacting via a Weeks-Chandler-Andersen potential in three spatial dimensions using Brownian dynamics simulations. Besides extracting the structure of the pair-distribution function from the simulations, we obtain an analytical representation for this function, parametrized by activity and concentration, which takes into account the symmetries of a homogeneous stationary state. Our results are useful as input to quantitative models of active Brownian particles and advance our understanding of the microstructure in dense active fluids.

Selecting active matter according to motility in an acoustofluidic setup: self-propelled particles and sperm cells

Active systems - including sperm cells, living organisms like bacteria, fish, birds, or active soft matter systems like synthetic "microswimmers" - are characterized by motility, i.e., the ability to propel using their own "engine". Motility is the key feature that distinguishes active systems from passive or externally driven systems. In a large ensemble, motility of individual species can vary in a wide range. Selecting active species according to their motility represents an exciting and challenging problem. We propose a new method for selecting active species based on their motility using an acoustofluidic setup where highly motile species escape from the acoustic trap. This is demonstrated in simulations and in experiments with self-propelled Janus particles and human sperm. The immediate application of this method is selecting highly motile sperm for medically assisted reproduction (MAR). Due to the tunable acoustic trap, the proposed method is more flexible than the existing passive microfluidic methods. The proposed selection method based on motility can also be applied to other active systems that require selecting highly motile species or removing immotile species.

Survival probabilities and first-passage distributions of self-propelled particles in spherical cavities

A model of self-propelled motion in a closed compartment containing simple or complex fluids is formulated in this paper in terms of the dynamics of a point particle moving in a spherical cavity under the action of random thermal forces and exponentially correlated noise. The particle's time evolution is governed by a generalized Langevin equation (GLE) in which the memory function, connected to the thermal forces by a fluctuation-dissipation relation, is described by Jeffrey's model of viscoelasticity (which reduces to a model of ordinary viscous dynamics in a suitable limit). The GLE is transformed exactly to a Fokker-Planck equation that in spherical polar coordinates is in turn found to admit of an exact solution for the particle's probability density function under absorbing boundary conditions at the surface of the sphere. The solution is used to derive an expression (that is also exact) for the survival probability of the particle in the sphere, starting from its center, which is then used to calculate the distribution of the particle's first-passage times to the boundary. The behavior of these quantities is investigated as a function of the Péclet number and the persistence time of the athermal forces, providing insight into the effects of nonequilibrium fluctuations on confined particle motion in three dimensions.

Free and enclosed inertial active gas

In this work, the free expansion of an inertial active gas in three dimensions made of spherical non-interactive active Brownian particles with both translational and rotational inertia (IABPs) is studied. After elucidating the active particles' orientational correlation in three dimensions by employing a Fokker-Planck formalism, the diffusion, mean-square speed, persistence length, reorientation time, Swim and Reynolds pressures and total pressure of this system, are obtained theoretically and corroborated by performing Langevin dynamics simulations. Afterwards, a numerical study on particles' distribution and the mechanical pressure exerted by the active gas enclosed in a cubic box and its dependence on inertia is also carried out. This experiment highlights two important observations: first, as inertia in the system grows while fixing activity, a more uniform particle distribution within the box is achieved. In other words, the classical accumulation of active particles at the walls is seen to be suppressed by inertia. Second, an active gas with translational and rotational inertiae and made of spherical particles still has a state equation which is offered here. This is supported by the fact that both the mechanical pressure definition and the bulk pressure definition as the trace of the swim and Reynolds stress tensors, coincide in the thermodynamic limit.

Effects of inertia on conformation and dynamics of tangentially driven active filaments

Active filamentlike systems propelling along their backbone exist across scales ranging from motor-driven biofilaments to worms and robotic chains. In macroscopic active filaments such as a chain of robots, in contrast to their microscopic counterparts, inertial effects on their motion cannot be ignored. Nonetheless, the consequences of the interplay between inertia and flexibility on the shape and dynamics of active filaments remain unexplored. Here we examine inertial effects on a flexible tangentially driven active polymer model pertinent to the above examples and we determine the conditions under which inertia becomes important. Performing Langevin dynamics simulations of active polymers with underdamped and overdamped dynamics for a wide range of contour lengths and activities, we uncover striking inertial effects on conformation and dynamics for high levels of activities. Inertial collisions increase the persistence length of active polymers and remarkably alter their scaling behavior. In stark contrast to passive polymers, inertia leaves its fingerprint at long times by an enhanced diffusion of the center of mass. We rationalize inertia-induced enhanced dynamics by analytical calculations of center-of-mass velocity correlations, applicable to any active polymer model, which reveal significant contributions from active force fluctuations convoluted by inertial relaxation.

Growth and aging in a few phase-separating active matter systems

Via computer simulations we study evolution dynamics in systems of continuously moving active Brownian particles. The obtained results are discussed against those from the passive 2D Ising case. Following sudden quenches of random configurations to state points lying within the miscibility gaps and to the critical points, we investigate the far-from-steady-state dynamics by calculating quantities associated with structure and characteristic length scales. We also study aging for quenches into the miscibility gap and provide a quantitative picture for the scaling behavior of the two-time order-parameter correlation function. The overall structure and dynamics are consistent with expectations from the Ising model. This remains true for certain active lattice models as well, for which we present results for quenches to the critical points.

Inverted Sedimentation of Active Particles in Unbiased ac Fields

Gaining control over the motion of active particles is crucial for applications ranging from targeted cargo delivery to nanomedicine. While much progress has been made recently to control active motion based on external forces, flows, or gradients in concentration or light intensity, which all have a well-defined direction or bias, little is known about how to steer active particles in situations where no permanent bias can be realized. Here, we show that ac fields with a vanishing time average provide an alternative route to steering active particles. We exemplify this route for inertial active particles in a gravitational field, observing that a substantial fraction of them persistently travels in the upward direction upon switching on the ac field, resulting in an inverted sedimentation profile at the top wall of a confining container. Our results offer a generic control principle that could be used in the future to steer active motion, direct collective behaviors, and purify mixtures.

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.

Light, Matter, Action: Shining Light on Active Matter

Light carries energy and momentum. It can therefore alter the motion of objects on the atomic to astronomical scales. Being widely available, readily controllable, and broadly biocompatible, light is also an ideal tool to propel microscopic particles, drive them out of thermodynamic equilibrium, and make them active. Thus, light-driven particles have become a recent focus of research in the field of soft active matter. In this Perspective, we discuss recent advances in the control of soft active matter with light, which has mainly been achieved using light intensity. We also highlight some first attempts to utilize light's additional properties, such as its wavelength, polarization, and momentum. We then argue that fully exploiting light with all of its properties will play a critical role in increasing the level of control over the actuation of active matter as well as the flow of light itself through it. This enabling step will advance the design of soft active matter systems, their functionalities, and their transfer toward technological applications.

3-D rotation tracking from 2-D images of spherical colloids with textured surfaces

Tracking the three-dimensional rotation of colloidal particles is essential to elucidate many open questions, e.g. concerning the contact interactions between particles under flow, or the way in which obstacles and neighboring particles affect self-propulsion in active suspensions. In order to achieve rotational tracking, optically anisotropic particles are required. We synthesise here rough spherical colloids that present randomly distributed fluorescent asperities and track their motion under different experimental conditions. Specifically, we propose a new algorithm based on a 3-D rotation registration, which enables us to track the 3-D rotation of our rough colloids at short time-scales, using time series of 2-D images acquired at high frame rates with a conventional wide-field microscope. The method is based on the image correlation between a reference image and rotated 3-D prospective images to identify the most likely angular displacements between frames. We first validate our approach against simulated data and then apply it to the cases of: particles flowing through a capillary, freely diffusing at solid-liquid and liquid-liquid interfaces, and self-propelling above a substrate. By demonstrating the applicability of our algorithm and sharing the code, we hope to encourage further investigations in the rotational dynamics of colloidal systems.

Inertial and geometrical effects of self-propelled elliptical Brownian particles

Active particles that self-propel by transforming energy into mechanical motion represent a growing area of research in mathematics, physics, and chemistry. Here we investigate the dynamics of nonspherical inertial active particles moving in a harmonic potential, introducing geometric parameters which take into account the role of eccentricity for nonspherical particles. A comparison between the overdamped and underdamped models for elliptical particles is performed. The model of overdamped active Brownian motion has been used to describe most of the basic aspects of micrometer-sized particles moving in a liquid ("microswimmers"). We consider active particles by extending the active Brownian motion model to incorporate translation and rotation inertia and account for the role of eccentricity. We show how the overdamped and the underdamped models behave in the same way for small values of activity (Brownian case) if eccentricity is equal to zero, but increasing eccentricity leads the two dynamics to substantially depart from each other-in particular, the action of a torque induced by external forces, induced a marked difference close to the walls of the domain if eccentricity is high. Effects induced by inertia include an inertial delay time of the self-propulsion direction from the particle velocity, and the differences between the overdamped and underdamped systems are particularly evident in the first and second moments of the particle velocities. Comparison with the experimental results of vibrated granular particles shows good agreement and corroborates the notion that self-propelling massive particles moving in gaseous media are dominated by inertial effects.

Inertial active ratchet: Simulation versus theory

We present the inertial active dynamics of an Ornstein-Uhlenbeck particle in a piecewise sawtooth ratchet potential. Using the Langevin simulation and matrix continued fraction method (MCFM), the particle transport, steady-state diffusion, and coherence in transport are investigated in different parameter regimes of the model. Spatial asymmetry is found to be a key criterion for the possibility of directed transport in the ratchet. The MCFM results for net particle current of overdamped dynamics of the particle agree well with the simulation results. The simulated particle trajectories for the inertial dynamics and the corresponding position and velocity distribution functions reveal that the system passes through an activity-induced transition in the transport from the running phase to the locked phase of the dynamics. This is further corroborated by the mean square displacement (MSD) calculations, where the MSD gets suppressed with increase in the persistent duration of activity or self-propulsion in the medium and finally approaches zero for a very large value of self propulsion time. The nonmonotonic behavior of the particle current and Péclet number with self-propulsion time confirms that the particle transport and its coherence can be enhanced or reduced by fine tuning the persistent duration of activity. Moreover, for intermediate ranges of self-propulsion time as well as mass of the particle, even though the particle current shows a pronounced unusual maximum with mass, there is no enhancement in the Péclet number, instead the Péclet number decreases with mass, confirming the degradation of coherence in transport.

Dynamics of active particles with translational and rotational inertia

Inertial effects affecting both the translational and rotational dynamics are inherent to a broad range of active systems at the macroscopic scale. Thus, there is a pivotal need for proper models in the framework of active matter to correctly reproduce experimental results, hopefully achieving theoretical insights. For this purpose, we propose an inertial version of the active Ornstein-Uhlenbeck particle (AOUP) model accounting for particle mass (translational inertia) as well as its moment of inertia (rotational inertia) and derive the full expression for its steady-state properties. The inertial AOUP dynamics introduced in this paper is designed to capture the basic features of the well-established inertial active Brownian particle model, i.e. the persistence time of the active motion and the long-time diffusion coefficient. For a small or moderate rotational inertia, these two models predict similar dynamics at all timescales and, in general, our inertial AOUP model consistently yields the same trend upon changing the moment of inertia for various dynamical correlation functions.

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.

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

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

From a microscopic inertial active matter model to the Schrödinger equation

Active field theories, such as the paradigmatic model known as 'active model B+', are simple yet very powerful tools for describing phenomena such as motility-induced phase separation. No comparable theory has been derived yet for the underdamped case. In this work, we introduce active model I+, an extension of active model B+ to particles with inertia. The governing equations of active model I+ are systematically derived from the microscopic Langevin equations. We show that, for underdamped active particles, thermodynamic and mechanical definitions of the velocity field no longer coincide and that the density-dependent swimming speed plays the role of an effective viscosity. Moreover, active model I+ contains an analog of the Schrödinger equation in Madelung form as a limiting case, allowing one to find analoga of the quantum-mechanical tunnel effect and of fuzzy dark matter in active fluids. We investigate the active tunnel effect analytically and via numerical continuation.

Tuning nonequilibrium phase transitions with inertia

In striking contrast to equilibrium systems, inertia can profoundly alter the structure of active systems. Here, we demonstrate that driven systems can exhibit effective equilibrium-like states with increasing particle inertia, despite rigorously violating the fluctuation-dissipation theorem. Increasing inertia progressively eliminates motility-induced phase separation and restores equilibrium crystallization for active Brownian spheres. This effect appears to be general for a wide class of active systems, including those driven by deterministic time-dependent external fields, whose nonequilibrium patterns ultimately disappear with increasing inertia. The path to this effective equilibrium limit can be complex, with finite inertia sometimes acting to accentuate nonequilibrium transitions. The restoration of near equilibrium statistics can be understood through the conversion of active momentum sources to passive-like stresses. Unlike truly equilibrium systems, the effective temperature is now density dependent, the only remnant of the nonequilibrium dynamics. This density-dependent temperature can in principle introduce departures from equilibrium expectations, particularly in response to strong gradients. Our results provide additional insight into the effective temperature ansatz while revealing a mechanism to tune nonequilibrium phase transitions.

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.

Attractor-driven matter

The state of a classical point-particle system may often be specified by giving the position and momentum for each constituent particle. For non-pointlike particles, the center-of-mass position may be augmented by an additional coordinate that specifies the internal state of each particle. The internal state space is typically topologically simple, in the sense that the particle's internal coordinate belongs to a suitable symmetry group. In this paper, we explore the idea of giving internal complexity to the particles, by attributing to each particle an internal state space that is represented by a point on a strange (or otherwise) attracting set. It is, of course, very well known that strange attractors arise in a variety of nonlinear dynamical systems. However, rather than considering strange attractors as emerging from complex dynamics, we may employ strange attractors to drive such dynamics. In particular, by using an attractor (strange or otherwise) to model each particle's internal state space, we present a class of matter coined "attractor-driven matter." We outline the general formalism for attractor-driven matter and explore several specific examples, some of which are reminiscent of active matter. Beyond the examples studied in this paper, our formalism for attractor-driven dynamics may be applicable more broadly, to model complex dynamical and emergent behaviors in a variety of contexts.

Persistent motion of a Brownian particle subject to repulsive feedback with time delay

Based on analytical and numerical calculations we study the dynamics of an overdamped colloidal particle moving in two dimensions under time-delayed, nonlinear feedback control. Specifically, the particle is subject to a force derived from a repulsive Gaussian potential depending on the difference between its instantaneous position, r(t), and its earlier position r(t-τ), where τ is the delay time. Considering first the deterministic case, we provide analytical results for both the case of small displacements and the dynamics at long times. In particular, at appropriate values of the feedback parameters, the particle approaches a steady state with a constant, nonzero velocity whose direction is constant as well. In the presence of noise, the direction of motion becomes randomized at long times, but the (numerically obtained) velocity autocorrelation still reveals some persistence of motion. Moreover, the mean-squared displacement (MSD) reveals a mixed regime at intermediate times with contributions of both ballistic motion and diffusive translational motion, allowing us to extract an estimate for the effective propulsion velocity in presence of noise. We then analyze the data in terms of exact, known results for the MSD of active Brownian particles. The comparison indeed indicates a strong similarity between the dynamics of the particle under repulsive delayed feedback and active motion. This relation carries over to the behavior of the long-time diffusion coefficient D_{eff} which, similarly to active motion, is strongly enhanced compared to the free case. Finally, we show that, for small delays, D_{eff} can be estimated analytically.

Two-temperature activity induces liquid-crystal phases inaccessible in equilibrium

In equilibrium hard-rod fluids, and in effective hard-rod descriptions of anisotropic soft-particle systems, the transition from the isotropic (I) phase to the nematic phase (N) is observed above the rod aspect ratio L/D=3.70 as predicted by Onsager. We examine the fate of this criterion in a molecular dynamics study of a system of soft repulsive spherocylinders rendered active by coupling half the particles to a heat bath at a higher temperature than that imposed on the other half. We show that the system phase-separates and self-organizes into various liquid-crystalline phases that are not observed in equilibrium for the respective aspect ratios. In particular, we find a nematic phase for L/D=3 and a smectic phase for L/D=2 above a critical activity.

Rotation and separation of chiral active particles in a ring-shaped channel

Transport of chiral active particles is numerically investigated in a two-dimensional ring-shaped channel. The ring-shaped channel is transversal asymmetric and can induce the directed transport (rotation) of chiral active particles. For the particles with small chirality, they slide along the outer boundary of the channel. For the particles with large chirality, the particles move along some small local circular orbits and can also exhibit directed rotation. Moreover, the rotation effect can be strongly enhanced by modifying the inner boundary geometry. Based on the study of particle rotation, we further study the separation of active particles with different chiralities. It is found that the particles with different chiralities may be distributed in different regions of the ring-shaped channel. Interestingly, these particles can be completely separated by shifting the channel's inner boundary or adding a blocking plate in the channel. Our results may be useful for understanding relevant experimental phenomena and provide a scheme for the separation of binary mixtures.

Bouncing dynamics of inertial self-propelled particles reveals directional asymmetry

This study aims to examine experimental conditions in which active particles are forced by their surroundings to move forward and backward in a continuous oscillatory manner. The experimental design is based on using a vibrating self-propelled toyrobot called hexbug, which is placed inside a narrow channel closed on one end by a rigid moving wall. Using the end-wall velocity as a controlling factor, the main forward mode of the hexbug movement can be turned to mostly rearward mode. We investigate the bouncing hexbug motion on both experimental and theoretical grounds. The Brownian model of active particles with inertia is employed in the theoretical framework. The model itself uses a pulsed Langevin equation in order to simulate abrupt changes in velocity that mimic hexbug propulsion in the moments when its legs make contact with the base plate. Significant directional asymmetry is caused by the legs bending backward. We demonstrate that the simulation successfully reproduces the experimental characteristics of hexbug motion after regressing the spatial and temporal statistical characteristics, especially when directional asymmetry is under consideration.

Accurate analytical calculation of the rate coefficient for the diffusion-controlled reactions due to hyperbolic diffusion

Using an approach based on the diffusion analog of the Cattaneo-Vernotte differential model, we find the exact analytical solution to the corresponding time-dependent linear hyperbolic initial boundary value problem, describing irreversible diffusion-controlled reactions under Smoluchowski's boundary condition on a spherical sink. By means of this solution, we extend exact analytical calculations for the time-dependent classical Smoluchowski rate coefficient to the case that includes the so-called inertial effects, occurring in the host media with finite relaxation times. We also present a brief survey of Smoluchowski's theory and its various subsequent refinements, including works devoted to the description of the short-time behavior of Brownian particles. In this paper, we managed to show that a known Rice's formula, commonly recognized earlier as an exact reaction rate coefficient for the case of hyperbolic diffusion, turned out to be only its approximation being a uniform upper bound of the exact value. Here, the obtained formula seems to be of great significance for bridging a known gap between an analytically estimated rate coefficient on the one hand and molecular dynamics simulations together with experimentally observed results for the short times regime on the other hand. A particular emphasis has been placed on the rigorous mathematical treatment and important properties of the relevant initial boundary value problems in parabolic and hyperbolic diffusion theories.

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.

Mode-coupling theory for the dynamics of dense underdamped active Brownian particle system

We present a theory to study the inertial effect on glassy dynamics of the underdamped active Brownian particle (UABP) system. Using the assumption of the nonequilibrium steady-state, we obtain an effective Fokker-Planck equation for the probability distribution function (PDF) as a function of positions and momentums. With this equation, we achieve the evolution equation of the intermediate scattering function through the Zwanzig-Mori projection operator method and the mode-coupling theory (MCT). Theoretical analysis shows that the inertia of the particle affects the memory function and corresponding glass transition by influencing the structure factor and a velocity correlation function. The theory provides theoretical support and guidance for subsequent simulation work.

Entropy production of active particles formulated for underdamped dynamics

The present work investigates the effect of inertia on the entropy production rate Π for all canonical models of active particles for different dimensions and the type of confinement. To calculate Π, the link between the entropy production and dissipation of heat rate is explored, resulting in a simple and intuitive expression. By analyzing the Kramers equation, alternative formulations of Π are obtained and the virial theorem for active particles is derived. Exact results are obtained for particles in an unconfined environment and in a harmonic trap. In both cases, Π is independent of temperature. For the case of a harmonic trap, Π attains a maximal value for τ=ω^{-1}, where τ is the persistence time and ω is the natural frequency of an oscillator. For active particles in one-dimensional box, or other nonharmonic potentials, thermal fluctuations are found to reduce Π.

Collective behavior of soft self-propelled disks with rotational inertia

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

Diffusion Coefficient of a Brownian Particle in Equilibrium and Nonequilibrium: Einstein Model and Beyond

The diffusion of small particles is omnipresent in many processes occurring in nature. As such, it is widely studied and exerted in almost all branches of sciences. It constitutes such a broad and often rather complex subject of exploration that we opt here to narrow our survey to the case of the diffusion coefficient for a Brownian particle that can be modeled in the framework of Langevin dynamics. Our main focus centers on the temperature dependence of the diffusion coefficient for several fundamental models of diverse physical systems. Starting out with diffusion in equilibrium for which the Einstein theory holds, we consider a number of physical situations outside of free Brownian motion and end by surveying nonequilibrium diffusion for a time-periodically driven Brownian particle dwelling randomly in a periodic potential. For this latter situation the diffusion coefficient exhibits an intriguingly non-monotonic dependence on temperature.

Active glassy dynamics is unaffected by the microscopic details of self-propulsion

Recent years have seen a rapid increase of interest in dense active materials, which, in the disordered state, share striking similarities with the conventional passive glass-forming matter. For such passive glassy materials, it is well established (at least in three dimensions) that the details of the microscopic dynamics, e.g., Newtonian or Brownian, do not influence the long-time glassy behavior. Here, we investigate whether this still holds true in the non-equilibrium active case by considering two simple and widely used active particle models, i.e., active Ornstein-Uhlenbeck particles (AOUPs) and active Brownian particles (ABPs). In particular, we seek to gain more insight into the role of the self-propulsion mechanism on the glassy dynamics by deriving a mode-coupling theory (MCT) for thermal AOUPs, which can be directly compared to a recently developed MCT for ABPs. Both theories explicitly take into account the active degrees of freedom. We solve the AOUP- and ABP-MCT equations in two dimensions and demonstrate that both models give almost identical results for the intermediate scattering function over a large variety of control parameters (packing fractions, active speeds, and persistence times). We also confirm this theoretical equivalence between the different self-propulsion mechanisms numerically via simulations of a polydisperse mixture of active quasi-hard spheres, thereby establishing that, at least for these model systems, the microscopic details of self-propulsion do not alter the active glassy behavior.

Circular motion subject to external alignment under active driving: Nonlinear dynamics and the circle map

Hardly any real self-propelling or actively driven object is perfect. Thus, undisturbed motion will generally not follow straight lines but rather bent or circular trajectories. We here address self-propelled or actively driven objects that move in discrete steps and additionally tend to migrate towards a certain direction by discrete angular adjustment. Overreaction in the angular alignment is possible. This competition implies pronounced nonlinear dynamics including period doubling and chaotic behavior in a broad parameter regime. Such behavior directly affects the appearance of the trajectories. Furthermore, we address collective motion and effects of spatial self-concentration.

Active matter in space

In the last 20 years, active matter has been a highly dynamic field of research, bridging fundamental aspects of non-equilibrium thermodynamics with applications to biology, robotics, and nano-medicine. Active matter systems are composed of units that can harvest and harness energy and information from their environment to generate complex collective behaviours and forms of self-organisation. On Earth, gravity-driven phenomena (such as sedimentation and convection) often dominate or conceal the emergence of these dynamics, especially for soft active matter systems where typical interactions are of the order of the thermal energy. In this review, we explore the ongoing and future efforts to study active matter in space, where low-gravity and microgravity conditions can lift some of these limitations. We envision that these studies will help unify our understanding of active matter systems and, more generally, of far-from-equilibrium physics both on Earth and in space. Furthermore, they will also provide guidance on how to use, process and manufacture active materials for space exploration and colonisation.