- Journal Home
- Volume 36 - 2024
- Volume 35 - 2024
- Volume 34 - 2023
- Volume 33 - 2023
- Volume 32 - 2022
- Volume 31 - 2022
- Volume 30 - 2021
- Volume 29 - 2021
- Volume 28 - 2020
- Volume 27 - 2020
- Volume 26 - 2019
- Volume 25 - 2019
- Volume 24 - 2018
- Volume 23 - 2018
- Volume 22 - 2017
- Volume 21 - 2017
- Volume 20 - 2016
- Volume 19 - 2016
- Volume 18 - 2015
- Volume 17 - 2015
- Volume 16 - 2014
- Volume 15 - 2014
- Volume 14 - 2013
- Volume 13 - 2013
- Volume 12 - 2012
- Volume 11 - 2012
- Volume 10 - 2011
- Volume 9 - 2011
- Volume 8 - 2010
- Volume 7 - 2010
- Volume 6 - 2009
- Volume 5 - 2009
- Volume 4 - 2008
- Volume 3 - 2008
- Volume 2 - 2007
- Volume 1 - 2006
Commun. Comput. Phys., 24 (2018), pp. 44-68.
Published online: 2018-03
Cited by
- BibTex
- RIS
- TXT
The minimum action method (MAM) is to calculate the most probable transition path in randomly perturbed stochastic dynamics, based on the idea of action minimization in the path space. The accuracy of the numerical path between different metastable states usually suffers from the "clustering problem" near fixed points. The adaptive minimum action method (aMAM) solves this problem by relocating image points equally along arc-length with the help of moving mesh strategy. However, when the time interval is large, the images on the path may still be locally trapped around the transition state in a tangle, due to the singularity of the relationship between arc-length and time at the transition state. Additionally, in most non-gradient dynamics, the tangent direction of the path is not continuous at the transition state so that a geometric corner forms, which brings extra challenges for the aMAM. In this note, we improve the aMAM by proposing a better monitor function that does not contain the numerical approximation of derivatives, and taking use of a generalized scheme of the Euler-Lagrange equation to solve the minimization problem, so that both the path-tangling problem and the non-smoothness in parametrizing the curve do not exist. To further improve the accuracy, we apply the Weighted Essentially non-oscillatory (WENO) method for the interpolation to achieve better performance. Numerical examples are presented to demonstrate the advantages of our new method.
}, issn = {1991-7120}, doi = {https://doi.org/10.4208/cicp.OA-2016-0230}, url = {http://global-sci.org/intro/article_detail/cicp/10927.html} }The minimum action method (MAM) is to calculate the most probable transition path in randomly perturbed stochastic dynamics, based on the idea of action minimization in the path space. The accuracy of the numerical path between different metastable states usually suffers from the "clustering problem" near fixed points. The adaptive minimum action method (aMAM) solves this problem by relocating image points equally along arc-length with the help of moving mesh strategy. However, when the time interval is large, the images on the path may still be locally trapped around the transition state in a tangle, due to the singularity of the relationship between arc-length and time at the transition state. Additionally, in most non-gradient dynamics, the tangent direction of the path is not continuous at the transition state so that a geometric corner forms, which brings extra challenges for the aMAM. In this note, we improve the aMAM by proposing a better monitor function that does not contain the numerical approximation of derivatives, and taking use of a generalized scheme of the Euler-Lagrange equation to solve the minimization problem, so that both the path-tangling problem and the non-smoothness in parametrizing the curve do not exist. To further improve the accuracy, we apply the Weighted Essentially non-oscillatory (WENO) method for the interpolation to achieve better performance. Numerical examples are presented to demonstrate the advantages of our new method.