Global Recursion Polynomial Expansions of the Green’s Function and Time Evolution Operator for the Schrödinger Equation with Absorbing Boundary Conditions

Author(s):  
Vladimir A. Mandelshtam
2017 ◽  
Vol 22 (3) ◽  
pp. 408-423 ◽  
Author(s):  
Andrej Bugajev ◽  
Raimondas Čiegis ◽  
Rima Kriauzienė ◽  
Teresė Leonavičienė ◽  
Julius Žilinskas

A detailed analysis of absorbing boundary conditions for the linear Schrodinger equation is presented in this paper. It is focused on absorbing boundary conditions that are obtained by using rational functions to approximate the exact transparent boundary conditions. Different strategies are investigated for the optimal selection of the coefficients of these rational functions, including the Pade approximation, the L2 norm approximations of the Fourier symbol, L2 minimization of a reflection coefficient, and two error minimization techniques for the chosen benchmark problems with known exact solutions. The results of computational experiments are given and a detailed comparison of the efficiency of these techniques is presented.


2010 ◽  
Vol 21 (11) ◽  
pp. 1391-1406 ◽  
Author(s):  
F. L. DUBEIBE

In recent decades a lot of research has been done on the numerical solution of the time-dependent Schrödinger equation. On the one hand, some of the proposed numerical methods do not need any kind of matrix inversion, but source terms cannot be easily implemented into these schemes; on the other, some methods involving matrix inversion can implement source terms in a natural way, but are not easy to implement into some computational software programs widely used by non-experts in programming (e.g. Mathematica). We present a simple method to solve the time-dependent Schrödinger equation by using a standard Crank–Nicholson method together with a Cayley's form for the finite-difference representation of evolution operator. Here, such standard numerical scheme has been simplified by inverting analytically the matrix of the evolution operator in position representation. The analytical inversion of the N × N matrix let us easily and fully implement the numerical method, with or without source terms, into Mathematica or even into any numerical computing language or computational software used for scientific computing.


2012 ◽  
Vol 22 (10) ◽  
pp. 1250026 ◽  
Author(s):  
XAVIER ANTOINE ◽  
CHRISTOPHE BESSE ◽  
PAULINE KLEIN

The aim of this paper is to construct some classes of absorbing boundary conditions for the two-dimensional Schrödinger equation with a time and space varying exterior potential and for general convex smooth boundaries. The construction is based on asymptotics of the inhomogeneous pseudodifferential operators defining the related Dirichlet-to-Neumann operator. Furthermore, a priori estimates are developed for the truncated problems with various increasing order boundary conditions. The effective numerical approximation will be treated in a second paper.


Sign in / Sign up

Export Citation Format

Share Document