ترغب بنشر مسار تعليمي؟ اضغط هنا

Strong rates of convergence of a splitting scheme for Schr{o}dinger equations with nonlocal interaction cubic nonlinearity and white noise dispersion

378   0   0.0 ( 0 )
 نشر من قبل Charles-Edouard Br\\'ehier
 تاريخ النشر 2020
  مجال البحث الهندسة المعلوماتية
والبحث باللغة English




اسأل ChatGPT حول البحث

We analyse a splitting integrator for the time discretization of the Schrodinger equation with nonlocal interaction cubic nonlinearity and white noise dispersion. We prove that this time integrator has order of convergence one in the $p$-th mean sense, for any $pgeq1$ in some Sobolev spaces. We prove that the splitting schemes preserves the $L^2$-norm, which is a crucial property for the proof of the strong convergence result. Finally, numerical experiments illustrate the performance of the proposed numerical scheme.



قيم البحث

اقرأ أيضاً

213 - Lu Zhang 2021
This paper proposes and analyzes an ultra-weak local discontinuous Galerkin scheme for one-dimensional nonlinear biharmonic Schr{o}dinger equations. We develop the paradigm of the local discontinuous Galerkin method by introducing the second-order sp atial derivative as an auxiliary variable instead of the conventional first-order derivative. The proposed semi-discrete scheme preserves a few physically relevant properties such as the conservation of mass and the conservation of Hamiltonian accompanied by its stability for the targeted nonlinear biharmonic Schr{o}dinger equations. We also derive optimal $L^2$-error estimates of the scheme that measure both the solution and the auxiliary variable. Several numerical studies demonstrate and support our theoretical findings.
We analyze the qualitative properties and the order of convergence of a splitting scheme for a class of nonlinear stochastic Schrodinger equations driven by additive It^o noise. The class of nonlinearities of interest includes nonlocal interaction cu bic nonlinearities. We show that the numerical solution is symplectic and preserves the expected mass for all times. On top of that, for the convergence analysis, some exponential moment bounds for the exact and numerical solutions are proved. This enables us to provide strong orders of convergence as well as orders of convergence in probability and almost surely. Finally, extensive numerical experiments illustrate the performance of the proposed numerical scheme.
Asymptotic reductions of a defocusing nonlocal nonlinear Schr{o}dinger model in $(3+1)$-dimensions, in both Cartesian and cylindrical geometry, are presented. First, at an intermediate stage, a Boussinesq equation is derived, and then its far-field, in the form of a variety of Kadomtsev-Petviashvilli (KP) equations for right- and left-going waves, is found. KP models includ
141 - Zhihui Liu , Zhonghua Qiao 2018
We establish a general theory of optimal strong error estimation for numerical approximations of a second-order parabolic stochastic partial differential equation with monotone drift driven by a multiplicative infinite-dimensional Wiener process. The equation is spatially discretized by Galerkin methods and temporally discretized by drift-implicit Euler and Milstein schemes. By the monotone and Lyapunov assumptions, we use both the variational and semigroup approaches to derive a spatial Sobolev regularity under the $L_omega^p L_t^infty dot H^{1+gamma}$-norm and a temporal Holder regularity under the $L_omega^p L_x^2$-norm for the solution of the proposed equation with an $dot H^{1+gamma}$-valued initial datum for $gammain [0,1]$. Then we make full use of the monotonicity of the equation and tools from stochastic calculus to derive the sharp strong convergence rates $O(h^{1+gamma}+tau^{1/2})$ and $O(h^{1+gamma}+tau^{(1+gamma)/2})$ for the Galerkin-based Euler and Milstein schemes, respectively.
In this paper, we introduce two new families of generalised Hermite polynomials/functions (GHPs/GHFs) in arbitrary dimensions, and develop efficient and accurate generalised Hermite spectral algorithms for PDEs with integral fractional Laplacian (IFL ) and/or Schr{o}dinger operators in $mathbb R^d.$ As a generalisation of the G. Szeg{o}s family in 1D (1939), the first family of GHPs (resp. GHFs) are orthogonal with respect to $|bx|^{2mu} e^{-|bx|^2}$ (resp. $|bx |^{2mu}$) in $mathbb R^d$. We further define adjoint generalised Hermite functions (A-GHFs) which have an interwoven connection with the corresponding GHFs through the Fourier transform, and which are orthogonal with respect to the inner product $[u,v]_{H^s(mathbb R^d)}=((-Delta)^{s/ 2}u, (-Delta)^{s/2} v )_{mathbb R^d}$ associated with the IFL of order $s>0$. Thus, the spectral-Galerkin method using A-GHFs as basis functions leads to a diagonal stiffness matrix for the IFL (which is known to be notoriously difficult and expensive to discretise). The new basis also finds efficient and accurate in solving PDEs with the fractional Schr{o}dinger operator: $(-Delta)^s +|bs x|^{2mu}$ with $sin (0,1]$ and $mu>-1/2.$ Following the same spirit, we construct the second family of GHFs, dubbed as Muntz-type generalised Hermite functions (M-GHFs), which are orthogonal with respect to an inner product associated with the underlying Schr{o}dinger operator, and are tailored to the singularity of the solution at the origin. We demonstrate that the Muntz-type GHF spectral method leads to sparse matrices and spectrally accurate to some Schr{o}dinger eigenvalue problems.
التعليقات
جاري جلب التعليقات جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
mircosoft-partner

هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا