“Solving the nonlinear Schrödinger equation using exponential integrators” Authors: Håvard Berland, Brynjulf Owren and Bård Skaflestad Affiliation: NTNU (Department of Mathematical Sciences) Reference: 2006, Vol. 27, No. 4, pp. 201-218.

Keywords: Exponential integrators, nonlinear Schro¨dinger equation, regularity requirements, stiff and non-stiff order conditions

Abstract: Using the notion of integrating factors, Lawson developed a class of numerical methods for solving stiff systems of ordinary differential equations. However, the performance of these "Generalized Runge - Kutta processes" was demonstrably poorer when compared to the ETD schemes of Certaine and Nørsett, recently rediscovered by Cox and Matthews. The deficit is particularly pronounced when the schemes are applied to parabolic problems. In this paper we compare a fourth order Lawson scheme and a fourth order ETD scheme due to Cox and Matthews, using the nonlinear Schro¨dinger equation as the test problem. The primary testing parameters are degree of regularity of the potential function and the initial condition, and numerical performance is heavily dependent upon these values. The Lawson and ETD schemes exhibit significant performance differences in our tests, and we present some analysis on this.

PDF (220 Kb)        DOI: 10.4173/mic.2006.4.1

DOI forward links to this article:
[1] Håvard Berland Alvaro L. Islas Constance M. Schober, (2007), doi:10.1016/j.jcp.2006.11.030
[2] Morten Dahlby Brynjulf Owren, (2009), doi:10.1051/m2an/2009022
[3] Brett N. Ryland Robert I. Mclachlan Jason Frank, (2007), doi:10.1080/00207160701458633


References:
[1] BERLAND, H., OWREN, B. and SKAFLESTAD, B. (2005). B-series and order conditions for exponential integrators, SIAM J. Numer. Anal., 43 pp. 1715 - 1727, doi:10.1137/040612683
[2] BERLAND, H., SKAFLESTAD, B. and WRIGHT, W. (2005). Expint - A Matlab package for exponential integrators, Tech. Rep. 4/05, Department of Mathematical Sciences, NTNU, Norway, 2005. http://www.math.ntnu.no/preprint/. To appear in ACM TOMS 2006.
[3] BRIDGES, T. J. and REICH, S. (2001). Multi-symplectic integrators: numerical schemes for Hamiltonian PDEs that conserve symplecticity, Phys. Lett. A, 284 pp. 184 - 193, doi:10.1016/S0375-9601(01)00294-8
[4] CAZENAVE, T. (1996). An introduction to nonlinear Schro¨dinger equations, no. 26 in Textos de Metodos Mathema´ticos, Universidade Federal do Rio de Janeiro, Instituto de Matema´tica, third ed., 1996.
[5] CERTAINE, J. (1960). The solution of ordinary differential equations with large time constants, in Mathematical methods for digital computers, Wiley, New York, pp. 128 - 132.
[6] COX, S. M. and MATTHEWS, P. C. (2002). Exponential time differencing for stiff systems, J. Comput. Phys., 176, pp. 430 - 455, doi:10.1006/jcph.2002.6995
[7] EHLE, B. L. and LAWSON, J. D. (1975). Generalized Runge - Kutta processes for stiff initial-value problems, J. Inst. Maths. Applics., 16, pp. 11 - 21, doi:10.1093/imamat/16.1.11
[8] HOCHBRUCK, M. and OSTERMANN, A. (2005). Explicit exponential Runge - Kutta methods for semilinear parabolic problems, SIAM J. Numer. Anal., 43, pp. 1069 - 1090, doi:10.1137/040611434
[9] HOCHBRUCK, M., LUBICH, C. and SELHOFER, H. (1998). Exponential integrators for large systems of differential equations, SIAM J. Sci. Comput., 19, pp. 1552 - 1574, doi:10.1137/S1064827595295337
[10] ISLAS, A. L., KARPEEV, D. A. and SCHOBER, C. M. (2001). Geometric integrators for the nonlinear Schrödinger equation, J. of Comp. Phys., 173, pp. 116 - 148, doi:10.1006/jcph.2001.6854
[11] JAHNKE, T. and LUBICH, C. (2000). Error bounds for exponential operator splittings, BIT, 40, pp. 735 - 744, doi:10.1023/A:1022396519656
[12] KASSAM, A.-K. and TREFETHEN, L. N. (2005). Fourth-order time-stepping for stiff PDEs, SIAM J. Sci. Comput., 26, pp. 1214 - 1233 (electronic), doi:10.1137/S1064827502410633
[13] KROGSTAD, S. (2005). Generalized integrating factor methods for stiff PDEs, J. of Comp. Phys., 203, pp. 72 - 88, doi:10.1016/j.jcp.2004.08.006
[14] LAWSON, J. D. (1967). Generalized Runge - Kutta processes for stable systems with large Lipschitz constants, SIAM J. Numer. Anal., 4, pp. 372 - 380, doi:10.1137/0704033
[15] MINCHEV, B. and WRIGHT, W. M. (2005). A review of exponential integrators for semilinear problems, Tech. Rep. 2/05, Department of Mathematical Sciences, NTNU, Norway, 2005. http://www.math.ntnu.no/preprint/.
[16] NØRSETT, S. P. (1969). An A-stable modification of the Adams - Bashforth methods, in Conf. on Numerical Solution of Differential Equations (Dundee, Scotland, 1969), Springer, Berlin, pp. 214 - 219, doi:10.1007/BFb0060031
[17] OSTERMANN, A., THALHAMMER, M. and WRIGHT, W. M. (2005). A class of explicit exponential general linear methods. In preparation.


BibTeX:
@article{MIC-2006-4-1,
  title={{Solving the nonlinear Schrödinger equation using exponential integrators}},
  author={Berland, H{\aa}vard and Owren, Brynjulf and Skaflestad, B{\aa}rd},
  journal={Modeling, Identification and Control},
  volume={27},
  number={4},
  pages={201--218},
  year={2006},
  doi={10.4173/mic.2006.4.1},
  publisher={Norwegian Society of Automatic Control}
};