A "converse" stability condition is necessary for a compact higher order scheme on non-uniform meshes for the time-dependent Schrödinger equation. (June 2018)

• Record Type: Journal Article
• Authors: Zlotnik, Alexander; Čiegis, Raimondas
• Is Part Of: Applied mathematics letters. Volume 80(2018)
• Page Start: 35
• Page End: 40
• Publication Date: 2018-06
• Subjects: Stability -- Non-uniform spatial mesh -- Compact scheme -- Numerov–Crank–Nicolson scheme -- Time-dependent Schrödinger equation; Applied mathematics -- Periodicals; 519.05
• DOI: 10.1016/j.aml.2018.01.005
• Languages: English
• ISSNs: 0893-9659

Abstract

Abstract: The stability bounds and error estimates for a compact higher order Numerov–Crank–Nicolson scheme on non-uniform spatial meshes for the 1D time-dependent Schrödinger 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 ω ≤ c 0 τ, where h ω is the mean spatial step, τ 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 spatial norm. Also computational results for zero and non-zero potentials show unacceptably wrong behavior of numerical solutions when τ decreases and this condition is violated.