Abstract
A second-order Crank-Nicolson finite difference method, integrating a fast approximation of an exact discrete absorbing boundary condition, is proposed for solving the one-dimensional Schrödinger equation in the whole space. The fast approximation is based on Gaussian quadrature approximation of the convolution coefficients in the discrete absorbing boundary conditions. It approximates the time convolution in the discrete absorbing boundary conditions by a system of O(log 2 N) ordinary differential equations at each time step, where N denotes the total number of time steps. Stability and an error estimate are presented for the numerical solutions given by the proposed fast algorithm. Numerical experiments are provided, which agree with the theoretical results and show the performance of the proposed numerical method.
Original language | English |
---|---|
Pages (from-to) | A4083-A4104 |
Number of pages | 22 |
Journal | SIAM Journal on Scientific Computing |
Volume | 40 |
Issue number | 6 |
DOIs | |
Publication status | Published - 13 Dec 2018 |
Keywords
- Absorbing boundary condition
- Error estimate
- Fast algorithm
- Gaussian quadrature
- Schrödinger equation
- Stability
ASJC Scopus subject areas
- Computational Mathematics
- Applied Mathematics