Escape of an active ring from an attractive surface: Behaving like a self-propelled Brownian particle.

Escape of active agents from metastable states is of great interest in statistical and biological physics. In this paper, we investigate the escape of a flexible active ring, composed of active Brownian particles, from a flat attractive surface using Brownian dynamics simulations. To systematically explore the effects of activity, persistence time, and the shape of attractive potentials, we calculate escape time τ_{e} and effective temperature T_{eff}. We observe two distinct escape mechanisms: Kramers-like thermal activation at small persistence times, where ln(τ_{e})∼1/(k_{B}T_{eff}), and the maximal force problem at large persistence time, where τ_{e} is determined by persistence time. The escape time explicitly depends on the shape of the potential barrier at high activity and large persistence time. Moreover, when the propulsion force is biased along the ring's contour, escape becomes more difficult and is primarily driven by thermal noise. Our findings highlight that, despite its intricate configuration, the active ring can be effectively modeled as a self-propelled Brownian particle when studying its escape from a smooth surface.

Deformation-induced phase separation of active vesicles.

Many active materials, such as bacteria and cells, are deformable. Deformability significantly affects their collective behaviors and movements in complex environments. Here, we introduce a two-dimensional deformable active vesicle (DAV) model to emulate cell-like deformable active matter, wherein the deformability can be continuously adjusted. We find that changes in deformability can induce phase separation of DAVs. The system can transition between a homogeneous gas state, a coexistence of gas and liquid, and a coexistence of gas and solid. The occurrence of deformation-induced phase separation is accompanied by nonmonotonic changes in effective concentration, particle size and shape. Moreover, the degree of deformability also impacts the motility and stress within the dense phase following phase separation. Our results offer new insights into the role of deformability in the collective behavior of active matter.

Constrained motion of self-propelling eccentric disks linked by a spring.

It has been supposed that the interplay of elasticity and activity plays a key role in triggering the non-equilibrium behaviors in biological systems. However, the experimental model system is missing to investigate the spatiotemporally dynamical phenomena. Here, a model system of an active chain, where active eccentric-disks are linked by a spring, is designed to study the interplay of activity, elasticity, and friction. Individual active chain exhibits longitudinal and transverse motions; however, it starts to self-rotate when pinning one end and self-beat when clamping one end. In addition, our eccentric-disk model can qualitatively reproduce such behaviors and explain the unusual self-rotation of the first disk around its geometric center. Furthermore, the structure and dynamics of long chains were studied via simulations without steric interactions. It was found that a hairpin conformation emerges in free motion, while in the constrained motions, the rotational and beating frequencies scale with the flexure number (the ratio of self-propelling force to bending rigidity), χ, as ∼(χ)4/3. Scaling analysis suggests that it results from the balance between activity and energy dissipation. Our findings show that topological constraints play a vital role in non-equilibrium synergy behaviors.

Configuration and dynamics of a self-propelled diblock copolymer chain.

Active polymers are slender or chain-like self-propelled objects. Synthetic chains of self-propelled colloidal particles are one of the examples, which provide a potential way to develop varied active polymers. Here, we study the configuration and dynamics of an active diblock copolymer chain. Our focus is on the competition and the cooperation between the equilibrium self-assembly due to chain heterogeneity and the dynamic self-assembly due to propulsion. Simulations show that an active diblock copolymer chain can form the spiral(+)/tadpole(+) states under forward propulsion and the spiral(-)/tadpole(-)/bean states under backward propulsion. Interestingly, it is easier for the backward-propelled chain to form a spiral. The transitions between the states can be analyzed in terms of work and energy. For forward propulsion, we found a key quantity, i.e. the chirality of the packed self-attractive A block, which determines the configuration of the whole chain and the dynamics. However, no such quantity is found for the backward propulsion. Our results set the foundation for further study of the self-assembly of multiple active copolymer chains and provide a reference for the design and application of polymeric active materials.

