Weak Ill-Posedness of Spatial Discretizations of Absorbing Boundary Conditions for Schrödinger-Type Equations

When we wish to solve numerically a differential problem defined on an infinite domain, it is necessary to consider a finite subdomain and to use artificial boundary conditions in such a way that the solutions in the finite subdomain approximate the original solution. These boundary conditions are called absorbing when small reflections to the interior domain are allowed. In this paper, we develop a general class of absorbing boundary conditions for Schrodinger-type equations by using rational approximations to the transparent boundary conditions. With this approach, previous absorbing boundary conditions in the literature are included in this class. We use the method of lines for the discretization of the initial boundary value problems obtained this way. We show that the ordinary differential systems that arise after the spatial discretization are weakly ill-posed, explaining a previous conjecture of Fevens and Jiang. The time discretization is carried out with A-stable Runge--Kutta methods, where the high order ones may be used to compensate for the possible troubles present in the problems semidiscretized in space.

[1]  Eric Dubach Artificial boundary conditions for diffusion equations: numerical study , 1996 .

[2]  Thomas Fevens,et al.  ABSORBING BOUNDARY CONDITIONS FOR THE SCHR ¨ ODINGER EQUATION , 1999 .

[3]  L. Trefethen,et al.  III‐posedness of absorbing boundary conditions for migration , 1988 .

[4]  David Yevick,et al.  Optimal absorbing boundary conditions , 1995 .

[5]  Frank Schmidt,et al.  An adaptive approach to the numerical solution of Fresnel's wave equation , 1993 .

[6]  A. V. Popov,et al.  Implementation of transparent boundaries for numerical solution of the Schrödinger equation , 1991 .

[7]  Laurence Halpern,et al.  Absorbing boundary conditions for diffusion equations , 1995 .

[8]  H. Kreiss Initial boundary value problems for hyperbolic systems , 1970 .

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

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

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

[12]  Isaías Alonso Mallo,et al.  Condiciones frontera transparentes y absorbentes para la ecuación de Schrödinger en una dimensión , 1999 .

[13]  Hong Jiang,et al.  Absorbing Boundary Conditions for the Schrödinger Equation , 1999, SIAM J. Sci. Comput..

[14]  L. Trefethen Group velocity interpretation of the stability theory of Gustafsson, Kreiss, and Sundström , 1983 .

[15]  B. Kågström Bounds and perturbation bounds for the matrix exponential , 1977 .

[16]  F. R. Gantmakher The Theory of Matrices , 1984 .

[17]  Kuska Absorbing boundary conditions for the Schrödinger equation on finite intervals. , 1992, Physical review. B, Condensed matter.

[18]  Shibata Absorbing boundary conditions for the finite-difference time-domain calculation of the one-dimensional Schrödinger equation. , 1991, Physical review. B, Condensed matter.

[19]  R. Higdon Absorbing boundary conditions for difference approximations to the multi-dimensional wave equation , 1986 .

[20]  A. Majda,et al.  Absorbing boundary conditions for the numerical simulation of waves , 1977 .

[21]  L. Trefethen,et al.  Well-Posedness of one-way wave equations and absorbing boundary conditions , 1986 .

[22]  L. Trefethen,et al.  Stability of the method of lines , 1992, Spectra and Pseudospectra.