Instability of the finite‐difference split‐step method applied to the nonlinear Schrödinger equation. I. standing soliton

We consider numerical instability that can be observed in simulations of solitons of the nonlinear Schrödinger equation (NLS) by a split‐step method (SSM) where the linear part of the evolution is solved by a finite‐difference discretization. The von Neumann analysis predicts that this method is unconditionally stable on the background of a constant‐amplitude plane wave. However, simulations show that the method can become unstable on the background of a soliton. We present an analysis explaining this instability. Both this analysis and the features and threshold of the instability are substantially different from those of the Fourier SSM, which computes the linear part of the NLS by a spectral discretization. For example, the modes responsible for the numerical instability are not similar to plane waves, as for the Fourier SSM or, more generally, in the von Neumann analysis. Instead, they are localized at the sides of the soliton. This also makes them different from “physical” (as opposed to numerical) unstable modes of nonlinear waves, which (the modes) are localized around the “core” of a solitary wave. Moreover, the instability threshold for thefinite‐difference split‐step method is considerably relaxed compared with that for the Fourier split‐step. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 1002–1023, 2016     

Instability of the finite‐difference split‐step method applied to the generalized nonlinear Schrödinger equation. III. external potential and oscillating pulse solutions

This is the final part of a series of articles where we have studied numerical instability (NI) of localized solutions of the generalized nonlinear Schrödinger equation (gNLS). It extends our earlier studies of this topic in two ways. First, it examines differences in the development of the NI between the case of the purely cubic NLS and the case where the gNLS has an external bounded potential. Second, it investigates how the NI is affected by the oscillatory dynamics of the simulated pulse. The latter situation is common when the initial condition is not an exact stationary soliton. We have found that in this case, the NI may remain weak when the time step exceeds the threshold quite significantly. This means that the corresponding numerical solution, while formally numerically unstable, can remain sufficiently accurate over long times, because the numerical noise will stay small. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 633–650, 2017     

Instability analysis of the split‐step Fourier method on the background of a soliton of the nonlinear Schrödinger equation

We analyze a numerical instability that occurs in the well‐known split‐step Fourier method on the background of a soliton. This instability is found to be very sensitive to small changes of the parameters of both the numerical grid and the soliton, unlike the instability of most finite‐difference schemes. Moreover, the principle of “frozen coefficients,” in which variable coefficients are treated as “locally constant” for the purpose of stability analysis, is strongly violated for the instability of the split‐step method on the soliton background. Our analysis quantitatively explains all these features. It is enabled by the fact that the period of oscillations of the unstable Fourier modes is much smaller than the width of the soliton. Our analysis is different from the von Neumann analysis in that it requires spatially growing or decaying harmonics (not localized near the boundaries) as opposed to purely oscillatory ones. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 28: 641–669, 2012     

Two‐grid mixed finite‐element methods for nonlinear Schrödinger equations

Two‐grid mixed finite element schemes are developed for solving both steady state and unsteady state nonlinear Schrödinger equations. The schemes use discretizations based on a mixed finite‐element method. The two‐grid approach yields iterative procedures for solving the nonlinear discrete equations. The idea is to relegate all of the Newton‐like iterations to grids much coarser than the final one, with no loss in order of accuracy. Numerical tests are performed. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 28: 63‐73, 2012     

On a discrete‐time collocation method for the nonlinear Schrödinger equation with wave operator

We consider a discrete‐time orthogonal spline collocation scheme for solving Schrödinger equation with wave operator. The scheme is proposed recently by Wang et al. (J Comput Appl Math 235 (2011), 1993–2005) and is showed to have high‐order convergence rate when a parameter θ in the scheme is not less than $\frac{1}{4}$. In this article, we show that the result can be extended to include $\theta\in(0,\frac{1}{4})$ under an assumption. Numerical example is given to justify the theoretical result. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

Extended Lagrangian approach for the defocusing nonlinear Schrödinger equation

We study the defocusing nonlinear Schrödinger (NLS) equation written in hydrodynamic form through the Madelung transform. From the mathematical point of view, the hydrodynamic form can be seen as the Euler–Lagrange equations for a Lagrangian submitted to a differential constraint corresponding to the mass conservation law. The dispersive nature of the NLS equation poses some major numerical challenges. The idea is to introduce a two‐parameter family of extended Lagrangians, depending on a greater number of variables, whose Euler–Lagrange equations are hyperbolic and accurately approximate NLS equation in a certain limit. The corresponding hyperbolic equations are studied and solved numerically using Godunov‐type methods. Comparison of exact and asymptotic solutions to the one‐dimensional cubic NLS equation (“gray” solitons and dispersive shocks) and the corresponding numerical solutions to the extended system was performed. A very good accuracy of such a hyperbolic approximation was observed.     

