Improved error estimates for splitting methods applied to highly-oscillatory nonlinear Schrödinger equations
HTML articles powered by AMS MathViewer
- by Philippe Chartier, Florian Méhats, Mechthild Thalhammer and Yong Zhang PDF
- Math. Comp. 85 (2016), 2863-2885 Request permission
Abstract:
In this work, the error behavior of operator splitting methods is analyzed for highly-oscillatory differential equations. The scope of applications includes time-dependent nonlinear Schrödinger equations, where the evolution operator associated with the principal linear part is highly-oscillatory and periodic in time. In a first step, a known convergence result for the second-order Strang splitting method applied to the cubic Schrödinger equation is adapted to a wider class of nonlinearities. In a second step, the dependence of the global error on the decisive parameter $0 < \varepsilon <\!\!< 1$, defining the length of the period, is examined. The main result states that, compared to established error estimates, the Strang splitting method is more accurate by a factor $\varepsilon$, provided that the time stepsize is chosen as an integer fraction of the period. This improved error behavior over a time interval of fixed length, which is independent of the period, is due to an averaging effect. The extension of the convergence result to higher-order splitting methods and numerical illustrations complement the investigations.References
- Serge Alinhac and Patrick Gérard, Pseudo-differential operators and the Nash-Moser theorem, Graduate Studies in Mathematics, vol. 82, American Mathematical Society, Providence, RI, 2007. Translated from the 1991 French original by Stephen S. Wilson. MR 2304160, DOI 10.1090/gsm/082
- F. Castella, Ph. Chartier, F. Méhats, and A. Murua, Stroboscopic averaging for the nonlinear Schrödinger equation, Found. Comput. Math. 15 (2015), no. 2, 519–559. MR 3320932, DOI 10.1007/s10208-014-9235-7
- Ludwig Gauckler, Convergence of a split-step Hermite method for the Gross-Pitaevskii equation, IMA J. Numer. Anal. 31 (2011), no. 2, 396–415. MR 2813177, DOI 10.1093/imanum/drp041
- Christian Lubich, On splitting methods for Schrödinger-Poisson and cubic nonlinear Schrödinger equations, Math. Comp. 77 (2008), no. 264, 2141–2153. MR 2429878, DOI 10.1090/S0025-5718-08-02101-7
- Camil Muscalu and Wilhelm Schlag, Classical and multilinear harmonic analysis. Vol. II, Cambridge Studies in Advanced Mathematics, vol. 138, Cambridge University Press, Cambridge, 2013. MR 3052499
- J. A. C. Weideman and B. M. Herbst, Split-step methods for the solution of the nonlinear Schrödinger equation, SIAM J. Numer. Anal. 23 (1986), no. 3, 485–507. MR 842641, DOI 10.1137/0723033
- Haruo Yoshida, Construction of higher order symplectic integrators, Phys. Lett. A 150 (1990), no. 5-7, 262–268. MR 1078768, DOI 10.1016/0375-9601(90)90092-3
Additional Information
- Philippe Chartier
- Affiliation: INRIA-Rennes, IRMAR, ENS Cachan Bretagne, IPSO Project Team, Campus de Beaulieu, 35042 Rennes Cedex, France.
- MR Author ID: 335517
- Email: Philippe.Chartier@inria.fr
- Florian Méhats
- Affiliation: IRMAR, Université de Rennes 1, INRIA-Rennes, IPSO Project Team, Campus de Beaulieu, 35042 Rennes Cedex, France.
- MR Author ID: 601414
- Email: Florian.Mehats@univ-rennes1.fr
- Mechthild Thalhammer
- Affiliation: Leopold-Franzens Universität Innsbruck, Institut für Mathematik, Technikerstraße 13/VII, 6020 Innsbruck, Austria.
- MR Author ID: 661917
- Email: Mechthild.Thalhammer@uibk.ac.at
- Yong Zhang
- Affiliation: Wolfgang Pauli Institut c/o Universität Wien, Fakultät für Mathematik, Oskar-Morgenstern-Platz 1, 1090 Vienna, Austria.
- Email: Yong.Zhang@univie.ac.at
- Received by editor(s): March 7, 2014
- Received by editor(s) in revised form: March 30, 2015
- Published electronically: February 16, 2016
- Additional Notes: Corresponding author: Philippe Chartier
The authors acknowledge financial support by the Agence nationale de la recherche (ANR) within the project LODIQUAS ANR-11-IS01-0003 and the project Moonrise ANR-14-CE23-0007-01, by the Austrian Science Fund (FWF) under SFP Vienna Computational Materials Laboratory (ViCoM) and project P21620-N13, and by the Austrian Ministry of Science and Research via its grant for the WPI. The presented numerical results have been achieved by using the Vienna Scientific Cluster. - © Copyright 2016 American Mathematical Society
- Journal: Math. Comp. 85 (2016), 2863-2885
- MSC (2010): Primary 34K33, 37L05, 35Q55
- DOI: https://doi.org/10.1090/mcom/3088
- MathSciNet review: 3522973