Monte Carlo approach to the decay rate of a metastable system with an arbitrarily shaped barrier.

A path integral Monte Carlo method based on the fast-Fourier transform technique combined with the important sampling method is proposed to calculate the decay rate of a metastable quantum system with an arbitrary shape of a potential barrier. The contribution of all fluctuation actions is included which can be used to check the accuracy of the usual steepest-descent approximation, namely, the perturbation expansion of potential. The analytical approximation is found to produce the decay rate of a particle in a cubic potential being about 20% larger than the Monte Carlo data at the crossover temperature. This disagreement increases with increasing complexity of the potential shape. We also demonstrate via Langevin simulation that the postsaddle potential influences strongly upon the classical escape rate. [ABSTRACT FROM AUTHOR]