Abstract:The transparent boundary condition for the free Schrödinger equation on a rectangular computational domain requires implementation of an operator of the form $\sqrt{\partial_t-i\triangle_{\Gamma}}$ where $\triangle_{\Gamma}$ is the Laplace-Beltrami operator. It is known that this operator is nonlocal in time as well as space which poses a significant challenge in developing an efficient numerical method of solution. The computational complexity of the existing methods scale with the number of time-steps which can be attributed to the nonlocal nature of the boundary operator. In this work, we report an effectively local approximation for the boundary operator such that the resulting complexity remains independent of number of time-steps. At the heart of this algorithm is a Padé approximant based rational approximation of certain fractional operators that handles corners of the domain adequately. For the spatial discretization, we use a Legendre-Galerkin spectral method with a new boundary adapted basis which ensures that the resulting linear system is banded. A compatible boundary-lifting procedure is also presented which accommodates the segments as well as the corners on the boundary. The proposed novel scheme can be implemented within the framework of any one-step time marching schemes. In particular, we demonstrate these ideas for two one-step methods, namely, the backward-differentiation formula of order 1 (BDF1) and the trapezoidal rule (TR). For the sake of comparison, we also present a convolution quadrature based scheme conforming to the one-step methods which is computationally expensive but serves as a golden standard. Finally, several numerical tests are presented to demonstrate the effectiveness of our novel method as well as to verify the order of convergence empirically.

Spectral Collocation Technique for Absorbing Boundary Conditions with Increasingly High Order Approximation

Efficient Chebyshev Spectral Method for Solving Linear Elliptic Pdes Using Quasi-Inverse Technique

High-order Symplectic FDTD Scheme for Solving a Time-Dependent Schrödinger Equation

Adaptive absorbing boundary conditions for Schrödinger-type equations: Application to nonlinear and multi-dimensional problems

Adaptive absorbing boundary conditions for Schrodinger-type equations: application to nonlinear and multi-dimensional problems

An Efficient Second-Order Finite Difference Method for the One-Dimensional Schrödinger Equation with Absorbing Boundary Conditions

Fast and Stable Evaluation of the Exact Absorbing Boundary Condition for the Semi-Discrete Linear Schrödinger Equation in Unbounded Domains

Efficient Spectral and Spectral Element Methods for Eigenvalue Problems of Schrödinger Equations with an Inverse Square Potential

A unified approach to split absorbing boundary conditions for nonlinear Schrödinger equations

Adaptive Artificial Boundary Condition for the Two-Level Schrödinger Equation with Conical Crossings

Stability and Error Analysis for a Second-Order Fast Approximation of the One-dimensional Schrödinger Equation under Absorbing Boundary Conditions

Exponential Collocation Methods Based on Continuous Finite Element Approximations for Efficiently Solving the Cubic Schrödinger Equation

Almost Exact Boundary Condition for One-Dimensional Schrödinger Equations

Transparent boundary condition and its effectively local approximation for the Schrödinger equation on a rectangular computational domain

Domain truncation, absorbing boundary conditions, Schur complements, and Padé approximation

Nonreflecting Boundary Condition for the free Schrödinger equation for hyperrectangular computational domains

Local absorbing boundary conditions on fixed domains give order-one errors for high-frequency waves

Absorbing Boundary Conditions for Variable Potential Schrödinger Equations via Titchmarsh-Weyl Theory

Transparent boundary condition and its high frequency approximation for the Schrödinger equation on a rectangular computational domain

Exact Boundary Condition for Semi-discretized Schrödinger Equation and Heat Equation in a Rectangular Domain

EFFICIENT SCALING AND MOVING TECHNIQUES FOR SPECTRAL METHODS IN UNBOUNDED DOMAINS