We present a simple method of calculation of the stability (monodromy) matrix that enters the widely used semiclassical propagator of Herman and Kluk and almost all other semiclassical propagators. The method is based on the unitarity of classical propagation and does not involve any approximations. The number of auxiliary differential equations per trajectory scales linearly rather than quadratically with the system size. Just the first derivatives of the potential surface are needed. The method is illustrated on the collinear H3 system.
Published in The Journal of Chemical Physics, Volume 113, Issue 21, 2000, pages 9390-9392.
Copyright 2000 American Institute of Physics. This article may be downloaded for personal use only. Any other use requires prior permission of the author and the American Institute of Physics.
The following article appeared in
Garashchuk, S. & Light, J. C. (2000). Simplified calculation of the stability matrix for semiclassical propagation. The Journal of Chemical Physics, 113(21), 9390-9392. http://dx.doi.org/10.1063/1.1321032
and may be found at