Symplectic integration approach for metastable systems

Nonadiabatic behavior of metastable systems modeled by anharmonic Hamiltonians is reproduced by the Fokker-Planck and imaginary time Schrödinger equation scheme with subsequent symplectic integration. Example solutions capture ergodicity breaking, reassure the H-theorem of global stability [M. Shiino, Phys. Rev. A <Emphasis Type="Bold">36, 2393 (1987)], and reproduce spatially extended response under alternate source fields. Copyright EDP Sciences/Società Italiana di Fisica/Springer-Verlag 2006