High Order Stable Finite Difference Methods for the Schrödinger Equation

Authors:
A. Nissen;G. Kreiss;M. Gerritsen
Affiliations:
Department of Energy Resources Engineering, Stanford University, Stanford, USA;Department of Information Technology, Uppsala University, Uppsala, Sweden;Department of Energy Resources Engineering, Stanford University, Stanford, USA
Venue:
Journal of Scientific Computing
Year:
2013

Citing 10
Cited 1

Time-stable boundary conditions for finite-difference schemes solving hyperbolic systems: methodology and application to high-order compact schemes

Journal of Computational Physics
Summation by parts for finite difference approximations for d/dx

Journal of Computational Physics
A stable and conservative interface treatment of arbitrary spatial accuracy

Journal of Computational Physics
Summation by parts operators for finite difference approximations of second derivatives

Journal of Computational Physics
On the order of accuracy for difference approximations of initial-boundary value problems

Journal of Computational Physics
A class of energy stable, high-order finite-difference interface schemes suitable for adaptive mesh refinement of hyperbolic problems

Journal of Computational Physics
Nondissipative and energy-stable high-order finite-difference interface schemes for 2-D patch-refined grids

Journal of Computational Physics
Stable Boundary Treatment for the Wave Equation on Second-Order Form

Journal of Scientific Computing
Stable and Accurate Interpolation Operators for High-Order Multiblock Finite Difference Methods

SIAM Journal on Scientific Computing
Stability at Nonconforming Grid Interfaces for a High Order Discretization of the Schrödinger Equation

Journal of Scientific Computing

High-fidelity numerical solution of the time-dependent Dirac equation

Journal of Computational Physics

Quantified Score

Hi-index	0.01

Visualization

Abstract

In this paper we extend the Summation-by-parts-simultaneous approximation term (SBP-SAT) technique to the Schrödinger equation. Stability estimates are derived and the accuracy of numerical approximations of interior order 2m, m=1,2,3, are analyzed in the case of Dirichlet boundary conditions. We show that a boundary closure of the numerical approximations of order m lead to global accuracy of order m+2. The results are supported by numerical simulations.