Escape of a passive particle from an activity-induced energy landscape: emergence of slow and fast effective diffusion

The effects of asymmetry in active noises on the efficiency of single colloidal Stirling engines with active noises

Molecular machines, which operate in highly fluctuating environments far from equilibrium, may benefit from their non-equilibrium environments. It is, however, a topic of controversy how the efficiency of the microscopic engines can be enhanced. Recent experiments showed that microscopic Stirling engines in bacterial reservoirs could show high performance beyond the equilibrium thermodynamics. In this work, we perform overdamped Langevin dynamics simulations for microscopic Stirling heat engines in bacterial reservoirs and show that the temperature dependence of the magnitude of active noises should be responsible for such high efficiency. Only when we introduce temperature-dependent active noises, the efficiency of the microscopic Stirling engines is enhanced significantly as in experiments.

Effect of confinement and topology: 2-TIPS vs. MIPS

2-TIPS (two temperature induced phase separation) refers to the phase separation phenomenon observed in mixtures of active and passive particles which are modelled using scalar activity. The active particles are connected to a thermostat at high temperature while the passive particles are connected to the thermostat at low temperature and the relative temperature difference between "hot" and "cold" particles is taken as the measure of the activity χ of the non-equilibrium system. The study of such binary mixtures of hot and cold particles under various kinds of confinement is an important problem in many physical and biological processes. The nature and extent of phase separation are heavily influenced by the geometry of confinement, activity, and density of the non-equilibrium binary mixture. Investigating such 3D binary mixtures confined by parallel walls, we observe that the active and passive particles phase separate, but the extent of phase separation is reduced compared to bulk phase separation at high densities and enhanced at low densities. However, when the binary mixture of active and passive particles is confined inside a spherical cavity, the phase separation is radial for small radii of the confining sphere and the extent of phase separation is higher compared to their bulk counterparts. Confinement leads to interesting properties in the passive (cold) region like enhanced layering and high compression in the direction parallel to the confining wall. In 2D, both the bulk and confined systems of the binary mixture show a significant decrement in the extent of phase separation at higher densities. This observation is attributed to the trapping of active particles inside the passive cluster, which increases with density. Thus the 2D systems show structures more akin to dense-dilute phase co-existence, which is observed in motility induced phase separation in 2D active systems. The binary mixture constrained on the spherical surface also shows similar phase co-existence. Our analyses reveal that the coexistent densities observed in 2-TIPS on the spherical surface agree with the findings of previous studies on MIPS in active systems on a sphere.

Escape dynamics of a self-propelled nanorod from circular confinements with narrow openings

We perform computer simulations to explore the escape dynamics of a self-propelled (active) nanorod from circular confinements with narrow opening(s). Our results clearly demonstrate how the persistent and directed motion of the nanorod helps it to escape. Such escape events are absent if the nanorod is passive. To quantify the escape dynamics, we compute the radial probability density function (RPDF) and mean first escape time (MFET) and show how the activity is responsible for the bimodality of RPDF, which is clearly absent if the nanorod is passive. Broadening of displacement distributions with activity has also been observed. The computed mean first escape time decreases with activity. In contrast, the fluctuations of the first escape times vary in a non-monotonic way. This results in high values of the coefficient of variation and indicates the presence of multiple timescales in first escape time distributions and multimodality in uniformity index distributions. We hope our study will help in differentiating activity-driven escape dynamics from purely thermal passive diffusion in confinement.

Method to probe the pronounced growth of correlation lengths in active glass-forming liquids using an elongated probe

The growth of correlation lengths in equilibrium glass-forming liquids near the glass transition is considered a critical finding in the quest to understand the physics of glass formation. These understandings helped us understand various dynamical phenomena observed in supercooled liquids. It is known that at least two different length scales exist; one is of thermodynamic origin, while the other is dynamical in nature. Recent observations of glassy dynamics in biological and synthetic systems where the external or internal driving source controls the dynamics, apart from the usual thermal noise, lead to the emergence of the field of active glassy matter. A question of whether the physics of glass formation in these active systems is also accompanied by growing dynamic and static lengths is indeed timely. In this article, we probe the growth of dynamic and static lengths in a model active glass system using rod-like elongated probe particles, an experimentally viable method. We show that the dynamic and static lengths in these nonequilibrium systems grow much more rapidly than their passive counterparts. We then offer an understanding of the violation of the Stokes-Einstein relation and Stokes-Einstein-Debye relation using these lengths via a scaling theory.

Phase behavior of active and passive dumbbells

We report phase separation in a mixture of "hot" and "cold" three-dimensional dumbbells which interact by Lennard-Jones potential. We also have studied the effect of asymmetry of dumbbells and the variation of ratio of "hot" and "cold" dumbbells on their phase separation. The ratio of the temperature difference between hot and cold dumbbells to the temperature of cold dumbbells is a measure of the activity χ of the system. From constant density simulations of symmetric dumbbells, we observe that the "hot" and "cold" dumbbells phase separate at higher activity ratio (χ>5.80) compared to that of a mixture of hot and cold Lennard-Jones monomers (χ>3.44). We find that, in the phase-separated system, the hot dumbbells have high effective volume and hence high entropy which is calculated by two-phase thermodynamic method. The high kinetic pressure of hot dumbbells forces the cold dumbbells to form dense clusters such that at the interface the high kinetic pressure of hot dumbbells is balanced by the virial pressure of cold dumbbells. We find that phase separation pushes the cluster of cold dumbbells to have solidlike ordering. Bond orientation order parameters reveal that the cold dumbbells form solidlike ordering consisting of predominantly face-centered cubic and hexagonal-close packing packing, but the individual dumbbells have random orientations. The simulation of the nonequilibrium system of symmetric dumbbells at different ratios of number of hot dumbbells to cold dumbbells reveals that the critical activity of phase separation decreases with increase in fraction of hot dumbbells. The simulation of equal mixture of hot and cold asymmetric dumbbells revealed that the critical activity of phase separation was independent of the asymmetry of dumbbells. We also observed that the clusters of cold asymmetric dumbbells showed both crystalline and noncrystalline order depending on the asymmetry of dumbbells.

Enhanced short time peak in four-point dynamic susceptibility in dense active glass-forming liquids

Active glassy systems are simple model systems that imitate complex biological processes. Sometimes, it becomes crucial to estimate the amount of activity present in such biological systems, such as predicting the progression rate of the cancer cells or the healing time of the wound, etc. In this work, we study a model active glassy system to quantify the degree of activity from the collective, long-wavelength fluctuations in the system. These long-wavelength fluctuations present themselves as an additional peak in the four-point dynamic susceptibility (χ₄(t)) apart from the usual peak at structural relaxation time. We then show how the degree of the activity at such a small timescale can be obtained by measuring the variation in χ₄(t) due to changing activity. A Detailed finite size analysis of the peak height of χ₄(t) suggests the existence of an intrinsic dynamic length scale that grows with increasing activity. Finally, we show that this peak height is a unique function of effective activity across all system sizes, serving as a possible parameter for characterizing the degree of activity in a system.

Chemically symmetric and asymmetric self-driven rigid dumbbells in a 2D polymer gel

We employ computer simulations to unveil the translational and rotational dynamics of self-driven chemically symmetric and asymmetric rigid dumbbells in a two-dimensional polymer gel. Our results show that the activity or the self-propulsion always enhances the dynamics of the dumbbells. Making the self-propelled dumbbell chemically asymmetric leads to further enhancement in dynamics. Additionally, the direction of self-propulsion is a key factor for chemically asymmetric dumbbells, where self-propulsion towards the non-sticky half of the dumbbell results in faster translational and rotational dynamics compared to the case with the self-propulsion towards the sticky half of the dumbbell. Our analyses show that both the symmetric and asymmetric passive rigid dumbbells get trapped inside the mesh of the polymer gel, but the chemical asymmetry always facilitates the mesh to mesh motion of the dumbbell and it is even more pronounced when the dumbbell is self-propelled. This results in multiple peaks in the van Hove function with increasing self-propulsion. In a nutshell, we believe that our in silico study can guide researchers to design efficient artificial microswimmers possessing potential applications in site-specific delivery.

Motion of an active particle with dynamical disorder

We propose a model for investigating the motion of a single active particle in a heterogeneous environment where the heterogeneity may arise due to crowding, conformational fluctuations and/or slow rearrangement of the surroundings. Describing the active particle in terms of the Ornstein-Uhlenbeck process (OUP) and incorporating heterogeneity in a thermal bath using two separate models, namely "diffusing diffusivity" and "switching diffusion", we explore the essential dynamical properties of the particle for its one-dimensional motion. In addition, we show how the dynamical behavior is controlled by dynamical variables associated with active noise such as strength and persistence time. Our model is relevant in the context of single particle dynamics in a crowded environment, driven by activity.

