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

Variationally Optimized Free Energy Flooding for Rate Calculation

82   0   0.0 ( 0 )
 نشر من قبل James McCarty
 تاريخ النشر 2015
  مجال البحث فيزياء
والبحث باللغة English




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

We propose a new method to obtain kinetic properties of infrequent events from molecular dynamics simulation. The procedure employs a recently introduced variational approach [Valsson and Parrinello, Phys. Rev. Lett. 113, 090601 (2014)] to construct a bias potential as a function of several collective variables that is designed to flood only the associated free energy surface up to a predefined level. The resulting bias potential effectively accelerates transitions between metastable free energy minima while ensuring bias-free transition states, thus allowing accurate kinetic rates to be obtained. We test the method on a few illustrative systems for which we obtain an order of magnitude improvement in efficiency relative to previous approaches, and several orders of magnitude relative to unbiased molecular dynamics. We expect an even larger improvement in more complex systems. This and the ability of the variational approach to deal efficiently with a large number of collective variables will greatly enhance the scope of these calculations. This work is a vindication of the potential that the variational principle has if applied in innovative ways

قيم البحث

اقرأ أيضاً

Free energy is energy that is available to do work. Maximizing the free energy gain and the gain in work that can be extracted from a system is important for a wide variety of physical and technological processes, from energy harvesting processes suc h as photosynthesis to energy storage systems such as fuels and batteries. This paper extends recent results from non-equilibrium thermodynamics and quantum resource theory to derive closed-form solutions for the maximum possible gain in free energy and extractable work that can be obtained by varying the initial states of classical and quantum stochastic processes. Simple formulae allow the comparison the free energy increase for the optimal procedure with that for a sub-optimal procedure. The problem of finding the optimal free-energy harvesting procedure is shown to be convex and solvable via gradient descent.
354 - A. Laio , M. Parrinello 2002
We introduce a novel and powerful method for exploring the properties of the multidimensional free energy surfaces of complex many-body systems by means of a coarse-grained non-Markovian dynamics in the space defined by a few collective coordinates.A characteristic feature of this dynamics is the presence of a history-dependent potential term that, in time, fills the minima in the free energy surface, allowing the efficient exploration and accurate determination of the free energy surface as a function of the collective coordinates. We demonstrate the usefulness of this approach in the case of the dissociation of a NaCl molecule in water and in the study of the conformational changes of a dialanine in solution.
170 - Jun Zhang , Yao-Kun Lei , Xing Che 2019
In this paper we first analyzed the inductive bias underlying the data scattered across complex free energy landscapes (FEL), and exploited it to train deep neural networks which yield reduced and clustered representation for the FEL. Our parametric method, called Information Distilling of Metastability (IDM), is end-to-end differentiable thus scalable to ultra-large dataset. IDM is also a clustering algorithm and is able to cluster the samples in the meantime of reducing the dimensions. Besides, as an unsupervised learning method, IDM differs from many existing dimensionality reduction and clustering methods in that it neither requires a cherry-picked distance metric nor the ground-true number of clusters, and that it can be used to unroll and zoom-in the hierarchical FEL with respect to different timescales. Through multiple experiments, we show that IDM can achieve physically meaningful representations which partition the FEL into well-defined metastable states hence are amenable for downstream tasks such as mechanism analysis and kinetic modeling.
71 - A. Hietanen , B. Lucini 2011
We propose a general method (based on the Wang-Landau algorithm) to compute numerically free energies that are obtained from the logarithm of the ratio of suitable partition functions. As an application, we determine with high accuracy the order-orde r interface tension of the four-state Potts model in three dimensions on cubic lattices of linear extension up to L=56. The infinite volume interface tension is then extracted at each beta from a fit of the finite volume interface tension to a known universal behavior. A comparison of the order-order and order-disorder interface tension at the critical value of beta provides a clear numerical evidence of perfect wetting.
Monte Carlo simulations of crystal nuclei coexisting with the fluid phase in thermal equilibrium in finite volumes are presented and analyzed, for fluid densities from dense melts to the vapor. Generalizing the lever-rule for two-phase coexistence in the canonical ensemble to finite volume, measurements of the nucleus volume together with the pressure and chemical potential of the surrounding fluid allows to extract the surface free energy of the nucleus. Neither the knowledge of the (in general non-spherical) nucleus shape nor of the angle-dependent interface tension is required for this task. The feasibility of the approach is demonstrated for a variant of the Asakura-Oosawa model for colloid-polymer mixtures, which form face-centered cubic colloidal crystals. For a polymer to colloid size ratio of $0.15$, the colloid packing fraction in the fluid phase can be varied from melt values to zero by the variation of an effective attractive potential between the colloids. It is found that the approximation of spherical crystal nuclei often underestimates actual nucleation barriers significantly. Nucleation barriers are found to scale as $Delta F^*=(4pi/3)^{1/3}bar{gamma}(V^*)^{2/3}+const.$ with the nucleus volume $V^*$, and the effective surface tension $bar{gamma}$ that accounts implicitly for the nonspherical shape can be precisely estimated.
التعليقات
جاري جلب التعليقات جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
mircosoft-partner

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