Driven anomalous diffusion: An example from polymer stretching

Transition Path Dynamics of Non-Markovian Systems across a Rough Potential Barrier

Transition paths refer to rare events in physics, chemistry, and biology where the molecules cross barriers separating stable molecular conformations. The conventional analysis of the transition path times employs a diffusive and memoryless transition over a smooth potential barrier. However, it is widely acknowledged that the free energy profile between two minima in biomolecular processes is inherently not smooth. In this article, we discuss a theoretical model with a parabolic rough potential barrier and obtain analytical results of the transition path distribution and mean transition path times by incorporating absorbing boundary conditions across the boundaries under the driving of Gaussian white noise. Further, the influence of anomalous dynamics in rough potential driven by a power-law memory kernel is analyzed by deriving a time-dependent scaled diffusion coefficient that coarse-grains the effects of roughness, and the system's dynamics is reduced to a scaled diffusion on a smooth potential. Our theoretical results are tested and validated against numerical simulations. The findings of our study show the influence of the boundary conditions, barrier height, barrier roughness, and memory effect on the transition path time distributions in a rough potential, and the validity of the scaling diffusion coefficient has been discussed.

Structural dynamics and optimal transport of an active polymer

We study the spontaneous configuration transitions of an active semi-flexible polymer between spiral and non-spiral states, and show that the configuration dynamics is fully described by a subcritical pitchfork bifurcation. Exploiting the fact that an active polymer barely moves in spiral states and exhibits net displacements in non-spiral states, we theoretically prove that the motion of the active polymer is consistent with a run-and-tumble-like dynamics. Moreover, we find that there exists an optimal self-propelling force that maximizes the diffusion coefficient.

Generalized Langevin dynamics for single beads in linear elastic networks

We derive generalized Langevin equations (GLEs) for single beads in linear elastic networks. In particular, the derivations of the GLEs are conducted without employing normal modes, resulting in two distinct representations in terms of resistance and mobility kernels. The fluctuation-dissipation relations are also confirmed for both GLEs. Subsequently, we demonstrate that these two representations are interconnected via Laplace transforms. Furthermore, another GLE is derived by utilizing a projection operator method, and it is shown that the equation obtained through the projection scheme is consistent with the GLE with the resistance kernel. As simple examples, the general theory is applied to the Rouse model and the ring polymer, where the GLEs with the resistance and mobility kernels are explicitly derived for arbitrary positions of the tagged bead in these models. Finally, the GLE with the mobility kernel is also derived for the elastic network with hydrodynamic interactions under the pre-averaging approximation.

Chromatin Remodeling Due to Transient-Link-and-Pass Activity Enhances Subnuclear Dynamics

Spatiotemporal coordination of chromatin and subnuclear compartments is crucial for cells. Numerous enzymes act inside nucleus-some of those transiently link and pass two chromatin segments. Here, we study how such an active perturbation affects fluctuating dynamics of an inclusion in the chromatic medium. Using numerical simulations and a versatile effective model, we categorize inclusion dynamics into three distinct modes. The transient-link-and-pass activity speeds up inclusion dynamics by affecting a slow mode related to chromatin remodeling, viz., size and shape of the chromatin meshes.

Langevin analogy between particle trajectories and polymer configurations

A diffusive trajectory drawn by the generalized Langevin equation (GLE) for a colloidal particle evokes a random fractal of a static polymer configuration. This article proposes a static GLE-like description that enables the generation of a single configuration of a polymer chain with the noise formulated to satisfy the static fluctuation-response relation (FRR) along a one-dimensional chain structure but not along a temporal coordinate. A remarkable point is qualitative differences and similarities in the FRR formulation between the static and the dynamic GLEs. Guided by the static FRR, we further make analogous arguments in light of stochastic energetics and the steady-state fluctuation theorem.

Preaveraging description of polymer nonequilibrium stretching

This article focuses on a preaveraging description of polymer nonequilibrium stretching, where a single polymer undergoes a transient process from equilibrium to nonequilibrium steady state by pulling one chain end. The preaveraging method combined with mode analysis reduces the original Langevin equation to a simplified form for both a stretched steady state and an equilibrium state, even in the presence of self-avoiding repulsive interactions spanning a long range. However, the transient stretching process exhibits evolution of a hierarchal regime structure, which means a qualitative temporal change in probabilistic distributions assumed in preaveraging. We investigate the preaveraging method for evolution of the regime structure with consideration of the nonequilibrium work relations and deviations from the fluctuation-dissipation relation.

Projection of strong coupling interaction with thermal bath in a polymer

We investigate modifications of a stochastic polymer picture through a shift in the boundary between the system and an external environment. A conventional bead-and-spring model serving as the coarse-graining model is given by the Langevin equation for all the monomers subject to white noise. However, stochastic motion for only a tagged monomer is observed to occur in the presence of colored noise. The qualitative change in the observations arises from the boundary shift decided by the observer. The Langevin dynamics analyses interpret the colored noise as the emergence of the polymeric elastic force, resulting in additional heat in the tagged monomer observation. Being distinguished from coarse-graining based on scale separation, the projection of comparable internal degrees of freedom is also discussed in light of the fluctuation theorem and the stochastic polymer thermodynamics.

