Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

The Journal of Chemical Physics

Classical mechanics and the spreading of localized wave packets in condensed phase molecular systems

Michael F. Herman, D. F. Coker

Published Aug 1, 1999• 43 citations
Publisher

Abstract

The relationship between the diverging of classical trajectories in chaotic many-body systems, the spreading of quantum wave packets, and the validity and use of classical molecular dynamics is explored. This analysis, which is based on the semiclassical description of wave function propagation in terms of a weighted integration over a traveling fixed width coherent state basis, suggests that the exponential divergence of nearby classical trajectories in chaotic many-body systems should result in the rapid delocalization of an initially localized quantum wave packet describing the state of the system. Thus the justification for the use of classical molecular dynamics procedures for these supposedly classical systems cannot be based on the picture of the system wave function remaining localized as its center follows a nearly classical trajectory. The quantum evolution of the system density, on the other hand, requires two propagators, and each of these propagators is represented as an integration over trajectories in the semiclassical picture. The interference between the contributions from these two integrations over classical trajectories focuses the analysis on the most important points in this trajectory pair space, which are shown to occur when both trajectories in the pair are the same. Given reasonable assumptions for the initial density for a system that is expected to be well described by classical molecular dynamics, and given an appropriate choice for the width of the coherent state basis which is employed in the semiclassical description, it is shown that the semiclassical expressions for time dependent observables and correlation functions reduce the purely classical expressions, despite the fact that an initially localized wave packet would rapidly delocalize for the same system.

DOI
10.1063/1.479449

Related papers

Structures in classical phase space and quantum chaotic dynamics

B. S. Helmkamp, D. A. Browne

We compare the quantum and classical dynamics of a particle moving in a cosine potential while subject to a time-dependent force. We concentrate here on the behavior of an initially well-localized wave packet at times before the classically chaotic motion is fully developed. We find that the quantum and classical dynamics are indistinguishable well beyond the Ehrenfest time where the wave packet delocalizes. The quantum and classical descriptions first differ precisely when the classical probability density is folded in the vicinity of a hyperbolic fixed point. At this point, the wave function acquires a nodal structure which we show to be the result of a simple beating phenomenon between paths in the semiclassical propagator.

Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics•Mar 1, 1994•21 citations

Phase Space Structures in Chaotic Quantum Dynamics.

B. S. Helmkamp

I compare the quantum and classical dynamics of a particle moving in a cosine potential while subject to a time-dependent force. I concentrate here on the behavior of an initially well-localized wave packet at times before the classically chaotic motion is fully developed. I find that the quantum and classical dynamics are indistinguishable well beyond the Ehrenfest time where the wave packet delocalizes. The quantum and classical descriptions first differ precisely when the classical probability density is folded in the vicinity of a hyperbolic fixed point. At this point, the wave function acquires a nodal structure which I show to be the result of a simple beating phenomenon between paths in the semiclassical propagator. When the interaction of the classical manifold with the hyperbolic fixed point leads to escape from the remnant separatrix, representing the onset of classical chaos, the interference associated with the tendril leads to a long-lived accumulation of quantum amplitude on top of the cosine barrier. This effect also has a semiclassical interpretation, meaning that it arises from the interference between classical paths, but one must expand to second order in Planck's constant to describe the behavior correctly. Both the nodal structure and this pinning of amplitude on the barrier are dynamic mechanisms for the quantum inhibition of mixing. I then couple the system to a bath of harmonic oscillators in order to study the effect of the environment on these mechanisms. Although an oscillator bath environment brings dissipation as well as noise to the problem, the noise effect dominates for the high temperature, weak coupling regime I study. When there is sufficient noise to render the interfering classical paths indistinguishable, I find that the quantum interference gets erased. This dephasing occurs at very early times, long before there is appreciable dissipation of energy to the environment. Consequently, one can argue that the presence of an environment, even if its effect would be negligible in a nonchaotic setting, allows for the possibility of quantum mixing.

Civil War Book Review•Jan 1, 1995•0 citations

Time-dependent approach to semiclassical dynamics

Eric J. Heller

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.

The Journal of Chemical Physics•Feb 15, 1975•1,619 citations

