Tanh Jacobi spectral collocation method for the numerical simulation of nonlinear Schrödinger equations on unbounded domain

We present a class of orthogonal functions on infinite domain based on Jacobi polynomials. These functions are generated by applying a tanh transformation to Jacobi polynomials. We construct interpolation and projection error estimates using weighted pseudo-derivatives tailored to the involved mapping. Then, using the nodes of the newly introduced tanh Jacobi functions, we develop an efficient spectral tanh Jacobi collocation method for the numerical simulation of nonlinear Schrödinger equations on the infinite domain without using artificial boundary conditions. The applicability and accuracy of the solution method are demonstrated by two numerical examples for solving the nonlinear Schrödinger equation and the nonlinear Ginzburg–Landau equation.     

Initial‐Boundary Value Problems for the Defocusing Nonlinear Schrödinger Equation in the Semiclassical Limit

Initial‐boundary value problems for integrable nonlinear partial differential equations have become tractable in recent years due to the development of so‐called unified transform techniques. The main obstruction to applying these methods in practice is that calculation of the spectral transforms of the initial and boundary data requires knowledge of too many boundary conditions, more than are required to make the problem well‐posed. The elimination of the unknown boundary values is frequently addressed in the spectral domain via the so‐called global relation, and types of boundary conditions for which the global relation can be solved are called linearizable. For the defocusing nonlinear Schrödinger equation, the global relation is only known to be explicitly solvable in rather restrictive situations, namely homogeneous boundary conditions of Dirichlet, Neumann, and Robin (mixed) type. General nonhomogeneous boundary conditions are not known to be linearizable. In this paper, we propose an explicit approximation for the nonlinear Dirichlet‐to‐Neumann map supplied by the defocusing nonlinear Schrödinger equation and use it to provide approximate solutions of general nonhomogeneous boundary value problems for this equation posed as an initial‐boundary value problem on the half‐line. Our method sidesteps entirely the solution of the global relation. The accuracy of our method is proven in the semiclassical limit, and we provide explicit asymptotics for the solution in the interior of the quarter‐plane space‐time domain.     

Conservative finite difference methods for fractional Schrödinger–Boussinesq equations and convergence analysis

In this paper, two conservative finite difference schemes for fractional Schrödinger–Boussinesq equations are formulated and investigated. The convergence of the nonlinear fully implicit scheme is established via discrete energy method, while the linear semi‐implicit scheme is analyzed by means of mathematical induction method. Our schemes are proved to preserve the total mass and energy in discrete level. The numerical results are given to confirm the theoretical analysis.     

Multidomain spectral method for Schrödinger equations

A multidomain spectral method with compactified exterior domains combined with stable second and fourth order time integrators is presented for Schrödinger equations. The numerical approach allows high precision numerical studies of solutions on the whole real line. At examples for the linear and cubic nonlinear Schrödinger equation, this code is compared to transparent boundary conditions and perfectly matched layers approaches. The code can deal with asymptotically non vanishing solutions as the Peregrine breather being discussed as a model for rogue waves. It is shown that the Peregrine breather can be numerically propagated with essentially machine precision, and that localized perturbations of this solution can be studied.     

Fast evaluation of exact transparent boundary condition for one-dimensional cubic nonlinear Schrödinger equation

Fast evaluation of the exact transparent boundary condition for the one-dimensional cubic nonlinear Schrödinger equation is considered in this paper. In [J. Comput. Math., 2007, 25(6): 730–745], the author proposed a fast evaluation method for the half-order time derivative operator. In this paper, we apply this method for the exact transparent boundary condition for the one-dimensional cubic nonlinear Schrödinger equation. Numerical tests demonstrate the effectiveness of the proposed method.     

Transparent boundary conditions based on the pole condition for time‐dependent,two‐dimensional problems

The pole condition approach for deriving transparent boundary conditions is extended to the time‐dependent, two‐dimensional case. Nonphysical modes of the solution are identified by the position of poles of the solution's spatial Laplace transform in the complex plane. By requiring the Laplace transform to be analytic on some problem‐dependent complex half‐plane, these modes can be suppressed. The resulting algorithm computes a finite number of coefficients of a series expansion of the Laplace transform, thereby providing an approximation to the exact boundary condition. The resulting error decays super‐algebraically with the number of coefficients, so relatively few additional degrees of freedom are sufficient to reduce the error to the level of the discretization error in the interior of the computational domain. The approach shows good results for the Schrödinger and the drift‐diffusion equation but, in contrast to the one‐dimensional case, exhibits instabilities for the wave and Klein–Gordon equation. Numerical examples are shown that demonstrate the good performance in the former and the instabilities in the latter case. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

Unconditional convergence and optimal error estimates of the Euler semi-implicit scheme for a generalized nonlinear Schrödinger equation

