Numerical methods for generalized nonlinear Schrödinger equations

We present and analyze different splitting algorithms for numerical solution of the both classical and generalized nonlinear Schr"odinger equations describing propagation of wave packets with special emphasis on applications to nonlinear fiber-optics. The considered generalizations take into account the higher-order corrections of the linear differential dispersion operator as well as the saturation of nonlinearity and the self-steepening of the field envelope function. For stabilization of the pseudo-spectral splitting schemes for generalized Schr"odinger equations a regularization based on the approximation of the derivatives by the low number of Fourier modes is proposed. To illustrate the theoretically predicted performance of these schemes several numerical experiments have been done

ABC Method and Fractional Momentum Layer for the FDTD Method to Solve the Schrödinger Equation on Unbounded Domains

The finite­difference time­domain (FDTD) method and its generalized variant (G­FDTD) are efficient numerical tools for solving the linear and nonlinear Schrödinger equations because not only are they explicit, allowing parallelization, but they also provide high­order accuracy with relatively inexpensive computational costs. In addition, the G­FDTD method has a relaxed stability condition when compared to the original FDTD method. It is important to note that the existing simulations of the G­FDTD scheme employed analytical solutions to obtain function values at the points along the boundary; however, in simulations for which the analytical solution is unknown, theoretical approximations for values at points along the boundary are desperately needed. Hence, the objective of this dissertation research is to develop absorbing boundary conditions (ABCs) so that the G­FDTD method can be used to solve the nonlinear Schrödinger equation when the analytical solution is unknown. To create the ABCs for the nonlinear Schrödinger equation, we initially determine the associated Engquist­Majda one­way wave equations and then proceed to develop a finite difference scheme for them. These ABCs are made to be adaptive using a windowed Fourier transform to estimate a value of the wavenumber of the carrier wave. These ABCs were tested using the nonlinear Schrödinger equation for 1D and 2D soliton propagation as well as Gaussian packet collision and dipole radiation. Results show that these ABCs perform well, but they have three key limitations. First, there are inherent reflections at the interface of the interior and boundary domains due to the different schemes used the two regions; second, to use the ABCs, one needs to estimate a value for the carrier wavenumber and poor estimates can cause even more reflection at the interface; and finally, the ABCs require different schemes in different regions of the boundary, and this domain decomposition makes the ABCs tedious both to develop and to implement. To address these limitations for the FDTD method, we employ the fractional­order derivative concept to unify the Schrödinger equation with its one­way wave equation over an interval where the fractional order is allowed to vary. Through careful construction of a variable­order fractional momentum operator, outgoing waves may enter the fractionalorder region with little to no reflection and, inside this region, any reflected portions of the wave will decay exponentially with time. The fractional momentum operator is then used to create a fractional­order FDTD scheme. Importantly, this single scheme can be used for the entire computational domain, and the scheme smooths the abrupt transition between the FDTD method and the ABCs. Furthermore, the fractional FDTD scheme relaxes the precision needed for the estimated carrier wavenumber. This fractional FDTD scheme is tested for both the linear and nonlinear Schrödinger equations. Example cases include a 1D Gaussian packet scattering off of a potential, a 1D soliton propagating to the right, as well as 2D soliton propagation, and the collision of Gaussian packets. Results show that the fractional FDTD method outperforms the FDTD method with ABCs

Numerical investigation on nonlocal problems with the fractional Laplacian

Nonlocal models have recently become a powerful tool for studying complex systems with long-range interactions or memory effects, which cannot be described properly by the traditional differential equations. So far, different nonlocal (or fractional differential) models have been proposed, among which models with the fractional Laplacian have been well applied. The fractional Laplacian (-Δ)α/2 represents the infinitesimal generator of a symmetric α-stable Lévy process. It has been used to describe anomalous diffusion, turbulent flows, stochastic dynamics, finance, and many other phenomena. However, the nonlocality of the fractional Laplacian introduces considerable challenges in its mathematical modeling, numerical simulations, and mathematical analysis. To advance the understanding of the fractional Laplacian, two novel and accurate finite difference methods -- the weighted trapezoidal method and the weighted linear interpolation method are developed for discretizing the fractional Laplacian. Numerical analysis is provided for the error estimates, and fast algorithms are developed for their efficient implementation. Compared to the current state-of-the-art methods, these two methods have higher accuracy but less computational complexity. As an application, the solution behaviors of the fractional Schördinger equation are investigated to understand the nonlocal effects of the fractional Laplacian. First, the eigenvalues and eigenfunctions of the fractional Schrödinger equation in an infinite potential well are studied, and the results provide insights into an open problem in the fractional quantum mechanics. Second, three Fourier spectral methods are developed and compared in solving the fractional nonlinear Schördinger equation (NLS), among which the SSFS method is more effective in the study of the plane wave dynamics. Sufficient conditions are provided to avoid the numerical instability of the SSFS method. In contrast to the standard NLS, the plane wave dynamics of the fractional NLS are more chaotic due to the long-range interactions --Abstract, page iii

SCALAR: an AMR code to simulate axion-like dark matter models

We present a new code, SCALAR, based on the high-resolution hydrodynamics and N-body code RAMSES, to solve the Schr\"odinger equation on adaptive refined meshes. The code is intended to be used to simulate axion or fuzzy dark matter models where the evolution of the dark matter component is determined by a coupled Schr\"odinger-Poisson equation, but it can also be used as a standalone solver for both linear and non-linear Schr\"odinger equations with any given external potential. This paper describes the numerical implementation of our solver and presents tests to demonstrate how accurately it operates.Comment: 17 pages, 11 figure

Generalization of splitting methods based on modified potentials to nonlinear evolution equations of parabolic and Schr\"odinger type

The present work is concerned with the extension of modified potential operator splitting methods to specific classes of nonlinear evolution equations. The considered partial differential equations of Schr{\"o}dinger and parabolic type comprise the Laplacian, a potential acting as multiplication operator, and a cubic nonlinearity. Moreover, an invariance principle is deduced that has a significant impact on the efficient realisation of the resulting modified operator splitting methods for the Schr{\"o}dinger case.} Numerical illustrations for the time-dependent Gross--Pitaevskii equation in the physically most relevant case of three space dimensions and for its parabolic counterpart related to ground state and excited state computations confirm the benefits of the proposed fourth-order modified operator splitting method in comparison with standard splitting methods. The presented results are novel and of particular interest from both, a theoretical perspective to inspire future investigations of modified operator splitting methods for other classes of nonlinear evolution equations and a practical perspective to advance the reliable and efficient simulation of Gross--Pitaevskii systems in real and imaginary time.Comment: 30 pages, 6 figure