Strong Approximation of Stochastic Semiclassical Schroedinger Equation with Multiplicative Noise

Lihai Ji,Zhihui Liu
2024-08-17
Abstract:We consider the stochastic nonlinear Schroedinger equation driven by a multiplicative noise in a semiclassical regime, where the Plank constant is small. In this regime, the solution of the equation exhibits high-frequency oscillations. We design an efficient numerical method combining the spectral Galerkin approximation and the midpoint scheme. This accurately approximates the solution, or at least of the associated physical observables. Furthermore, the strong convergence rate for the proposed scheme is derived, which explicitly depends on the Planck constant. This conclusion implies the semiclassical regime's admissible meshing strategies for obtaining "correct" physical observables.
Numerical Analysis
What problem does this paper attempt to address?
The problem that this paper attempts to solve is: in the semi - classical regime, the strong approximation problem of the stochastic nonlinear Schrödinger equation (SNLSE) with multiplicative noise. Specifically, the article focuses on how to efficiently and accurately approximate the solution of this equation or the related physical observables numerically when the Planck constant \(\varepsilon\) is small. ### Specific description of the problem 1. **Equation form**: The paper studied the stochastic semi - classical Schrödinger equation in the following form: \[ i\varepsilon du^\varepsilon_t+\left(\frac{\varepsilon^2}{2}\Delta u^\varepsilon_t + F(u^\varepsilon_t)\right)dt = G(u^\varepsilon_t)dW_t,\quad(t,x)\in(0,T]\times T, \] where \(T:=[0,1]^d\), \(W_t\) is an \(L^2(T;\mathbb{R})\)-valued \(Q\)-Wiener process, and \(\varepsilon\) is the scaled Planck constant, representing the ratio of the micro - and macro - scales. 2. **Challenges**: - When \(\varepsilon\ll1\), the solution of the equation will exhibit high - frequency oscillations, which poses a great challenge to numerical calculations. - It is necessary to design efficient numerical methods to approximate the solution or at least the related physical observables. 3. **Objectives**: - Design a numerical method that combines spectral Galerkin approximation and midpoint scheme to efficiently and accurately approximate the solution of the equation. - Derive the strong convergence rate of the proposed scheme and clarify its dependence on the Planck constant \(\varepsilon\). - Propose an acceptable mesh strategy applicable to the semi - classical regime to obtain the "correct" physical observables. ### Main contributions - **Numerical method**: A method that combines spectral Galerkin approximation and midpoint scheme was proposed, which can efficiently and accurately approximate the solution of the stochastic semi - classical Schrödinger equation with multiplicative noise. - **Theoretical analysis**: The strong convergence rate of the proposed numerical method was derived, and its dependence on the Planck constant \(\varepsilon\) was clarified, providing a theoretical basis for numerical simulations. - **Mesh strategy**: An acceptable mesh strategy applicable to the semi - classical regime was proposed to ensure the accuracy of numerical results. Through these works, the paper aims to deepen the understanding of the highly oscillatory behavior of the semi - classical NLS equation perturbed by Wiener processes and provide effective numerical tools to deal with such problems.