Gaussian wave packet transform based numerical scheme for the semi-classical Schrödinger equation with random inputs

DOI10.1016/J.JCP.2019.109015MaRDI QIDQ2222699zbMATH OpenOpenAlexFDO

Authors Liu Liu, Zhennan Zhou, Shi Jin, Giovanni Russo

Publication date 27 January 2021

Published in Journal of Computational Physics (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/1903.08740

zbMATH Keywords

Gaussian wave packetSchrödinger equationstochastic collocationuncertainty quantification

Mathematics Subject Classification ID

PDEs with randomness, stochastic partial differential equations (35R60) Numerical solutions to stochastic differential and integral equations (65C30) Spectral, collocation and related methods for initial value and initial-boundary value problems involving PDEs (65M70) Semiclassical techniques, including WKB and Maslov methods applied to problems in quantum theory (81Q20)

Abstract: In this work, we study the semi-classical limit of the Schr"odinger equation with random inputs, and show that the semi-classical Schr"odinger equation produces $O(varepsilon)$ oscillations in the random variable space. With the Gaussian wave packet transform, the original Schr"odinger equation is mapped to an ODE system for the wave packet parameters coupled with a PDE for the quantity $w$ in rescaled variables. Further, we show that the $w$ equation does not produce $varepsilon$ dependent oscillations, and thus it is more amenable for numerical simulations. We propose multi-level sampling strategy in implementing the Gaussian wave packet transform, where in the most costly part, i.e. simulating the $w$ equation, it is sufficient to use $varepsilon$ independent samples. We also provide extensive numerical tests as well as meaningful numerical experiments to justify the properties of the numerical algorithm, and hopefully shed light on possible future directions.

Recommendations

• The Gaussian wave packet transform for the semi-classical Schrödinger equation with vector potentials
• The Gaussian wave packet transform: efficient computation of the semi-classical limit of the Schrödinger equation. I: Formulation and the one dimensional case
• The Gaussian wave packet transform: efficient computation of the semi-classical limit of the Schrödinger equation. II: Multidimensional case
• Fast Gaussian wavepacket transforms and Gaussian beams for the Schrödinger equation
• Convergence of a semiclassical wavepacket based time-splitting for the Schrödinger equation

Cites work

• A semi-Lagrangian time splitting method for the Schrödinger equation with vector potentials
• A Sparse Grid Stochastic Collocation Method for Partial Differential Equations with Random Input Data
• An improved semi-Lagrangian time splitting spectral method for the semi-classical Schrödinger equation with vector potentials using NUFFT
• Computing Semiclassical Quantum Dynamics with Hagedorn Wavepackets
• Convergence of frozen Gaussian approximation for high-frequency wave propagation
• Fast multiscale Gaussian wavepacket transforms and multiscale Gaussian beams for the wave equation
• Frozen Gaussian approximation with surface hopping for mixed quantum-classical dynamics: a mathematical justification of fewest switches surface hopping algorithms
• High-Order Collocation Methods for Differential Equations with Random Inputs
• scientific article; zbMATH DE number 3846088 (Why is no real title available?)
• scientific article; zbMATH DE number 1284216 (Why is no real title available?)
• Mathematical and computational methods for semiclassical Schrödinger equations
• Numerical approximation of the Schrödinger equation with the electromagnetic field by the Hagedorn wave packets
• Numerical methods for stochastic computations. A spectral method approach.
• On a Schrödinger-Landau-Lifshitz system: variational structure and numerical methods
• On time-splitting spectral approximations for the Schrödinger equation in the semiclassical regime
• On Wigner measures
• Radiative transport in a periodic structure
• Raising and lowering operators for semiclassical wave packets
• Semi-Eulerian and High Order Gaussian Beam Methods for the Schrödinger Equation in the Semiclassical Regime
• Spectral Methods for Uncertainty Quantification
• Stochastic finite element methods for partial differential equations with random input data
• The Gaussian wave packet transform for the semi-classical Schrödinger equation with vector potentials
• The Gaussian wave packet transform: efficient computation of the semi-classical limit of the Schrödinger equation. I: Formulation and the one dimensional case
• The Gaussian wave packet transform: efficient computation of the semi-classical limit of the Schrödinger equation. II: Multidimensional case
• The Wiener--Askey Polynomial Chaos for Stochastic Differential Equations

Cited in

(11)

• Fourier spectral method for a class of nonlinear Schrödinger models
• The Gaussian wave packet transform: efficient computation of the semi-classical limit of the Schrödinger equation. I: Formulation and the one dimensional case
• The Gaussian wave packet transform: efficient computation of the semi-classical limit of the Schrödinger equation. II: Multidimensional case
• The Gaussian wave packet transform for the semi-classical Schrödinger equation with vector potentials
• On a neural network approach for solving potential control problem of the semiclassical Schrödinger equation
• Efficient numerical method for the Schrödinger equation with high-contrast potentials
• On a class of multifidelity methods for the semiclassical Schrödinger equation with uncertainties
• Linear canonical wavelet transform and linear canonical wave packet transform on the Schwartz type spaces
• Strong approximation of stochastic semiclassical Schrödinger equation with multiplicative noise
• Structure-preserving conservative phase-amplitude PINNs for semiclassical Schrödinger-type equations
• A deep learning WKB-PINNs framework for forward and inverse problems of semiclassical nonlinear Schrödinger equations

This page was built for publication: Gaussian wave packet transform based numerical scheme for the semi-classical Schrödinger equation with random inputs

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2222699)