Waiting time distribution for continuous stochastic systems

Analytical refractory period distribution for a class of time-variant biochemical systems with second-order reactions

Refractory period (RP), the waiting time between signals, can induce complex signaling dynamics, such as acceleration, adaptation, and oscillation, within many cellular biochemical networks. However, its underlying molecular mechanisms are still unclear. Rigorously estimating the RP distribution may be essential to identify its causal regulatory mechanisms. Traditional methods of estimating the RP distribution depend on solving the underlying Chemical Master Equations (CMEs), the dominant modeling formalism of biochemical systems. However, exact solutions of the CME are only known for simple reaction systems with zero- and first-order reactions or specific systems with second-order reactions. General solutions still need to be derived for systems with bimolecular reactions. It is even more challenging if large state-space and nonconstant reaction rates are involved. Here, we developed a direct method to gain the analytical RP distribution for a class of second-order reaction systems with nonconstant reaction rates and large state space. Instead of using the CME, we used an equivalent path-wise representation, which is the solution to a transformed martingale problem of the CME. This allowed us to bypass solving a CME. We then applied the method to derive the analytical RP distribution of a real complex biochemical network with second-order reactions, the Drosophila phototransduction cascade. Our approach provides an alternative to the CMEs in deriving the analytical RP distributions of a class of second-order reaction systems. Since the bimolecular reactions are common in biological systems, our approach could enhance understanding real-world biochemical processes.

Data-Driven Langevin Modeling of Nonequilibrium Processes

Given nonstationary data from molecular dynamics simulations, a Markovian Langevin model is constructed that aims to reproduce the time evolution of the underlying process. While at equilibrium the free energy landscape is sampled, nonequilibrium processes can be associated with a biased energy landscape, which accounts for finite sampling effects and external driving. When the data-driven Langevin equation (dLE) approach [Phys. Rev. Lett. 2015, 115, 050602] is extended to the modeling of nonequilibrium processes, an efficient way to calculate multidimensional Langevin fields is outlined. The dLE is shown to correctly account for various nonequilibrium processes, including the enforced dissociation of sodium chloride in water, the pressure-jump induced nucleation of a liquid of hard spheres, and the conformational dynamics of a helical peptide sampled from nonstationary short trajectories.

Dynamical modes of sheared confined microscale matter

Based on (overdamped) Stokesian dynamics simulations and video microscopy experiments, we study the non equilibrium dynamics of a sheared colloidal cluster, which is confined to a two-dimensional disk. The experimental system is composed of a mixture of paramagnetic and non magnetic polystyrene particles, which are held in the disk by time shared optical tweezers. The paramagnetic particles are located at the center of the disk and are actuated by an external, rotating magnetic field that induces a magnetic torque. We identify two different steady states by monitoring the mean angular velocities per ring. The first one is characterized by rare slip events, where the inner rings momentarily depin from the outer ring, which is kept static by the set of optical traps. For the second state, we find a bistability of the mean angular velocities, which can be understood from the analysis of the slip events in the particle trajectories. We calculate the particle waiting- and jumping time distributions and estimate a time scale between slips, which is also reflected by a plateau in the mean squared azimuthal displacement. The dynamical transition is further reflected by the components of the stress tensor, revealing a shear-thinning behavior as well as shear stress overshoots. Finally, we briefly discuss the observed transition in the context of stochastic thermodynamics and how it may open future directions in this field.

Metadynamics Enhanced Markov Modeling of Protein Dynamics

Enhanced sampling techniques represent a versatile approach to account for rare conformational transitions in biomolecules. A particularly promising strategy is to combine massive parallel computing of short molecular dynamics (MD) trajectories (to sample the free energy landscape of the system) with Markov state modeling (to rebuild the kinetics from the sampled data). To obtain well-distributed initial structures for the short trajectories, it is proposed to employ metadynamics MD, which quickly sweeps through the entire free energy landscape of interest. Being only used to generate initial conformations, the implementation of metadynamics can be simple and fast. The conformational dynamics of helical peptide Aib₉ is adopted to discuss various technical issues of the approach, including metadynamics settings, minimal number and length of short MD trajectories, and the validation of the resulting Markov models. Using metadynamics to launch some thousands of nanosecond trajectories, several Markov state models are constructed that reveal that previous unbiased MD simulations of in total 16 μs length cannot provide correct equilibrium populations or qualitative features of the pathway distribution of the short peptide.

