It is shown that application of the time-dependent variational principle to trial wave functions of a particular form results in a description of the quantum system as a classical system consisting of two mutually interacting parts: the ordinary classical analog of the quantum system and a (generally finite) number of additional degrees of freedom. This holds out the prospect of simulating quantum dynamics by coupling the corresponding classical system to additional degrees of freedom, and it is pointed out that such calculations have, in effect, already been performed. Since the trial wave functions that qualify for this treatment may be arbitrarily accurate, there is no limitation to the accuracy of the resulting classical description. In a particular case, the canonical coordinates and momenta for all degrees of freedom of the overall classical system are identified and investigated. In general, certain parameters appearing in the trial wave functions treated in this work may become redundant at certain times, leading to singularities in the Hamiltonian equations of motion. It is shown, however, that in certain such cases, a classical formalism for the quantum dynamics remains valid if the Hamiltonian function is suitably modified. Calculations applying the present Hamiltonian formalism are carried out that demonstrate that it has better convergence properties than an alternative, nonvariational, Hamiltonian technique suggested by previous work. These numerical studies also suggest that the observed spreading of wave packets and development of complicated structure in wave functions may be interpreted classically in terms of energy transfer among the various degrees of freedom in the overall classical system.
It is shown that application of the time-dependent variational principle to trial wave functions of a particular form results in a description of the quantum system as a classical system consisting of two mutually interacting parts: the ordinary classical analog of the quantum system and a (generally finite) number of additional degrees of freedom. This holds out the prospect of simulating quantum dynamics by coupling the corresponding classical system to additional degrees of freedom, and it is pointed out that such calculations have, in effect, already been performed. Since the trial wave functions that qualify for this treatment may be arbitrarily accurate, there is no limitation to the accuracy of the resulting classical description. In a particular case, the canonical coordinates and momenta for all degrees of freedom of the overall classical system are identified and investigated. In general, certain parameters appearing in the trial wave functions treated in this work may become redundant at certain times, leading to singularities in the Hamiltonian equations of motion. It is shown, however, that in certain such cases, a classical formalism for the quantum dynamics remains valid if the Hamiltonian function is suitably modified. Calculations applying the present Hamiltonian formalism are carried out that demonstrate that it has better convergence properties than an alternative, nonvariational, Hamiltonian technique suggested by previous work. These numerical studies also suggest that the observed spreading of wave packets and development of complicated structure in wave functions may be interpreted classically in terms of energy transfer among the various degrees of freedom in the overall classical system.
An exact Hamiltonian formulation of nonrelativistic quantum mechanics that simplifies to ordinary classical mechanics as \ensuremath{\Elzxh}\ensuremath{\rightarrow}0 is presented. The theory allows one to treat a quantum system as the corresponding classical system mechanically connected to an infinite number of additional classical degrees of freedom. The onset of quantum effects such as tunneling and interference results from the mechanical interactions between the different parts of the overall infinite system. The formalism is applied to the one-dimensional motion of a particle in a symmetric double well. It is shown that the tunneling in that system involves multidimensional classical motion that circumvents the barrier.
In this paper we develop a new approach to semiclassical dynamics which exploits the fact that extended wavefunctions for heavy particles (or particles in harmonic potentials) may be decomposed into time−dependent wave packets, which spread minimally and which execute classical or nearly classical trajectories. A Gaussian form for the wave packets is assumed and equations of motion are derived for the parameters characterizing the Gaussians. If the potential (which may be nonseparable in many coordinates) is expanded in a Taylor series about the instantaneous center of the (many−particle) wave packet, and up to quadratic terms are kept, we find the classical parameters of the wave packet (positions, momenta) obey Hamilton’s equation of motion. Quantum parameters (wave packet spread, phase factor, correlation terms, etc.) obey similar first order quantum equations. The center of the wave packet is shown to acquire a phase equal to the action integral along the classical path. State−specific quantum information is obtained from the wave packet trajectories by use of the superposition principle and projection techniques. Successful numerical application is made to the collinear He + H2 system widely used as a test case. Classically forbidden transitions are accounted for and obtained in the same manner as the classically allowed transitions; turning points present no difficulties and flux is very nearly conserved.
In classical mechanics the state of a dynamical system at any particular time can be described by the values of a set of coordinates and their conjugate momenta, thus, if the system has n degrees of freedom, by 2n numbers. In quantum mechanics, on the other hand, we have to describe a state of the system by a wave function involving a set of coordinates, thus by a function of n variables. The quantum description is, therefore, much more complicated than the classical one. Let us consider, however, an ensemble of systems in Gibbs' sense, i.e. not a large number of actual systems which could, perhaps, interact with one another, but a large number of hypothetical systems which are introduced to describe one actual system of which our knowledge is only of a statistical nature. The basis of the quantum treatment of such an ensemble has been given by Neumann. The description obtained by Neumann of an ensemble on the quantum theory is no more complicated than the corresponding classical description. Thus the quantum theory, which appears to such a disadvantage on the score of complication when applied to individual systems, recovers its own when applied to an ensemble. It is the object of the present note to examine this question more closely and to show how complete the analogy is between the quantum and classical treatments of an ensemble.
The general, N state matrix representation of the time-dependent Schrödinger equation is equivalent to an N degree of freedom classical Hamiltonian system. We describe how classical mechanical methods and ideas can be applied towards understanding and modeling exact quantum dynamics. Two applications are presented. First, we illustrate how qualitative insights may be gained by treating the two state problem with a time-dependent coupling. In the case of periodic coupling, Poincaré surfaces of section are used to view the quantum dynamics, and features such as the Floquet modes take on interesting interpretations. The second application illustrates computational implications by showing how Liouville’s theorem, or more generally the symplectic nature of classical Hamiltonian dynamics, provides a new perspective for carrying out numerical wave packet propagation. We show how certain simple and explicit symplectic integrators can be used to numerically propagate wave packets. The approach is illustrated with an application to the problem of a diatomic molecule interacting with a laser, although it and related approaches may be useful for describing a variety of problems.
A path-integral formulation of quantum mechanics is investigated which is closely related to that of Feynman. It differs from Feynman's formulation in that it involves the Hamiltonian function of the canonically conjugate coordinates and momenta. The classical limit yields the variational principle: $\ensuremath{\delta}ff(p\ifmmode\ddagger\else\textdaggerdbl\fi{}\ensuremath{-}N)dt=0$. A path-integral formula is also obtained for the energy eigenstate projection operator associated with the time-independent Schr\"odinger equation. The classical limit of the projection operator formula yields a modified form of the well-known variational principle for the phase-space orbit of given energy. Relativistically covariant Hamiltonian variational principles are analyzed and lead naturally to a relativistic scalar wave equation which involves a proper time variable which is canonically conjugate to the mass in the same manner as the ordinary time variable is conjugate to the energy in nonrelativistic quantum theory.
Bose-Einstein condensates are bosonic many-body systems in which the quantum mechanical ground state is macroscopically occupied. If the interparticle interaction between the single bosons is attractive, this ground state is metastable and separated from the collapse by a barrier. A sufficient thermal excitation of the condensate can enable the crossing of this barrier and induce its collapse. In this thesis, the corresponding decay mechanism of the thermally induced coherent collapse of Bose-Einstein condensates with additional long-range interaction is investigated for the first time. For the description of the condensate’s dynamics a variational approach to the quantum wave function is used, which allows for the precise determination of the system’s transition state that governs the reaction dynamics. In order to quantitatively calculate the decay rate and lifetime of the condensate, respectively, a novel transition state theory for quantum wave packets is developed. This is based on a variational approach as well as normal form expansions of the dynamical equations and the energy functional in the vicinity of the transition state. The dynamical equations obtained from the variational approach induce a symplectic structure onto the space of variational parameters. Furthermore, the normal form expansions serve as a tool to locally construct canonical coordinates on this symplectic manifold. For this purpose, the dynamical equations as well as the energy functional are expanded and diagonalized. The eigenvectors of the linearized equations of motion thereby define a symplectic basis of the system. Afterwards, the transformations of the higher-order terms are carried out in two steps via nonlinear near-identity transformations: In the first step, the nonresonant terms of the expansions are eliminated and in the second step the remaining resonant terms are transformed in such a way that they fulfill canonical equations. The combination of these methods finally allows for the application of the well-established classical transition state theory. The procedure is valid for a wide variety of trial wave functions, for arbitrary dimensions of the variational space, in arbitrary order of the expansions, and it is independent of the precise form of the Hamilton operator. In addition, a procedure to handle bifurcations in the transition state is developed using a uniform rate formula in connection with an appropriate normal form of the potential energy surface. In contrast to the usual treatment which breaks down at the point of the bifurcation, the new method allows to cross the bifurcation smoothly, because it is valid on both sides of the bifurcation as well as arbitrarily close to it. At the same time, one can calculate the flux over the whole configuration of saddles which emerge in the bifurcation. The applicability of the newly developed transition state theory for wave packets is demonstrated for a cubic potential well. The thermal decay rates of a point particle and a frozen Gaussian wave packet, respectively, are calculated within the classical and quantum normal form approaches as well as the variational one. It is shown that the variational approach is capable to perfectly reproduce both the classical and quantum normal form results of the point particle in different limits. In the limit of a narrow trial wave function it agrees perfectly with the results using the classical normal form. The results with a broader trial wave function go even beyond the classical approach, i.e. they coincide with those of the quantum normal form. In the field of Bose-Einstein condensates with long-range interaction the application of transition state theory for quantum wave packets is extended to the field of nonlinear Schrödinger systems. For condensates with attractive 1/r-interaction, the thermal decay rates are calculated for different numbers of coupled Gaussian wave functions as well as different normal form orders. Significant corrections of the decay rate are obtained with increasing order of the normal forms. In the low-temperature regime convergence is observed in eighth order of the expansions. Finally, the thermally induced coherent collapse of dipolar condensates is investigated. These exhibit the feature of a biconcave ground state density distribution for certain values of the physical parameters. It is shown that such condensates reveal the astonishing effect of a symmetry-breaking collapse dynamics. This effect is directly related to the occurrence of the structured ground state, and it depends on the geometry of the external trapping potential. In case of a “conventional” density distribution, the collapse possesses an s-wave symmetry, while it has a d-wave symmetry in case of a blood-cell shaped ground state. The thermally induced coherent collapse plays a major role close to the stability threshold where the energy barrier is small and the decay rates are high. In case of Cr atoms, dipolar condensate have been first realized with, the region of interest is on the order of a Bohr radius above the critical scattering length. There, the lifetime of the condensate can be reduced to a millisecond. For the elements Dy and Er which have also been condensed in the laboratory, the effects become even stronger, i.e. the parameter region of importance becomes larger, the energy barriers become smaller, and the decay rates are higher.