Volume 24, Issue 1
An Improved Adaptive Minimum Action Method for the Calculation of Transition Path in Non-Gradient Systems

Yiqun Sun & Xiang Zhou

Commun. Comput. Phys., 24 (2018), pp. 44-68.

Published online: 2018-03

Preview Full PDF 515 1936
Export citation
  • Abstract

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.

  • Keywords

Rare event, transition path, minimum action method, moving mesh, WENO.

  • AMS Subject Headings

82C26, 60H30, 34F05

  • Copyright

COPYRIGHT: © Global Science Press

  • Email address
  • BibTex
  • RIS
  • TXT
@Article{CiCP-24-44, author = {}, title = {An Improved Adaptive Minimum Action Method for the Calculation of Transition Path in Non-Gradient Systems}, journal = {Communications in Computational Physics}, year = {2018}, volume = {24}, number = {1}, pages = {44--68}, abstract = {

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} }
TY - JOUR T1 - An Improved Adaptive Minimum Action Method for the Calculation of Transition Path in Non-Gradient Systems JO - Communications in Computational Physics VL - 1 SP - 44 EP - 68 PY - 2018 DA - 2018/03 SN - 24 DO - http://doi.org/10.4208/cicp.OA-2016-0230 UR - https://global-sci.org/intro/article_detail/cicp/10927.html KW - Rare event, transition path, minimum action method, moving mesh, WENO. AB -

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.

Yiqun Sun & Xiang Zhou. (2020). An Improved Adaptive Minimum Action Method for the Calculation of Transition Path in Non-Gradient Systems. Communications in Computational Physics. 24 (1). 44-68. doi:10.4208/cicp.OA-2016-0230
Copy to clipboard
The citation has been copied to your clipboard