Scaling behavior for the detachment of a self-propelling filament from an attractive surface.

Desorption of a self-propelling filament from an attractive surface is studied by computer simulations and the influence of activity, chain length, and chain rigidity is explored. For the flexible filament, we find three scaling regimes of desorption time vs activity with various scaling exponents. At low activity, the scaling law results from the spiral-like detachment kinetics. And at high activity, by theoretical analysis, the desorption is reminiscent of the escaping mechanism of a super-diffusive blob from a potential well at a short time scale. Additionally, the desorption time decreases first and then increases with chain length at low activity, since it is hard to form a spiral for short filaments due to the limited volume repulsion. For high activities, the desorption time approximately scales with chain length, with a scaling exponent ∼0.5, which can be explained by the theory and numerically fitting scaling law between the end-to-end distance of the "globule-like" filament and chain length. Furthermore, a non-monotonic behavior is observed between the desorption time and the chain stiffness. Desorption time slightly decreases first and then rapidly increases with stiffness due to the opposed effects of increasing rigidity on headiing-up time and leaving-away time. In contrast to traditional polymers, the scaling behavior suggests unique desorption characteristics of active polymers.

Insights into aggregation dynamics of NACore peptides from coarse-grained simulations.

Alpha(α)-synuclein is closely related to the pathogenesis of Parkinson's disease (PD). The NACore, a fragment of α-synuclein, is considered to be the key region of α-synuclein that causes PD. The aggregation dynamics of NACores are studied via coarse-grained molecular dynamics simulations. We find that NACores can self-assemble into a large cluster at high concentrations. The aggregation dynamics can be divided into three stages. The growth kinetics for the first and second stages follows the power law, Smax ~ tγ , with the second stage faster than the first one. The characteristic lifetime for the high concentration is 40 times larger than that for the low concentration, implying the low fluidity. Understanding the aggregation dynamics of NACores is helpful to develop drugs for therapeutic prevention and intervention.

Obstacle-induced giant jammed aggregation of active semiflexible filaments.

Filaments driven by bound motor proteins and chains of self-propelled colloidal particles are a typical example of active polymers (APs). Due to deformability, APs exhibit very rich dynamic behaviors and collective assembling structures. Here, we are concerned with a basic question: how APs behave near a single obstacle? We find that, in the presence of a big single obstacle, the assembly of APs becomes a two-state system, i.e. APs either gather nearly completely together into a giant jammed aggregate (GJA) on the surface of the obstacle or distribute freely in space. No partial aggregation is observed. Such a complete aggregation/collection is unexpected since it happens on a smooth convex surface instead of, e.g., a concave wedge. We find that the formation of a GJA experiences a process of nucleation and the curves of the transition between the GJA and the non-aggregate state form hysteresis-like loops. Statistical analysis of massive data on the growing time, chirality and angular velocity of both the GJAs and the corresponding nuclei shows the strong random nature of the phenomenon. Our results provide new insights into the behavior of APs in contact with porous media and also a reference for the design and application of polymeric active materials.

Kick effect of enzymes causes filament compression.

We investigate the influence of enzymes on the structure and dynamics of a filament by dissipative particle dynamics simulations. Enzyme exerts a kick force on the filament monomer. We pay particular attention to two factors: the magnitude of kick force and enzyme concentration. Large kick force as well as high enzyme concentration prefers a remarkable compression of the filament reminiscent of the effective depletion interaction owing to an effective increase in enzyme size and the reduction of solvent quality. Additionally, the kick effect gives rise to an increase of enzyme density from the center-of-mass of the filament to its periphery. Moreover, the increase of enzyme concentration and kick force also causes a decrease in relaxation time. Our finding is helpful to understand the role of catalytic force in chemo-mechano-biological function and the filament behavior under chemical reaction via kick-induced change of solvent quality.

Structure and dynamics of an active polymer adsorbed on the surface of a cylinder.

