Computing the chemical reaction path with a ray-based fast marching technique for solving the Hamilton-Jacobi equation in a general coordinate system

This article presents a ray-based fast marching approach for solving the static Hamilton-Jacobi equation. The approach is very general and can be used for both orthogonal and non-orthogonal coordinate system. The method is unconditionally stable, algorithmatically simple and highly accurate. As an application, we use the method to compute different types of reaction path. Specifically, we consider the path for which the change in action or time is less than that of all other conceivable paths connecting two states. Such reaction paths are efficiently evaluated by back-tracing on the least-action or least-time surfaces. The method is illustrated by applying it to the collinear reactions, F + H2 →HF + H and HF + H→H + FH.