Symplectic integration approach for metastable systems

Abstract

Nonadiabatic behavior of metastable systems modeled by anharmonic Hamiltonians is reproduced by the Fokker-Planck and imaginary time Schrodinger equation scheme with subsequent symplectic integration. Example solutions capture ergodicity breaking, reassure the H-theorem of global stability [M. Shiino, Phys. Rev. A, Vol 36, pp. 2393-2411 (1987)], and reproduce spatially extended response under alternate source fields.