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

We present a fast algorithm for the evaluation of exact, nonreflecting boundary conditions for the time-dependent Schrodinger 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 Kν(√is). © 2007 Wiley Periodicals, Inc.

[1]  L. Milne‐Thomson.  A Treatise on the Theory of Bessel Functions , 1945, Nature.

[2]  A. Arnold,et al.  On absorbing boundary conditions for quantum transport equations , 1994 .

[3]  Matthias Ehrhardt,et al.  Discrete transparent boundary conditions for the Schrödinger equation: fast calculation, approximation, and stability , 2003 .

[4]  Laurent Di Menza.  Transparent and absorbing boundary conditions for the schrödinger equation in a bounded domain , 1997 .

[5]  M. Abramowitz,et al.  Handbook of Mathematical Functions With Formulas, Graphs and Mathematical Tables (National Bureau of Standards Applied Mathematics Series No. 55) , 1965 .

[6]  Philip Rabinowitz,et al.  Methods of Numerical Integration , 1985 .

[7]  Matthias Ehrhardt,et al.  Discrete Transparent Boundary Conditions for Wide Angle Parabolic Equations in Underwater Acoustics , 1998 .

[8]  Leslie Greengard,et al.  A fast algorithm for particle simulations , 1987 .

[9]  Achim Schädle,et al.  Non-reflecting boundary conditions for the two-dimensional Schrödinger equation , 2002 .

[10]  Leslie Greengard,et al.  A Fast Adaptive Numerical Method for Stiff Two-Point Boundary Value Problems , 1997, SIAM J. Sci. Comput..

[11]  F. W. J. Olver,et al.  The asymptotic expansion of bessel functions of large order , 1954, Philosophical Transactions of the Royal Society of London. Series A, Mathematical and Physical Sciences.

[12]  T. Hagstrom.  Radiation boundary conditions for the numerical simulation of waves , 1999, Acta Numerica.

[13]  L. Register,et al.  Numerical simulation of mesoscopic systems with open boundaries using the multidimensional time‐dependent Schrödinger equation , 1991 .

[14]  Vladimir Rokhlin,et al.  Generalized Gaussian quadrature rules for systems of arbitrary functions , 1996 .

[15]  L. Greengard,et al.  A fast algorithm for the evaluation of heat potentials , 1990 .

[16]  Frank Schmidt,et al.  A comparison of transparent boundary conditions for the Fresnel equation , 2001 .

[17]  Leslie Greengard,et al.  Fast evaluation of nonreflecting boundary conditions for the Schrödinger equation in one dimension , 2004 .

[18]  Leslie Greengard,et al.  Rapid Evaluation of Nonreflecting Boundary Kernels for Time-Domain Wave Propagation , 2000, SIAM J. Numer. Anal..

[19]  Frank Schmidt,et al.  Discrete transparent boundary conditions for the numerical solution of Fresnel's equation , 1995 .

[20]  Christian Lubich,et al.  Fast Convolution for Nonreflecting Boundary Conditions , 2002, SIAM J. Sci. Comput..

[21]  Matthias Ehrhardt,et al.  Discrete transparent boundary conditions for the Schrödinger equation , 2001 .

[22]  Leslie Greengard,et al.  Spectral Approximation of the Free-Space Heat Kernel , 2000 .

[23]  Frank Schmidt,et al.  Discrete transparent boundary conditions for Schrödinger-type equations , 1997 .

[24]  R. Kosloff,et al.  Absorbing boundaries for wave propagation problems , 1986 .

[25]  D. Horner,et al.  Time-dependent approach to collisional ionization using exterior complex scaling , 2002 .