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

Mathematical Models and Computer Simulations

Simulation of quantum dynamics of a wave packet in a semiclassical limit

K. S. Arakelov

Published Sep 16, 2011• 0 citations
Publisher

Abstract

A new algorithm for modeling quantum dynamics of wave packets in a semiclassical limit was formulated and justified. The algorithm is valid in the case of two-particle systems and in the case of a single particle in any arbitrary central potential. The profile of the wave packet depending on the time for a single wave packet evolving in a one-dimensional (1D) Morse potential was obtained by direct computer simulations using the new algorithm. The results were compared with the calculation which was previously performed using the Initial Value Representation (IVR) method.

DOI
10.1134/s2070048211050048

Related papers

A note on the application of the semiclassical wave packet dynamics

Michal Kolář, M. K. Ali, S. F. O’Shea

The optimal choice of the initial parameters of the Gaussian wave packets in the time dependent quantum mechanical calculations using a Gaussian wave packet representation of the wave function is discussed. The proposed procedure is illustrated with the simple case of one dimensional scattering and the propagation of the ground state of the Morse potential.

The Journal of Chemical Physics•Jan 15, 1989•2 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

Wave packet propagation across barriers by semiclassical initial value methods

Jakob Petersen, Kenneth G. Kay

Semiclassical initial value representation (IVR) formulas for the propagator have difficulty describing tunneling through barriers. A key reason is that these formulas do not automatically reduce, in the classical limit, to the version of the Van Vleck-Gutzwiller (VVG) propagator required to treat barrier tunneling, which involves trajectories that have complex initial conditions and that follow paths in complex time. In this work, a simple IVR expression, that has the correct tunneling form in the classical limit, is derived for the propagator in the case of one-dimensional barrier transmission. Similarly, an IVR formula, that reduces to the Generalized Gaussian Wave Packet Dynamics (GGWPD) expression [D. Huber, E. J. Heller, and R. Littlejohn, J. Chem. Phys. 89, 2003 (1988)] in the classical limit, is derived for the transmitted wave packet. Uniform semiclassical versions of the IVR formulas are presented and simplified expressions in terms of real trajectories and WKB penetration factors are described. Numerical tests show that the uniform IVR treatment gives good results for wave packet transmission through the Eckart and Gaussian barriers in all cases examined. In contrast, even when applied with the proper complex trajectories, the VVG and GGWPD treatments are inaccurate when the mean energy of the wave packet is near the classical transmission threshold. The IVR expressions for the propagator and wave packet are cast as contour integrals in the complex space of initial conditions and these are generalized to potentially allow treatment of a larger variety of systems. A steepest descent analysis of the contour integral formula for the wave packet in the present cases confirms its relationship to the GGWPD method, verifies its semiclassical validity, and explains results of numerical calculations.

The Journal of Chemical Physics•Jul 2, 2015•8 citations

Hamiltonian dynamics of semiclassical Gaussian wave packets in electromagnetic potentials

Nolan King, Tomoki Ohsawa

Abstract We extend our previous work on symplectic semiclassical Gaussian wave packet dynamics to incorporate electromagnetic interactions by including a vector potential. The main advantage of our formulation is that the equations of motion derived are naturally Hamiltonian. We obtain an asymptotic expansion of our equations in terms of and show that our equations have corrections to those presented by Zhou, whereas ours also recover the equations of Zhou in the case of a linear vector potential and quadratic scalar potential. One and two dimensional examples of a particle in a magnetic field are given and numerical solutions are presented and compared with the classical solutions and the expectation values of the corresponding observables as calculated by the Egorov or initial value representation (IVR) method. We numerically demonstrate that the correction terms improve the accuracy of the classical or Zhou’s equations for short times in the sense that our solutions converge to the expectation values calculated using the Egorov/IVR method faster than the classical solutions or those of Zhou as .

arXiv:1908.11454•Journal of Physics A Mathematical and Theoretical•Jan 27, 2020•3 citations

A new method for wave packet dynamics: Derivative propagation along quantum trajectories

Corey Jason Trahan, Keith H. Hughes, Róbert E. Wyatt

A new method is proposed for computing the time evolution of quantum mechanical wave packets. Equations of motion for the real-valued functions C and S in the complex action S=C(r,t)+iS(r,t)/ℏ, with ψ(r,t)=exp(S), involve gradients and curvatures of C and S. In previous implementations of the hydrodynamic formulation, various time-consuming fitting techniques of limited accuracy were used to evaluate these derivatives around each fluid element in an evolving ensemble. In this study, equations of motion are developed for the spatial derivatives themselves and a small set of these are integrated along quantum trajectories concurrently with the equations for C and S. Significantly, quantum effects can be included at various orders of approximation, no spatial fitting is involved, there are no basis set expansions, and single quantum trajectories (rather than correlated ensembles) may be propagated, one at a time. Excellent results are obtained when the derivative propagation method is applied to anharmonic potentials involving barrier transmission.

The Journal of Chemical Physics•May 22, 2003•85 citations

