Numerical solution of coupled nonlinear Schrödinger equation by Galerkin method

The coupled nonlinear Schrödinger equation models several interesting physical phenomena presents a model equation for optical fiber with linear birefringence. In this paper we derive a finite element scheme to solve this equation, we test this method for stability and accuracy, many numerical tests have been conducted. The scheme is quite accurate and describe the interaction picture clearly.     

Multi-symplectic splitting method for the coupled nonlinear Schrödinger equation

In this paper, we propose a multi-symplectic splitting method to solve the coupled nonlinear Schrödinger (CNLS) equation by using the idea of splitting the multi-symplectic partial differential equation (PDE). Numerical experiments show that the proposed method can simulate the propagation and collision of solitons well. The corresponding errors in global energy and momentum are also presented to show the good preservation property of the proposed method during long-time numerical calculation.     

A Chebyshev pseudospectral multidomain method for the soliton solution of coupled nonlinear Schrödinger equations

The collision of solitary waves is an important problem in both physics and applied mathematics. In this paper, we study the solution of coupled nonlinear Schrödinger equations based on pseudospectral collocation method with domain decomposition algorithm for approximating the spatial variable. The problem is converted to a system of nonlinear ordinary differential equations which will be integrated in time by explicit Runge–Kutta method of order four. The multidomain scheme has much better stability properties than the single domain. Thus this permits using much larger step size for the time integration which fulfills stability restrictions. The proposed scheme reduces the effects of round-of-error for the Chebyshev collocation and also uses less memory without sacrificing the accuracy. The numerical experiments are presented which show the multidomain pseudospectral method has excellent long-time numerical behavior and preserves energy conservation property.     

Multi-symplectic methods for the coupled 1D nonlinear Schrödinger system

In the paper, the multi-symplectic formulation of the coupled 1D nonlinear Schrödinger system (CNLS) is considered. For the multi-symplectic formulation, a new six point scheme, which is equivalent to the multi-symplectic Preissman integrator, is derived. We also present numerical experiments, which show that the multi-symplectic scheme has excellent long-time numerical behaviour and energy conservation property.     

The global solution for a class of systems of fractional nonlinear Schrödinger equations with periodic boundary condition

In this paper, we consider a class of systems of fractional nonlinear Schrödinger equations. We prove the existence and uniqueness of the global solution to the periodic boundary value problem by using the Faedo-Galërkin method.     

Numerical simulation of a nonlinearly coupled Schrödinger system: A linearly uncoupled finite difference scheme

This article describes a finite difference scheme which is linearly uncoupled in computation for a nonlinearly coupled Schrödinger system. This numerical scheme is proved to preserve the original conservative properties. Using the discrete energy analysis method, we also prove that the scheme is unconditionally stable and second-order convergent in discrete _L2$L_2$-norm based on some preliminary estimations. The results show that the new scheme is efficiency.     

Time-splitting pseudo-spectral domain decomposition method for the soliton solutions of the one- and multi-dimensional nonlinear Schrödinger equations

In this paper, we study the simulation of nonlinear Schrödinger equation in one, two and three dimensions. The proposed method is based on a time-splitting method that decomposes the original problem into two parts, a linear equation and a nonlinear equation. The linear equation in one dimension is approximated with the Chebyshev pseudo-spectral collocation method in space variable and the Crank–Nicolson method in time; while the nonlinear equation with constant coefficients can be solved exactly. As the goal of the present paper is to study the nonlinear Schrödinger equation in the large finite domain, we propose a domain decomposition method. In comparison with the single-domain, the multi-domain methods can produce a sparse differentiation matrix with fewer memory space and less computations. In this study, we choose an overlapping multi-domain scheme. By applying the alternating direction implicit technique, we extend this efficient method to solve the nonlinear Schrödinger equation both in two and three dimensions, while for the solution at each time step, it only needs to solve a sequence of linear partial differential equations in one dimension, respectively. Several examples for one- and multi-dimensional nonlinear Schrödinger equations are presented to demonstrate high accuracy and capability of the proposed method. Some numerical experiments are reported which show that this scheme preserves the conservation laws of charge and energy.     