Negative friction memory induces persistent motion

We investigate the mean-square displacement (MSD) for random motion governed by the generalized Langevin equation for memory functions that contain two different time scales: In the first model, the memory kernel consists of a delta peak and a single-exponential and in the second model of the sum of two exponentials. In particular, we investigate the scenario where the long-time exponential kernel contribution is negative. The competition between positive and negative friction memory contributions produces an enhanced transient persistent regime in the MSD, which is relevant for biological motility and active matter systems.

Approximate dynamical eigenmodes of the Ising model with local spin-exchange moves

We establish that the Fourier modes of the magnetization serve as the dynamical eigenmodes for the two-dimensional Ising model at the critical temperature with local spin-exchange moves, i.e., Kawasaki dynamics. We obtain the dynamical scaling properties for these modes and use them to calculate the time evolution of two dynamical quantities for the system, namely, the autocorrelation function and the mean-square deviation of the line magnetizations. At intermediate times 1≲t≲L^{z_{c}}, where z_{c}=4-η=15/4 is the dynamical critical exponent of the model, we find that the line magnetization undergoes anomalous diffusion. Following our recent work on anomalous diffusion in spin models, we demonstrate that the generalized Langevin equation with a memory kernel consistently describes the anomalous diffusion, verifying the corresponding fluctuation-dissipation theorem with the calculation of the force autocorrelation function.

Inferring Active Noise Characteristics from the Paired Observations of Anomalous Diffusion

Anomalous diffusion has been most often argued in terms of a position fluctuation of a tracer. We here propose the other fluctuating observable, i.e., momentum transfer defined as the time integral of applied force to hold a tracer’s position. Being a conjugated variable, the momentum transfer is thought of as generating the anomalous diffusion paired with the position’s one. By putting together the paired anomalous diffusions, we aim to extract useful information in complex systems, which can be applied to experiments like tagged monomer observations in chromatin. The polymer being in the equilibrium, the mean square displacement (or variance) of position displacement or momentum transfer exhibits the sub- or superdiffusion, respectively, in which the sum of the anomalous diffusion indices is conserved quite generally, but the nonequilibrium media that generate the active noise may manifest the derivations from the equilibrium relation. We discuss the deviations that reflect the characteristics of the active noise.

Translocation of Charged Polymers through a Nanopore in Monovalent and Divalent Salt Solutions: A Scaling Study Exploring over the Entire Driving Force Regimes

Langevin dynamics simulations are performed to study polyelectrolytes driven through a nanopore in monovalent and divalent salt solutions. The driving electric field E is applied inside the pore, and the strength is varied to cover the four characteristic force regimes depicted by a rederived scaling theory, namely the unbiased (UB) regime, the weakly-driven (WD) regime, the strongly-driven trumpet (SD(T)) regime and the strongly-driven isoflux (SD(I)) regime. By changing the chain length N, the mean translocation time is studied under the scaling form 〈 τ 〉 ∼ N α E - δ . The exponents α and δ are calculated in each force regime for the two studied salt cases. Both of them are found to vary with E and N and, hence, are not universal in the parameter's space. We further investigate the diffusion behavior of translocation. The subdiffusion exponent γ p is extracted. The three essential exponents ν s , q, z p are then obtained from the simulations. Together with γ p , the validness of the scaling theory is verified. Through a comparison with experiments, the location of a usual experimental condition on the scaling plot is pinpointed.

Two-tag correlations and nonequilibrium fluctuation-response relation in ageing single-file diffusion

Spatiotemporally correlated motions of interacting Brownian particles, confined in a narrow channel of infinite length, are studied in terms of statistical quantities involving two particles. A theoretical framework that allows analytical calculation of two-tag correlations is presented on the basis of the Dean-Kawasaki equation describing density fluctuations in colloidal systems. In the equilibrium case, the time-dependent Einstein relation holds between the two-tag displacement correlation and the response function corresponding to it, which is a manifestation of the fluctuation-dissipation theorem for the correlation of density fluctuations. While the standard procedure of closure approximation for nonlinear density fluctuations is known to be obstructed by inconsistency with the fluctuation-dissipation theorem, this difficulty is naturally avoided by switching from the standard Fourier representation of the density field to the label-based Fourier representation of the vacancy field. In the case of ageing dynamics started from equidistant lattice configuration, the time-dependent Einstein relation is violated, as the two-tag correlation depends on the waiting time for equilibration while the response function is not sensitive to it. Within linear approximation, however, there is a simple relation between the density (or vacancy) fluctuations and the corresponding response function, which is valid even if the system is out of equilibrium. This non-equilibrium fluctuation-response relation can be extended to the case of nonlinear fluctuations by means of closure approximation for the vacancy field.

Generalized Langevin equation formulation for anomalous diffusion in the Ising model at the critical temperature

