The Herman−Kluk (HK) or coherent state initial value representation (IVR) of semiclassical (SC) theory expresses a time-evolving wave function as a sum (integral) over the initial conditions of classically evolved coherent states (aka minimum uncertainty wave packets). An exact expansion of the wave function can also be made in terms of these same coherent states, thinking of them simply as a (time-dependent) basis set. This paper shows the precise set of approximations that are necessary to degenerate the exact expansion into the SC-IVR result. As a byproduct, we also obtain a generalization of the HK “prefactor”.
The Herman−Kluk (HK) or coherent state initial value representation (IVR) of semiclassical (SC) theory expresses a time-evolving wave function as a sum (integral) over the initial conditions of classically evolved coherent states (aka minimum uncertainty wave packets). An exact expansion of the wave function can also be made in terms of these same coherent states, thinking of them simply as a (time-dependent) basis set. This paper shows the precise set of approximations that are necessary to degenerate the exact expansion into the SC-IVR result. As a byproduct, we also obtain a generalization of the HK “prefactor”.
The coherent path integral is generalized such that the identity operator represented in a complete (actually overcomplete) set of the coherent states with the ``time-variable'' exponents are inserted between two consecutive short-time propagators. Since such a complete set of any given exponent can constitute the identity operator, the exponent may be varied from time to time without loss of generality as long as it is set common to all the Gaussians. However, a finite truncation of the coherent state expansion should result in different values of the propagator depending on the choice of the exponents. Furthermore, approximation methodology to treat with the exact propagator can also depend on this choice, and thereby many different semiclassical propagators may emerge from these combinations. Indeed, we show that the well-known semiclassical propagators such as those of Van Vleck, Herman-Kluk, Heller's thawed Gaussian, and many others can be derived in a systematic manner, which enables one to comprehend these semiclassical propagators from a unified point of view. We are particularly interested in our generalized form of the Herman-Kluk propagator, since the relative accuracy of this propagator has been well established by Kay, and since, nevertheless, its derivation was not necessarily clear. Thus our generalized Herman-Kluk propagator replaces the classical Hamiltonian with a Gaussian averaged quantum Hamiltonian, generating non-Newtonian trajectories. We perform a numerical test to assess the quality of such a family of generalized Herman-Kluk propagators and find that the original Herman-Kluk gives an accurate result. The reason why this has come about is also discussed.
Semiclassical time-dependent solutions to the Schroedinger equation can be obtained by use of integral formulas known as initial value representations (IVRs). Over the last few decades IVR formulas have established a sucessful reputation, most notably in the field of physical chemistry, for it is often found that the quantum behavior of molecular systems is very well approximated by semiclassical wavefunctions. The IVR input consists in a swarm of independent classical trajectories, each of which obeys the dynamics dictated by the system’s unquantized Hamiltonian. In a recent work [Aguiar et al, Chem. Phys. 370 (2010)] an IVR expression based on complex classical trajectories was formulated. The complex nature of the classical trajectories in this particular formula is a consequence of employing minimum uncertainty wavepackets (coherent states) to represent the time evolution operator upon which the formula is constructed. Here, an improved and simpler version of this complex IVR formula is presented. In order to test its accuracy, we apply it to a simple model system.
Abstract Recent experimental progress in monitoring the dynamics of ultracold gases in optical lattices necessitates a quantitative theoretical description for a significant number of bosons. In the present paper, we investigate if time-dependent semiclassical initial value methodology, with propagators expressed as integrals over phase space and using classical trajectories, is suitable to describe interacting bosons, concentrating on a single mode. Despite the nonlinear contribution from the self-interaction, the corresponding classical dynamics allows for a largely analytical treatment of the semiclassical propagator. We find that application of the Herman–Kluk (HK) propagator conserves unitarity in the semiclassical limit, but a decay of the norm is seen for low particle numbers n . The frozen Gaussian approximation (FGA) (HK with unit prefactor) is explicitly shown to violate unitarity in the present system for non-vanishing interaction strength, even in the semiclassical limit. Furthermore, we show by evaluating the phase space integral in steepest descent approximation, that the HK propagator reproduces the exact spectrum correctly in the semiclassical limit ( ). An error is, however, incurred in next-to-next-to-leading order (small parameter ), as seen upon numerical evaluation of the integral and confirmed analytically by considering finite n corrections to the steepest descent calculations. The FGA, in contrast, is only accurate to lowest order, and an erroneous next-to-leading order term in the energy spectrum was found analytically. Finally, as an example application, we study the dynamics of wave packets by computing the time evolution of the Wigner function. While the often-used truncated Wigner approximation cannot capture any interferences present in the exact quantum mechanical solution (known analytically), we find that the HK approach, despite also using classical information only, reproduces the salient features of the exact solution correctly.
Using phase-space complexification, an Initial Value Representation (IVR) for the semiclassical propagator in position space is obtained as a composition of inverse Segal-Bargmann (S-B) transforms of the semiclassical coherent state propagator. The result is shown to be free of caustic singularities and identical to the Herman-Kluk (H-K) propagator, found ubiquitously in physical and chemical applications. We contrast the theoretical aspects of this particular IVR with the van Vleck-Gutzwiller (vV-G) propagator and one of its IVRs, often employed in order to evade the non-linear "root-search" for trajectories required by vV-G. We demonstrate that bypassing the root-search comes at the price of serious numerical instability for all IVRs except the H-K propagator. We back up our theoretical arguments with comprehensive numerical calculations performed using the homogeneous Kerr system, about which we also unveil some unexpected new phenomena, namely: (1) the observation of a clear mark of half the Ehrenfest's time in semiclassical dynamics; and (2) the accumulation of trajectories around caustics as a function of increasing time (dubbed "caustic stickiness"). We expect these phenomena to be more general than for the Kerr system alone.
The Herman–Kluk (HK) formula was shown in (Baranger et al J. Phys. A: Math. Gen. 34 7227 ) not to be a correct semiclassical limit of an exact quantum mechanical formula. Two previous attempts to derive it using semiclassical arguments contain serious errors. These statements are left totally untouched by Herman and Grossmann's comment. They argue that the formula which we found to be at fault is not the one that should be called the HK formula. However, the formula we criticized is definitely one of the steps, in fact the main step, in these two published derivations of the HK formula. Very recently, a new derivation was published by Miller. It is interesting, but it is not semiclassical.
The semiclassical Herman–Kluk initial value approximation of the quantum-mechanical propagator, which is heavily used in dynamical studies of atomic and molecular systems, takes into account nonlinearity of the underlying classical dynamics by using multiple classical initial value solutions. In applications to the propagation of Gaussian wavepackets non-Gaussian distortions in the course of time and therefore a realistic description of the quantum dynamics are possible. This is in contrast to the single-trajectory approximations investigated in a recent paper by Baranger et al , one of which was erroneously termed Herman–Kluk approximation.