In this paper, we focus on a linearized backward Euler scheme with a Galerkin finite element approximation for the time-dependent nonlinear Schrödinger equation. By splitting an error estimate into two parts, one from the spatial discretization and the other from the temporal discretization, we obtain unconditionally optimal error estimates of the fully-discrete backward Euler method for a generalized nonlinear Schrödinger equation. Numerical results are provided to support our theoretical analysis and efficiency of this method.     

N solutions for a derivative nonlinear Schrödinger‐type equation via Riemann‐Hilbert approach

The inverse scattering transform for the derivative nonlinear Schrödinger‐type equation is studied via the Riemann‐Hilbert approach. In the direct scattering process, the spectral analysis of the Lax pair is performed, from which a Riemann‐Hilbert problem is established for the derivative nonlinear Schrödinger‐type equation. In the inverse scattering process, N‐soliton solutions of the derivative nonlinear Schrödinger‐type equation are obtained by solving Riemann‐Hilbert problems corresponding to the reflectionless cases. Moreover, the dynamics of the exact solutions are discussed.     

An efficient spectral‐Galerkin method based on a dimension reduction scheme for eigenvalue problems of Schrödinger equations

In this paper, we propose an efficient spectral‐Galerkin method based on a dimension reduction scheme for eigenvalue problems of Schrödinger equations. Firstly, we carry out a truncation from a three‐dimensional unbounded domain to a bounded spherical domain. By using spherical coordinate transformation and spherical harmonic expansion, we transform the original problem into a series of one‐dimensional eigenvalue problem that can be solved effectively. Secondly, we introduce a weighted Sobolev space to treat the singularity in the effective potential. Using the property of orthogonal polynomials in weighted Sobolev space, the error estimate for the approximate eigenvalues and corresponding eigenfunctions are proved. Error estimates show that our numerical method can achieve spectral accuracy for approximate eigenvalues and eigenfunctions. Finally, we give some numerical examples to demonstrate the efficiency of our algorithms and the correctness of the theoretical results.     

Approximation of Invariant Measure for Damped Stochastic Nonlinear Schrödinger Equation via an Ergodic Numerical Scheme

In order to inherit numerically the ergodicity of the damped stochastic nonlinear Schrödinger equation with additive noise, we propose a fully discrete scheme, whose spatial direction is based on spectral Galerkin method and temporal direction is based on a modification of the implicit Euler scheme. We not only prove the unique ergodicity of the numerical solutions of both spatial semi-discretization and full discretization, but also present error estimations on invariant measures, which gives order 2 in spatial direction and order \({\frac 12}\) in temporal direction under certain hypotheses.     

Error and stability analysis of numerical solution for the time fractional nonlinear Schrödinger equation on scattered data of general‐shaped domains

In present work, a kind of spectral meshless radial point interpolation (SMRPI) technique is applied to the time fractional nonlinear Schrödinger equation in regular and irregular domains. The applied approach is based on erudite combination of meshless methods and spectral collocation techniques. The point interpolation method with the help of radial basis functions is used to construct shape functions which play as basis functions in the frame of SMRPI. It is proved the scheme is unconditionally stable with respect to the time variable in and also convergent by the order of convergence , . In the current work, the thin plate spline are used as the basis functions and to eliminate the nonlinearity, a simple predictor‐corrector (P‐C) scheme is performed. It is shown that the SMRPI solution, as a complex function, is suitable one for the time fractional nonlinear Schrödinger equation. The results of numerical experiments are compared to analytical solutions to confirm the reliable treatment of these stable solutions. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1043–1069, 2017     

Galerkin-Legendre spectral schemes for nonlinear space fractional Schrödinger equation

In the paper, we first propose a Crank-Nicolson Galerkin-Legendre (CN-GL) spectral scheme for the one-dimensional nonlinear space fractional Schrödinger equation. Convergence with spectral accuracy is proved for the spectral approximation. Further, a Crank-Nicolson ADI Galerkin-Legendre spectral method for the two-dimensional nonlinear space fractional Schrödinger equation is developed. The proposed schemes are shown to be efficient with second-order accuracy in time and spectral accuracy in space which are higher than some recently studied methods. Moreover, some numerical results are demonstrated to justify the theoretical analysis.     

The one‐dimensional Schrödinger operator on bounded time scales

In this paper, we consider the one‐dimensional Schrödinger operator on bounded time scales. We construct a space of boundary values of the minimal operator and describe all maximal dissipative, maximal accretive, self‐adjoint, and other extensions of the dissipative Schrödinger operators in terms of boundary conditions. In particular, using Lidskii's theorem, we prove a theorem on completeness of the system of root vectors of the dissipative Schrödinger operators on bounded time scales. Copyright © 2016 John Wiley & Sons, Ltd.     

Splitting Integrators for Nonlinear Schrödinger Equations Over Long Times