Numerical Study of Quantum Wave Packet in Meshless Particle Method

Fumiaki HIRONO, Misako IWASAWA, Satoru S. Kano +1 more

The particle method has no restriction on the particle arrangement, where the calculation is performed. The purpose is, by using this feature, to develop the method to analyze the dynamics of the electronic state. In order to describe the time evolution of the electronic state, we have developed a new method to solve the time-dependent wave equation, using the Bohmian that is very compatible with the particle method. In this form, there is also numerical instability in the region, where the value of the wave function is very small. We have applied our technique to the wave packet dynamics on harmonic potential and the analysis of interference by double slits as a simple system with analytical solution and easy to compare accuracy, and realized that both results are in good agreement.

Journal of Computer Chemistry Japan•Jan 1, 2019•2 citations
A new algorithm for modeling quantum dynamics of wave packets in a semiclassical limit was formulated and justified. The algorithm is valid in the case of two-particle systems and in the case of a single particle in any arbitrary central potential. The profile of the wave packet depending on the time for a single wave packet evolving in a one-dimensional (1D) Morse potential was obtained by direct computer simulations using the new algorithm. The results were compared with the calculation which was previously performed using the Initial Value Representation (IVR) method.
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.
Semiclassical initial value representation (IVR) formulas for the propagator have difficulty describing tunneling through barriers. A key reason is that these formulas do not automatically reduce, in the classical limit, to the version of the Van Vleck-Gutzwiller (VVG) propagator required to treat barrier tunneling, which involves trajectories that have complex initial conditions and that follow paths in complex time. In this work, a simple IVR expression, that has the correct tunneling form in the classical limit, is derived for the propagator in the case of one-dimensional barrier transmission. Similarly, an IVR formula, that reduces to the Generalized Gaussian Wave Packet Dynamics (GGWPD) expression [D. Huber, E. J. Heller, and R. Littlejohn, J. Chem. Phys. 89, 2003 (1988)] in the classical limit, is derived for the transmitted wave packet. Uniform semiclassical versions of the IVR formulas are presented and simplified expressions in terms of real trajectories and WKB penetration factors are described. Numerical tests show that the uniform IVR treatment gives good results for wave packet transmission through the Eckart and Gaussian barriers in all cases examined. In contrast, even when applied with the proper complex trajectories, the VVG and GGWPD treatments are inaccurate when the mean energy of the wave packet is near the classical transmission threshold. The IVR expressions for the propagator and wave packet are cast as contour integrals in the complex space of initial conditions and these are generalized to potentially allow treatment of a larger variety of systems. A steepest descent analysis of the contour integral formula for the wave packet in the present cases confirms its relationship to the GGWPD method, verifies its semiclassical validity, and explains results of numerical calculations.
Abstract We extend our previous work on symplectic semiclassical Gaussian wave packet dynamics to incorporate electromagnetic interactions by including a vector potential. The main advantage of our formulation is that the equations of motion derived are naturally Hamiltonian. We obtain an asymptotic expansion of our equations in terms of and show that our equations have corrections to those presented by Zhou, whereas ours also recover the equations of Zhou in the case of a linear vector potential and quadratic scalar potential. One and two dimensional examples of a particle in a magnetic field are given and numerical solutions are presented and compared with the classical solutions and the expectation values of the corresponding observables as calculated by the Egorov or initial value representation (IVR) method. We numerically demonstrate that the correction terms improve the accuracy of the classical or Zhou’s equations for short times in the sense that our solutions converge to the expectation values calculated using the Egorov/IVR method faster than the classical solutions or those of Zhou as .
A new method is proposed for computing the time evolution of quantum mechanical wave packets. Equations of motion for the real-valued functions C and S in the complex action S=C(r,t)+iS(r,t)/ℏ, with ψ(r,t)=exp(S), involve gradients and curvatures of C and S. In previous implementations of the hydrodynamic formulation, various time-consuming fitting techniques of limited accuracy were used to evaluate these derivatives around each fluid element in an evolving ensemble. In this study, equations of motion are developed for the spatial derivatives themselves and a small set of these are integrated along quantum trajectories concurrently with the equations for C and S. Significantly, quantum effects can be included at various orders of approximation, no spatial fitting is involved, there are no basis set expansions, and single quantum trajectories (rather than correlated ensembles) may be propagated, one at a time. Excellent results are obtained when the derivative propagation method is applied to anharmonic potentials involving barrier transmission.
The particle method has no restriction on the particle arrangement, where the calculation is performed. The purpose is, by using this feature, to develop the method to analyze the dynamics of the electronic state. In order to describe the time evolution of the electronic state, we have developed a new method to solve the time-dependent wave equation, using the Bohmian that is very compatible with the particle method. In this form, there is also numerical instability in the region, where the value of the wave function is very small. We have applied our technique to the wave packet dynamics on harmonic potential and the analysis of interference by double slits as a simple system with analytical solution and easy to compare accuracy, and realized that both results are in good agreement.