A semi-explicit multi-symplectic splitting scheme for a 3-coupled nonlinear Schrödinger equation

In this paper, we mainly propose an efficient semi-explicit multi-symplectic splitting scheme to solve a 3-coupled nonlinear Schrödinger (3-CNLS) equation. Based on its multi-symplectic formulation, the 3-CNLS equation can be split into one linear multi-symplectic subsystem and one nonlinear infinite-dimensional Hamiltonian subsystem. For the linear subsystem, the multi-symplectic Fourier pseudospectral method and symplectic Euler method are employed in spatial and temporal discretizations, respectively. For the nonlinear subsystem, the mid-point symplectic scheme is used. Numerical experiments for the unstable plane waves show the effectiveness of the proposed method during long-time numerical calculation.     

A compact split-step finite difference method for solving the nonlinear Schrödinger equations with constant and variable coefficients

We propose a compact split-step finite difference method to solve the nonlinear Schrödinger equations with constant and variable coefficients. This method improves the accuracy of split-step finite difference method by introducing a compact scheme for discretization of space variable while this improvement does not reduce the stability range and does not increase the computational cost. This method also preserves some conservation laws. Numerical tests are presented to confirm the theoretical results for the new numerical method by using the cubic nonlinear Schrödinger equation with constant and variable coefficients and Gross-Pitaevskii equation.     

Numerical solution to coupled nonlinear Schrödinger equations on unbounded domains

The numerical simulation of coupled nonlinear Schrödinger equations on unbounded domains is considered in this paper. By using the operator splitting technique, the original problem is decomposed into linear and nonlinear subproblems in a small time step. The linear subproblem turns out to be two decoupled linear Schrödinger equations on unbounded domains, where artificial boundaries are introduced to truncate the unbounded physical domains into finite ones. Local absorbing boundary conditions are imposed on the artificial boundaries. On the other hand, the coupled nonlinear subproblem is an ODE system, which can be solved exactly. To demonstrate the effectiveness of our method, some comparisons in terms of accuracy and computational cost are made between the PML approach and our method in numerical examples.     

The accurate numerical solution of the Schrödinger equation with an explicitly time-dependent Hamiltonian

We show how the highly accurate and efficient Constant Perturbation (CP) technique for steady-state Schrödinger problems can be used in the solution of time-dependent Schrödinger problems with explicitly time-dependent Hamiltonians, following a technique suggested by Ixaru (2010). By introducing a sectorwise spatial discretization using bases of accurately CP-computed eigenfunctions of carefully-chosen stationary problems, we deal with the possible highly oscillatory behavior of the wave function while keeping the dimension of the resulting ODE system low. Also for the time-integration of the ODE system a very effective CP-based approach can be used.     

Solitary wave solutions for a higher order nonlinear Schrödinger equation

We consider a higher order nonlinear Schrödinger equation with third- and fourth-order dispersions, cubic–quintic nonlinearities, self steepening, and self-frequency shift effects. This model governs the propagation of femtosecond light pulses in optical fibers. In this paper, we investigate general analytic solitary wave solutions and derive explicit bright and dark solitons for the considered model. The derived analytical dark and bright wave solutions are expressed in terms of the model coefficients. These exact solutions are useful to understand the mechanism of the complicated nonlinear physical phenomena which are related to wave propagation in a higher-order nonlinear and dispersive Schrödinger system.     

CP methods of higher order for Sturm-Liouville and Schrödinger equations

The algorithm upon which the code SLCPM12, described in Computer Physics Communications 118 (1999) 259-277, is based, is extended to higher order. The implementation of the original algorithm, which was of order {12,10} (meaning order 12 at low energies and order 10 at high energies), was more efficient than the well-established codes SL02F, SLEDGE and SLEIGN. In the new algorithm the orders {14,12}, {16,14} and {18,16} are introduced. Besides regular Sturm-Liouville and one-dimensional Schrödinger problems also radial Schrödinger equations are considered with potentials of the form V(r)=S(r)/r+R(r), where S(r) and R(r) are well behaved functions which tend to some (not necessarily equal) constants when r→0 and r→∞. Numerical illustrations are given showing the accuracy, the robustness and the CPU-time gain of the proposed algorithms.     