The structure and dynamics of an active polymer on a smooth cylindrical surface are studied by Brownian dynamics simulations. The effect of an active force on the polymer adsorption behavior and the combined effect of chain mobility, length N, rigidity κ, and cylinder radius, R, on the phase diagrams are systemically investigated. We find that complete adsorption is replaced by the irregular alternative adsorption/desorption process at a large driving force. Three typical (spiral, helix-like, and rod-like) conformations of the active polymer are observed, dependent on N, κ, and R. Dynamically, the polymer shows rotational motion in the spiral state, snake-like motion in the intermediate state, and straight translational motion without turning back in the rod-like state. In the spiral state, we find that the rotation velocity ω and the chain length follow a power-law relation ω ∼ N-0.42, consistent with the torque-balance theory of general Archimedean spirals. And the polymer shows super-diffusive behavior along the cylinder for a long time in the helix-like and rod-like states. Our results highlight that the mobility, rigidity, and curvature of surface can be used to regulate the polymer behavior.

Extracting free energies of counterion binding to polyelectrolytes by molecular dynamics simulations.

We use all-atom molecular dynamics simulations to extract ΔGeff, the free energy of binding of potassium ions K+ to the partially charged polyelectrolyte poly(acrylic acid), or PAA, in dilute regimes. Upon increasing the charge fraction of PAA, the chains adopt more extended conformations, and simultaneously, potassium ions bind more strongly (i.e., with more negative ΔGeff) to the highly charged chains to relieve electrostatic repulsions between charged monomers along the chains. We compare the simulation results with the predictions of a model that describes potassium binding to PAA chains as a reversible reaction whose binding free energy (ΔGeff) is adjusted from its intrinsic value (ΔG) by electrostatic correlations, captured by a random phase approximation. The bare or intrinsic binding free energy ΔG, which is an input in the model, depends on the binding species and is obtained from the radial distribution function of K+ around the charged monomer of a singly charged, short PAA chain in dilute solutions. We find that the model yields semi-quantitative predictions for ΔGeff and the degree of potassium binding to PAA chains, α, as a function of PAA charge fraction without using fitting parameters.

Absorption of self-propelled particles into a dense porous medium.

We study the absorption of self-propelled particles into a finite-size dense porous medium, which is mimicked by an obstacle array. We find that, depending on the competition of the propelling strength versus the repulsive barrier formed by obstacles and the contrast between the characteristic time scales of permeation and propelling persistence, the absorption process exhibits three distinct types of behavior. In Type I and II behavior, the propelling strength is not large enough to surmount the barrier, and hence particles transport in the medium by barrier-hopping dynamics. The initial permeation of particles toward the medium center is phenomenologically similar to a normal slow diffusion process. But, surprisingly, after the initial permeation process, a concentrated nucleus of particle aggregates forms and grows at the medium center in Type I, due to the long propelling persistence. Such an abnormal "nucleation" phenomenon does not appear in Type II, in which the propelling persistence is low. When the propelling strength is very high (Type III), particles transport smoothly in the medium, hence the initial slow diffusion process disappears and small particle clusters form and merge randomly in the medium. Our results provide a foundation for applications of active objects in a complex environment and also suggest the possible usage of a porous medium, for example, in the selection or sorting of active matter.

Vortex formation of spherical self-propelled particles around a circular obstacle.

A vortex is a common ratchet phenomenon in active systems. The spatial symmetry is usually broken by introducing asymmetric shapes or spontaneously by collective motion in the presence of hydrodynamic interactions or other alignment effects. Unexpectedly, we observe, by simulations, the formation of a vortex in the simplest model of a circular obstacle immersed in a bath of spherical self-propelled particles. No symmetry-breaking factors mentioned above are included in this model. The vortex forms only when the particle activity is high, i.e. large persistence. The obstacle size is also a key factor and the vortex only forms in a limited range of obstacle sizes. The sustainment of the vortex originates from the bias of the rotating particle cluster around the obstacle in accepting the incoming particles based on their propelling directions. Our results provide new understanding of and insights into the spontaneous symmetry-breaking and ratchet phenomena in active matter.

