2000 character limit reached
WKB-based third order method for the highly oscillatory 1D stationary Schrödinger equation (2402.18406v1)
Published 28 Feb 2024 in math.NA and cs.NA
Abstract: This paper introduces an efficient high-order numerical method for solving the 1D stationary Schr\"odinger equation in the highly oscillatory regime. Building upon the ideas from [2], we first analytically transform the given equation into a smoother (i.e. less oscillatory) equation. By developing sufficiently accurate quadratures for several (iterated) oscillatory integrals occurring in the Picard approximation of the solution, we obtain a one-step method that is third order w.r.t. the step size. The accuracy and efficiency of the method are illustrated through several numerical examples.
- Advanpix LLC. Multiprecision Computing Toolbox for Matlab version 5.1.0.15432. Advanpix LLC, 2023.
- WKB-Based Schemes for the Oscillatory 1D Schrödinger Equation in the Semiclassical Limit. SIAM J. Numer. Anal., 49:1436–1460, 2011.
- Optimally truncated WKB approximation for the 1D stationary Schrödinger equation in the highly oscillatory regime. arXiv, abs/2401.10141, 2024.
- WKB-method for the 1D Schrödinger equation in the semi-classical limit: enhanced phase treatment. BIT Numerical Mathematics, 62:1–22, 2022.
- Barycentric lagrange interpolation. SIAM Rev., 46:501–517, 2004.
- A method for numerical integration on an automatic computer. Numerische Mathematik, 2:197–205, 1960.
- Theory of the alternating-gradient synchrotron. Annals of Physics, 281:360–408, 1958.
- F. Ihlenburg and I. Babuška. Finite element solution of the Helmholtz equation with high wave number Part I: The hℎhitalic_h-version of the FEM. Computers and Mathematics with Applications, 30(9):9–37, 1995.
- Highly oscillatory quadrature: the story so far. Proceedings of ENUMATH 2005, Santiago de Compostela, pages 97–118, 2006.
- T. Jahnke and C. Lubich. Numerical integrators for quantum dynamics close to the adiabatic limit. Numerische Mathematik, 94:289–314, 2003.
- WKB-based scheme with adaptive step size control for the Schrödinger equation in the highly oscillatory regime. J. Comput. Appl. Math., 404:113905, 2021.
- Quantenmechanik. Akademie-Verlag, 1985.
- H. R. Lewis. Motion of a Time-Dependent Harmonic Oscillator, and of a Charged Particle in a Class of Time-Dependent, Axially Symmetric Electromagnetic Fields. Phys. Rev., 172:1313–1315, 1968.
- Adiabatic Integrators for Highly Oscillatory Second-Order Linear Differential Equations with Time-Varying Eigendecomposition. BIT Numerical Mathematics, 45:91–115, 2005.
- J. Martin and D. J. Schwarz. WKB approximation for inflationary cosmological perturbations. Physical Review D, 67:083512, 2003.
- Transient Schrödinger-Poisson simulations of a high-frequency resonant tunneling diode oscillator. J. Comput. Phys., 239:187–205, 2013.
- C. Negulescu. Asymptotical models and numerical schemes for quantum systems. PhD thesis, Université Paul Sabatier, Toulouse, 2005.
- NIST Handbook of Mathematical Functions. Cambridge University Press, 2010.
- Resonant tunneling diodes: models and properties. Proc. IEEE, 86:641–660, 1998.
- S. Winitzki. Cosmological particle production and the precision of the WKB approximation. Physical Review D, 72:104011, 2005.