Passive and active tracer dynamics in polymer solutions with isotropic-to-nematic phase transition

Macromolecular crowding plays a crucial role in determining the dynamics in a living cell. We adopt Langevin dynamics simulations to investigate the anomalous diffusion dynamics of passive and active particles in a solution of polymer chains with tunable stiffness. The solution's anisotropic feature is modulated by changing both the polymer stiffness and volume fraction, where isotropic-to-nematic phase transition is involved. Our results demonstrate the significant impact of polymer flexibility on the dynamics of both passive and active probes. The distinct diffusion mechanism for an active particle is clarified by the interplay between polymer stiffness, crowdedness and activity. Polymer stiffness leads to a global inhibition effect on passive particle diffusion. The diffusion coefficient exhibits an intriguing non-monotonic variation at increasing polymer stiffness, which is due to the fact that the alignment of polymer chains is beneficial for diffusion along the nematic direction but unfavorable for that in the direction perpendicular to it. In sharp contrast, polymer stiffness plays a dominant role in facilitating active particle diffusion. Self-propulsion of the particle can utilize stiffness-induced elastic interactions more efficiently, which promotes its mobility in both directions. Meanwhile, an active particle might have a stronger ability to take advantage of the polymer alignment, contributing substantially enhanced diffusivity. In addition, the diffusion coefficient of an active particle is subject to a tendency of degeneration against varying volume fraction. This counter-intuitive behavior is due to the contrasting factors that increasing crowdedness induces a lower particle speed but a longer persistent motion time.

Correlated escape of active particles across a potential barrier

We study the dynamics of one-dimensional active particles confined in a double-well potential, focusing on the escape properties of the system, such as the mean escape time from a well. We first consider a single-particle both in near and far-from-equilibrium regimes by varying the persistence time of the active force and the swim velocity. A non-monotonic behavior of the mean escape time is observed with the persistence time of the activity, revealing the existence of an optimal choice of the parameters favoring the escape process. For small persistence times, a Kramers-like formula with an effective potential obtained within the unified colored noise approximation is shown to hold. Instead, for large persistence times, we developed a simple theoretical argument based on the first passage theory, which explains the linear dependence of the escape time with the persistence of the active force. In the second part of the work, we consider the escape on two active particles mutually repelling. Interestingly, the subtle interplay of active and repulsive forces may lead to a correlation between particles, favoring the simultaneous jump across the barrier. This mechanism cannot be observed in the escape process of two passive particles. Finally, we find that in the small persistence regime, the repulsion favors the escape, such as in passive systems, in agreement with our theoretical predictions, while for large persistence times, the repulsive and active forces produce an effective attraction, which hinders the barrier crossing.

Polymer looping kinetics in active heterogeneous environments

A typical biological environment is usually featured by crowding and heterogeneity, leading to complex reaction kinetics of the immersed macromolecules. In the present work, we adopt Langevin dynamics simulations to systematically investigate polymer looping kinetics in active heterogeneous media crowded with a mixture of mobile active particles and immobile obstacles. For comparison, a parallel study is also performed in the passive heterogeneous media. We explicitly analyze the change of looping time and looping probability with the variation of obstacle ratio, volume fraction and crowder size. We reveal the novel phenomena of inhibition-facilitation transition of the looping rate induced by heterogeneity, crowdedness and activity. In addition, our results demonstrate a very non-trivial crowder size effect on the looping kinetics. The underlying mechanism is rationalized by the interplay of polymer diffusion, conformational change and looping free-energy barrier. The competing effect arising from active particles and obstacles on structural and dynamical properties of the polymer yields a consistent scenario for our observations. Lastly, the non-exponential kinetics of the looping process is also analyzed. We find that both activity and crowding can strengthen the heterogeneity degree of the looping kinetics.

Breakdown of the Stokes-Einstein relation in supercooled water: the jump-diffusion perspective

