Global geometrical optics method for vector-valued Schrödinger problems

Abstract

We extend the theory of global geometrical optics method, proposed originally for the linear scalar high-frequency wave-like equations in [Commun. Math. Sci., 2013, 11(1): 105-140], to the more general vectorvalued Schrödinger problems in the semi-classical regime. The key ingredient in the global geometrical optics method is a moving frame technique in the phase space. The governing equation is transformed into a new equation but of the same type when expressed in any moving frame induced by the underlying Hamiltonian ow. The classical Wentzel-Kramers-Brillouin (WKB) analysis benefits from this treatment as it maintains valid for arbitrary but fixed evolutionary time. It turns out that a WKB-type function defined merely on the underlying Lagrangian submanifold can be obtained with the help of this moving frame technique, and from which a uniform first-order approximation of the wave field can be derived, even around caustics. The general theory is exemplified by two specific instances. One is the two-level Schrödinger system and the other is the periodic Schrödinger equation. Numerical tests validate the theoretical results.

Keywords

Global geometrical optics method / Hamiltonian system / unitary representation / caustics / semiclassical approximation

References

[1] Bao W, Jin S, Markowich P A. On time-splitting spectral approximations for the Schrödinger equation in the semiclassical regime. J Comput Phys, 2002, 175: 487–524

[2] Carles R, Sparber C. Semiclassical wave packet dynamics in Schrödinger equations with periodic potentials. Discrete Contin Dyn Syst Ser B, 2012, 17(3): 757–774

[3] Chai L, Jin S, Li Q. Semi-classical models for the Schrödinger equation with periodic potentials and band crossings. Kinet Relat Models, 2013, 6: 505–532

[4] Chai L, Jin S, Li Q, Morandi O. A multiband semiclassical model for surface hopping quantum dynamics. Multiscale Model Simul, 2015, 13(1): 205–230

[5] Chai L, Jin S, Markowich P A. A hybrid method for computing the Schrödinger equations with periodic potential with band-crossings in the momentum space. Commun Comput Phys, Special Issue in Honor of the 80th Birthday of Prof. Houde Han (to appear)

[6] Faraj A, Jin S. The Landau-Zener transition and the surface hopping method for the 2D Dirac equation for graphene. Commun Comput Phys, 2017, 21(2): 313–357

[7] Folland G B. Harmonic Analysis in Phase Space.Princeton: Princeton Univ Press, 1989

[8] Gerard P, Markowich P A, Mauser N, Poupaud F. Homogenization limits and Wigner transforms. Comm Pure Appl Math, 1997, 50: 323–379

[9] Hagedorn G, Robinson S. Bohr-Sommerfeld quantization rules in the semiclassical limit. J Phys A, 1998, 31(50): 10113–10130

[10] Holstein B. Topics in Advanced Quantum Mechanics.Redwood City: Addison-Wesley, 1992

[11] Huang Z, Jin S, Markowich P A, Sparber C. A Bloch decomposition based split-step pseudo spectral method for quantum dynamics with periodic potentials. SIAM J Sci Comput, 2007, 29(2): 515–538

[12] Huang Z, Jin S, Markowich P A, Sparber C, Zheng C. A time-splitting spectral scheme for the Maxwell-Dirac system. J Comput Phys, 2005, 208(2): 761–789

[13] Jin S, Markowich P A, Sparber C. Mathematical and computational methods for semi-classical Schrödinger equations. Acta Numer, 2011, 20: 211–289

[14] Karasev M V. Connections on Lagrangian submanifolds and some quasiclassical approximation problems I. J Soviet Math, 1992, 59(5): 1053–1062

[15] Maslov V P, Fedoryuk M V. Semi-Classical Approximation in Quantum Mechanics.Dordrecht: D Reidel Pub Co, 1981

[16] Voros A. Wentzel-Kramers-Brillouin method in the Bargmann representation. Phys Rev A, 1989, 40(3): 6814–6825

[17] Wu H, Huang Z, Jin S, Yin D. Gaussian beam methods for the Dirac equation in the semi-classical regime. Commun Math Sci, 2012, 10: 1301–1315

[18] Zheng C. Global geometrical optics method. Commun Math Sci, 2013, 11(1): 105–140

RIGHTS & PERMISSIONS

Higher Education Press and Springer-Verlag GmbH Germany, part of Springer Nature