ترغب بنشر مسار تعليمي؟ اضغط هنا

Markov models from the Square Root Approximation of the Fokker-Planck equation: calculating the grid-dependent flux

202   0   0.0 ( 0 )
 نشر من قبل Luca Donati Dr.
 تاريخ النشر 2020
  مجال البحث فيزياء
والبحث باللغة English




اسأل ChatGPT حول البحث

Molecular dynamics are extremely complex, yet understanding the slow components of their dynamics is essential to understanding their macroscopic properties. To achieve this, one models the molecular dynamics as a stochastic process and analyses the dominant eigenfunctions of the associated Fokker-Planck operator, or of closely related transfer operators. So far, the calculation of the discretized operators requires extensive molecular dynamics simulations. The Square-root approximation of the Fokker-Planck equation is a method to calculate transition rates as a ratio of the Boltzmann densities of neighboring grid cells times a flux, and can in principle be calculated without a simulation. In a previous work we still used molecular dynamics simulations to determine the flux. Here, we propose several methods to calculate the exact or approximate flux for various grid types, and thus estimate the rate matrix without a simulation. Using model potentials we test computational efficiency of the methods, and the accuracy with which they reproduce the dominant eigenfunctions and eigenvalues. For these model potentials, rate matrices with up to $mathcal{O}(10^6)$ states can be obtained within seconds on a single high-performance compute server if regular grids are used.



قيم البحث

اقرأ أيضاً

61 - C.-L. Ho 2016
We consider the solvability of the Fokker-Planck equation with both time-dependent drift and diffusion coefficients by means of the similarity method. By the introduction of the similarity variable, the Fokker-Planck equation is reduced to an ordinar y differential equation. Adopting the natural requirement that the probability current density vanishes at the boundary, the resulting ordinary differential equation turns out to be integrable, and the probability density function can be given in closed form. New examples of exactly solvable Fokker-Planck equations are presented.
76 - Choon-Lin Ho 2021
A procedure is presented for solving the Fokker-Planck equation with constant diffusion but non-stationary drift. It is based on the correspondence between the Fokker-Planck equation and the non-stationary Schrodinger equation. The formalism of super symmetric quantum mechanics is extended by applying the Darboux transformation directly to the non-stationary Schrodinger equation. From a solution of a Fokker-Planck equation a solution of the Darboux partner is obtained. For drift coefficients satisfying the condition of shape invariance, a supersymmetric hierarchy of Fokker-Planck equations with solutions related by the Darboux transformation is obtained.
For the analysis of molecular processes, the estimation of time-scales, i.e., transition rates, is very important. Estimating the transition rates between molecular conformations is -- from a mathematical point of view -- an invariant subspace projec tion problem. A certain infinitesimal generator acting on function space is projected to a low-dimensional rate matrix. This projection can be performed in two steps. First, the infinitesimal generator is discretized, then the invariant subspace is approxi-mated and used for the subspace projection. In our approach, the discretization will be based on a Voronoi tessellation of the conformational space. We will show that the discretized infinitesimal generator can simply be approximated by the geometric average of the Boltzmann weights of the Voronoi cells. Thus, there is a direct correla-tion between the potential energy surface of molecular structures and the transition rates of conformational changes. We present results for a 2d-diffusion process and Alanine dipeptide.
We consider extinction times for a class of birth-death processes commonly found in applications, where there is a control parameter which determines whether the population quickly becomes extinct, or rather persists for a long time. We give an exact expression for the discrete case and its asymptotic expansion for large values of the population. We have results below the threshold, at the threshold, and above the threshold (where there is a quasi-stationary state and the extinction time is very long.) We show that the Fokker-Planck approximation is valid only quite near the threshold. We compare our analytical results to numerical simulations for the SIS epidemic model, which is in the class that we treat. This is an interesting example of the delicate relationship between discrete and continuum treatments of the same problem.
161 - S. I. Denisov 2009
We study the connection between the parameters of the fractional Fokker-Planck equation, which is associated with the overdamped Langevin equation driven by noise with heavy-tailed increments, and the transition probability density of the noise gener ating process. Explicit expressions for these parameters are derived both for finite and infinite variance of the rescaled transition probability density.
التعليقات
جاري جلب التعليقات جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
mircosoft-partner

هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا