The time-dependent solution of the Fokker-Planck equation with bistable potentials is considered in terms of the eigenfunctions and eigenvalues of the linear Fokker-Planck operator. The Fokker-Planck equation is the high friction limit of the corresponding Kramers' equation. Two different potentials are considered defined with a constant diffusion coefficient, $\epsilon$, and position-dependent drift coefficients. The smallest nonzero eigenvalue of the Fokker-Planck operator, ${\lambda }_1$, provides the long-time rate coefficient for the transformation of the different species in the two stable states. A novel pseudospectral method with nonclassical polynomials is applied to this class of systems. The convergence of the eigenvalues and eigenfunctions of the Fokker-Planck operator versus the number of basis functions is studied and compared with previous results. The results are consistent with Kramers' theory, and a linear relationship between $ln{\lambda }_1$ and $1/\epsilon$ for sufficiently small $\epsilon$ values is verified. A comparison with analytic approximations to ${\lambda }_1$ is provided.

• Received 10 December 2018

DOI:https://doi.org/10.1103/PhysRevE.99.042101