Conservation properties of a full discretization via a spectral semi-discretization in space and a Lie–Trotter splitting in time for cubic Schrödinger equations with small initial data (or small nonlinearity) are studied. The approximate conservation of the actions of the linear Schrödinger equation, energy, and momentum over long times is shown using modulated Fourier expansions. The results are valid in arbitrary spatial dimension.     

Numerical solution of the heat equation with nonlinear boundary conditions in unbounded domains

The numerical solution of the heat equation in unbounded domains (for a 1D problem‐semi‐infinite line and for a 2D one semi‐infinite strip) is considered. The artificial boundaries are introduced and the exact artificial boundary conditions are derived. The original problems are transformed into problems on finite domains. The space semi‐discretization by finite element method and the full approximation by the implicit‐explicit Euler's method are presented. The solvability of the full discretization schemes is analyzed. Computational examples demonstrate the accuracy and the efficiency of the algorithms. Also, the behavior of blowing up solutions is examined numerically. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 23: 379–399, 2007     

Spectral problems with mixed Dirichlet–Neumann boundary conditions: Isospectrality and beyond

In this study, an implicit semi-discrete higher order compact (HOC) scheme, with an averaged time discretization, has been presented for the numerical solution of unsteady two-dimensional (2D) Schrödinger equation. The scheme is second order accurate in time and fourth order accurate in space. The results of numerical experiments are presented, and are compared with analytical solutions and well established numerical results of some other finite difference schemes. In all cases, the present scheme produces highly accurate results with much better computational efficiency.     

Pathfollowing for essentially singular boundary value problems with application to the complex Ginzburg-Landau equation

We present a pathfollowing strategy based on pseudo-arclength parametrization for the solution of parameter-dependent boundary value problems for ordinary differential equations. We formulate criteria which ensure the successful application of this method for the computation of solution branches with turning points for problems with an essential singularity. The advantages of our approach result from the possibility to use efficient mesh selection, and a favorable conditioning even for problems posed on a semi-infinite interval and subsequently transformed to an essentially singular problem. This is demonstrated by a Matlab implementation of the solution method based on an adaptive collocation scheme which is well suited to solve problems of practical relevance. As one example, we compute solution branches for the complex Ginzburg-Landau equation which start from non-monotone ‘multi-bump’ solutions of the nonlinear Schrödinger equation. Following the branches around turning points, real-valued solutions of the nonlinear Schrödinger equation can easily be computed.     

A Novel Time‐Domain BIEM for Wave Propagation Analysis in Anisotropic Solids With Cracks

Analysis of elastic wave propagation in anisotropic solids with cracks is of particular interest to quantitative non‐destructive testing and fracture mechanics. For this purpose, a novel time‐domain boundary integral equation method (BIEM) is presented in this paper. A finite crack in an unbounded elastic solid of general anisotropy subjected to transient elastic wave loading is considered. Two‐dimensional plane strain or plane stress condition is assumed. The initial‐boundary value problem is formulated as a set of hypersingular time‐domain traction boundary integral equations (BIEs) with the crack‐opening‐displacements (CODs) as unknown quantities. A time‐stepping scheme is developed for solving the hypersingular time‐domain BIEs. The scheme uses the convolution quadrature formula of Lubich [1] for temporal convolution and a Galerkin method for spatial discretization of the BIEs. An important feature of the present time‐domain BIEM is that it uses the Laplace‐domain instead of the more complicated time‐domain Green's functions. Fourier integral representations of Laplace‐domain Green's functions are applied. No special technique is needed in the present time‐domain BIEM for evaluating hypersingular integrals.     

Symplectic and multi-symplectic wavelet collocation methods for two-dimensional Schrödinger equations

In this paper, we develop symplectic and multi-symplectic wavelet collocation methods to solve the two-dimensional nonlinear Schrödinger equation in wave propagation problems and the two-dimensional time-dependent linear Schrödinger equation in quantum physics. The Hamiltonian and the multi-symplectic formulations of each equation are considered. For both formulations, wavelet collocation method based on the autocorrelation function of Daubechies scaling functions is applied for spatial discretization and symplectic method is used for time integration. The conservation of energy and total norm is investigated. Combined with splitting scheme, splitting symplectic and multi-symplectic wavelet collocation methods are also constructed. Numerical experiments show the effectiveness of the proposed methods.     

From the AKNS system to the matrix Schrödinger equation with vanishing potentials: Direct and inverse problems

We relate the scattering theory of the focusing AKNS system with vanishing boundary conditions to that of the matrix Schrödinger equation. The corresponding Miura transformation, which allows this connection, converts the focusing matrix nonlinear Schrödinger (NLS) equation into a new nonlocal integrable equation. We apply the matrix triplet method to derive the multisoliton solutions of the nonlocal integrable equation, thus proposing a new method to solve the matrix NLS equation.