We develop methodology for performing time dependent quantum mechanical calculations by representing the wave function as a sum of Gaussian wave packets (GWP) each characterized by a set of parameters such as width, position, momentum, and phase. The problem of computing the time evolution of the wave function is thus reduced to that of finding the time evolution of the parameters in the Gaussians. This parameter motion is determined by minimizing the error made by replacing the exact wave function in the time dependent Schrodinger equation with its Gaussian representation approximant. This leads to first order differential equations for the time dependence of the parameters, and those describing the packet position and the momentum of each packet have some resemblance with the classical equations of motion. The paper studies numerically the strategy needed to achieve the best GWP representation of time dependent processes. The issues discussed are the representation of the initial wave function, the numerical stability and the solution of the differential equations giving the evolution of the parameters, and the analysis of the final wave function. Extensive comparisons are made with an approximate method which assumes that the Gaussians are independent and their width is smaller than the length scale over which the potential changes. This approximation greatly simplifies the calculations and has the advantage of a greater resemblance to classical mechanics, thus being more intuitive. We find, however, that its range of applications is limited to problems involving localized degrees of freedom that participate in the dynamic process for a very short time. Finally we give particular attention to the notion that the GWP representation of the wave function reduces the dynamics of one quantum degree of freedom to that of a set of pseudoparticles (each represented by one packet) moving according to a pseudoclassical (i.e., classical‐like) mechanics whose ‘‘phase space’’ is described by a position and momentum as well as a complex phase and width.
[1]
H. Metiu,et al.
An examination of the use of wave packets for the calculation of atom diffraction by surfaces
,
1985
.
[2]
R. T. Skodje,et al.
Localized Gaussian wave packet methods for inelastic collisions involving anharmonic oscillators
,
1984
.
[3]
Gerhard Drolshagen,et al.
A time dependent wave packet approach to three‐dimensional gas–surface scattering
,
1983
.
[4]
E. Heller,et al.
Angular momentum projection operators and molecular bound states
,
1983
.
[5]
J. L. Kinsey,et al.
Spectroscopic characterization of repulsive potential energy surfaces: fluorescence spectrum of ozone
,
1982
.
[6]
E. Heller,et al.
Polyatomic Raman scattering for general harmonic potentials
,
1982
.
[7]
E. Heller.
The semiclassical way to molecular spectroscopy
,
1981
.
[8]
Eric J. Heller,et al.
Frozen Gaussians: A very simple semiclassical approximation
,
1981
.
[9]
Michael J Davis,et al.
Semiclassical Gaussian basis set method for molecular vibrational wave functions
,
1979
.
[10]
E. Heller,et al.
Time dependent formulation of polyatomic photofragmentation: Application to H3+
,
1978
.
[11]
E. Heller.
Phase space interpretation of semiclassical theory
,
1977
.
[12]
E. Heller.
Classical S‐matrix limit of wave packet dynamics
,
1976
.
[13]
E. Heller.
Time‐dependent approach to semiclassical dynamics
,
1975
.
[14]
Edmund Taylor Whittaker,et al.
The Calculus of Observations
,
1969,
The Mathematical Gazette.
[15]
P. Ehrenfest.
Bemerkung über die angenäherte Gültigkeit der klassischen Mechanik innerhalb der Quantenmechanik
,
1927
.
[16]
E. Heller.
Time dependent variational approach to semiclassical dynamics
,
1976
.
[17]
A. Messiah.
Quantum Mechanics
,
1961
.