**Page description appears here**

“Solving the nonlinear Schrödinger equation using exponential integrators”

Authors: Håvard Berland, Brynjulf Owren and Bård Skaflestad,
Affiliation: NTNU
Reference: 2006, Vol 27, No 4, pp. 201-218.

     Valid XHTML 1.0 Strict


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 PDF (220 Kb)        DOI: 10.4173/mic.2006.4.1

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


References:
[1] BERLAND, H., OWREN, B. 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. WRIGHT, W. (2005). Expint - A Matlab package for exponential integrators, Tech. Rep. 4/05, Department of Mathematical Sciences, NTNU, Norway. http://www.math.ntnu.no/preprint/, To appear in ACM TOMS 2006.
[3] BRIDGES, T. J. 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 Schrödinger equations, no, 26 in Textos de Metodos Mathemáticos. Universidade Federal do Rio de Janeiro, Instituto de Matema´tica, third ed.
[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. 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. 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. 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. 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. 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. LUBICH, C. (2000). Error bounds for exponential operator splittings, BIT, 40, pp. 735 - 744 doi:10.1023/A:1022396519656
[12] KASSAM, A.-K. 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. WRIGHT, W. M. (2005). A review of exponential integrators for semilinear problems, Tech. Rep. 2/05, Department of Mathematical Sciences, NTNU, Norway. 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, Springer, Berlin, pp. 214 - 219 doi:10.1007/BFb0060031
[17] OSTERMANN, A., THALHAMMER, M. 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åvard and Owren, Brynjulf and Skaflestad, Bå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}
};

News

Aug 2014: For the 3rd year in a row MIC's impact factor increases. It is now 0.826.


Dec 2013: New database-driven web-design enabling extended statistics. Article number 500 is published and MIC reaches 1000 DOI Forward Links.


Jan 2012: Follow MIC on your smartphone by using the RSS feed.

Smartphone


July 2011: MIC passes 1000 ISI Web of Science citations.


Mar 2010: MIC is now indexed by DOAJ and has received the Sparc Seal seal for open access journals.


Dec 2009: A MIC group is created at LinkedIn and Twitter.


Oct 2009: MIC is now fully updated in ISI Web of Knowledge.