A modified splitting method for the cubic nonlinear Schrödinger equation
share
As a classical time-stepping method, it is well-known that the Strang splitting method reaches the first-order accuracy by losing two spatial derivatives. In this paper, we propose a modified splitting method for the 1D cubic nonlinear Schrödinger equation: u^n+1=e^iτ/2∂_x^2𝒩_τ[e^iτ/2∂_x^2(Π_τ +e^-2π iλ M_0τΠ^τ)u^n], with 𝒩_t(ϕ)=e^-iλ t|Π_τϕ|^2ϕ, and M_0 is the mass of the initial data. Suitably choosing the filters Π_τ and Π^τ, it is shown rigorously that it reaches the first-order accuracy by only losing 3/2-spatial derivatives. Moreover, if γ∈ (0,1), the new method presents the convergence rate of τ^4γ/4+γ in L^2-norm for the H^γ-data; if γ∈ [1,2], it presents the convergence rate of τ^2/5(1+γ)- in L^2-norm for the H^γ-data. for the first-order convergence in L^2-norm is only H^3/2+. These results are better than the expected ones for the standard (filtered) Strang splitting methods. Moreover, the mass is conserved: 1/2π∫_𝕋 |u^n(x)|^2 d x≡ M_0, n=0,1,…, L . The key idea is based on the observation that the low frequency and high frequency components of solutions are almost separated (up to some smooth components). Then the algorithm is constructed by tracking the solution behavior at the low and high frequency components separately.
READ FULL TEXT