The discreteness of the spectrum of the Schrödinger operator equation and some properties of the s‐numbers of the inverse Schrödinger operator

In this article, we investigate the discreteness and some other properties of the spectrum for the Schrödinger operator L defined by the formula on the space L₂(H, [0, ∞)) , where H is a Hilbert space. For the first time, an estimate is obtained for sum of the s‐numbers of the inverse Schrödinger operator. The obtained results were applied to the Laplace's equation in an angular region.     

Exact soliton solution for the fourth‐order nonlinear Schrödinger equation with generalized cubic‐quintic nonlinearity

In this paper, we investigate the fourth‐order nonlinear Schrödinger equation with parameterized nonlinearity that is generalized from regular cubic‐quintic formulation in optics and ultracold physics scenario. We find the exact solution of the fourth‐order generalized cubic‐quintic nonlinear Schrödinger equation through modified F‐expansion method, identifying the particular bright soliton behavior under certain external experimental setting, with the system's particular nonlinear features demonstrated. Copyright © 2016 John Wiley & Sons, Ltd.     

Finite‐dimensional global attractor for a semi‐discrete fractional nonlinear Schrödinger equation

We consider a semi‐discrete in time Crank–Nicolson scheme to discretize a weakly damped forced nonlinear fractional Schrödinger equation u _t ?i (?Δ)^α u +i |u |²u +γ u =f for considered in the the whole space . We prove that such semi‐discrete equation provides a discrete infinite‐dimensional dynamical system in that possesses a global attractor in . We show also that if the external force is in a suitable weighted Lebesgue space, then this global attractor has a finite fractal dimension. Copyright © 2017 John Wiley & Sons, Ltd.     

An application of exp‐function method to approximate general and explicit solutions for nonlinear Schrödinger equations

In this work, we implement a relatively new analytical technique, the exp‐function method, for solving nonlinear equations and absolutely a special form of generalized nonlinear Schrödinger equations, which may contain high‐nonlinear terms. This method can be used as an alternative to obtain analytical and approximate solutions of different types of fractional differential equations, which is applied in engineering mathematics. Some numerical examples are presented to illustrate the efficiency and the reliability of exp‐function method. It is predicted that exp‐function method can be found widely applicable in engineering. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 27: 1016–1025, 2011     

Nonvanishing at spatial extremity solutions of the defocusing nonlinear Schrödinger equation

We show local existence of certain type of solutions for the Cauchy problem of the defocusing nonlinear Schrödinger equation with pure power nonlinearity, in various cases of open sets, unbounded or bounded. These solutions do not vanish at the boundary or at infinity. We also show, in certain cases, that these solutions are unique and global.     

Asymptotics for the third‐order nonlinear Schrödinger equation in the critical case

We consider the Cauchy problem for the third‐order nonlinear Schrödinger equation where and is the Fourier transform. Our purpose in this paper is to prove the large time asymptoitic behavior of solutions for the defocusing case λ > 0 with a logarithmic correction under the non zero mass condition Copyright © 2016 John Wiley & Sons, Ltd.     

Efficient energy‐preserving scheme of the three‐coupled nonlinear Schrödinger equation

An energy‐preserving scheme is proposed for the three‐coupled nonlinear Schrödinger (T‐CNLS) equation. The T‐CNLS equation is rewritten into the classical Hamiltonian form. Then the spatial variable is discretized by using high‐order compact method to convert it into a finite‐dimensional Hamiltonian system. Next, a second‐order averaged vector field (AVF) method is employed in time which results in an energy‐preserving scheme. Some theoretical results such as convergence are investigated. In addition, it provides some numerical examples to illustrate the robustness and reliability of the theoretical results. It also explores the role of the parameters in the model and initial condition on the wave propagation.     

A complex tanh-function method applied to nonlinear equations of Schrödinger type

A complex tanh-function method is introduced for constructing exact travelling wave solutions of nonlinear partial differential equations with complex phases and solutions. The scheme is implemented for obtaining multiple soliton solutions to the nonlinear cubic Schrödinger equation and a generalized Schrödinger-like equation. In additon. an ansätze is proposed to obtain stationary soliton solutions of the cubic Schrödinger equation.