Semiclassical dynamics with quantum trajectories: Formulation and comparison with the semiclassical initial value representation propagator

Sophya V. Garashchuk, Vitaly A. Rassolov

We present a time-dependent semiclassical method based on quantum trajectories. Quantum-mechanical effects are described via the quantum potential computed from the wave function density approximated as a linear combination of Gaussian fitting functions. The number of the fitting functions determines the accuracy of the approximate quantum potential (AQP). One Gaussian fit reproduces time-evolution of a Gaussian wave packet in a parabolic potential. The limit of the large number of fitting Gaussians and trajectories gives the full quantum-mechanical result. The method is systematically improvable from classical to fully quantum. The fitting procedure is implemented as a gradient minimization. We also compare AQP method to the widely used semiclassical propagator of Herman and Kluk by computing energy-resolved transmission probabilities for the Eckart barrier from the wave packet time-correlation functions. We find the results obtained with the Herman–Kluk propagator to be essentially equivalent to those of AQP method with a one-Gaussian density fit for several barrier widths.

The Journal of Chemical Physics•Jan 30, 2003•54 citations

Quantum manifestations of classical stochasticity. II. Dynamics of wave packets of bound states

Yitzhak Weissman, Joshua Jortner

In this paper we consider some of the consequences of classical stochasticity phenomena in the quantum mechanical Henon–Heiles and Barbanis systems. To explore classical–quantum analogies we have introduced a quantum mechanical phase space based on the coherent state representation. Quantum Poincaré maps (QPMs) were constructed from contour plots of the stationary phase-space densities. We have observed a correlation between the topological features of a QPM at an eigenenergy E and the sensitivity ‖dE/dε‖ of that energy to the strength ε of the nonlinear coupling. States with extreme values of ‖dE/dε‖, corresponding to both high and to low values of ‖dE/dε‖, show regular regions in the quantum phase space, while states with intermediate values of ‖dE/dε‖ span large parts of the quantum phase-space regions. These general features of the QPMs apply for almost each energy multiplet, and no qualitative change in the character of the QPMs was observed above the classical critical energy Ec for the stochastic transition. To investigate the quantum analog of classical trajectories, we have studied the dynamics of initially coherent Gaussian wave packets. We have discovered two limiting types of time evolution of intially coherent wave packets, which exhibit quasiperiodic time evolution or rapid dephasing. Quasiperiodic quantum dynamics is analogous to quasiperiodic classical trajectory dynamics, while rapid dephasing of the wave packet reflects the effects of spreading of the wave packet or/and stochastic classical dynamics. We were able to establish a correspondence between the topology of the quantum phase-space density and the wave packet dynamics. Quasiperiodic time evolution and rapid dephasing are exhibited by wave packets initially located in the regular and in the irregular regions of the quantum phase space, respectively.

The Journal of Chemical Physics•Aug 1, 1982•65 citations

Semiclassical approximations based on complex trajectories

A. D. Ribeiro, Marcus A. M. de Aguiar, Michel Baranger