Multiresolution scheme for Time-Dependent Schrödinger Equation

This paper is devoted to a multiresolution approach for solving laser-molecule Time-Dependent Schrödinger Equations (TDSE) in strong and high frequency fields. It is well known that short and intense laser-molecule interactions lead to complex nonlinear phenomena that necessitate an accurate numerical approximation of the TDSE. In particular, intense electric fields rapidly delocalize molecule wavefunctions so that their support can vary a lot during the interaction. In this kind of physical configurations, mesh adaption is a usual compromise between precision and computational efficiency. We then propose to explore numerically mesh adaptation for TDSE using a multiresolution analysis coupled with a Crank-Nicolson-like scheme. We then discuss the efficiency and the drawbacks of such a strategy.     

Domain wall arrays, fronts, and bright and dark solitons in a generalized derivative nonlinear Schrödinger equation

We obtain some exact solutions of a generalized derivative nonlinear Schrödinger equation, including domain wall arrays (periodic solutions in terms of elliptic functions), fronts, and bright and dark solitons. In certain parameter domains, fundamental bright and dark solitons are chiral, and the propagation direction is determined by the sign of the self-steepening parameter. Moreover, we also find the chirping reversal phenomena of fronts, and bright and dark solitons, and discuss two different ways to produce the chirping reversal.     

Soliton dynamics in linearly coupled discrete nonlinear Schrödinger equations

We study soliton dynamics in a system of two linearly coupled discrete nonlinear Schrödinger equations, which describe the dynamics of a two-component Bose gas, coupled by an electromagnetic field, and confined in a strong optical lattice. When the nonlinear coupling strengths are equal, we use a unitary transformation to remove the linear coupling terms, and show that the existing soliton solutions oscillate from one species to the other. When the nonlinear coupling strengths are different, the soliton dynamics is numerically investigated and the findings are compared to the results of an effective two-mode model. The case of two linearly coupled Ablowitz–Ladik equations is also briefly discussed.     

An algebraic method for Schrödinger equations in quaternionic quantum mechanics

In the study of theory and numerical computations of quaternionic quantum mechanics and quantum chemistry, one of the most important tasks is to solve the Schrödinger equation with A an anti-self-adjoint real quaternion matrix, and |f〉 an eigenstate to A. The quaternionic Schrödinger equation plays an important role in quaternionic quantum mechanics, and it is known that the study of the quaternionic Schrödinger equation is reduced to the study of quaternionic eigen-equation Aα=αλ with A an anti-self-adjoint real quaternion matrix (time-independent). This paper, by means of complex representation of quaternion matrices, introduces concepts of norms of quaternion matrices, studies the problems of quaternionic Least Squares eigenproblem, and give a practical algebraic technique of computing approximate eigenvalues and eigenvectors of a quaternion matrix in quaternionic quantum mechanics.     

A dissipative exponentially-fitted method for the numerical solution of the Schrödinger equation and related problems

In this paper a dissipative exponentially-fitted method for the numerical integration of the Schrödinger equation and related problems is developed. The method is called dissipative since is a nonsymmetric multistep method. An application to the the resonance problem of the radial Schrödinger equation and to other well known related problems indicates that the new method is more efficient than the corresponding classical dissipative method and other well known methods. Based on the new method and the method of Raptis and Cash a new variable-step method is obtained. The application of the new variable-step method to the coupled differential equations arising from the Schrödinger equation indicates the power of the new approach.     

Multi-symplectic wavelet collocation method for the nonlinear Schrödinger equation and the Camassa–Holm equation

In this paper, we develop a novel multi-symplectic wavelet collocation method for solving multi-symplectic Hamiltonian system with periodic boundary conditions. Based on the autocorrelation function of Daubechies scaling functions, collocation method is conducted for the spatial discretization. The obtained semi-discrete system is proved to have semi-discrete multi-symplectic conservation laws and semi-discrete energy conservation laws. Then, appropriate symplectic scheme is applied for time integration, which leads to full-discrete multi-symplectic conservation laws. Numerical experiments for the nonlinear Schrödinger equation and Camassa–Holm equation show the high accuracy, effectiveness and good conservation properties of the proposed method.