Force-linearization closure for non-Markovian Langevin systems with time delay

This paper is concerned with the Fokker-Planck (FP) description of classical stochastic systems with discrete time delay. The non-Markovian character of the corresponding Langevin dynamics naturally leads to a coupled infinite hierarchy of FP equations for the various n-time joint distribution functions. Here, we present an approach to close the hierarchy at the one-time level based on a linearization of the deterministic forces in all members of the hierarchy starting from the second one. This leads to a closed equation for the one-time probability density in the steady state. Considering two generic nonlinear systems, a colloidal particle in a sinusoidal or bistable potential supplemented by a linear delay force, we demonstrate that our approach yields a very accurate representation of the density as compared to quasiexact numerical results from direct solution of the Langevin equation. Moreover, the results are significantly improved against those from a small-delay approximation and a perturbation-theoretical approach. We also discuss the possibility of accessing transport-related quantities, such as escape times, based on an additional Kramers approximation. Our approach applies to a wide class of models with nonlinear deterministic forces.

Depinning and heterogeneous dynamics of colloidal crystal layers under shear flow

Using Brownian dynamics (BD) simulations and an analytical approach we investigate the shear-induced, nonequilibrium dynamics of dense colloidal suspensions confined to a narrow slit-pore. Focusing on situations where the colloids arrange in well-defined layers with solidlike in-plane structure, the confined films display complex, nonlinear behavior such as collective depinning and local transport via density excitations. These phenomena are reminiscent of colloidal monolayers driven over a periodic substrate potential. In order to deepen this connection, we present an effective model that maps the dynamics of the shear-driven colloidal layers to the motion of a single particle driven over an effective substrate potential. This model allows us to estimate the critical shear rate of the depinning transition based on the equilibrium configuration, revealing the impact of important parameters, such as the slit-pore width and the interaction strength. We then turn to heterogeneous systems where a layer of small colloids is sheared with respect to bottom layers of large particles. For these incommensurate systems we find that the particle transport is dominated by density excitations resembling the so-called "kink" solutions of the Frenkel-Kontorova (FK) model. In contrast to the FK model, however, the corresponding "antikinks" do not move.

Time- and ensemble-averages in evolving systems: the case of Brownian particles in random potentials

Anomalous diffusion is a ubiquitous phenomenon in complex systems. It is often quantified using time- and ensemble-averages to improve statistics, although time averages represent a non-local measure in time and hence can be difficult to interpret. We present a detailed analysis of the influence of time- and ensemble-averages on dynamical quantities by investigating Brownian particles in a rough potential energy landscape (PEL). Initially, the particle ensemble is randomly distributed, but the occupancy of energy values evolves towards the equilibrium distribution. This relaxation manifests itself in the time evolution of time- and ensemble-averaged dynamical measures. We use Monte Carlo simulations to study particle dynamics in a potential with a Gaussian distribution of energy values, where the long-time limit of the diffusion coefficient is known from theory. In our experiments, individual colloidal particles are exposed to a laser speckle pattern inducing a non-Gaussian roughness and are followed by optical microscopy. The relaxation depends on the kind and degree of roughness of the PEL. It can be followed and quantified by the time- and ensemble-averaged mean squared displacement. Moreover, the heterogeneity of the dynamics is characterized using single-trajectory analysis. The results of this work are relevant for the correct interpretation of single-particle tracking experiments in general.

Enhancement of mobility in an interacting colloidal system under feedback control

Feedback control schemes are a promising way to manipulate transport properties of driven colloidal suspensions. In the present article, we suggest a feedback scheme to enhance the collective transport of colloidal particles with repulsive interactions through a one-dimensional tilted washboard potential. The control is modeled by a harmonic confining potential, mimicking an optical "trap," with the center of this trap moving with the (instantaneous) mean particle position. Our theoretical analysis is based on the Smoluchowski equation combined with dynamical density functional theory for systems with hard-core or ultrasoft (Gaussian) interactions. For either type of interaction, we find that the feedback control can lead to an enhancement of the mobility by several orders of magnitude relative to the uncontrolled case. The largest effects occur for intermediate stiffness of the trap and large particle numbers. Moreover, in some regions of the parameter space the feedback control induces oscillations of the mean velocity. Finally, we show that the enhancement of mobility is robust against a small time delay in implementing the feedback control.