A High-Order Integrator for the Schrödinger Equation with Time-Dependent, Homogeneous Magnetic Field
The SMAI Journal of computational mathematics, Tome 6 (2020), pp. 253-271.

We construct a family of numerical methods for the Pauli equation of charged particles in a time-dependent, homogeneous magnetic field. These methods are described in a general setting comprising systems of multiple particles and extend the usual splitting and Fourier grid approach. The issue is that the magnetic field causes charged particles to rotate. The corresponding rotations of the wave function are highly incompatible with the Fourier grid approach used for the standard Schrödinger equation. Motivated by the theory of Lie algebras and their representations, our new approach approximates the exact flow map in terms of rotated potentials and rotated initial data, and thereby avoids this issue. Finally, we provide numerical examples to examine convergence and preservation of norm and energy.

Publié le :
DOI : 10.5802/smai-jcm.69
Classification : 81Q05, 65M12, 65M22
Mots clés : quantum mechanics, time-dependent Schrödinger equation, magnetic field, discrete Fourier transform, numerical method, splitting, convergence rate, Lie group, Lie algebra
Gradinaru, Vasile 1 ; Rietmann, Oliver 1

1 Seminar for Applied Mathematics, ETH Zurich, Switzerland
@article{SMAI-JCM_2020__6__253_0,
     author = {Gradinaru, Vasile and Rietmann, Oliver},
     title = {A {High-Order} {Integrator} for the {Schr\"odinger} {Equation} with {Time-Dependent,} {Homogeneous} {Magnetic} {Field}},
     journal = {The SMAI Journal of computational mathematics},
     pages = {253--271},
     publisher = {Soci\'et\'e de Math\'ematiques Appliqu\'ees et Industrielles},
     volume = {6},
     year = {2020},
     doi = {10.5802/smai-jcm.69},
     language = {en},
     url = {http://archive.numdam.org/articles/10.5802/smai-jcm.69/}
}
TY  - JOUR
AU  - Gradinaru, Vasile
AU  - Rietmann, Oliver
TI  - A High-Order Integrator for the Schrödinger Equation with Time-Dependent, Homogeneous Magnetic Field
JO  - The SMAI Journal of computational mathematics
PY  - 2020
SP  - 253
EP  - 271
VL  - 6
PB  - Société de Mathématiques Appliquées et Industrielles
UR  - http://archive.numdam.org/articles/10.5802/smai-jcm.69/
DO  - 10.5802/smai-jcm.69
LA  - en
ID  - SMAI-JCM_2020__6__253_0
ER  - 
%0 Journal Article
%A Gradinaru, Vasile
%A Rietmann, Oliver
%T A High-Order Integrator for the Schrödinger Equation with Time-Dependent, Homogeneous Magnetic Field
%J The SMAI Journal of computational mathematics
%D 2020
%P 253-271
%V 6
%I Société de Mathématiques Appliquées et Industrielles
%U http://archive.numdam.org/articles/10.5802/smai-jcm.69/
%R 10.5802/smai-jcm.69
%G en
%F SMAI-JCM_2020__6__253_0
Gradinaru, Vasile; Rietmann, Oliver. A High-Order Integrator for the Schrödinger Equation with Time-Dependent, Homogeneous Magnetic Field. The SMAI Journal of computational mathematics, Tome 6 (2020), pp. 253-271. doi : 10.5802/smai-jcm.69. http://archive.numdam.org/articles/10.5802/smai-jcm.69/

[1] Blanes, S.; Moan, P. C. Splitting methods for the time-dependent Schrödinger equation, Phys. Lett., A, Volume 265 (2000) no. 1, pp. 35-42 http://www.sciencedirect.com/science/article/pii/s037596019900866x | DOI | Zbl

[2] Blanes, S.; Moan, P. C. Practical symplectic partitioned Runge–Kutta and Runge–Kutta–Nyström methods, J. Comput. Appl. Math., Volume 142 (2002) no. 2, pp. 313-330 http://www.sciencedirect.com/science/article/pii/s0377042701004927 | DOI | Zbl

[3] Blanes, S.; Moan, P. C. Fourth- and sixth-order commutator-free Magnus integrators for linear and non-linear dynamical systems, Appl. Numer. Math., Volume 56 (2006) no. 12, pp. 1519-1537 http://www.sciencedirect.com/science/article/pii/s0168927405002163 | DOI | MR | Zbl

[4] Caliari, M.; Ostermann, A.; Piazzola, C. A splitting approach for the magnetic Schrödinger equation, J. Comput. Appl. Math., Volume 316 (2017), pp. 74-85 http://www.sciencedirect.com/science/article/pii/s0377042716304101 (Selected Papers from NUMDIFF-14) | DOI | Zbl

[5] Descombes, S.; Thalhammer, M. An exact local error representation of exponential operator splitting methods for evolutionary problems and applications to linear Schrödinger equations in the semi-classical regime, BIT Numer. Math., Volume 50 (2010) no. 4, pp. 729-749 | DOI | Zbl

[6] Dong, S.-H. Wave equations in higher dimensions, Springer, 2011 | MR

[7] Faou, E.; Gradinaru, V.; Lubich, C. Computing Semiclassical Quantum Dynamics with Hagedorn Wavepackets, SIAM J. Sci. Comput., Volume 31 (2009), pp. 3027-3041 | DOI | MR | Zbl

[8] Hairer, E.; Gheorghend Nørsett, S. P.; Wanner, G. Solving Ordinary Differential Equations I (2nd Revised. Ed.): Nonstiff Problems, Springer, 1993

[9] Hairer, E.; Lubich, C.; Wanner, G. Geometric numerical integration: Structure-preserving algorithms for ordinary differential equations, Springer Series in Computational Mathematics, 31, Springer, 2010 | Zbl

[10] Hall, B. C. Quantum theory for mathematicians, Graduate Texts in Mathematics, 267, Springer, 2013 | MR | Zbl

[11] Jin, S.; Zhou, Z. A semi-Lagrangian time splitting method for the Schrödinger equation with vector potentials, Commun. Inf. Syst., Volume 13 (2013), pp. 247-289 | DOI | Zbl

[12] Kahan, W.; Li, R.-C. Composition constants for raising the orders of unconventional schemes for ordinary differential equations, Math. Comp., Volume 66 (1997), pp. 1089-1099 | DOI | MR | Zbl

[13] Louck, J. D. Generalized orbital angular momentum and the n-fold degenerate quantum-mechanical oscillator: Part II. The n-fold degenerate oscillator, Journal of Molecular Spectroscopy, Volume 4 (1960) no. 1, pp. 298-333 http://www.sciencedirect.com/science/article/pii/0022285260900916 | DOI

[14] Ma, Z.; Zhang, Y.; Zhou, Z. An improved semi-Lagrangian time splitting spectral method for the semi-classical Schrödinger equation with vector potentials using NUFFT, Appl. Numer. Math., Volume 111 (2017), pp. 144-159 http://www.sciencedirect.com/science/article/pii/s0168927416301659 | DOI | Zbl

[15] Magnus, W. On the exponential solution of differential equations for a linear operator, Commun. Pure Appl. Math., Volume 7 (1954) no. 4, pp. 649-673 | DOI | MR

[16] Mie, G. Zur kinetischen Theorie der einatomigen Körper, Annalen der Physik, Volume 316 (1903) no. 8, pp. 657-697 | arXiv | DOI | Zbl

[17] Pauli, W. Zur Quantenmechanik des magnetischen Elektrons, Zeitschrift für Physik, Volume 43 (1927) no. 9, pp. 601-623 | DOI | Zbl

[18] Simon, B.; Reed, M. Fourier analysis, self-adjointness, Methods of modern mathematical physics, 2, Academic Press Inc., 1975 | Zbl

[19] Strang, G. On the construction and comparison of difference schemes, SIAM J. Numer. Anal., Volume 5 (1968), pp. 506-517 | DOI | MR | Zbl

[20] Werth, G.; Gheorghe, V. N.; Major, F. G. Charged Particle Traps II: Applications, Springer Series on Atomic, Optical, and Plasma Physics, Springer, 2009 https://books.google.ch/books?id=-wfwb45ntuwc

[21] Yoshida, H. Construction of higher order symplectic integrators, Phys. Lett., A, Volume 150 (1990) no. 5, pp. 262-268 http://www.sciencedirect.com/science/article/pii/0375960190900923 | DOI | MR

[22] Zhou, Z. Numerical approximation of the Schrödinger equation with the electromagnetic field by the Hagedorn wave packets, J. Comput. Phys., Volume 272 (2014), pp. 386-407 | DOI | Zbl

Cité par Sources :