Emerging activity in bilayered dispersions with wake-mediated interactions

Impact of non-reciprocal interactions on colloidal self-assembly with tunable anisotropy

Non-reciprocal (NR) effective interactions violating Newton's third law occur in many biological systems, but can also be engineered in synthetic, colloidal systems. Recent research has shown that such NR interactions can have tremendous effects on the overall collective behavior and pattern formation, but can also influence aggregation processes on the particle scale. Here, we focus on the impact of non-reciprocity on the self-assembly of a colloidal system (originally passive) with anisotropic interactions whose character is tunable by external fields. In the absence of non-reciprocity, that is, under equilibrium conditions, the colloids form square-like and hexagonal aggregates with extremely long lifetimes yet no large-scale phase separation [Kogler et al., Soft Matter 11, 7356 (2015)], indicating kinetic trapping. Here, we study, based on Brownian dynamics simulations in 2D, an NR version of this model consisting of two species with reciprocal isotropic, but NR anisotropic interactions. We find that NR induces an effective propulsion of particle pairs and small aggregates ("active colloidal molecules") forming at the initial stages of self-assembly, an indication of the NR-induced non-equilibrium. The shape and stability of these initial clusters strongly depend on the degree of anisotropy. At longer times, we find, for weak NR interactions, large (even system-spanning) clusters where single particles can escape and enter at the boundaries, in stark contrast to the small rigid aggregates appearing at the same time in the passive case. In this sense, weak NR shortcuts the aggregation. Increasing the degree of NR (and thus, propulsion), we even observe large-scale phase separation if the interactions are weakly anisotropic. In contrast, systems with strong NR and anisotropy remain essentially disordered. Overall, the NR interactions are shown to destabilize the rigid aggregates interrupting self-assembly and phase separation in the passive case, thereby helping the system to overcome kinetic barriers.

Robustness of traveling states in generic nonreciprocal mixtures

Emergent nonreciprocal interactions violating Newton's third law are widespread in out-of-equilibrium systems. Phase separating mixtures with such interactions exhibit traveling states with no equilibrium counterpart. Using extensive Brownian dynamics simulations, we investigate the existence and stability of such traveling states in a generic nonreciprocal particle system. By varying a broad range of parameters including aggregate state of mixture components, diffusivity, degree of nonreciprocity, effective spatial dimension and density, we determine that traveling states do exist below the predator-prey regime, but nonetheless are only found in a narrow region of the parameter space. Our work also sheds light on the physical mechanisms for the disappearance of traveling states when relevant parameters are being varied, and has implications for a range of nonequilibrium systems including nonreciprocal phase separating mixtures, nonequilibrium pattern formation and predator-prey models.

Method of spectral response to stochastic processes for measuring the nonreciprocal effective interactions

The theoretical background of the nonperturbative method of spectral response to stochastic processes (SRSP) for measuring the nonreciprocal interparticle effective interactions in strongly coupled underdamped systems is described. Analytical expressions for vibrational spectral density of confined Brownian particles with a nonreciprocal effective interaction are presented. The changes in the vibrational spectral density with varying different parameters of the system (nonreciprocity, viscosity, ratios of particle sizes, and intensities of random processes acting on each particle) are discussed using the example of a pair of nonidentical particles in a harmonic trap. The SRSP method is compared to three other nonperturbative methods. The SRSP method demonstrates an undeniable advantage when processing particle trajectories with errors in particle tracking.

Heating leads to liquid-crystal and crystalline order in a two-temperature active fluid of rods

We report phase separation and liquid-crystal ordering induced by scalar activity in a system of soft repulsive spherocylinders (SRSs) of shape anisotropy L/D=5 using molecular dynamics (MD) simulations. Activity is introduced by increasing the temperature of half of the SRSs (labeled hot) while maintaining the temperature of the other half constant at a lower value (labeled cold). The difference between the two temperatures scaled by the lower temperature provides a measure of the activity. Starting from different equilibrium initial phases, we find that activity leads to segregation of the hot and cold particles. Activity also drives the cold particles through a phase transition to a more ordered state and the hot particles to a state of less order compared to the initial equilibrium state. The cold components of a homogeneous isotropic structure acquire nematic and, at higher activity, crystalline order. Similarly, the cold zone of a nematic initial state undergoes smectic and crystal ordering above a critical value of activity while the hot component turns isotropic. We find that the hot particles occupy a larger volume and exert an extra kinetic pressure, confining, compressing, and provoking an ordering transition of the cold-particle domains.

Active Ornstein-Uhlenbeck model for self-propelled particles with inertia

