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

Two-dimensional quantum propagation using wavelets in space and time

Douglas K. Sparks, Bruce R. Johnson

Published Sep 20, 2006• 8 citations
Publisher

Abstract

A recent method for solving the time-dependent Schrodinger equation has been developed using expansions in compact-support wavelet bases in both space and time [H. Wang et al., J. Chem. Phys. 121, 7647 (2004)]. This method represents an exact quantum mixed time-frequency approach, with special initial temporal wavelets used to solve the initial value problem. The present work is a first extension of the method to multiple spatial dimensions applied to a simple two-dimensional (2D) coupled anharmonic oscillator problem. A wavelet-discretized version of norm preservation for time-independent Hamiltonians discovered in the earlier one-dimensional investigation is verified to hold as well in 2D and, by implication, in higher numbers of spatial dimensions. The wavelet bases are not restricted to rectangular domains, a fact which is exploited here in a 2D adaptive version of the algorithm.

DOI
10.1063/1.2338318

Related papers

Multiscale quantum propagation using compact-support wavelets in space and time

Haixiang Wang, Ramiro Acevedo, Heather Mollé +3 more

Orthogonal compact-support Daubechies wavelets are employed as bases for both space and time variables in the solution of the time-dependent Schrodinger equation. Initial value conditions are enforced using special early-time wavelets analogous to edge wavelets used in boundary-value problems. It is shown that the quantum equations may be solved directly and accurately in the discrete wavelet representation, an important finding for the eventual goal of highly adaptive multiresolution Schrodinger equation solvers. While the temporal part of the basis is not sharp in either time or frequency, the Chebyshev method used for pure time-domain propagations is adapted to use in the mixed domain and is able to take advantage of Hamiltonian matrix sparseness. The orthogonal separation into different time scales is determined theoretically to persist throughout the evolution and is demonstrated numerically in a partially adaptive treatment of scattering from an asymmetric Eckart barrier.

The Journal of Chemical Physics•Oct 12, 2004•7 citations

Numerical method for solving Schrödinger equations using wavelet basis functions

Myoung Won Cho

Abstract This paper presents a numerical method for solving Schrödinger equations with arbitrary potentials with the adoption of the qunatum diagonalization scheme and basis functions in the discrete wavelet transform. The quantum diagonalization scheme, developed to solve the Hubbard model and so on, is known to allow the diagonalization of a Hamiltonian matrix of a much large size because it does not require to maintain the entire matrix in memory. And, the adoption of wavelet basis functions enables representing eigenstates efficiently using a number of bases much smaller than the number of grid points to integral and reducing the computation time to perform the diagonalization process. With the help of another additional techniques, the method allows to find solutions of a Schrödinger equation much exactly and efficiently. The validity of the method and its efficiency depending on several conditions are proven and surveyed by applying it to several problems for which the exact solutions are known.

•Jun 30, 2025•0 citations

Application of moment wavelet transform to quantum mechanics

Asmerom Ogbazghi