Transport of self-propelled particles across a porous medium: trapping, clogging, and the Matthew effect.

We study the transport of self-propelled particles from one free chamber to another across two stripe-like areas of dense porous medium. The medium is mimicked by arrays of obstacles. We find that active motion could greatly speed up the transport of particles. However, more and more particles become trapped in the obstacle arrays with the enhancement of activity. At high persistence (low rotational diffusion rate) and moderate particle concentration, we observe the Matthew effect in the aggregation of particles in the two obstacle arrays. This effect is weakened by introduction of randomness or deformability into the obstacle arrays. Moreover, the dependence on deformability shows the characteristics of first-order phase transition. In rare situations, the system could be stuck in a dynamic unstable state, e.g. the particles alternatively gather more in one of the two obstacle arrays, exhibiting oscillation of particle number between the arrays. Our results reveal new features in the transport of active objects in a complex medium and have implications for manipulating their collective assembly.

Assembly structures and dynamics of active colloidal cells.

Many types of active matter are deformable, such as epithelial cells and bacteria. To mimic the feature of deformability, we built a model called an active colloidal cell (ACC), i.e. a vesicle enclosed with self-propelled particles (SPPs), which as a whole can move actively. Based on the model, we then study the role of deformability in the assembly structures and dynamics of ACCs by Langevin dynamics simulation. We find that deformability weakens the self-trapping effect and hence suppresses the clustering and phase separation of the deformable soft ACCs (sACCs). Instead of forming a large compact cluster like ordinary SPPs, sACCs pack into a loose network or porous structure in the phase-separation region. The condensed phase is liquid-like, in which sACCs are strongly compressed and deformed but still keep high motility. The interface between the gas and the condensed phases is blurry and unstable, and the effective interfacial energy is very low. Our work gives new insights into the role of deformability in the assembly of active matter and also provides a reference for further studies on different types of deformable active matter.

Unfolding of a diblock chain and its anomalous diffusion induced by active particles.

We study the structural and dynamical behavior of an A-B diblock chain in the bath of active Brownian particles (ABPs) by Brownian dynamics simulations in two dimensions. We are interested in the situation that the effective interaction between the A segments is attractive, while that between the B segments is repulsive. Therefore, in thermal (nonactive) equilibrium, the A block "folds" into a compact globule, while the B block is in the expanded coil state. Interestingly, we find that the A block could "unfold" sequentially like unknitting a sweater, driven by the surrounding ABPs when the propelling strength on them is beyond a certain value. This threshold value decreases and then levels off as the length of the B block increases. We also find a simple power-law relation between the unfolding time of the A block and the self-propelling strength and an exponential relation between the unfolding time and the length of the B block. Finally, we probe the translational and rotational diffusion of the chain and find that both of them show "super-diffusivity" in a large time window, especially when the self-propelling strength is small and the A block is in the folded state. Such super-diffusivity is due to the strong asymmetric distribution of ABPs around the chain. Our work provides new insights into the behavior of a polymer chain in the environment of active objects.

Shape transformation and manipulation of a vesicle by active particles.

Langevin dynamics simulations are employed to study the shape transformation of a two-dimensional vesicle induced by active particles both inside and outside. We find that the shape of the vesicle changes from circle, to capsule, and eventually to dumbbell with the enhancement of the particle activity. Under the cooperation between the inside and the outside active particles, such significant shape transformation is realized by tuning the activity in a small range. And unexpectedly, the fluctuations of the capsule and the dumbbell shapes are not completely random but mostly along the direction of the short axis. In the situation of strong activity, the inside of the dumbbell vesicle is analogous to a system of two chambers, which are connected by a narrow channel. Intriguingly, we observe the vibration of the channel width, accompanied with the exchange of active particles between the two chambers. We also find that dynamical manipulation of the vesicle shape is possible through tuning the particle activity dynamically. This work provides new ideas to the control of the vesicle morphology and new insights into the dynamics in the vesicle's shape transformation.

