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

Exact wave packet propagation using time-dependent basis sets

Jasna Kučar, Hans‐Dieter Meyer

Published May 15, 1989• 46 citations
Publisher

Abstract

The use of a time-dependent basis set for solving the time-dependent Schrödinger equation is proposed. The basis set is a generalization of the harmonic oscillator functions. The proposed method is similar to the approaches formulated by Lee and Heller and by Coalson and Karplus. However, in the present method the time dependence of the basis functions is determined differently. Numerical studies on the quartic oscillator and on the Morse potential are carried out to investigate the performance of the various methods.

DOI
10.1063/1.456410

Related papers

Time-dependent treatment of scattering. II. Novel integral equation approach to quantum wave packets

Omar A. Sharafeddin, Richard Judson, Donald Jack Kouri +1 more

The time-dependent form of the Lippmann–Schwinger integral equation is used as the basis for a novel wave-packet propagation scheme. The method has the advantage over a previous integral equation treatment in that it does not require extensive matrix inversions involving the potential. This feature will be important when applications are made to systems where in some degrees of freedom the potential is expressed in a basis expansion. As was the case for the previous treatment, noniterated and iterated versions of the equations are given; the iterated equations, which are much simpler in the present new scheme than in the old, eliminate a matrix inversion that is required for solving the earlier noniterated equations. In the present noniterated equations, the matrix to be inverted is a function of the kinetic energy operator and thus is diagonal in a Bessel function basis set (or a sine basis set, if the centrifugal potential operator is incorporated into an effective potential). Transition amplitudes for various orbital angular momentum quantum numbers can be obtained from: (1) Fourier transform of the amplitude density from the time to the energy domain, and (2) direct analysis of the scattered wave packet. The approach is illustrated by an application to a standard potential scattering model problem.

The Journal of Chemical Physics•Oct 15, 1990•21 citations

Analytic time-dependent solutions of the one-dimensional Schrödinger equation

W. van Dijk, F. Masafumi Toyama, Sjirk Jan Prins +1 more

We show that one can obtain analytic solutions of the time-dependent Schrödinger equation that are more complex than the well-known oscillating coherent wave packet. Such Hermite-Gaussian or initially square wave packets exist for a free particle or for one subject to the harmonic oscillator potential. In either case, the Hermite-Gaussian packets retain their nodal structure even after long times. There is a single class of exact solutions for the system with oscillator constant K > 0, K = 0, or K < 0, leading to wave functions for the harmonic oscillator, the free particle, and the inverted oscillator, respectively.

American Journal of Physics•Sep 24, 2014•16 citations

Two-dimensional oscillator in time-dependent fields: Comparison of some exact and approximate calculations

Ochbadrakh Chuluunbaatar, Alexander A. Gusev, Sergue Il'ich Vinitsky +3 more

Operator-difference multilayer schemes for solving the time-dependent Schrödinger equation up to sixth order of accuracy in the time step are presented. Reduced schemes for solving a set of coupled time-dependent Schrödinger equations with respect to the hyper-radial variable are devised using expansion of a wave packet over the set of appropriate basis angular functions. Further discretization of the resulting problem is realized by means of the finite-element method. The convergence of the expansion with respect to the number of basis functions and the efficiency of the numerical schemes are demonstrated in the exactly solvable model of an electric-field-driven two-dimensional oscillator (or a charged particle in a constant uniform magnetic field), in which we explicitly observed an effect of the periodical focusing and defocusing of the probability density flux.

Physical Review E•Jul 22, 2008•2 citations

A strategy for time dependent quantum mechanical calculations using a Gaussian wave packet representation of the wave function

Shin‐ichi Sawada, Robert Heather, Bret Jackson +1 more

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 Schrödinger 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.

The Journal of Chemical Physics•Sep 15, 1985•205 citations

Quantum dynamics via mobile basis sets: The Dirac variational principle

David Hsu, D. F. Coker

We propose a variation on the wave-packet interpretation of quantum mechanics introduced by Heller, here using the Dirac variational principle to optimize the time evolution of all the parameters that enter our basis set. We solve the general one-dimensional problem, and test our method on the Morse oscillator and the symmetric quartic double well.

The Journal of Chemical Physics•Mar 15, 1992•20 citations

Time Dependent Schrödinger Equation

Pawel Scharoch, Maciej P. Polak, Radosław Szymon +1 more

