A novel spectral method for the semi-classical Schrödinger equation based on the Gaussian wave-packet transform

02/20/2021
by   Borui Miao, et al.
0

In this article, we develop and analyse a new spectral method to solve the semi-classical Schrödinger equation based on the Gaussian wave-packet transform (GWPT) and Hagedorn's semi-classical wave-packets (HWP). The GWPT equivalently recasts the highly oscillatory wave equation as a much less oscillatory one (the w equation) coupled with a set of ordinary differential equations governing the dynamics of the so-called GWPT parameters. The Hamiltonian of the w equation consists of a quadratic part and a small non-quadratic perturbation, which is of order 𝒪(√(ε)), where ε≪ 1 is the rescaled Planck's constant. By expanding the solution of the w equation as a superposition of Hagedorn's wave-packets, we construct a spectral method while the 𝒪(√(ε)) perturbation part is treated by the Galerkin approximation. This numerical implementation of the GWPT avoids imposing artificial boundary conditions and facilitates rigorous numerical analysis. For arbitrary dimensional cases, we establish how the error of solving the semi-classical Schrödinger equation with the GWPT is determined by the errors of solving the w equation and the GWPT parameters. We prove that this scheme has the spectral convergence with respect to the number of Hagedorn's wave-packets in one dimension. Extensive numerical tests are provided to demonstrate the properties of the proposed method.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
01/07/2021

Approximation of wave packets on the real line

In this paper we compare three different orthogonal systems in L_2(ℝ) wh...
research
10/07/2020

The Gaussian Wave Packet Transform via Quadrature Rules

We study variants of the Gaussian wave packet transform and investigate ...
research
11/09/2022

Analysis and Hermite spectral approximation of diffusive-viscous wave equations in unbounded domains arising in geophysics

The diffusive-viscous wave equation (DVWE) is widely used in seismic exp...
research
04/14/2023

Numerical approximation of the boundary control for the wave equation with a spectral collocation method

We propose a spectral collocation method to approximate the exact bounda...
research
07/15/2016

Spectral Echolocation via the Wave Embedding

Spectral embedding uses eigenfunctions of the discrete Laplacian on a we...
research
06/05/2022

Efficient Frozen Gaussian Sampling Algorithms for Nonadiabatic Quantum Dynamics at Metal Surfaces

In this article, we propose a Frozen Gaussian Sampling (FGS) algorithm f...
research
12/07/2021

Computing spectral properties of topological insulators without artificial truncation or supercell approximation

Topological insulators (TIs) are renowned for their remarkable electroni...

Please sign up or login with your details

Forgot password? Click here to reset