Although water is the most ubiquitous liquid it shows many thermodynamic and dynamic anomalies. Some of the anomalies further intensify in the supercooled regime. While many experimental and theoretical studies have focused on the thermodynamic anomalies of supercooled water, fewer studies explored the dynamical anomalies very extensively. This is due to the intricacy of the experimental measurement of the dynamical properties of supercooled water. Violation of the Stokes-Einstein relation (SER), an important relation connecting the diffusion of particles with the viscosity of the medium, is one of the major dynamical anomalies. In absence of experimentally measured viscosity, researchers used to check the validity of SER indirectly using average translational relaxation time or α-relaxation time. Very recently, the viscosity of supercooled water was accurately measured at a wide range of temperatures and pressures. This allowed direct verification of the SER at different temperature-pressure thermodynamic state points. An increasing breakdown of the SER was observed with decreasing temperature. Increasing pressure reduces the extent of breakdown. Although some well-known theories explained the above breakdown, a detailed molecular mechanism was still elusive. Recently, a translational jump-diffusion (TJD) approach has been able to quantitatively explain the breakdown of the SER in pure supercooled water and an aqueous solution of methanol. The objective of this article is to present a detailed and state-of-the-art analysis of the past and present works on the breakdown of SER in supercooled water with a specific focus on the new TJD approach for explaining the breakdown of the SER.

Solution to the Kramers barrier crossing problem caused by two noises: Thermal noise and Poisson white noise

We consider the escape of a particle trapped in a metastable potential well and acted upon by two noises. One of the noises is thermal and the other is Poisson white noise, which is non-Gaussian. Using path integral techniques, we find an analytic solution to this generalization of the classic Kramers barrier crossing problem. Using the "barrier climbing" path, we calculate the activation exponent. We also derive an approximate expression for the prefactor. The calculated results are compared with the simulations, and a good agreement between the two is found. Our results show that, unlike in the case of thermal noise, the rate depends not just on the barrier height but also on the shape of the whole barrier. A comparison between the simulations and the theory also shows that a better approximation for the prefactor is needed for agreement for all values of the parameters.

Comparative study of polymer looping kinetics in passive and active environments

Intra-chain looping in complex environments is significant in advancing our understanding of biological processes in life. We adopt Langevin dynamics simulations to perform a comparative study of polymer looping kinetics in passive and active environments. From the analysis of looping quantities, including looping-unlooping times and looping probabilities, we unraveled the intriguing effects of active crowder size, activity and crowding. Firstly, we figured out the phase diagram involving a novel facilitation-inhibition transition in the parameter space of active crowder size and active force, and the two-fold roles of activity are clarified. In particular, we find that active particles of a size comparable to the polymer monomer are most favorable for facilitated looping, while those with a similar size to the polymer gyration radius impede the looping most seriously. Secondly, the underlying looping mechanisms in different active crowder size regimes are rationalized by the interplay among diffusion, polymer conformational change and the free-energy barrier. For small active crowders, activity significantly promotes end-to-end distance diffusion, which dominantly facilitates both looping and unlooping processes. In the case of moderate active crowders, the polymer chain suffers from prominent swelling, and thus inevitable inhibited looping will occur. For large active crowders, activity induces a counterintuitive non-cage effect on the looping kinetics, through yielding a higher effective temperature and larger unlooping free-energy barrier. This is in sharp contrast to the caging phenomena observed in passive media. Lastly, the volume-fraction dependence of the looping quantities in an active bath demonstrates dramatic discrepancies from that in a passive bath, which highlights the contrasting effects of activity and crowding.

Fluid Viscoelasticity Triggers Fast Transitions of a Brownian Particle in a Double Well Optical Potential

Thermally activated transitions are ubiquitous in nature, occurring in complex environments which are typically conceived as ideal viscous fluids. We report the first direct observations of a Brownian bead transiting between the wells of a bistable optical potential in a viscoelastic fluid with a single long relaxation time. We precisely characterize both the potential and the fluid, thus enabling a neat comparison between our experimental results and a theoretical model based on the generalized Langevin equation. Our findings reveal a drastic amplification of the transition rates compared to those in a Newtonian fluid, stemming from the relaxation of the fluid during the particle crossing events.

Effects of active crowder size and activity-crowding coupling on polymer translocation

Polymer translocation in complex environments is crucially important to many biological processes in life. In the present work, we adopted two-dimensional Langevin dynamics simulations to study the forced and unbiased polymer translocation dynamics in active and crowded media. The translocation time and probability are analyzed in terms of active force F_a, volume fraction φ and also the crowder size. The non-trivial active crowder size effect and activity-crowding coupling effect as well as the novel mechanism of unbiased translocation between two active environments with different active particle sizes are clarified. Firstly, for forced translocation, we reveal an intriguing non-monotonic dependence of the translocation time on the crowder size in the case of large activity. In particular, crowders of intermediate size similar to the polymer segment are proven to be the most favorable for translocation. Moreover, a facilitation-inhibition crossover of the translocation time with increasing volume fraction is observed, indicating a crucial activity-crowding coupling effect. Secondly, for unbiased translocation driven by different active crowder sizes, the translocation probability demonstrates a novel turnover phenomenon, implying the appearance of an opposite directional preference as the active force exceeds a critical value. The translocation time in both directions decreases monotonically with the active force. The asymmetric activity effect together with the entropic driving scenario provides a reasonable picture for the peculiar behavior observed in unbiased translocation.