This project looks into the time evolution of a wave function within a two-dimensional quantum well. We start by solving the time-dependent Schrödinger equation for stationary states in a quantum well. Next, we express any wave function as a linear combination of stationary states, allowing us to understand their time evolution. Two methods are presented: one relies on decomposing the wave function into a basis of stationary states and the other on discretisation of the time-dependent Schrödinger equation, incorporating three-point formulas for derivatives. These approaches necessitate confronting intricate boundary conditions and require maintaining energy conservation for numerical accuracy. We further demonstrate the methods using a wave packet, revealing fundamental phenomena in quantum physics. Our results demonstrate the utility of these methods in understanding quantum systems, despite the challenges in determining stationary states for a given potential. This study enhances our comprehension of the dynamics of quantum states in constrained systems, essential for fields like quantum computing and nanotechnology.

Cambridge University Press eBooks•Jan 31, 2024•0 citations
The use of a time-dependent basis set for solving the time-dependent Schrödinger equation is proposed. The basis set is a generalization of the harmonic oscillator functions. The proposed method is similar to the approaches formulated by Lee and Heller and by Coalson and Karplus. However, in the present method the time dependence of the basis functions is determined differently. Numerical studies on the quartic oscillator and on the Morse potential are carried out to investigate the performance of the various methods.
The time-dependent form of the Lippmann–Schwinger integral equation is used as the basis for a novel wave-packet propagation scheme. The method has the advantage over a previous integral equation treatment in that it does not require extensive matrix inversions involving the potential. This feature will be important when applications are made to systems where in some degrees of freedom the potential is expressed in a basis expansion. As was the case for the previous treatment, noniterated and iterated versions of the equations are given; the iterated equations, which are much simpler in the present new scheme than in the old, eliminate a matrix inversion that is required for solving the earlier noniterated equations. In the present noniterated equations, the matrix to be inverted is a function of the kinetic energy operator and thus is diagonal in a Bessel function basis set (or a sine basis set, if the centrifugal potential operator is incorporated into an effective potential). Transition amplitudes for various orbital angular momentum quantum numbers can be obtained from: (1) Fourier transform of the amplitude density from the time to the energy domain, and (2) direct analysis of the scattered wave packet. The approach is illustrated by an application to a standard potential scattering model problem.
We show that one can obtain analytic solutions of the time-dependent Schrödinger equation that are more complex than the well-known oscillating coherent wave packet. Such Hermite-Gaussian or initially square wave packets exist for a free particle or for one subject to the harmonic oscillator potential. In either case, the Hermite-Gaussian packets retain their nodal structure even after long times. There is a single class of exact solutions for the system with oscillator constant K > 0, K = 0, or K < 0, leading to wave functions for the harmonic oscillator, the free particle, and the inverted oscillator, respectively.
Operator-difference multilayer schemes for solving the time-dependent Schrödinger equation up to sixth order of accuracy in the time step are presented. Reduced schemes for solving a set of coupled time-dependent Schrödinger equations with respect to the hyper-radial variable are devised using expansion of a wave packet over the set of appropriate basis angular functions. Further discretization of the resulting problem is realized by means of the finite-element method. The convergence of the expansion with respect to the number of basis functions and the efficiency of the numerical schemes are demonstrated in the exactly solvable model of an electric-field-driven two-dimensional oscillator (or a charged particle in a constant uniform magnetic field), in which we explicitly observed an effect of the periodical focusing and defocusing of the probability density flux.
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 Schrödinger 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.
This project looks into the time evolution of a wave function within a two-dimensional quantum well. We start by solving the time-dependent Schrödinger equation for stationary states in a quantum well. Next, we express any wave function as a linear combination of stationary states, allowing us to understand their time evolution. Two methods are presented: one relies on decomposing the wave function into a basis of stationary states and the other on discretisation of the time-dependent Schrödinger equation, incorporating three-point formulas for derivatives. These approaches necessitate confronting intricate boundary conditions and require maintaining energy conservation for numerical accuracy. We further demonstrate the methods using a wave packet, revealing fundamental phenomena in quantum physics. Our results demonstrate the utility of these methods in understanding quantum systems, despite the challenges in determining stationary states for a given potential. This study enhances our comprehension of the dynamics of quantum states in constrained systems, essential for fields like quantum computing and nanotechnology.