Numerical Stability of Explicit Runge-Kutta Finite-Difference Schemes for the Nonlinear Schrödinger Equation
arXiv:1107.4810
Abstract
Linearized numerical stability bounds for solving the nonlinear time-dependent Schrödinger equation (NLSE) using explicit finite-differencing are shown. The bounds are computed for the fourth-order Runge-Kutta scheme in time and both second-order and fourth-order central differencing in space. Results are given for Dirichlet, modulus-squared Dirichlet, Laplacian-zero, and periodic boundary conditions for one, two, and three dimensions. Our approach is to use standard Runge-Kutta linear stability theory, treating the nonlinearity of the NLSE as a constant. The required bounds on the eigenvalues of the scheme matrices are found analytically when possible, and otherwise estimated using the Gershgorin circle theorem.
21 pages, 6 figures
References in corpus (4)
- Nonlinear Waves in Bose-Einstein Condensates: Physical Relevance and Mathematical Techniques
- A Modulus-Squared Dirichlet Boundary Condition for Time-Dependent Complex Partial Differential Equations and its Application to the Nonlinear Schrödinger Equation
- NLSEmagic: Nonlinear Schrödinger Equation Multidimensional Matlab-based GPU-accelerated Integrators using Compact High-order Schemes
- A Two-Step High-Order Compact Scheme for the Laplacian Operator and its Implementation in an Explicit Method for Integrating the Nonlinear Schrödinger Equation