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 H[sub 3] system.
Published in Journal of Chemical Physics, Volume 113, Issue 21, 2000, pages 9390-9392.
© 2000 by American Institute of Physics