Self-propelled particles, which convert energy into mechanical motion, exhibit inertia if they have a macroscopic size or move inside a gaseous medium, in contrast to micron-sized overdamped particles immersed in a viscous fluid. Here we study an extension of the active Ornstein-Uhlenbeck model, in which self-propulsion is described by colored noise, to access these inertial effects. We summarize and discuss analytical solutions of the particle's mean-squared displacement and velocity autocorrelation function for several settings ranging from a free particle to various external influences, like a linear or harmonic potential and coupling to another particle via a harmonic spring. Taking into account the particular role of the initial particle velocity in a nonstationary setup, we observe all dynamical exponents between zero and four. After the typical inertial time, determined by the particle's mass, the results inherently revert to the behavior of an overdamped particle with the exception of the harmonically confined systems, in which the overall displacement is enhanced by inertia. We further consider an underdamped model for an active particle with a time-dependent mass, which critically affects the displacement in the intermediate time-regime. Most strikingly, for a sufficiently large rate of mass accumulation, the particle's motion is completely governed by inertial effects as it remains superdiffusive for all times.

Asymmetric acoustic wave scattering by a nonreciprocal and position-dependent mass defect

We investigate the asymmetric wave scattering in a phononic one-dimensional lattice with a nonreciprocal defect and position dependent masses coupled by the defect spring. The nonreciprocal interaction is characterized by a single parameter Δ while the nonlinear contribution due to position-dependent masses are controlled by a parameterχ. The transmission and reflection coefficients are analytically computed and the effects of the nonreciprocity and nonlinearity are detailed. We show that, in opposite with the linear case, the rectification factor has a frequency dependence, which leads to a more efficient diode-like action at large wavevectors. Further, the nonlinearity leads to an asymmetry of the reflected component, absent in the linear regime. We extend our analysis to a system with frictional forces which suppresses the multistability window promoted by the nonlinear mass contribution without compromising the rectification action.

Active Brownian particle in homogeneous media of different viscosities: numerical simulations

Self-propelled colloids, active polymers and membranes, driven (vibrated) granular layers and hybrid synthetic-biological systems are striking examples of systems containing synthetic active Brownian particles. Such particles autonomously convert the available energy of the environment into their own directed mechanical motion. In most studies the self-propelled Brownian particles move in overdamped media. Recently, experiments with Janus particles in a low-pressure plasma have appeared. A distinctive feature of such a medium is an extremely low viscosity at which the inertial effects play a significant role, resulting in underdamped Brownian motion. At present, there is a lack of statistical theory describing the underdamped Brownian motion of self-propelled particles at all time scales. This paper presents the numerical simulation results of active Brownian motion in homogeneous media of different viscosities. The calculations are performed using a mathematical model of a self-propelled Brownian sphere with translational and rotational inertia. The time-dependent mean square displacement and mean linear displacement (the noise-averaged trajectory) of the particle are investigated as a function of medium viscosity, self-propulsion velocity and moment of inertia. Our simulation reveals that the dynamics of a self-propelled spherical particle significantly depends on two independent dimensionless parameters of the particle: the ratio of the self-propulsion velocity to the characteristic thermal velocity and the ratio of the friction coefficient to the rotational diffusion coefficient. The obtained statistical characteristics of active Brownian motion are compared with the known theoretical models in a wide range of medium viscosities. We propose simple corrections to the basic theory of overdamped active Brownian motion, which allow one to calculate the effective diffusion coefficient and the persistence length of a self-propelled Brownian particle in a medium with any dynamic viscosity. The results obtained are discussed in relation to active particles in a colloidal plasma and superfluid helium.

Buckling of two-dimensional plasma crystals with nonreciprocal interactions

Laboratory realizations of two-dimensional (2D) plasma crystals typically involve monodisperse microparticles confined into horizontal monolayers in radio-frequency (rf) plasma sheaths. This gives rise to the so-called plasma wakes beneath the microparticles. The presence of wakes renders the interactions in such systems nonreciprocal, a fact that can lead to a quite different behavior from the one expected for their reciprocal counterparts. Here we examine the buckling of a hexagonal 2D plasma crystal, occurring as the confinement strength is decreased, taking explicitly into account the nonreciprocity of the system via a well-established point-wake model. We observe that for a finite wake charge, the monolayer hexagonal crystal undergoes a transition first to a bilayer hexagonal structure, unrealizable in harmonically confined reciprocal Yukawa systems, and subsequently to a bilayer square structure. Our theoretical results are confirmed by molecular dynamics simulations for experimentally relevant parameters, indicating the potential of their observation in state-of-the-art experiments with 2D complex plasmas.

Experimental study of the nonreciprocal effective interactions between microparticles in an anisotropic plasma

