اشترك بالحزمة الذهبية واحصل على وصول غير محدود شمرا أكاديميا
تسجيل مستخدم جديدExtremely rare ultra-fast non-equilibrium processes can be close to equilibrium: RNA unfolding and refolding
103
0
0.0
(
0
)
اسأل ChatGPT حول البحث
ﻻ يوجد ملخص باللغة العربية
We study numerically the behavior of RNA secondary structures under influence of a varying external force. This allows to measure the work $W$ during the resulting fast unfolding and refolding processes. Here, we investigate a medium-size hairpin structure. Using a sophisticated large-deviation algorithm, we are able to measure work distributions with high precision down to probabilities as small as $10^{-46}$. Due to this precision and by comparison with exact free-energy calculations we are able to verify the theorems of Crooks and Jarzynski. Furthermore, we analyze force-extension curves and the configurations of the secondary structures during unfolding and refolding for typical equilibrium processes and non-equilibrium processes, conditioned to selected values of the measured work $W$, typical and rare ones. We find that the non-equilibrium processes where the work values are close to those which are most relevant for applying Crooks and Jarzynski theorems, respectively, are most and quite similar to the equilibrium processes. Thus, a similarity of equilibrium and non-equilibrium behavior with respect to a mere scalar variable, which occurs with a very small probability but can be generated in a controlled but non-targeted way, is related to a high similarity for the set of configurations sampled along the full dynamical trajectory.
قيم البحث
اقرأ أيضاً
In this paper we propose a new formalism to map history-dependent metadynamics in a Markovian process. We apply this formalism to a model Langevin dynamics and determine the equilibrium distribution of a collection of simulations. We demonstrate that
the reconstructed free energy is an unbiased estimate of the underlying free energy and analytically derive an expression for the error. The present results can be applied to other history-dependent stochastic processes such as Wang-Landau sampling.
The mechanical unfolding of a simple RNA hairpin and of a 236--bases portion of the Tetrahymena thermophila ribozyme is studied by means of an Ising--like model. Phase diagrams and free energy landscapes are computed exactly and suggest a simple two-
-state behaviour for the hairpin and the presence of intermediate states for the ribozyme. Nonequilibrium simulations give the possible unfolding pathways for the ribozyme, and the dominant pathway corresponds to the experimentally observed one.
Estimating the expected value of an observable appearing in a non-equilibrium stochastic process usually involves sampling. If the observables variance is high, many samples are required. In contrast, we show that performing the same task without sam
pling, using tensor network compression, efficiently captures high variances in systems of various geometries and dimensions. We provide examples for which matching the accuracy of our efficient method would require a sample size scaling exponentially with system size. In particular, the high variance observable $mathrm{e}^{-beta W}$, motivated by Jarzynskis equality, with $W$ the work done quenching from equilibrium at inverse temperature $beta$, is exactly and efficiently captured by tensor networks.
The Jarzynski identity can be applied to instances when a microscopic system is pulled repeatedly but quickly along some coordinate, allowing the calculation of an equilibrium free energy profile along the pulling coordinate from a set of independent
non-equilibrium trajectories. Using the formalism of Wiener stochastic path integrals in which we assign temperature-dependent weights to Langevin trajectories, we derive exact formulae for the temperature derivatives of the free energy profile. This leads naturally to analytical expressions for decomposing a free energy profile into equilibrium entropy and internal energy profiles from non-equilibrium pulling. This decomposition can be done from trajectories evolved at a unique temperature without repeating the measurement as done in finite-difference decompositions. Three distinct analytical expressions for the entropy-energy decomposition are derived: using a time-dependent generalization of the weighted histogram analysis method, a quasi harmonic spring limit, and a Feynman-Kac formula. The three novel formulae of reconstructing the pair of entropy-energy profiles are exemplified by Langevin simulations of a two-dimensional model system prototypical for force-induced biomolecular conformational changes. Connections to single-molecule experimental means to probe the functionals needed in the decomposition are suggested.
We demonstrate that the clustering statistics and the corresponding phase transition to non-equilibrium clustering found in many experiments and simulation studies with self-propelled particles (SPPs) with alignment can be obtained from a simple kine
tic model. The key elements of this approach are the scaling of the cluster cross-section with the cluster mass -- characterized by an exponent $alpha$ -- and the scaling of the cluster perimeter with the cluster mass -- described by an exponent $beta$. The analysis of the kinetic approach reveals that the SPPs exhibit two phases: i) an individual phase, where the cluster size distribution (CSD) is dominated by an exponential tail that defines a characteristic cluster size, and ii) a collective phase characterized by the presence of non-monotonic CSD with a local maximum at large cluster sizes. At the transition between these two phases the CSD is well described by a power-law with a critical exponent $gamma$, which is a function of $alpha$ and $beta$ only. The critical exponent is found to be in the range $0.8 < gamma < 1.5$ in line with observations in experiments and simulations.
سجل دخول لتتمكن من نشر تعليقات
التعليقات
جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
