Renormalization of the frozen Gaussian approximation to the quantum propagator

The frozen Gaussian approximation to the quantum propagator may be a viable method for obtaining "on the fly" quantum dynamical information on systems with many degrees of freedom. However, it has two severe limitations, it rapidly loses normalization and one needs to know the Gaussian averaged potential, hence it is not a purely local theory in the force field. These limitations are in principle remedied by using the Herman-Kluk (HK) form for the semiclassical propagator. The HK propagator approximately conserves unitarity for relatively long times and depends only locally on the bare potential and its second derivatives. However, the HK propagator involves a much more expensive computation due to the need for evaluating the monodromy matrix elements. In this paper, we (a) derive a new formula for the normalization integral based on a prefactor free HK propagator which is amenable to "on the fly" computations; (b) show that a frozen Gaussian version of the normalization integral is not readily computable "on the fly"; (c) provide a new insight into how the HK prefactor leads to approximate unitarity; and (d) how one may construct a prefactor free approximation which combines the advantages of the frozen Gaussian and the HK propagators. The theoretical developments are backed by numerical examples on a Morse oscillator and a quartic double well potential.