Abstract:Coupled nonlinear Schrödinger equations model various physical phenomena, such as wave propagation in nonlinear optics, multi-component Bose-Einstein condensates, and shallow water waves. Despite their extensive applications, analytical solutions of coupled nonlinear Schrödinger equations are widely either unknown or challenging to compute, prompting the need for stable and efficient numerical methods to understand the nonlinear phenomenon and complex dynamics of systems governed by coupled nonlinear Schrödinger equations. This paper explores the use of the fourth-order Runge-Kutta based exponential time-differencing and integrating factor methods combined with the Fourier spectral method to simulate multi-dimensional M-coupled nonlinear Schrödinger equations.
The theoretical derivation and stability of the methods, as well as the runtime complexity of the algorithms used for their implementation, are examined. Numerical experiments are performed on systems of two and four multi-dimensional coupled nonlinear Schrödinger equations. It is demonstrated by the results that both methods effectively conserve mass and energy while maintaining fourth-order temporal and spectral spatial convergence. Overall, it is shown by the numerical results that the exponential time-differencing method outperforms the integrating factor method in this application, and both may be considered further in modeling more nonlinear dynamics in future work.
What problem does this paper attempt to address?
This paper attempts to solve the problems of accuracy and efficiency in the numerical simulation of multi - dimensional M - coupled nonlinear Schrödinger equations (CNLSE). Specifically, the article focuses on how to use high - order time - integration methods and Fourier spectral methods to simulate these equations efficiently and accurately.
### Core problems of the paper
1. **Limitations of analytical solutions**: Coupled nonlinear Schrödinger equations are widely used in modeling physical phenomena, such as wave propagation in nonlinear optics, multi - component Bose - Einstein condensates, and shallow water waves. However, in most cases, the analytical solutions of these equations are either unknown or difficult to calculate. Therefore, it is necessary to develop stable and efficient numerical methods to understand and simulate the complex dynamical systems described by these equations.
2. **Deficiencies of existing numerical methods**: Although there are various numerical methods for solving CNLSE, they still have limitations when dealing with multi - dimensional, non - integrable or chaotic systems. For example, finite difference methods, finite element methods, operator splitting methods, exponential time - difference methods (ETD), integrating factor methods (IF), spectral methods and pseudo - spectral methods each have their own advantages and disadvantages, but may not be efficient or stable enough under certain specific conditions.
3. **Proposed new methods**: This paper proposes and evaluates two time - stepping methods based on the fourth - order Runge - Kutta method, combined with the Fourier spectral method to solve multi - dimensional M - coupled nonlinear Schrödinger equations. These two methods are:
- **Krogstad - P22 method**: An improved fourth - order exponential time - difference Runge - Kutta method (ETDRK4 - B), which uses Padé-(2, 2) approximation to handle matrix exponential terms.
- **IFRK4 - P13 method**: An improved fourth - order integrating factor Runge - Kutta method, which uses L - acceptable (1, 3)-Padé approximation to handle matrix exponential terms.
### Research objectives
- **Verify the effectiveness of the methods**: Through theoretical derivation and stability analysis, verify the effectiveness of the two proposed methods in solving multi - dimensional M - coupled nonlinear Schrödinger equations.
- **Evaluate the computational performance**: Through numerical experiments, evaluate the accuracy and computational efficiency of these two methods under different boundary conditions, and compare their advantages and disadvantages with existing methods.
- **Conservative properties**: Check whether these methods can effectively maintain mass conservation and energy conservation, which are important properties for simulating nonlinear waves.
### Conclusions
The numerical experiment results show that both methods can effectively maintain mass and energy conservation, and have fourth - order time convergence and spectral - space convergence. In particular, the exponential time - difference method (Krogstad - P22) shows better performance than the integrating factor method (IFRK4 - P13) in this application. These results provide strong support for the future modeling of more complex nonlinear dynamical systems.
In summary, this paper aims to improve the simulation accuracy and efficiency of multi - dimensional M - coupled nonlinear Schrödinger equations by introducing and evaluating new numerical methods, so as to better understand the physical phenomena described by these equations.