In this work we reconstruct Quantum Mechanical wave functions for some confining potentials, using the moment-wavelet method of Handy and Murenzi. This method consists in transforming the Schrodinger equation into an equivalent continuous wavelet transform (CWT) representation involving scaled and translated moments, µ1/ α,b (p)=∫χpe-Q\\( Ψ(x +b), where ϐ -Q becomes the mother wavelet. The discrete bound states are determined through a multiscale process involving the integration of a finite number of coupled linear first order differential equation in the moments µ1/ α,b (p). The underlying initial value problem depends on moment quantization methods to determine the infinite scale (a = ∞) moments and energy. Using this method we calculate the energies and wavefunctions for the first two quantum states of quartic and dectic anharmonic oscillator potentials, V(x) = mχ2 + gχ4, V(χ) = χ 2+χ10 respectively.

•Jul 1, 1998•0 citations

Hamilton‐jacobi dynamics for the solution of time dependent quantum problems I. Formalism and wave packet propagation in one dimension

Ersin Yurtsever, Jürgen Brickmann

Abstract Two methods for the numerical integration of the time‐dependent Schrödinger equation with given initial conditions (initial wave packet) are presented. The first method (method A) is based on the Schrödinger representation of the quantum‐dynamical system while the second one (method B) is based upon the intermediate representation. In both cases the quantum dynamical equation is transformed into a system of Hamilton‐Jacobi type equations of motion as occurring in multi particle classical dynamics, i.e. standard molecular dynamics techniques can be applied for the integration. The dynamics of a minimum uncertainty Gaussian wave packet in a strongly anharmonic oscillator is taken as an example.

Berichte der Bunsengesellschaft für physikalische Chemie•Apr 1, 1994•4 citations

Sources of quantum waves

Andrés D. Baute, I. L. Egusquiza, Juan Gonzalo Muga

Due to the space and time dependence of the wavefunction in the time-dependent Schrödinger equation, different boundary conditions are possible. The equation is usually solved as an `initial-value problem', by fixing the value of the wavefunction in all space at a given instant. We compare this standard approach to `source boundary conditions' that fix the wave at all times in a given region, in particular at a point in one dimension. In contrast to the well known physical interpretation of the initial-value-problem approach, the interpretation of the source approach has remained unclear, since it introduces negative energy components, even for `free motion', and a time-dependent norm. This work provides physical meaning to the source method by finding the link with equivalent initial-value problems.

arXiv:quant-ph/0007066•Journal of Physics A Mathematical and General•May 9, 2001•29 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 recent method for solving the time-dependent Schrodinger equation has been developed using expansions in compact-support wavelet bases in both space and time [H. Wang et al., J. Chem. Phys. 121, 7647 (2004)]. This method represents an exact quantum mixed time-frequency approach, with special initial temporal wavelets used to solve the initial value problem. The present work is a first extension of the method to multiple spatial dimensions applied to a simple two-dimensional (2D) coupled anharmonic oscillator problem. A wavelet-discretized version of norm preservation for time-independent Hamiltonians discovered in the earlier one-dimensional investigation is verified to hold as well in 2D and, by implication, in higher numbers of spatial dimensions. The wavelet bases are not restricted to rectangular domains, a fact which is exploited here in a 2D adaptive version of the algorithm.
Orthogonal compact-support Daubechies wavelets are employed as bases for both space and time variables in the solution of the time-dependent Schrodinger equation. Initial value conditions are enforced using special early-time wavelets analogous to edge wavelets used in boundary-value problems. It is shown that the quantum equations may be solved directly and accurately in the discrete wavelet representation, an important finding for the eventual goal of highly adaptive multiresolution Schrodinger equation solvers. While the temporal part of the basis is not sharp in either time or frequency, the Chebyshev method used for pure time-domain propagations is adapted to use in the mixed domain and is able to take advantage of Hamiltonian matrix sparseness. The orthogonal separation into different time scales is determined theoretically to persist throughout the evolution and is demonstrated numerically in a partially adaptive treatment of scattering from an asymmetric Eckart barrier.
Abstract This paper presents a numerical method for solving Schrödinger equations with arbitrary potentials with the adoption of the qunatum diagonalization scheme and basis functions in the discrete wavelet transform. The quantum diagonalization scheme, developed to solve the Hubbard model and so on, is known to allow the diagonalization of a Hamiltonian matrix of a much large size because it does not require to maintain the entire matrix in memory. And, the adoption of wavelet basis functions enables representing eigenstates efficiently using a number of bases much smaller than the number of grid points to integral and reducing the computation time to perform the diagonalization process. With the help of another additional techniques, the method allows to find solutions of a Schrödinger equation much exactly and efficiently. The validity of the method and its efficiency depending on several conditions are proven and surveyed by applying it to several problems for which the exact solutions are known.
In this work we reconstruct Quantum Mechanical wave functions for some confining potentials, using the moment-wavelet method of Handy and Murenzi. This method consists in transforming the Schrodinger equation into an equivalent continuous wavelet transform (CWT) representation involving scaled and translated moments, µ1/ α,b (p)=∫χpe-Q\\( Ψ(x +b), where ϐ -Q becomes the mother wavelet. The discrete bound states are determined through a multiscale process involving the integration of a finite number of coupled linear first order differential equation in the moments µ1/ α,b (p). The underlying initial value problem depends on moment quantization methods to determine the infinite scale (a = ∞) moments and energy. Using this method we calculate the energies and wavefunctions for the first two quantum states of quartic and dectic anharmonic oscillator potentials, V(x) = mχ2 + gχ4, V(χ) = χ 2+χ10 respectively.
Abstract Two methods for the numerical integration of the time‐dependent Schrödinger equation with given initial conditions (initial wave packet) are presented. The first method (method A) is based on the Schrödinger representation of the quantum‐dynamical system while the second one (method B) is based upon the intermediate representation. In both cases the quantum dynamical equation is transformed into a system of Hamilton‐Jacobi type equations of motion as occurring in multi particle classical dynamics, i.e. standard molecular dynamics techniques can be applied for the integration. The dynamics of a minimum uncertainty Gaussian wave packet in a strongly anharmonic oscillator is taken as an example.
Due to the space and time dependence of the wavefunction in the time-dependent Schrödinger equation, different boundary conditions are possible. The equation is usually solved as an `initial-value problem', by fixing the value of the wavefunction in all space at a given instant. We compare this standard approach to `source boundary conditions' that fix the wave at all times in a given region, in particular at a point in one dimension. In contrast to the well known physical interpretation of the initial-value-problem approach, the interpretation of the source approach has remained unclear, since it introduces negative energy components, even for `free motion', and a time-dependent norm. This work provides physical meaning to the source method by finding the link with equivalent initial-value problems.
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.