Efficient representation of nonreflecting boundary conditions for the time-dependent Schrödinger equation in two dimensions

Document Type

Article

Publication Date

1-1-2008

Abstract

We present a fast algorithm for the evaluation of exact, nonreflecting boundary conditions for the time-dependent Schrödinger equation in two dimensions on the unit circle. After separation of variables, the exact outgoing condition for each Fourier mode contains a nonlocal term that is a convolution integral in time. The kernel for that convolution is the inverse Laplace transform of the logarithmic derivative of a modified Bessel function, and the convolution integral can be split into two parts: a local part and a history part, which can be treated separately. The local part is easily handled by an appropriate quadrature. For the history part, we show that the convolution kernel can be well approximated by a sum of exponentials. Once such a representation is available, the convolution integrals can be evaluated recursively, reducing the cost from O (N2) work to O(N), where N is the number of time steps. The main technical development lies in the uniform rational approximation of the logarithmic derivative of the modified Bessel function Kv(√is). © 2007 Wiley Periodicals, Inc.

Identifier

38049144584 (Scopus)

Publication Title

Communications on Pure and Applied Mathematics

External Full Text Location

https://doi.org/10.1002/cpa.20200

e-ISSN

00103640

ISSN

00103640

First Page

261

Last Page

288

Issue

2

Volume

61

This document is currently not available here.

Share

COinS