Nonequilibrium diffusion of active particles bound to a semi-flexible polymer network: simulations and fractional Langevin equation

Abstract

In a viscoelastic environment, the diffusion of a particle becomes non-Markovian due to the memory effect. An open question is to quantitatively explain how self-propulsion particles with directional memory diffuse in such a medium. Based on simulations and analytic theory, we address this issue with active viscoelastic systems where an active particle is connected with multiple semi-flexible filaments. Our Langevin dynamics simulations show that the active cross-linker displays super- and sub-diffusive athermal motion with a time-dependent anomalous exponent $\alpha$. In such viscoelastic feedback, the active particle always has superdiffusion with $\alpha =3/2$ at times shorter than the self-propulsion time (${\tau }_A$). At times greater than ${\tau }_A$, the subdiffusion emerges with $\alpha$ bounded between $1/2$ and $3/4$. Remarkably, the active subdiffusion is reinforced as the active propulsion (Pe) is more vigorous. In the high-Pe limit, the athermal fluctuation in the stiff filament eventually leads to $\alpha =1/2$, which can be misinterpreted with the thermal Rouse motion in a flexible chain. We demonstrate that the motion of active particles cross-linking a network of semi-flexible filaments can be governed by a fractional Langevin equation combined with fractional Gaussian noise and an Ornstein-Uhlenbeck noise. We analytically derive the velocity autocorrelation function and mean-squared displacement of the model, explaining their scaling relations as well as the prefactors. We find that there exist the threshold Pe (${\mathrm{Pe}}^{\ast }$) and cross-over times (${\tau }^{\ast }$ and ${\tau }^{†}$) above which the active viscoelastic dynamics emerge on the timescales of ${\tau }^{\ast }\lesssim t\lesssim {\tau }^{†}$. Our study may provide a theoretical insight into various nonequilibrium active dynamics in intracellular viscoelastic environments.