We consider the two- (2D) and three-dimensional (3D) Ising models on a square lattice at the critical temperature T_{c}, under Monte Carlo spin flip dynamics. The bulk magnetization and the magnetization of a tagged line in the 2D Ising model, and the bulk magnetization and the magnetization of a tagged plane in the 3D Ising model, exhibit anomalous diffusion. Specifically, their mean-square displacements increase as power laws in time, collectively denoted as ∼t^{c}, where c is the anomalous exponent. We argue that the anomalous diffusion in all these quantities for the Ising model stems from time-dependent restoring forces, decaying as power laws in time-also with exponent c -in striking similarity to anomalous diffusion in polymeric systems. Prompted by our previous work that has established a memory-kernel based generalized Langevin equation (GLE) formulation for polymeric systems, we show that a closely analogous GLE formulation holds for the Ising model as well. We obtain the memory kernels from spin-spin correlation functions, and the formulation allows us to consistently explain anomalous diffusion as well as anomalous response of the Ising model to an externally applied magnetic field in a consistent manner.

Deviation from fluctuation-dissipation relation for driven superdiffusion: Polymer stretching as an example

We discuss a deviation of the fluctuation-dissipation relation (FDR) in a driven superdiffusive system as exemplified by polymer stretching. The superdiffusion is found by monitoring momentum transfer to a tracer, which is a conjugate observable with the position. Molecular-dynamics simulation demonstrates that the FDR deviates during the nonequilibrium transient process. We then propose nonequilibrium mode analysis for superdiffusion, which is a counterpart to that for driven subdiffusion. The mode analysis yields results that are in qualitative agreement with the simulation results, suggesting that the fluctuations of the stiffness in the system from initial equilibrium to stretching account for the FDR deviation.

Non-Markovian dynamics of reaction coordinate in polymer folding

We develop a theoretical description of the critical zipping dynamics of a self-folding polymer. We use tension propagation theory and the formalism of the generalized Langevin equation applied to a polymer that contains two complementary parts which can bind to each other. At the critical temperature, the (un)zipping is unbiased and the two strands open and close as a zipper. The number of broken base pairs n(t) displays a subdiffusive motion characterized by a variance growing as 〈Δn²(t)〉 ∼ t^α with α < 1 at long times. Our theory provides an estimate of both the asymptotic anomalous exponent α and of the subleading correction term, which are both in excellent agreement with numerical simulations. The results indicate that the tension propagation theory captures the relevant features of the dynamics and shed some new insights on related polymer problems characterized by anomalous dynamical behavior.

Complementary mode analyses between sub- and superdiffusion

Several subdiffusive stochastic processes in nature, e.g., the motion of a tagged monomer in polymers, the height fluctuation of interfaces, particle dynamics in single-file diffusion, etc., can be described rigorously or approximately by the superposition of various modes whose relaxation times are broadly distributed. In this paper, we propose a mode analysis generating superdiffusion, which is paired with or complementary to subdiffusion. The key point in our discussion lies in the identification of a pair of conjugated variables, which undergo sub- and superdiffusion, respectively. We provide a simple interpretation for the sub- and superdiffusion duality for these variables using the language of polymer physics. The analysis also suggests the usefulness of looking at the force fluctuation in experiments, where a polymer is driven by a constant velocity.

Active diffusion of model chromosomal loci driven by athermal noise

Active diffusion, i.e., fluctuating dynamics driven by athermal noise, is found in various out-of-equilibrium systems. Here we discuss the nature of the active diffusion of tagged monomers in a flexible polymer. A scaling argument based on the notion of tension propagation clarifies how the polymeric effect is reflected in the anomalous diffusion exponent, which may be of relevance to the dynamics of chromosomal loci in living cells.

Dynamics of Polymer Translocation: A Short Review with an Introduction of Weakly-Driven Regime

As emphasized in a recent review (by V.V. Palyulin, T. Ala-Nissila, R. Metzler), theoretical understanding of the unbiased polymer translocation lags behind that of the (strongly) driven translocation. Here, we suggest the introduction of a weakly-driven regime, as described by the linear response theory to the unbiased regime, which is followed by the strongly-driven regime beyond the onset of nonlinear response. This provides a concise crossover scenario, bridging the unbiased to strongly-driven regimes.

Conformation Change, Tension Propagation and Drift-Diffusion Properties of Polyelectrolyte in Nanopore Translocation

Using Langevin dynamics simulations, conformational, mechanical and dynamical properties of charged polymers threading through a nanopore are investigated. The shape descriptors display different variation behaviors for the cis- and trans-side sub-chains, which reflects a strong cis-trans dynamical asymmetry, especially when the driving field is strong. The calculation of bond stretching shows how the bond tension propagates on the chain backbone, and the chain section straightened by the tension force is determined by the ratio of the direct to the contour distances of the monomer to the pore. With the study of the waiting time function, the threading process is divided into the tension-propagation stage and the tail-retraction stage. At the end, the drift velocity, diffusive property and probability density distribution are explored. Owing to the non-equilibrium nature, translocation is not a simple drift-diffusion process, but exhibits several intermediate behaviors, such as ballistic motion, normal diffusion and super diffusion, before ending with the last, negative-diffusion behavior.