There is a variety of cases in nature when the action–reaction symmetry is broken. In particular, suitable conditions for this are realized in colloidal suspensions and complex plasmas. Since the first theories and simulations of the nonreciprocal effective interactions between microparticles in complex plasmas were published in 1995–1996, there have been hundreds of studies in the theoretical development of this theme. However, despite such a rich theoretical background, one of the important unsolved problems is a direct experimental determination of the nonreciprocal interparticle interaction forces. Here, we studied experimentally in detail the forces of the nonreciprocal effective interaction between microparticles suspended a radio-frequency produced plasma sheath. For this purpose, an experimental method based on an analysis of the spectral density of random processes in an open dissipative two-particle system was developed. In contrast to previous investigations, the proposed method takes into account random and dissipative processes in the system, does not require a special design of the experimental setup and any external perturbations, pre-measurements of external fields and any assumptions about the type of interaction. We found that even small charge changes of one particle, caused by its thermal motion in a wake field of another particle, can lead to a significant change in the effective (measurable) interaction between the particles.

Strange attractors induced by melting in systems with nonreciprocal effective interactions

Newton's third law-the action-reaction symmetry-can be violated for effective interbody forces in open and nonequilibrium systems that are ubiquitous in areas as diverse as complex plasmas, colloidal suspensions, active and living soft matter, and social behavior. While studying monolayer complex plasma (confined charged particles in an ionized gas) with nonreciprocal interactions mediated by plasma flows, in silico we found that an interplay between melting and thermal activation drastically transforms the collective dynamics: the order-disorder transition modifies the system's thermal steady state so that the crystal tends to melt, whereas the fluid tends to freeze, jumping chaotically between the two states. We identified this collective chaotic behavior as strange attractors formed in a monolayer complex plasma and link the strange attractor behavior to the specifics of interparticle interactions.

Dissipative phase transitions in systems with nonreciprocal effective interactions

The reciprocity of effective interparticle forces can be violated in various open and nonequilibrium systems, in particular, in colloidal suspensions and complex (dusty) plasmas. Here, we obtain a criterion under which a nonreciprocal system can be strictly reduced to a pseudo-Hamiltonian system with a detailed dynamic equilibrium. In particular, the criterion is satisfied for catalytically active colloids interacting via nonreciprocal diffusiophoretic forces. However, in the general case, when this criterion is not satisfied, the steady state is determined by the interplay between dissipation and the energy source due to the nonreciprocity of interactions. The results indicate the realization of bistability and dissipative spinodal decomposition in a broad class of systems with nonreciprocal effective interactions.

Effects of collectively induced scattering of gas stream by impurity ensembles: Shock-wave enhancement and disorder-stimulated nonlinear screening

We report on specific effects of collective scattering for a cloud of heavy impurities exposed to a gas stream. Formation is presented of a common density perturbation and shock waves, both generated collectively by a system of scatterers at sudden application of the stream-inducing external field. Our results demonstrate that (i) the scattering of gas stream can be essentially amplified, due to nonlinear collective effects, upon fragmentation of a solid obstacle into a cluster of impurities (heterogeneously fractured obstacle); (ii) a cluster of disordered impurities can produce considerably stronger scattering accompanied by enhanced and accelerated shock wave, as compared to a regularly ordered cluster. We also show that the final steady-state density distribution is formed as a residual perturbation left after the shock front passage. In particular, a kinklike steady distribution profile can be formed as a result of shock front stopping effect. The possibility of the onset of solitary diffusive density waves, reminiscent of precursor solitons, is shown and briefly discussed.

Cooperative behavior of biased probes in crowded interacting systems

We study, via extensive numerical simulations, dynamics of a crowded mixture of mutually interacting (with a short-range repulsive potential) colloidal particles immersed in a suspending solvent, acting as a heat bath. The mixture consists of a majority component - neutrally buoyant colloids subject to internal stimuli only, and a minority component - biased probes (BPs) also subject to a constant force. In such a system each of the BPs alters the distribution of the colloidal particles in its vicinity, driving their spatial distribution out of equilibrium. This induces effective long-range interactions and multi-tag correlations between the BPs, mediated by an out-of-equilibrium majority component, and prompts the BPs to move collectively assembling in clusters. We analyse the size-distribution of the self-assembling clusters in the steady-state, their specific force-velocity relations and also properties of the effective interactions emerging between the BPs.

Effects of gas interparticle interaction on dissipative wake-mediated forces

We examine how the short-range repulsive interaction in a gas of Brownian particles affects behavior of the nonequilibrium depletion forces between obstacles embedded into the gas flow. It is shown that for an ensemble of small and widely separated obstacles the dissipative wake-mediated interaction belongs to the type of induced dipole-dipole interaction governed by an anisotropic screened Coulomb-like potential. For closely located obstacles, formation of a common density perturbation "coat" around them leads to enhancement of dissipative interaction, manifested by characteristic peaks in its dependence on both the bath fraction and the external driving field. Moreover, additional screening of the gas flow due to nonlinear blockade effect gives rise to generation of a pronounced step-like profile of gas density distribution around the obstacles. This can lead to additional enhancement of dissipative interaction between obstacles. The possibility of the dissipative pairing effect and dissipative interaction switching provoked by wake inversion is briefly discussed. All the results are obtained within the classical lattice-gas model.