Globule-stretch transition of a self-attracting chain in the repulsive active particle bath.

Folding and unfolding of a chain structure are often manipulated in experiments by tuning the pH, temperature, single-molecule forces or shear fields. Here, we carry out Brownian dynamics simulations to explore the behavior of a single self-attracting chain in a suspension of self-propelling particles (SPPs). As the propelling force increases, the globule-stretch (G-S) transition of the chain occurs due to the enhanced disturbance from the SPPs. Two distinct mechanisms of the transition in the limits of low and high rotational diffusion rates of SPPs have been observed: shear-induced stretching at a low rate and collision-induced melting at a high rate. The G-S and S-G (stretch-globule) curves form a hysteresis loop at the low rate, while they merge at the high rate. Besides, we find that two competing effects result in a non-monotonic dependence of the G-S transition on SPP density at the low rate. Our results suggest an alternative approach to manipulating the folding and unfolding of (bio)polymers by utilizing active agents.

Microrotor of a chain-grafted colloidal disk immersed in the active bath: The impact of particle concentration, grafting density, and chain rigidity.

In an earlier work, we discussed the possibility to realize a microrotor by immersing a chain-grafted colloidal disk in a thin film of active-particle suspension. Under certain conditions, the colloidal disk rotates unidirectionally driven by the bath active particles. Here we systematically study the role of active-particle concentration, grafting density, and chain rigidity in the phenomenon of the spontaneous symmetry breaking of the chain configurations and the unidirectional rotation of the disk. We find that high chain rigidity can help stabilize both the collective asymmetric chain configurations and the rotation of the disk, while it has a weak impact on the rotational speed/efficiency. Increasing the number of grafted chains can also stabilize the rotation but has a negative impact on the rotational speed/efficiency. Active particles power the rotation of the colloidal disk, yet their contribution saturates beyond a certain concentration. Our work provides new insights into the active systems with chain-structured objects and the design of soft/deformable micromachines.

Molecular dynamics simulation of G-actin interacting with PAMAM dendrimers.

Understanding the interactions of dendrimers as drug/gene delivery vectors with proteins is important for functional optimization. Here, atomistic molecular dynamics simulations are employed to study the interactions between six positively-charged polyamidoamine dendrimers of the second generation (G2 PAMAM) and G-actin. We find that the structure of G-actin is relatively stable after dendrimers' binding. PAMAM dendrimers also do not significantly change the secondary structure of G-actin. Furthermore, we find the formation of dendrimer-actin complex is mainly driven by electrostatic interactions. Moreover, we suggest the secondary structure change of local domains of G-actin is probably responsible to the inhibition of actin polymerization.

Beating of grafted chains induced by active Brownian particles.

We study the interplay between active Brownian particles (ABPs) and a "hairy" surface in two-dimensional geometry. We find that the increase of propelling force leads to and enhances inhomogeneous accumulation of ABPs inside the brush region. Oscillation of chain bundles (beating like cilia) is found in company with the formation and disassembly of a dynamic cluster of ABPs at large propelling forces. Meanwhile chains are stretched and pushed down due to the effective shear force by ABPs. The decrease of the average brush thickness with propelling force reflects the growth of the beating amplitude of chain bundles. Furthermore, the beating phenomenon is investigated in a simple single-chain system. We find that the chain swings regularly with a major oscillatory period, which increases with chain length and decreases with the increase of propelling force. We build a theory to describe the phenomenon and the predictions on the relationship between the period and amplitude for various chain lengths, and propelling forces agree very well with simulation data.