The semiclassical limit of the coherent state propagator $⟨{\mathbf{z}}^{\ensuremath{''}}\ensuremath{\mid}{e}^{\ensuremath{-}i\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{H}T∕\ensuremath{\hbar}}\ensuremath{\mid}{\mathbf{z}}^{\ensuremath{'}}⟩$ involves complex classical trajectories of the Hamiltonian $\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{H}(\mathbf{u},\mathbf{v})=⟨\mathbf{v}\ensuremath{\mid}\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{H}\ensuremath{\mid}\mathbf{u}⟩$ satisfying $\mathbf{u}(0)={\mathbf{z}}^{\ensuremath{'}}$ and $\mathbf{v}(T)={\mathbf{z}}^{\ensuremath{''}}*$. In this work we study mostly the case ${\mathbf{z}}^{\ensuremath{'}}={\mathbf{z}}^{\ensuremath{''}}$. The propagator is then the return probability amplitude of a wave packet. We show that a plot of the exact return probability brings out the quantal images of the classical periodic orbits. Then we compare the exact return probability with its semiclassical approximation for a soft chaotic system with two degrees of freedom. We find two situations where classical trajectories satisfying the correct boundary conditions must be excluded from the semiclassical formula. The first occurs when the contribution of the trajectory to the propagator becomes exponentially large as $\ensuremath{\hbar}$ goes to zero. The second occurs when the contributing trajectories undergo bifurcations. Close to the bifurcation the semiclassical formula diverges. More interestingly, in the example studied, after the bifurcation, where more than one trajectory satisfying the boundary conditions exist, only one of them in fact contributes to the semiclassical formula, a phenomenon closely related to Stokes lines. When the contributions of these trajectories are filtered out, the semiclassical results show excellent agreement with the exact calculations.

Physical Review E•Jun 4, 2004•48 citations
The relationship between the diverging of classical trajectories in chaotic many-body systems, the spreading of quantum wave packets, and the validity and use of classical molecular dynamics is explored. This analysis, which is based on the semiclassical description of wave function propagation in terms of a weighted integration over a traveling fixed width coherent state basis, suggests that the exponential divergence of nearby classical trajectories in chaotic many-body systems should result in the rapid delocalization of an initially localized quantum wave packet describing the state of the system. Thus the justification for the use of classical molecular dynamics procedures for these supposedly classical systems cannot be based on the picture of the system wave function remaining localized as its center follows a nearly classical trajectory. The quantum evolution of the system density, on the other hand, requires two propagators, and each of these propagators is represented as an integration over trajectories in the semiclassical picture. The interference between the contributions from these two integrations over classical trajectories focuses the analysis on the most important points in this trajectory pair space, which are shown to occur when both trajectories in the pair are the same. Given reasonable assumptions for the initial density for a system that is expected to be well described by classical molecular dynamics, and given an appropriate choice for the width of the coherent state basis which is employed in the semiclassical description, it is shown that the semiclassical expressions for time dependent observables and correlation functions reduce the purely classical expressions, despite the fact that an initially localized wave packet would rapidly delocalize for the same system.
We compare the quantum and classical dynamics of a particle moving in a cosine potential while subject to a time-dependent force. We concentrate here on the behavior of an initially well-localized wave packet at times before the classically chaotic motion is fully developed. We find that the quantum and classical dynamics are indistinguishable well beyond the Ehrenfest time where the wave packet delocalizes. The quantum and classical descriptions first differ precisely when the classical probability density is folded in the vicinity of a hyperbolic fixed point. At this point, the wave function acquires a nodal structure which we show to be the result of a simple beating phenomenon between paths in the semiclassical propagator.
I compare the quantum and classical dynamics of a particle moving in a cosine potential while subject to a time-dependent force. I concentrate here on the behavior of an initially well-localized wave packet at times before the classically chaotic motion is fully developed. I find that the quantum and classical dynamics are indistinguishable well beyond the Ehrenfest time where the wave packet delocalizes. The quantum and classical descriptions first differ precisely when the classical probability density is folded in the vicinity of a hyperbolic fixed point. At this point, the wave function acquires a nodal structure which I show to be the result of a simple beating phenomenon between paths in the semiclassical propagator. When the interaction of the classical manifold with the hyperbolic fixed point leads to escape from the remnant separatrix, representing the onset of classical chaos, the interference associated with the tendril leads to a long-lived accumulation of quantum amplitude on top of the cosine barrier. This effect also has a semiclassical interpretation, meaning that it arises from the interference between classical paths, but one must expand to second order in Planck's constant to describe the behavior correctly. Both the nodal structure and this pinning of amplitude on the barrier are dynamic mechanisms for the quantum inhibition of mixing. I then couple the system to a bath of harmonic oscillators in order to study the effect of the environment on these mechanisms. Although an oscillator bath environment brings dissipation as well as noise to the problem, the noise effect dominates for the high temperature, weak coupling regime I study. When there is sufficient noise to render the interfering classical paths indistinguishable, I find that the quantum interference gets erased. This dephasing occurs at very early times, long before there is appreciable dissipation of energy to the environment. Consequently, one can argue that the presence of an environment, even if its effect would be negligible in a nonchaotic setting, allows for the possibility of quantum mixing.
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.
We present a time-dependent semiclassical method based on quantum trajectories. Quantum-mechanical effects are described via the quantum potential computed from the wave function density approximated as a linear combination of Gaussian fitting functions. The number of the fitting functions determines the accuracy of the approximate quantum potential (AQP). One Gaussian fit reproduces time-evolution of a Gaussian wave packet in a parabolic potential. The limit of the large number of fitting Gaussians and trajectories gives the full quantum-mechanical result. The method is systematically improvable from classical to fully quantum. The fitting procedure is implemented as a gradient minimization. We also compare AQP method to the widely used semiclassical propagator of Herman and Kluk by computing energy-resolved transmission probabilities for the Eckart barrier from the wave packet time-correlation functions. We find the results obtained with the Herman–Kluk propagator to be essentially equivalent to those of AQP method with a one-Gaussian density fit for several barrier widths.
In this paper we consider some of the consequences of classical stochasticity phenomena in the quantum mechanical Henon–Heiles and Barbanis systems. To explore classical–quantum analogies we have introduced a quantum mechanical phase space based on the coherent state representation. Quantum Poincaré maps (QPMs) were constructed from contour plots of the stationary phase-space densities. We have observed a correlation between the topological features of a QPM at an eigenenergy E and the sensitivity ‖dE/dε‖ of that energy to the strength ε of the nonlinear coupling. States with extreme values of ‖dE/dε‖, corresponding to both high and to low values of ‖dE/dε‖, show regular regions in the quantum phase space, while states with intermediate values of ‖dE/dε‖ span large parts of the quantum phase-space regions. These general features of the QPMs apply for almost each energy multiplet, and no qualitative change in the character of the QPMs was observed above the classical critical energy Ec for the stochastic transition. To investigate the quantum analog of classical trajectories, we have studied the dynamics of initially coherent Gaussian wave packets. We have discovered two limiting types of time evolution of intially coherent wave packets, which exhibit quasiperiodic time evolution or rapid dephasing. Quasiperiodic quantum dynamics is analogous to quasiperiodic classical trajectory dynamics, while rapid dephasing of the wave packet reflects the effects of spreading of the wave packet or/and stochastic classical dynamics. We were able to establish a correspondence between the topology of the quantum phase-space density and the wave packet dynamics. Quasiperiodic time evolution and rapid dephasing are exhibited by wave packets initially located in the regular and in the irregular regions of the quantum phase space, respectively.
The semiclassical limit of the coherent state propagator $⟨{\mathbf{z}}^{\ensuremath{''}}\ensuremath{\mid}{e}^{\ensuremath{-}i\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{H}T∕\ensuremath{\hbar}}\ensuremath{\mid}{\mathbf{z}}^{\ensuremath{'}}⟩$ involves complex classical trajectories of the Hamiltonian $\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{H}(\mathbf{u},\mathbf{v})=⟨\mathbf{v}\ensuremath{\mid}\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{H}\ensuremath{\mid}\mathbf{u}⟩$ satisfying $\mathbf{u}(0)={\mathbf{z}}^{\ensuremath{'}}$ and $\mathbf{v}(T)={\mathbf{z}}^{\ensuremath{''}}*$. In this work we study mostly the case ${\mathbf{z}}^{\ensuremath{'}}={\mathbf{z}}^{\ensuremath{''}}$. The propagator is then the return probability amplitude of a wave packet. We show that a plot of the exact return probability brings out the quantal images of the classical periodic orbits. Then we compare the exact return probability with its semiclassical approximation for a soft chaotic system with two degrees of freedom. We find two situations where classical trajectories satisfying the correct boundary conditions must be excluded from the semiclassical formula. The first occurs when the contribution of the trajectory to the propagator becomes exponentially large as $\ensuremath{\hbar}$ goes to zero. The second occurs when the contributing trajectories undergo bifurcations. Close to the bifurcation the semiclassical formula diverges. More interestingly, in the example studied, after the bifurcation, where more than one trajectory satisfying the boundary conditions exist, only one of them in fact contributes to the semiclassical formula, a phenomenon closely related to Stokes lines. When the contributions of these trajectories are filtered out, the semiclassical results show excellent agreement with the exact calculations.