A new fourth-order Fourier-Bessel split-step method for the extended nonlinear Schrödinger equation

DOI10.1016/J.JCP.2007.10.012MaRDI QIDQ2476871zbMATH OpenOpenAlexFDO

Authors Patrick L. Nash

Publication date 12 March 2008

Published in Journal of Computational Physics (Search for Journal in Brave)

Full work available at URL https://doi.org/10.1016/j.jcp.2007.10.012

zbMATH Keywords

associated Bessel functionsfinite differenceFourier-Bessel algorithmnonlinear Schrödinger equationnormal mode condensationnumerical examplesshort laser pulsessoliton-like numerical solutionssplit-step methodstability

Mathematics Subject Classification ID

Spectral, collocation and related methods for initial value and initial-boundary value problems involving PDEs (65M70)

Review text

Fourier split-step techniques are often used to compute soliton-like numerical solutions of the nonlinear Schrödinger equation. Here, a new fourth-order implementation of the Fourier split-step algorithm is described for problems possessing azimuthal symmetry in \(3 + 1\)-dimensions. This implementation is based, in part, on a finite difference approximation \(\Delta^{\text{FDA}}_\perp\) of \(\frac1r\,\frac{\partial}{\partial r}\,r\,\frac{\partial}{\partial r}\) that possesses an associated exact unitary representation of \(e^{\frac i2\lambda\Delta^{\text{FDA}}_\perp}\). The matrix elements of this unitary matrix are given by special functions known as the associated Bessel functions. Hence the attribute Fourier-Bessel for the method. The Fourier-Bessel algorithm is shown to be unitary and unconditionally stable. The Fourier-Bessel algorithm is employed to simulate the propagation of a periodic series of short laser pulses through a nonlinear medium. This numerical simulation calculates waveform intensity profiles in a sequence of planes that are transverse to the general propagation direction, and labeled by the cylindrical coordinate \(z\). These profiles exhibit a series of isolated pulses that are offset from the time origin by characteristic times, and provide evidence for a physical effect that may be loosely termed normal mode condensation. Normal mode condensation is consistent with experimentally observed pulse filamentation into a packet of short bursts, which may occur as a result of short, intense irradiation of a medium.

Recommendations

• Split-Step Methods for the Solution of the Nonlinear Schrödinger Equation
• Instability of the finite-difference split-step method applied to the nonlinear Schrödinger equation. I. Standing soliton
• Higher-order split-step Fourier schemes for the generalized nonlinear Schrödinger equation
• Instability analysis of the Split-step Fourier method on the background of a soliton of the nonlinear Schrödinger equation
• scientific article; zbMATH DE number 2110740

Cites work

• Associated Bessel functions and the discrete approximation of the free-particle time evolution operator in cylindrical coordinates
• Dynamics of Rotating Bose--Einstein Condensates and its Efficient and Accurate Numerical Computation
• Ground States and Dynamics of Multicomponent Bose--Einstein Condensates
• Higher-order split-step Fourier schemes for the generalized nonlinear Schrödinger equation
• scientific article; zbMATH DE number 3457423 (Why is no real title available?)
• Self-Focusing in the Perturbed and Unperturbed Nonlinear Schrödinger Equation in Critical Dimension
• Self-guiding of femtosecond light pulses in condensed media: plasma generation versus chromatic dispersion
• Solution of the Schrödinger equation by a spectral method

Cited in

(1)

• A fourth-order explicit schemes for the coupled nonlinear Schrödinger equation

Describes a project that uses

Uses Software

• FFTW

This page was built for publication: A new fourth-order Fourier-Bessel split-step method for the extended nonlinear Schrödinger equation

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2476871)