Semiclassical wavefunction propagation for chaotic systems suffer from numerical difficulties due to the chaotic nature of classical trajectories, resulting in reduced accuracy for standard initial value representation (IVR) methods. We compare four recent IVR methods developed to overcome these difficulties (Herman-Kluk with trajectory truncation; stationary-phase Herman-Kluk (SPHK); cellularized frozen Gaussian approximation; stationary-phase Monte Carlo) by computing the Franck-Condon spectrum of the 2-dimensional Henon-Heiles and quartic oscillator systems. The SPHK is found to be the most successful of the four methods. The SPHK is then used to determine the spectrum for collinear CO2 photodissociation.
Bibliographical noteFunding Information:
The authors would like to thank Prof. Ken Kay (Bar-Ilan University) for a helpful communication regarding symmetry adaption for the quartic oscillator potential. This work was supported by the US Office of Naval Research and by the Natural Sciences and Engineering Research Council of Canada.