Tunable depletion force in active and crowded environments

We adopt two-dimensional Langevin dynamics simulations to study the effective interactions between two passive colloids in a bath crowded with active particles. We mainly pay attention to the significant effects of active particle size, crowding-activity coupling, and chirality. First, a transition of depletion force from repulsion to attraction is revealed by varying particle size. Moreover, larger active crowders with sufficient activity can generate strong attractive force, which is in contrast to the cage effect in passive media. It is interesting that the attraction induced by large active crowders follows a linear scaling with the persistence length of active particles. Second, the effective force also experiences a transition from repulsion to attraction as volume fraction increases, as a consequence of the competition between the two contrastive factors of activity and crowding. As bath volume fraction is relatively small, activity generates a dominant repulsion force, while as the bath becomes concentrated, crowding-induced attraction becomes overwhelming. Lastly, in a chiral bath, we observe a very surprising oscillation phenomenon of active depletion force, showing an evident quasiperiodic variation with increasing chirality. Aggregation of active particles in the vicinity of the colloids is carefully examined, which serves as a reasonable picture for our observations. Our findings provide an inspiring strategy for the tunable active depletion force by crowding, activity, and chirality.

Active Brownian and inertial particles in disordered environments: Short-time expansion of the mean-square displacement

We consider an active Brownian particle moving in a disordered two-dimensional energy or motility landscape. The averaged mean-square displacement (MSD) of the particle is calculated analytically within a systematic short-time expansion. As a result, for overdamped particles, both an external random force field and disorder in the self-propulsion speed induce ballistic behavior adding to the ballistic regime of an active particle with sharp self-propulsion speed. Spatial correlations in the force and motility landscape contribute only to the cubic and higher-order powers in time for the MSD. Finally, for inertial particles two superballistic regimes are found where the scaling exponent of the MSD with time is α=3 and α=4. We confirm our theoretical predictions by computer simulations. Moreover, they are verifiable in experiments on self-propelled colloids in random environments.

Narrow-escape time and sorting of active particles in circular domains

It is now well established that microswimmers can be sorted or segregated fabricating suitable microfluidic devices or using external fields. A natural question is how these techniques can be employed for dividing swimmers of different motility. In this paper, using numerical simulations in the dilute limit, we investigate how motility parameters (time of persistence and velocity) impact the narrow-escape time of active particles from circular domains. We show that the escape time undergoes a crossover between two asymptotic regimes. The control parameters of the crossover is the ratio between the persistence length of the active motion and the typical length scale of the circular domain. We explore the possibility of taking advantage of this finding for sorting active particles by motility parameters.

Translational and rotational dynamics of a self-propelled Janus probe in crowded environments

We computationally investigate the dynamics of a self-propelled Janus probe in crowded environments. The crowding is caused by the presence of viscoelastic polymers or non-viscoelastic disconnected monomers. Our simulations show that the translational as well as rotational mean square displacements have a distinctive three-step growth for fixed values of self-propulsion force, and steadily increase with self-propulsion, irrespective of the nature of the crowder. On the other hand, in the absence of crowders, the rotational dynamics of the Janus probe is independent of self-propulsion force. On replacing the repulsive polymers with sticky ones, translational and rotational mean square displacements of the Janus probe show a sharp drop. Since different faces of a Janus particle interact differently with the environment, we show that the direction of self-propulsion also affects its dynamics. The ratio of long-time translational and rotational diffusivities of the self-propelled probe with a fixed self-propulsion, when plotted against the area fraction of the crowders, passes through a minimum and at higher area fraction merges to its value in the absence of the crowder. This points towards the decoupling of the translational and rotational dynamics of the self-propelled probe at an intermediate area fraction of the crowders. However, such translational-rotational decoupling is absent for passive probes.

Field-driven tracer diffusion through curved bottlenecks: fine structure of first passage events

Using scaling arguments and extensive numerical simulations, we study the dynamics of a tracer particle in a corrugated channel represented by a periodic sequence of broad chambers and narrow funnel-like bottlenecks enclosed by a hard-wall boundary.