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

The Barrier Method: A Technique for Calculating Very Long Transition Times

278   0   0.0 ( 0 )
 نشر من قبل David Adams
 تاريخ النشر 2010
  مجال البحث فيزياء
والبحث باللغة English




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

In many dynamical systems there is a large separation of time scales between typical events and rare events which can be the cases of interest. Rare-event rates are quite difficult to compute numerically, but they are of considerable practical importance in many fields: for example transition times in chemical physics and extinction times in epidemiology can be very long, but are quite important. We present a very fast numerical technique that can be used to find long transition times (very small rates) in low-dimensional systems, even if they lack detailed balance. We illustrate the method for a bistable non-equilibrium system introduced by Maier and Stein and a two-dimensional (in parameter space) epidemiology model.


قيم البحث

اقرأ أيضاً

We show that in certain one-dimensional spin chains with open boundary conditions, the edge spins retain memory of their initial state for very long times. The long coherence times do not require disorder, only an ordered phase. In the integrable Isi ng and XYZ chains, the presence of a strong zero mode means the coherence time is infinite, even at infinite temperature. When Ising is perturbed by interactions breaking the integrability, the coherence time remains exponentially long in the perturbing couplings. We show that this is a consequence of an edge almost strong zero mode that almost commutes with the Hamiltonian. We compute this operator explicitly, allowing us to estimate accurately the plateau value of edge spin autocorrelator.
Biomolecular conformational transitions are usually modeled as barrier crossings in a free energy landscape. The transition paths connect two local free energy minima and transition path times (TPT) are the actual durations of the crossing events. Th e simplest model employed to analyze TPT and to fit empirical data is that of a stochastic particle crossing a parabolic barrier. Motivated by some disagreement between the value of the barrier height obtained from the TPT distributions as compared to the value obtained from kinetic and thermodynamic analyses, we investigate here TPT for barriers which deviate from the symmetric parabolic shape. We introduce a continuous set of potentials, that starting from a parabolic shape, can be made increasingly asymmetric by tuning a single parameter. The TPT distributions obtained in the asymmetric case are very well-fitted by distributions generated by parabolic barriers. The fits, however, provide values for the barrier heights and diffusion coefficients which deviate from the original input values. We show how these findings can be understood from the analysis of the eigenvalues spectrum of the Fokker-Planck equation and highlight connections with experimental results.
We test an improved finite-size scaling method for reliably extracting the critical temperature $T_{rm BKT}$ of a Berezinskii-Kosterlitz-Thouless (BKT) transition. Using known single-parameter logarithmic corrections to the spin stiffness $rho_s$ at $T_{rm BKT}$ in combination with the Kosterlitz-Nelson relation between the transition temperature and the stiffness, $rho_s(T_{rm BKT})=2T_{rm BKT}/pi$, we define a size dependent transition temperature $T_{rm BKT}(L_1,L_2)$ based on a pair of system sizes $L_1,L_2$, e.g., $L_2=2L_1$. We use Monte Carlo data for the standard two-dimensional classical XY model to demonstrate that this quantity is well behaved and can be reliably extrapolated to the thermodynamic limit using the next expected logarithmic correction beyond the ones included in defining $T_{rm BKT}(L_1,L_2)$. For the Monte Carlo calculations we use GPU (graphical processing unit) computing to obtain high-precision data for $L$ up to 512. We find that the sub-leading logarithmic corrections have significant effects on the extrapolation. Our result $T_{rm BKT}=0.8935(1)$ is several error bars above the previously best estimates of the transition temperature; $T_{rm BKT} approx 0.8929$. If only the leading log-correction is used, the result is, however, consistent with the lower value, suggesting that previous works have underestimated $T_{rm BKT}$ because of neglect of sub-leading logarithms. Our method is easy to implement in practice and should be applicable to generic BKT transitions.
The original perturbative Kramers method (starting from the phase space coordinates) (Kramers, 1940) of determining the energy-controlled-diffusion equation for Newtonian particles with separable and additive Hamiltonians is generalized to yield the energy-controlled diffusion equation and thus the very low damping (VLD) escape rate including spin-transfer torque for classical giant magnetic spins with two degrees of freedom. These have dynamics governed by the magnetic Langevin and Fokker-Planck equations and thus are generally based on non-separable and non-additive Hamiltonians. The derivation of the VLD escape rate directly from the (magnetic) Fokker-Planck equation for the surface distribution of magnetization orientations in the configuration space of the polar and azimuthal angles $(vartheta, varphi)$ is much simpler than those previously used.
292 - E. Carlon , H. Orland , T. Sakaue 2018
An analytical expression is derived for the transition path time distribution for a one-dimensional particle crossing of a parabolic barrier. Two cases are analyzed: (i) A non-Markovian process described by a generalized Langevin equation with a powe r-law memory kernel and (ii) a Markovian process with a noise violating the fluctuation-dissipation theorem, modeling the stochastic dynamics generated by active forces. In the case (i) we show that the anomalous dynamics strongly affecting the short time behavior of the distributions, but this happens only for very rare events not influencing the overall statistics. At long times the decay is always exponential, in disagreement with a recent study suggesting a stretched exponential decay. In the case (ii) the active forces do not substantially modify the short time behavior of the distribution, but lead to an overall decrease of the average transition path time. These findings offer some novel insights, useful for the analysis of experiments of transition path times in (bio)molecular systems.
التعليقات
جاري جلب التعليقات جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
mircosoft-partner

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