Error estimate of a quasi-Monte Carlo time-splitting pseudospectral method for nonlinear Schrodinger equation with random potentials

Zhizhang Wu,Zhiwen Zhang,Xiaofei Zhao
2023-11-19
Abstract:In this paper, we consider the numerical solution of a nonlinear Schrodinger equation with spatial random potential. The randomly shifted quasi-Monte Carlo (QMC) lattice rule combined with the time-splitting pseudospectral discretization is applied and analyzed. The nonlinearity in the equation induces difficulties in estimating the regularity of the solution in random space. By the technique of weighted Sobolev space, we identify the possible weights and show the existence of QMC that converges optimally at the almost-linear rate without dependence on dimensions. The full error estimate of the scheme is established. We present numerical results to verify the accuracy and investigate the wave propagation.
Numerical Analysis
What problem does this paper attempt to address?
### Problems Addressed by the Paper The paper primarily focuses on the numerical solution and error estimation of the one-dimensional nonlinear Schrödinger equation (NLS) with a random potential. Specifically, the paper attempts to address the following core issues: 1. **Parameterization of the Random Potential**: - The random potential is parameterized using the Karhunen-Loève expansion. 2. **Numerical Discretization Method**: - By combining the quasi-Monte Carlo (QMC) method and the time-splitting pseudospectral method (TSFP), an efficient numerical solution scheme is proposed. - The convergence properties of this scheme in both time and space are studied, and strict error estimates are provided. 3. **Calculation of Expectation Values of Physical Quantities**: - The method for calculating the expectation values of physical observables (such as the modulus squared of the density wave function) is analyzed, and its regularity in random space is discussed. 4. **Truncation Error in Random Dimensions**: - The impact of truncating random dimensions on the solution is studied, and corresponding error estimates are given. ### Summary The paper aims to develop an efficient and accurate method to solve the one-dimensional nonlinear Schrödinger equation with a random potential. By combining the QMC method and the TSFP method, the researchers not only propose an effective numerical solution scheme but also provide a detailed analysis of its convergence and error. This offers important mathematical tools for understanding the behavior of the nonlinear Schrödinger equation under random potentials.