A ``converse stability condition is necessary for a compact higher order scheme on non-uniform meshes for the time-dependent Schrödinger equation

DOI10.1016/J.AML.2018.01.005MaRDI QIDQ1748595zbMATH OpenOpenAlexFDO

Authors Alexander Zlotnik, Raimondas Čiegis

Publication date 11 May 2018

Published in Applied Mathematics Letters (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/1707.09943

zbMATH Keywords

stabilitycompact schemenon-uniform spatial meshNumerov-Crank-Nicolson schemetime-dependent Schrödinger equation

Mathematics Subject Classification ID

Time-dependent Schrödinger equations and Dirac equations (35Q41) Finite difference methods for initial value and initial-boundary value problems involving PDEs (65M06) Stability and convergence of numerical methods for initial value and initial-boundary value problems involving PDEs (65M12) Error bounds for initial value and initial-boundary value problems involving PDEs (65M15)

Abstract: The stability bounds and error estimates for a compact higher order Numerov-Crank-Nicolson scheme on non-uniform space meshes for the 1D time-dependent Schr"odinger equation have been recently derived. This analysis has been done in $L^2$ and $H^1$ mesh norms and used the non-standard "converse" condition $h_{o}megaleq{c}_{0}au$, where $h_{o}mega$ is the mean space step, $au$ is the time step and $c_0>0$. Now we prove that such condition is necessary for some families of non-uniform meshes and any space norm. Also numerical results show unacceptably wrong behavior of numerical solutions (their dramatic mass non-conservation) when this condition is violated.

Recommendations

• The Numerov-Crank-Nicolson scheme on a non-uniform mesh for the time-dependent Schrödinger equation on the half-axis
• A compact higher-order finite-difference scheme for the wave equation can be strongly non-dissipative on non-uniform meshes
• Stability of a Numerov type finite-difference scheme with approximate transparent boundary conditions for the nonstationary Schrödinger equation on the half-axis
• On compact high order finite difference schemes for linear Schrödinger problem on non-uniform meshes
• High order compact finite difference schemes on nonuniform grids

Cites work

• scientific article; zbMATH DE number 3250158 (Why is no real title available?)
• A fast direct algorithm for implementing a high-order finite element method on rectangles as applied to boundary value problems for the Poisson equation
• A new Numerov-type method for the numerical solution of the Schrödinger equation
• Basic concepts in computational physics
• Finite element method with discrete transparent boundary conditions for the time-dependent 1D Schrödinger equation
• On a family of finite-difference schemes with approximate transparent boundary conditions for a generalized 1D Schrödinger equation
• On compact high order finite difference schemes for linear Schrödinger problem on non-uniform meshes
• The Numerov-Crank-Nicolson scheme on a non-uniform mesh for the time-dependent Schrödinger equation on the half-axis
• Variable mesh methods for the numerical solution of two-point singular perturbation problems

Cited in

(5)

• PRACTICAL ERROR ANALYSIS FOR THE THREE-LEVEL BILINEAR FEM AND FINITE-DIFFERENCE SCHEME FOR THE 1D WAVE EQUATION WITH NON-SMOOTH DATA
• A compact higher-order finite-difference scheme for the wave equation can be strongly non-dissipative on non-uniform meshes
• On compact 4th order finite-difference schemes for the wave equation
• On numerical simulation of electromagnetic field effects in the combustion process
• High order compact finite difference schemes on nonuniform grids

This page was built for publication: A ``converse stability condition is necessary for a compact higher order scheme on non-uniform meshes for the time-dependent Schrödinger equation

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q1748595)